Search

Molecular switch proteins whose cycling between states is controlled by opposing regulators1,2 are central to biological signal transduction. As switch proteins function within highly connected interaction networks3, the fundamental question arises of how functional specificity is achieved when different processes share common regulators. Here we show that functional specificity of the small GTPase switch protein Gsp1 in Saccharomyces cerevisiae (the homologue of the human protein RAN)4 is linked to differential sensitivity of biological processes to different kinetics of the Gsp1 (RAN) switch cycle. We make 55 targeted point mutations to individual protein interaction interfaces of Gsp1 (RAN) and show through quantitative genetic5 and physical interaction mapping that Gsp1 (RAN) interface perturbations have widespread cellular consequences. Contrary to expectation, the cellular effects of the interface mutations group by their biophysical effects on kinetic parameters of the GTPase switch cycle and not by the targeted interfaces. Instead, we show that interface mutations allosterically tune the GTPase cycle kinetics. These results suggest a model in which protein partner binding, or post-translational modifications at distal sites, could act as allosteric regulators of GTPase switching. Similar mechanisms may underlie regulation by other GTPases, and other biological switches. Furthermore, our integrative platform to determine the quantitative consequences of molecular perturbations may help to explain the effects of disease mutations that target central molecular switches. Interface mutations in the GTPase switch protein Gsp1 (the yeast homologue of human RAN) allosterically affect the kinetics of the switch cycle, revealing a systems-level mechanism of multi-specificity.

Allosteric regulation is central to protein function in cellular networks1. However, despite technological advances2,3 most studies of allosteric effects on function are conducted in heterologous environments2,4,5, limiting the discovery of allosteric mechanisms that rely on endogenous binding partners or posttranslational modifications to modulate activity. Here we report an approach that enables probing of new sites of allosteric regulation at residue-level resolution in essential eukaryotic proteins in their native biological context by comprehensive mutational scanning. We apply our approach to the central GTPase Gsp1/Ran. GTPases are highly regulated molecular switches that control signaling, with switching occurring via catalyzed GTP hydrolysis and nucleotide exchange. We find that 28% of 4,315 assayed mutations in Gsp1/Ran are highly deleterious, showing a toxic response identified by our assay as gain-of-function (GOF). Remarkably, a third of all positions enriched for GOF mutations (20/60) are outside the GTPase active site. Kinetic analysis shows that these distal sites are allosterically coupled to the active site, including a novel cluster of sites that alter the nucleotide preference of Gsp1 from GDP to GTP. We describe multiple distinct mechanisms by which allosteric mutations alter Gsp1/Ran cellular function by modulating GTPase switching. Our systematic discovery of new regulatory sites provides a functional map relevant to other GTPases such as Ras that could be exploited for targeting and reprogramming critical biological processes.

A newly described coronavirus named severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which is the causative agent of coronavirus disease 2019 (COVID-19), has infected over 2.3 million people, led to the death of more than 160,000 individuals and caused worldwide social and economic disruption1,2. There are no antiviral drugs with proven clinical efficacy for the treatment of COVID-19, nor are there any vaccines that prevent infection with SARS-CoV-2, and efforts to develop drugs and vaccines are hampered by the limited knowledge of the molecular details of how SARS-CoV-2 infects cells. Here we cloned, tagged and expressed 26 of the 29 SARS-CoV-2 proteins in human cells and identified the human proteins that physically associated with each of the SARS-CoV-2 proteins using affinity-purification mass spectrometry, identifying 332 high-confidence protein–protein interactions between SARS-CoV-2 and human proteins. Among these, we identify 66 druggable human proteins or host factors targeted by 69 compounds (of which, 29 drugs are approved by the US Food and Drug Administration, 12 are in clinical trials and 28 are preclinical compounds). We screened a subset of these in multiple viral assays and found two sets of pharmacological agents that displayed antiviral activity: inhibitors of mRNA translation and predicted regulators of the sigma-1 and sigma-2 receptors. Further studies of these host-factor-targeting agents, including their combination with drugs that directly target viral enzymes, could lead to a therapeutic regimen to treat COVID-19. A human–SARS-CoV-2 protein interaction map highlights cellular processes that are hijacked by the virus and that can be targeted by existing drugs, including inhibitors of mRNA translation and predicted regulators of the sigma receptors.

