Systematic profiling of dominant ubiquitin variants reveals key functional nodes contributing to evolutionary selection.

Dominant-negative mutations can help to investigate the biological mechanisms and to understand the selective pressures for multifunctional proteins. However, most studies have focused on recessive mutant effects that occur in the absence of a second functional gene copy, which overlooks the fact that most eukaryotic genomes contain more than one copy of many genes. We have identified dominant effects on yeast growth rate among all possible point mutations in ubiquitin expressed alongside a wild-type allele. Our results reveal more than 400 dominant-negative mutations, indicating that dominant-negative effects make a sizable contribution to selection acting on ubiquitin. Cellular and biochemical analyses of individual ubiquitin variants show that dominant-negative effects are explained by varied accumulation of polyubiquitinated cellular proteins and/or defects in conjugation of ubiquitin variants to ubiquitin ligases. Our approach to identify dominant-negative mutations is general and can be applied to other proteins of interest.

A complete allosteric map of a GTPase switch in its native cellular network.

Allosteric regulation is central to protein function in cellular networks. A fundamental open question is whether cellular regulation of allosteric proteins occurs only at a few defined positions or at many sites distributed throughout the structure. Here, we probe the regulation of GTPases-protein switches that control signaling through regulated conformational cycling-at residue-level resolution by deep mutagenesis in the native biological network. For the GTPase Gsp1/Ran, we find that 28% of the 4,315 assayed mutations show pronounced gain-of-function responses. Twenty of the sixty positions enriched for gain-of-function mutations are outside the canonical GTPase active site switch regions. Kinetic analysis shows that these distal sites are allosterically coupled to the active site. We conclude that the GTPase switch mechanism is broadly sensitive to cellular allosteric regulation. Our systematic discovery of new regulatory sites provides a functional map to interrogate and target GTPases controlling many essential biological processes.

The Adaptive Potential of the Middle Domain of Yeast Hsp90.

The distribution of fitness effects (DFEs) of new mutations across different environments quantifies the potential for adaptation in a given environment and its cost in others. So far, results regarding the cost of adaptation across environments have been mixed, and most studies have sampled random mutations across different genes. Here, we quantify systematically how costs of adaptation vary along a large stretch of protein sequence by studying the distribution of fitness effects of the same ≈2,300 amino-acid changing mutations obtained from deep mutational scanning of 119 amino acids in the middle domain of the heat shock protein Hsp90 in five environments. This region is known to be important for client binding, stabilization of the Hsp90 dimer, stabilization of the N-terminal-Middle and Middle-C-terminal interdomains, and regulation of ATPase-chaperone activity. Interestingly, we find that fitness correlates well across diverse stressful environments, with the exception of one environment, diamide. Consistent with this result, we find little cost of adaptation; on average only one in seven beneficial mutations is deleterious in another environment. We identify a hotspot of beneficial mutations in a region of the protein that is located within an allosteric center. The identified protein regions that are enriched in beneficial, deleterious, and costly mutations coincide with residues that are involved in the stabilization of Hsp90 interdomains and stabilization of client-binding interfaces, or residues that are involved in ATPase-chaperone activity of Hsp90. Thus, our study yields information regarding the role and adaptive potential of a protein sequence that complements and extends known structural information.

Analyses of HIV proteases variants at the threshold of viability reveals relationships between processing efficiency and fitness.

Investigating the relationships between protein function and fitness provides keys for understanding biochemical mechanisms that underly evolution. Mutations with partial fitness defects can delineate the threshold of biochemical function required for viability. We utilized a previous deep mutational scan of HIV-1 protease (PR) to identify variants with 15-45 per cent defects in replication and analysed the biochemical function of eight variants (L10M, L10S, V32C, V32I, A71V, A71S, Q92I, Q92N). We purified each variant and assessed the efficiency of peptide cleavage for three cut sites (MA-CA, TF-PR, and PR-RT) as well as gel-based analyses of processing of purified Gag. The cutting activity of at least one site was perturbed relative to WT protease for all variants, consistent with cutting activity being a primary determinant of fitness effects. We examined the correlation of fitness defects with cutting activity of different sites. MA-CA showed the weakest correlation (R 2 = 0.02) with fitness, suggesting relatively weak coupling with viral replication. In contrast, cutting of the TF-PR site showed the strongest correlation with fitness (R 2 = 0.53). Cutting at the TF-PR site creates a new PR protein with a free N-terminus that is critical for activity. Our findings indicate that increasing the pool of active PR is rate limiting for viral replication, making this an ideal step to target with inhibitors.

Genetic interaction mapping informs integrative structure determination of protein complexes.