Author(s): Gordon, David E; Jang, Gwendolyn M; Bouhaddou, Mehdi; Xu, Jiewei; Obernier, Kirsten; O'Meara, Matthew J; Guo, Jeffrey Z; Swaney, Danielle L; Tummino, Tia A; Huttenhain, Ruth; Kaake, Robyn M; Richards, Alicia L; Tutuncuoglu, Beril; Foussard, Helene; Batra, Jyoti; Haas, Kelsey; Modak, Maya; Kim, Minkyu; Haas, Paige; Polacco, Benjamin J; Braberg, Hannes; Fabius, Jacqueline M; Eckhardt, Manon; Soucheray, Margaret; Bennett, Melanie J; Cakir, Merve; McGregor, Michael J; Li, Qiongyu; Naing, Zun Zar Chi; Zhou, Yuan; Peng, Shiming; Kirby, Ilsa T; Melnyk, James E; Chorba, John S; Lou, Kevin; Dai, Shizhong A; Shen, Wenqi; Shi, Ying; Zhang, Ziyang; Barrio-Hernandez, Inigo; Memon, Danish; Hernandez-Armenta, Claudia; Mathy, Christopher JP; Perica, Tina; Pilla, Kala B; Ganesan, Sai J; Saltzberg, Daniel J; Ramachandran, Rakesh; Liu, Xi; Rosenthal, Sara B; Calviello, Lorenzo; Venkataramanan, Srivats; Lin, Yizhu; Wankowicz, Stephanie A; Bohn, Markus; Trenker, Raphael; Young, Janet M; Cavero, Devin; Hiatt, Joe; Roth, Theo; Rathore, Ujjwal; Subramanian, Advait; Noack, Julia; Hubert, Mathieu; Roesch, Ferdinand; Vallet, Thomas; Meyer, Bjorn; White, Kris M; Miorin, Lisa; Agard, David; Emerman, Michael; Ruggero, Davide; Garcia-Sastre, Adolfo; Jura, Natalia; von Zastrow, Mark; Taunton, Jack; Schwartz, Olivier; Vignuzzi, Marco; d'Enfert, Christophe; Mukherjee, Shaeri; Jacobson, Matt; Malik, Harmit S; Fujimori, Danica G; Ideker, Trey; Craik, Charles S | Abstract: An outbreak of the novel coronavirus SARS-CoV-2, the causative agent of COVID-19 respiratory disease, has infected over 290,000 people since the end of 2019, killed over 12,000, and caused worldwide social and economic disruption1,2. There are currently no antiviral drugs with proven efficacy nor are there vaccines for its prevention. Unfortunately, the scientific community has little knowledge of the molecular details of SARS-CoV-2 infection. To illuminate this, we cloned, tagged and expressed 26 of the 29 viral proteins in human cells and identified the human proteins physically associated with each using affinity- purification mass spectrometry (AP-MS), which identified 332 high confidence SARS-CoV-2-human protein-protein interactions (PPIs). Among these, we identify 66 druggable human proteins or host factors targeted by 69 existing FDA-approved drugs, drugs in clinical trials and/or preclinical compounds, that we are currently evaluating for efficacy in live SARS-CoV-2 infection assays. The identification of host dependency factors mediating virus infection may provide key insights into effective molecular targets for developing broadly acting antiviral therapeutics against SARS-CoV-2 and other deadly coronavirus strains.

Molecular switches are central to signal transduction in protein interaction networks. One switch protein can independently regulate distinct cellular processes, but the molecular mechanisms enabling this functional multi-specificity remain unclear. Here we integrate system-scale cellular and biophysical measurements to study how a paradigm switch, the small GTPase Ran/Gsp1, achieves its functional multi-specificity. We make 55 targeted point mutations to individual interactions of Ran/Gsp1 and show through quantitative, systematic genetic and physical interaction mapping that Ran/Gsp1 interface perturbations have widespread cellular consequences that cluster by biological processes but, unexpectedly, not by the targeted interactions. Instead, the cellular consequences of the interface mutations group by their biophysical effects on kinetic parameters of the GTPase switch cycle, and cycle kinetics are allosterically tuned by distal interface mutations. We propose that the functional multi-specificity of Ran/Gsp1 is encoded by a differential sensitivity of biological processes to different kinetic parameters of the Gsp1 switch cycle, and that Gsp1 partners binding to the sites of distal mutations act as allosteric regulators of the switch. Similar mechanisms may underlie biological regulation by other GTPases and biological switches. Finally, our integrative platform to determine the quantitative consequences of cellular perturbations may help explain the effects of disease mutations targeting central switches.