Determining structures of protein complexes is crucial for understanding cellular functions. Here, we describe an integrative structure determination approach that relies on in vivo measurements of genetic interactions. We construct phenotypic profiles for point mutations crossed against gene deletions or exposed to environmental perturbations, followed by converting similarities between two profiles into an upper bound on the distance between the mutated residues. We determine the structure of the yeast histone H3-H4 complex based on ~500,000 genetic interactions of 350 mutants. We then apply the method to subunits Rpb1-Rpb2 of yeast RNA polymerase II and subunits RpoB-RpoC of bacterial RNA polymerase. The accuracy is comparable to that based on chemical cross-links; using restraints from both genetic interactions and cross-links further improves model accuracy and precision. The approach provides an efficient means to augment integrative structure determination with in vivo observations.

Comprehensive fitness maps of Hsp90 show widespread environmental dependence.

Gene-environment interactions have long been theorized to influence molecular evolution. However, the environmental dependence of most mutations remains unknown. Using deep mutational scanning, we engineered yeast with all 44,604 single codon changes encoding 14,160 amino acid variants in Hsp90 and quantified growth effects under standard conditions and under five stress conditions. To our knowledge, these are the largest determined comprehensive fitness maps of point mutants. The growth of many variants differed between conditions, indicating that environment can have a large impact on Hsp90 evolution. Multiple variants provided growth advantages under individual conditions; however, these variants tended to exhibit growth defects in other environments. The diversity of Hsp90 sequences observed in extant eukaryotes preferentially contains variants that supported robust growth under all tested conditions. Rather than favoring substitutions in individual conditions, the long-term selective pressure on Hsp90 may have been that of fluctuating environments, leading to robustness under a variety of conditions.

Extending chemical perturbations of the ubiquitin fitness landscape in a classroom setting reveals new constraints on sequence tolerance.

Although the primary protein sequence of ubiquitin (Ub) is extremely stable over evolutionary time, it is highly tolerant to mutation during selection experiments performed in the laboratory. We have proposed that this discrepancy results from the difference between fitness under laboratory culture conditions and the selective pressures in changing environments over evolutionary timescales. Building on our previous work (Mavor et al., 2016), we used deep mutational scanning to determine how twelve new chemicals (3-Amino-1,2,4-triazole, 5-fluorocytosine, Amphotericin B, CaCl2, Cerulenin, Cobalt Acetate, Menadione, Nickel Chloride, p-Fluorophenylalanine, Rapamycin, Tamoxifen, and Tunicamycin) reveal novel mutational sensitivities of ubiquitin residues. Collectively, our experiments have identified eight new sensitizing conditions for Lys63 and uncovered a sensitizing condition for every position in Ub except Ser57 and Gln62. By determining the ubiquitin fitness landscape under different chemical constraints, our work helps to resolve the inconsistencies between deep mutational scanning experiments and sequence conservation over evolutionary timescales.

Rescue of conformational dynamics in enzyme catalysis by directed evolution.

Rational design and directed evolution have proved to be successful approaches to increase catalytic efficiencies of both natural and artificial enzymes. Protein dynamics is recognized as important, but due to the inherent flexibility of biological macromolecules it is often difficult to distinguish which conformational changes are directly related to function. Here, we use directed evolution on an impaired mutant of the proline isomerase CypA and identify two second-shell mutations that partially restore its catalytic activity. We show both kinetically, using NMR spectroscopy, and structurally, by room-temperature X-ray crystallography, how local perturbations propagate through a large allosteric network to facilitate conformational dynamics. The increased catalysis selected for in the evolutionary screen is correlated with an accelerated interconversion between the two catalytically essential conformational sub-states, which are both captured in the high-resolution X-ray ensembles. Our data provide a glimpse of an evolutionary trajectory and show how subtle changes can fine-tune enzyme function.

Determination of ubiquitin fitness landscapes under different chemical stresses in a classroom setting.

Ubiquitin is essential for eukaryotic life and varies in only 3 amino acid positions between yeast and humans. However, recent deep sequencing studies indicate that ubiquitin is highly tolerant to single mutations. We hypothesized that this tolerance would be reduced by chemically induced physiologic perturbations. To test this hypothesis, a class of first year UCSF graduate students employed deep mutational scanning to determine the fitness landscape of all possible single residue mutations in the presence of five different small molecule perturbations. These perturbations uncover 'shared sensitized positions' localized to areas around the hydrophobic patch and the C-terminus. In addition, we identified perturbation specific effects such as a sensitization of His68 in HU and a tolerance to mutation at Lys63 in DTT. Our data show how chemical stresses can reduce buffering effects in the ubiquitin proteasome system. Finally, this study demonstrates the potential of lab-based interdisciplinary graduate curriculum.

Thermodynamics of Rev-RNA interactions in HIV-1 Rev-RRE assembly.

The HIV-1 protein Rev facilitates the nuclear export of intron-containing viral mRNAs by recognizing a structured RNA site, the Rev-response-element (RRE), contained in an intron. Rev assembles as a homo-oligomer on the RRE using its α-helical arginine-rich-motif (ARM) for RNA recognition. One unique feature of this assembly is the repeated use of the ARM from individual Rev subunits to contact distinct parts of the RRE in different binding modes. How the individual interactions differ and how they contribute toward forming a functional complex is poorly understood. Here we examine the thermodynamics of Rev-ARM peptide binding to two sites, RRE stem IIB, the high-affinity site that nucleates Rev assembly, and stem IA, a potential intermediate site during assembly, using NMR spectroscopy and isothermal titration calorimetry (ITC). NMR data indicate that the Rev-IIB complex forms a stable interface, whereas the Rev-IA interface is highly dynamic. ITC studies show that both interactions are enthalpy-driven, with binding to IIB being 20-30 fold tighter than to IA. Salt-dependent decreases in affinity were similar at both sites and predominantly enthalpic in nature, reflecting the roles of electrostatic interactions with arginines. However, the two interactions display strikingly different partitioning between enthalpy and entropy components, correlating well with the NMR observations. Our results illustrate how the variation in binding modes to different RRE target sites may influence the stability or order of Rev-RRE assembly and disassembly, and consequently its function.

RNA-directed remodeling of the HIV-1 protein Rev orchestrates assembly of the Rev-Rev response element complex.

The HIV-1 protein Rev controls a critical step in viral replication by mediating the nuclear export of unspliced and singly-spliced viral mRNAs. Multiple Rev subunits assemble on the Rev Response Element (RRE), a structured region present in these RNAs, and direct their export through the Crm1 pathway. Rev-RRE assembly occurs via several Rev oligomerization and RNA-binding steps, but how these steps are coordinated to form an export-competent complex is unclear. Here, we report the first crystal structure of a Rev dimer-RRE complex, revealing a dramatic rearrangement of the Rev-dimer upon RRE binding through re-packing of its hydrophobic protein-protein interface. Rev-RNA recognition relies on sequence-specific contacts at the well-characterized IIB site and local RNA architecture at the second site. The structure supports a model in which the RRE utilizes the inherent plasticity of Rev subunit interfaces to guide the formation of a functional complex.

Ler of pathogenic Escherichia coli forms toroidal protein-DNA complexes.

Enteropathogenic and enterohaemorrhagic Escherichia coli are related pathotypes of bacteria that cause acute watery diarrhoea and haemorrhagic colitis, respectively, and enterohaemorrhagic E. coli can lead to a serious complication known as haemolytic uraemic syndrome. In both bacteria the global regulatory protein Ler controls virulence. The ler gene is found within the locus of enterocyte effacement, or LEE, encoding a type III secretion system necessary for injecting effector proteins into intestinal epithelial cells and causing net secretory diarrhoea. The nucleoid-associated protein H-NS silences, whereas Ler serves as an anti-silencer of, multiple LEE operons. Although Ler has a higher affinity for DNA than does H-NS, the precise molecular mechanism by which Ler increases LEE transcription remains to be determined. In this report we investigate the oligomerization activity of Ler. In solution, Ler forms dimers and soluble aggregates of up to 5000 kDa molecular mass, and appears to oligomerize more readily than the related protein H-NS. An insertional mutation into the Ler linker region diminished oligomerization activity. Despite being proteins of similar mass and having homologous DNA-binding domains, Ler and H-NS complexed to DNA migrated to distinct locations, as determined by an electrophoretic mobility shift assay, implying that the related proteins form different 3D shapes in the presence of DNA. Lastly, we present electron microscopy images of toroidal Ler-DNA structures that are predicted to be involved in stimulating gene expression.

Ler interdomain linker is essential for anti-silencing activity in enteropathogenic Escherichia coli.

Enteropathogenic Escherichia coli (EPEC) expresses a type III secretion system (T3SS) required for pathogenesis. Regulation of the genes encoding the T3SS is complex; two major regulators control transcription, the silencer H-NS, and the related H-NS-like protein Ler. Our laboratory is interested in understanding the molecular differences that distinguish the anti-silencer Ler from H-NS, and how Ler differentially regulates EPEC virulence genes. Here, we demonstrate that mutated Ler proteins either containing H-NS alpha-helices 1 and 2, missing from Ler, or truncated for the 11 aa C-terminal extension compared with the related H-NS protein, did not appreciably alter Ler function. In contrast, mutating the proline at position 92 of Ler, in the conserved C-terminal DNA binding motif, eliminated Ler activity. Inserting 11 H-NS-specific amino acids, 11 alanines or 6 alanines into the Ler linker severely impaired the ability of Ler to increase LEE5 transcription. To extend our analysis, we constructed six chimeric proteins containing the N terminus, linker region or C terminus of Ler in different combinations with the complementary domains of H-NS, and monitored their in vivo activities. Replacing the Ler linker domain with that of H-NS, or replacing the Ler C-terminal, DNA binding domain with that of H-NS eliminated the ability of Ler to increase transcription at the LEE5 promoter. Thus, the linker and C-terminal domains of Ler and H-NS are not functionally equivalent. Conversely, replacing the H-NS linker region with that of Ler caused increased transcription at LEE5 in a strain deleted for hns. In summary, the interdomain linker specific to Ler is necessary for anti-silencing activity in EPEC.