GTPase switches are hubs for multiple distinct cell signaling inputs and outputs. In a new study combining genetic and biochemical methods, Perica & Mathy et al. identify an unexpected connection between the kinetics of a GTPase switch cycle and functional specificity.

Here we review mechanisms of protein evolution leading to structural changes in protein complexes. These mechanisms include mutations directly within protein interfaces, as well as the effects of mutations that propagate from distant regions of the protein. We also discuss the constraints protein complex structures impose on sequence evolution. We interpret, wherever possible, these mechanisms using amino acid residue contact networks. Many insights into protein evolution come from studies of monomers, and these results facilitate our understanding of evolution of protein complexes. Finally, we highlight the potential of formalizing a phylogenetic framework to integrate residue evolution, structure evolution, and to quantify changes in residue contact networks in protein families.

All proteins require physical interactions with other proteins in order to perform their functions. Most of them oligomerize into homomers, and a vast majority of these homomers interact with other proteins, at least part of the time, forming transient or obligate heteromers. In the present paper, we review the structural, biophysical and evolutionary aspects of these protein interactions. We discuss how protein function and stability benefit from oligomerization, as well as evolutionary pathways by which oligomers emerge, mostly from the perspective of homomers. Finally, we emphasize the specificities of heteromeric complexes and their structure and evolution. We also discuss two analytical approaches increasingly being used to study protein structures as well as their interactions. First, we review the use of the biological networks and graph theory for analysis of protein interactions and structure. Secondly, we discuss recent advances in techniques for detecting correlated mutations, with the emphasis on their role in identifying pathways of allosteric communication.

Oligomerization plays an important role in the function of many proteins. Thus, understanding, predicting, and, ultimately, engineering oligomerization presents a long-standing interest. From the perspective of structural biology, protein–protein interactions have mainly been analyzed in terms of the biophysical nature and evolution of protein interfaces. Here, our aim is to quantify the importance of the larger structural context of protein interfaces in protein interaction evolution. Specifically, we ask to what extent intersubunit geometry affects oligomerization state. We define a set of structural parameters describing the overall geometry and relative positions of interfaces of homomeric complexes with different oligomeric states. This allows us to quantify the contribution of direct sequence changes in interfaces versus indirect changes outside the interface that affect intersubunit geometry. We find that such indirect, or allosteric mutations affecting intersubunit geometry via indirect mechanisms are as important as interface sequence changes for evolution of oligomeric states.

Introduction Evolution and design of protein complexes are frequently viewed through the lens of amino acid mutations at protein interfaces, but we showed previously that residues distant from interfaces are also commonly involved in the evolution of alternative quaternary structures. We hypothesized that in these protein families, the difference in oligomeric state is due to a change in intersubunit geometry. The indirect mutations would act by changing protein conformation and dynamics, similar to the way in which allosteric small molecules introduce functional conformational change. We refer to these substitutions as “allosteric mutations.” Rationale In this work, we investigate the mechanism of action of allosteric mutations on oligomeric state in the PyrR family of pyrimidine operon attenuators. In this family, an entirely sequence-conserved helix that forms a tetrameric interface in the thermophilic ortholog (BcPyrR) switches to being solvent-exposed in the mesophilic ortholog (BsPyrR). This results in a homodimeric structure in which the two subunits are clearly rotated relative to their orientation in the tetramer. What is the origin of this rotation and the change in quaternary structure? To dissect the role of the 49 substitutions between BsPyrR and BcPyrR, we used ancestral sequence reconstruction in combination with structural and biophysical methods to identify a set of allosteric mutations that are responsible for this shift in conformation. We compared the conformational changes introduced by the mutations to the protein motion during allosteric regulation by guanosine monophosphate (GMP). Results We identified 11 key mutations controlling oligomeric state, all distant from the interfaces and outside ligand-binding pockets. We confirmed the role of these allosteric mutations by engineering a shift in oligomeric state in an inferred ancestral PyrR protein (intermediate in sequence between the extant orthologs). We further used the inferred ancestral states and their mutants to show that the allosteric mutations are part of a downhill adaptation of the PyrR proteins to lower temperatures. We compared the x-ray crystal structures of ancestral and engineered PyrR proteins to the free and GMP-bound structure of the mesophilic BsPyrR, which shifts its equilibrium from dimer to tetramer upon ligand binding. Binding of the allosteric molecule introduces a change in intersubunit geometry that is equivalent to the evolutionary difference in intersubunit geometry between the dimeric and tetrameric homologs. We further find that the difference in oligomeric state is coupled to the difference in intrinsic dynamics of the dimers. Finally, we used the residue-residue contact network approach to show that the residues corresponding to the allosteric mutations undergo large contact rewiring when the intersubunit geometry and, in turn, oligomeric state change, either by GMP binding or by the introduction of allosteric mutations. Conclusion We show that evolution employs the intrinsic dynamics of this protein to toggle a conformational switch in a manner similar to that of small molecules. Shifting the relative populations of different states by subtle modifications is a process central to protein function and, as shown here, also to protein evolution. This suggests that we can learn from evolution and design proteins with multiple conformational states.

Oligomerisation is essential for the function of some proteins. The PyrR family of proteins are involved in pyrimidine operon attenuation. They are regulated by the presence of certain nucleotides such as guanosine monophosphate (GMP), which stabilises the tetrameric state. Notably, some members of this family can adopt a tetrameric oligomerisation state regardless of the presence of the GMP. Perica et al. (2012) previously found that this family, among several others, have differences in sequence outside the oligomeric interfaces and these are sufficient to explain the changes to their subunit geometry and oligomerisation states. Here, we compared the differences in structural flexibility linked to the changes in oligomeric state caused by a) specific mutations and, b) by the presence of bound GMP. We calculated the coarse-grained normal modes of dimeric units of the PyrR dimers and tetramers using an Elastic Network Model implemented in the Molecular Modelling ToolKit. We conducted several analyses including comparing the normal modes of these proteins with each other using the Bhattacharyya Coefficient similarity measure, correlations matrices, and the conformational overlap analysis, by which the contribution of these modes to the transition from one state to another can be quantified. Firstly, we found that while the dynamics were very similar between all the structures, there was a noticeable difference between the dimeric and tetrameric units. We also show that both sets of proteins transition from tetrameric to dimeric states similarly, indicating some overlap between the effects of allostery and evolution on oligomerisation in PyrR.

Most proteins can interact with other proteins to form complexes. Is the dynamic assembly of a protein complex connected to the intrinsic flexibility of its component subunits? Can the order of this assembly process be related to evolutionary changes in quaternary structure? Two methodological advances have allowed us to address these questions. First, we introduced a simple method for predicting the free-state flexibility of protein subunits from the structures of protein complexes based on accessible surface area calculations. Second, we showed that we can predict assembly pathways through analysis of intersubunit interface areas, and validated this method by experimentally characterizing the assembly of several heteromers with macromolecular mass spectrometry. using these new methods, we demonstrate that the intrinsic flexibility of protein subunits is crucial for assembly, as it facilitates the conformational changes that enable both the packing together of multiple distinct subunits and the formation of asymmetric interfaces required for cyclic symmetries. We also find that protein complex subunits have a strong tendency to assemble in an order corresponding to their flexibility and magnitude of conformational change, thus suggesting a possible mechanism for control of ordered assembly. Moreover, we reveal that changes in quaternary structure driven by evolutionary gene fusion events tend to occur in such a way so that they mimic the existing order of subunit assembly. This shows that there is selective pressure in the evolution of protein complexes to conserve assembly pathways. Finally, we observe that evolutionarily more recent subunits tend to be more flexible than more ancient subunits within a complex. Together, these results reveal the influence of both evolution and intrinsic structural dynamics on protein assembly and demonstrate, for the first time on a genome-wide scale, the biological importance of ordered assembly pathways.

Microbial communities represent the largest portion of the Earth’s biomass. Metagenomics projects use high-throughput sequencing to survey these communities and shed light on genetic capabilities that enable microbes to inhabit every corner of the biosphere. Metagenome studies are generally based on (i) classifying and ranking functions of identified genes ; and (ii) estimating the phyletic distribution of constituent microbial species. To understand microbial communities at the systems level, it is necessary to extend these studies beyond the species’ boundaries and capture higher levels of metabolic complexity. We evaluated 11 metagenome samples and demonstrated that microbes inhabiting the same ecological niche share common preferences for synonymous codons, regardless of their phylogeny. By exploring concepts of translational optimization through codon usage adaptation, we demonstrated that community-wide bias in codon usage can be used as a prediction tool for lifestyle-specific genes across the entire microbial community, effectively considering microbial communities as meta-genomes. These findings set up a ‘functional metagenomics’ platform for the identification of genes relevant for adaptations of entire microbial communities to environments. Our results provide valuable arguments in defining the concept of microbial species through the context of their interactions within the community.

The interface of protein structural biology, protein biophysics, molecular evolution, and molecular population genetics forms the foundations for a mechanistic understanding of many aspects of protein biochemistry. Current efforts in interdisciplinary protein modeling are in their infancy and the state-of-the art of such models is described. Beyond the relationship between amino acid substitution and static protein structure, protein function, and corresponding organismal fitness, other considerations are also discussed. More complex mutational processes such as insertion and deletion and domain rearrangements and even circular permutations should be evaluated. The role of intrinsically disordered proteins is still controversial, but may be increasingly important to consider. Protein geometry and protein dynamics as a deviation from static considerations of protein structure are also important. Protein expression level is known to be a major determinant of evolutionary rate and several considerations including selection at the mRNA level and the role of interaction specificity are discussed. Lastly, the relationship between modeling and needed high-throughput experimental data as well as experimental examination of protein evolution using ancestral sequence resurrection and in vitro biochemistry are presented, towards an aim of ultimately generating better models for biological inference and prediction. © 2012 The Protein Society., NSF EF; Grant number: 0905606

The role of ubiquitin in many of the known cellular processes, not just protein degradation, is based on its unique ability to bind a range of proteins that are structurally and functionally different. To understand how ubiquitin can bind to proteins with different structures, we review the extent of the conservation and variation that occur in the structures of two free ubiquitins and ubiquitins in 16 complexes that have been determined at high resolution (1.2–2 A). Around 80% of the atomic groups in these structures have positions that differ less than 1 A. This conserved core provides a rigid platform for flexible loop regions, 39 residues with side chains that can take up different conformations, and a flexible six-residue region at the C-terminus. In most cases the ability of ubiquitin to bind different structures is limited in part by a central set of residues that largely conserve their conformations. The accommodation of differences in binding proteins is enabled by changes in the flexible surface side chains, loop movements, different specific interactions, water molecules in the interface and the flexible C-terminus.

Wrnip1 (Werner helicase-interacting protein 1) has been implicated in the bypass of stalled replication forks in bakers' yeast. However, the function(s) of human Wrnip1 has remained elusive so far. Here we report that Wrnip1 is distributed inside heterogeneous structures detectable in nondamaged cells throughout the cell cycle. In an attempt to characterize these structures, we found that Wrnip1 resides in DNA replication factories. Upon treatments that stall replication forks, such as UVC light, the amount of chromatin-bound Wrnip1 and the number of foci significantly increase, further implicating Wrnip1 in DNA replication. Interestingly, the nuclear pattern of Wrnip1 appears to extend to a broader landscape, as it can be detected in promyelocytic leukemia nuclear bodies. The presence of Wrnip1 into these heterogeneous subnuclear structures requires its ubiquitin-binding zinc finger (UBZ) domain, which is able to interact with different ubiquitin (Ub) signals, including mono-Ub and chains linked via lysine 48 and 63. Moreover, the oligomerization of Wrnip1 mediated by its C terminus is also important for proper subnuclear localization. Our study is the first to reveal the composite and regulated topography of Wrnip1 in the human nucleus, highlighting its potential role in replication and other nuclear transactions.

Summary Is the order in which proteins assemble into complexes important for biological function? Here, we seek to address this by searching for evidence of evolutionary selection for ordered protein complex assembly. First, we experimentally characterize the assembly pathways of several heteromeric complexes and show that they can be simply predicted from their three-dimensional structures. Then, by mapping gene fusion events identified from fully sequenced genomes onto protein complex assembly pathways, we demonstrate evolutionary selection for conservation of assembly order. Furthermore, using structural and high-throughput interaction data, we show that fusion tends to optimize assembly by simplifying protein complex topologies. Finally, we observe protein structural constraints on the gene order of fusion that impact the potential for fusion to affect assembly. Together, these results reveal the intimate relationships among protein assembly, quaternary structure, and evolution and demonstrate on a genome-wide scale the biological importance of ordered assembly pathways., Graphical Abstract Highlights ► 3D structures predict assembly pathways of protein complexes in solution ► Assembly pathways tend to be conserved in evolution, based on gene fusion data ► Gene fusion tends to optimize assembly by simplifying protein complex topologies ► Quaternary structure influences selection for fusion events with close N/C termini, The assembly pathway of most protein complexes can be predicted based on analysis of their subunit interfaces in 3D structures. Large-scale analyses of protein sequence, structure, and interaction data suggest that gene fusion events during evolution tend to conserve and optimize these ordered assembly pathways.