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, CaCl₂, 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.

The Combined Effect of Oseltamivir and Favipiravir on Influenza A Virus Evolution

Influenza virus inflicts a heavy death toll annually and resistance to existing antiviral drugs has generated interest in the development of agents with novel mechanisms of action. Favipiravir is an antiviral drug that acts by increasing the genome-wide mutation rate of influenza A virus (IAV). Potential synergistic benefits of combining oseltamivir and favipiravir have been demonstrated in animal models of influenza, but the population-level effects of combining the drugs are unknown. In order to elucidate the underlying evolutionary processes at play, we performed genome-wide sequencing of IAV experimental populations subjected to serial passaging in vitro under a combined protocol of oseltamivir and favipiravir. We describe the interplay between mutation, selection, and genetic drift that ultimately culminates in population extinction. In particular, selective sweeps around oseltamivir resistance mutations reduce genome-wide variation while deleterious mutations hitchhike to fixation given the increased mutational load generated by favipiravir. This latter effect reduces viral fitness and accelerates extinction compared with IAV populations treated with favipiravir alone, but risks spreading both established and newly emerging mutations, including possible drug resistance mutations, if transmission occurs before the viral populations are eradicated.

Molecular Basis for Differential Patterns of Drug Resistance in Influenza N1 and N2 Neuraminidase

Neuraminidase (NA) inhibitors are used for the prevention and treatment of influenza A virus infections. Two subtypes of NA, N1 and N2, predominate in viruses that infect humans, but differential patterns of drug resistance have emerged in each subtype despite highly homologous active sites. To understand the molecular basis for the selection of these drug resistance mutations, structural and dynamic analyses on complexes of N1 and N2 NA with substrates and inhibitors were performed. Comparison of dynamic substrate and inhibitor envelopes and interactions at the active site revealed how differential patterns of drug resistance have emerged for specific drug resistance mutations, at residues I222, S246, and H274 in N1 and E119 in N2. Our results show that the differences in intermolecular interactions, especially van der Waals contacts, of the inhibitors versus substrates at the NA active site effectively explain the selection of resistance mutations in the two subtypes. Avoiding such contacts that render inhibitors vulnerable to resistance by better mimicking the dynamics and intermolecular interactions of substrates can lead to the development of novel inhibitors that avoid drug resistance in both subtypes.

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.

DOI:http://dx.doi.org/10.7554/eLife.15802.001

The ability of an organism to grow and reproduce, that is, it’s “fitness”, is determined by how its genes interact with the environment. Yeast is a model organism in which researchers can control the exact mutations present in the yeast’s genes (its genotype) and the conditions in which the yeast cells live (their environment). This allows researchers to measure how a yeast cell’s genotype and environment affect its fitness.

Ubiquitin is a protein that many organisms depend on to manage cell stress by acting as a tag that targets other proteins for degradation. Essential proteins such as ubiquitin often remain unchanged by mutation over long periods of time. As a result, these proteins evolve very slowly. Like all proteins, ubiquitin is built from a chain of amino acid molecules linked together, and the ubiquitin proteins of yeast and humans are made of almost identical sequences of amino acids.

Although ubiquitin has barely changed its sequence over evolution, previous studies have shown that – under normal growth conditions in the laboratory – most amino acids in ubiquitin can be mutated without any loss of cell fitness. This led Mavor et al. to hypothesize that treating the yeast cells with chemicals that cause cell stress might lead to amino acids in ubiquitin becoming more sensitive to mutation.

To test this idea, a class of graduate students at the University of California, San Francisco grew yeast cells with different ubiquitin mutations together, and with different chemicals that induce cell stress, and measured their growth rates. Sequencing the ubiquitin gene in the thousands of tested yeast cells revealed that three of the chemicals cause a shared set of amino acids in ubiquitin to become more sensitive to mutation.

This result suggests that these amino acids are important for the stress response, possibly by altering the ability of yeast cells to target certain proteins for degradation. Conversely, another chemical causes yeast to become more tolerant to changes in the ubiquitin sequence. The experiments also link changes in particular amino acids in ubiquitin to specific stress responses.

Mavor et al. show that many of ubquitin’s amino acids are sensitive to mutation under different stress conditions, while others can be mutated to form different amino acids without effecting fitness. By testing the effects of other chemicals, future experiments could further characterize how the yeast’s genotype and environment interact.

DOI:http://dx.doi.org/10.7554/eLife.15802.002

Positive Selection Drives Preferred Segment Combinations during Influenza Virus Reassortment

Influenza A virus (IAV) has a segmented genome that allows for the exchange of genome segments between different strains. This reassortment accelerates evolution by breaking linkage, helping IAV cross species barriers to potentially create highly virulent strains. Challenges associated with monitoring the process of reassortment in molecular detail have limited our understanding of its evolutionary implications. We applied a novel deep sequencing approach with quantitative analysis to assess the in vitro temporal evolution of genomic reassortment in IAV. The combination of H1N1 and H3N2 strains reproducibly generated a new H1N2 strain with the hemagglutinin and nucleoprotein segments originating from H1N1 and the remaining six segments from H3N2. By deep sequencing the entire viral genome, we monitored the evolution of reassortment, quantifying the relative abundance of all IAV genome segments from the two parent strains over time and measuring the selection coefficients of the reassorting segments. Additionally, we observed several mutations coemerging with reassortment that were not found during passaging of pure parental IAV strains. Our results demonstrate how reassortment of the segmented genome can accelerate viral evolution in IAV, potentially enabled by the emergence of a small number of individual mutations.

Influenza virus drug resistance: a time-sampled population genetics perspective

The challenge of distinguishing genetic drift from selection remains a central focus of population genetics. Time-sampled data may provide a powerful tool for distinguishing these processes, and we here propose approximate Bayesian, maximum likelihood, and analytical methods for the inference of demography and selection from time course data. Utilizing these novel statistical and computational tools, we evaluate whole-genome datasets of an influenza A H1N1 strain in the presence and absence of oseltamivir (an inhibitor of neuraminidase) collected at thirteen time points. Results reveal a striking consistency amongst the three estimation procedures developed, showing strongly increased selection pressure in the presence of drug treatment. Importantly, these approaches re-identify the known oseltamivir resistance site, successfully validating the approaches used. Enticingly, a number of previously unknown variants have also been identified as being positively selected. Results are interpreted in the light of Fisher's Geometric Model, allowing for a quantification of the increased distance to optimum exerted by the presence of drug, and theoretical predictions regarding the distribution of beneficial fitness effects of contending mutations are empirically tested. Further, given the fit to expectations of the Geometric Model, results suggest the ability to predict certain aspects of viral evolution in response to changing host environments and novel selective pressures.

Resistance to vemurafenib resulting from a novel mutation in the BRAFV600E kinase domain

Resistance to the BRAF inhibitor vemurafenib poses a significant problem for the treatment of BRAFV600E-positive melanomas. It is therefore critical to prospectively identify all vemurafenib resistance mechanisms prior to their emergence in the clinic. The vemurafenib resistance mechanisms described to date do not result from secondary mutations within BRAFV600E. To search for possible mutations within BRAFV600E that can confer drug resistance, we developed a systematic experimental approach involving targeted saturation mutagenesis, selection of drug-resistant variants, and deep sequencing. We identified a single nucleotide substitution (T1514A, encoding L505H) that greatly increased drug resistance in cultured cells and mouse xenografts. The kinase activity of BRAFV600E/L505H was higher than that of BRAFV600E, resulting in cross-resistance to a MEK inhibitor. However, BRAFV600E/L505H was less resistant to several other BRAF inhibitors whose binding sites were further from L505 than that of PLX4720. Our results identify a novel vemurafenib-resistant mutant and provide insights into the treatment for melanomas bearing this mutation.

Enforced N-domain proximity stimulates Hsp90 ATPase activity and is compatible with function in vivo

Hsp90 populates distinct open and closed conformations mediated by transient N-terminal dimerization. To investigate the mechanistic role of these large conformational changes, we designed Hsp90 with an N-terminal coiled-coil to clamp the termini together and enforce N-domain proximity. Biophysical analyses demonstrate that the coiled-coil effectively maintains N-domain proximity in the absence of ATP, a condition that favors the open state of Hsp90. Enforcing N-domain proximity results in increased ATPase activity, indicating that N-terminal dimerization is a rate-limiting step that is sped-up with the coiled-coil due to increased effective N-domain concentration. The relative difference in ATPase activity between coil-Hsp90 and wt was reduced in the presence of both an ATPase activating (Aha1) and an inhibiting (Sba1) co-chaperone. As both of these co-chaperones bind preferentially to N-terminally dimerized Hsp90, the buffering effect of these co-chaperones demonstrates the biochemical relevance of Hsp90 conformational properties in addition to N-terminal dimerization. Enforcing N-domain proximity is compatible with viability in yeast, underlining the mechanistic relevance of Hsp90 conformational changes that are less dramatic than the transition between fully open and closed.

Charge-rich regions modulate the anti-aggregation activity of Hsp90

Protein aggregation can have dramatic effects on cellular function and plays a causative role in many human diseases. In all cells, molecular chaperones bind to aggregation-prone proteins and hinder aggregation. The ability of a protein to resist aggregation and remain soluble in aqueous solution is linked to the physical properties of the protein. Numerous physical studies demonstrate that charged atoms favor solubility. We note that many molecular chaperones possess a substantial negative charge that may allow them to impart solubility on aggregation-prone proteins. Hsp90 is one such negatively charged molecular chaperone. The charge on Hsp90 is largely concentrated in two highly acidic regions. To investigate the relationship between chaperone charge and protein solubility, we deleted these charge-rich regions and analyzed the resulting Hsp90 constructs for anti-aggregation activity. We found that deletion of both charge-rich regions dramatically impaired Hsp90 anti-aggregation activity. The anti-aggregation role of the deleted charge-rich regions could be due to net charge or sequence-specific features. To distinguish these possibilities, we attached an acid-rich region with a distinct amino acid sequence to our double-deleted Hsp90 construct. This charge rescue construct displayed effective anti-aggregation activity indicating that the net charge of Hsp90 contributes to its anti-aggregation activity.

Modular control of cross-oligomerization: analysis of superstabilized Hsp90 homodimers in vivo

Homo-oligomeric proteins fulfill numerous functions in all cells. The ability to co-express subunits of these proteins that preferentially self-assemble without cross-oligomerizing provides for controlled experiments to analyze the function of mutant homo-oligomers in vivo. Hsp90 is a dimeric chaperone involved in the maturation of many kinases and steroid hormone receptors. We observed that co-expression of different Hsp90 subunits in Saccharomyces cerevisiae caused unpredictable synthetic growth defects due to cross-dimerization. We engineered superstabilized Hsp90 dimers that resisted cross-dimerization with endogenous Hsp90 and alleviated the synthetic growth defect. Superstabilized Hsp90 dimers supported robust growth of S. cerevisiae, indicating that dissociation of Hsp90 dimers could be hindered without compromising essential function. We utilized superstabilized dimers to analyze the activity of ATPase mutant homodimers in a temperature-sensitive yeast background where elevated temperature inactivated all other Hsp90 species. We found that ATP binding and hydrolysis by Hsp90 are both required for the efficient maturation of glucocorticoid receptor and v-Src, confirming the critical role of ATP hydrolysis in the maturation of steroid hormone receptors and kinases in vivo.

Dimerization of Hsp90 is required for in vivo function. Design and analysis of monomers and dimers

Heat shock protein 90 (Hsp90) plays a central role in signal transduction and has emerged as a promising target for anti-cancer therapeutics, but its molecular mechanism is poorly understood. At physiological concentration, Hsp90 predominantly forms dimers, but the function of full-length monomers in cells is not clear. Hsp90 contains three domains: the N-terminal and middle domains contribute directly to ATP binding and hydrolysis and the C domain mediates dimerization. To study the function of Hsp90 monomers, we used a single-chain strategy that duplicated the C-terminal dimerization domain. This novel monomerization strategy had the dual effect of stabilizing the C domain to denaturation and hindering intermolecular association of the ATPase domain. The resulting construct was predominantly monomeric at physiological concentration and did not function to support yeast viability as the sole Hsp90. The monomeric construct was also defective at ATP hydrolysis and the activation of a kinase and steroid receptor substrate in yeast cells. The ability to support yeast growth was rescued by the addition of a coiled-coil dimerization domain, indicating that the parental single-chain construct is functionally defective because it is monomeric.

Altered tethering of the SspB adaptor to the ClpXP protease causes changes in substrate delivery

SspB is a dimeric adaptor protein that increases the rate at which ssrA-tagged substrates are degraded by tethering them to the ClpXP protease. Each SspB subunit consists of a folded domain that forms the dimer interface and a flexible C-terminal tail. Ternary delivery complexes are stabilized by three sets of tethering interactions. The C-terminal XB peptide of each SspB subunit binds ClpX, the body of SspB binds one part of the ssrA-tag sequence, and ClpX binds another part of the tag. To test the functional importance of these tethering interactions, we engineered monomeric SspB variants and dimeric variants with different length linkers between the SspB body and the XB peptide and employed substrates with degradation tags that bind ClpX weakly and/or contain extensions between the binding sites for SspB and ClpX. We find that monomeric SspB variants can enhance ClpXP degradation of a subset of substrates, that doubling the number of tethering interactions stimulates degradation via changes in Km and Vmax, and that major alterations in the length of the 48-residue SspB linker cause only small changes in the efficiency of substrate delivery. These results indicate that the properties of the degradation tag and the number of SspB.ClpX tethering interactions are the major factors that determine the extent to which the substrate and ClpX are engaged in ternary delivery complexes.

Specificity versus stability in computational protein design

Protein-protein interactions can be designed computationally by using positive strategies that maximize the stability of the desired structure and/or by negative strategies that seek to destabilize competing states. Here, we compare the efficacy of these methods in reengineering a protein homodimer into a heterodimer. The stability-design protein (positive design only) was experimentally more stable than the specificity-design heterodimer (positive and negative design). By contrast, only the specificity-design protein assembled as a homogenous heterodimer in solution, whereas the stability-design protein formed a mixture of homodimer and heterodimer species. The experimental stabilities of the engineered proteins correlated roughly with their calculated stabilities, and the crystal structure of the specificity-design heterodimer showed most of the predicted side-chain packing interactions and a main-chain conformation indistinguishable from the wild-type structure. These results indicate that the design simulations capture important features of both stability and structure and demonstrate that negative design can be critical for attaining specificity when competing states are close in structure space.

Asymmetric interactions of ATP with the AAA+ ClpX6 unfoldase: allosteric control of a protein machine

ATP hydrolysis by AAA+ ClpX hexamers powers protein unfolding and translocation during ClpXP degradation. Although ClpX is a homohexamer, positive and negative allosteric interactions partition six potential nucleotide binding sites into three classes with asymmetric properties. Some sites release ATP rapidly, others release ATP slowly, and at least two sites remain nucleotide free. Recognition of the degradation tag of protein substrates requires ATP binding to one set of sites and ATP or ADP binding to a second set of sites, suggesting a mechanism that allows repeated unfolding attempts without substrate release over multiple ATPase cycles. Our results rule out concerted hydrolysis models involving ClpX(6)*ATP(6) or ClpX(6)*ADP(6) and highlight structures of hexameric AAA+ machines with three or four nucleotides as likely functional states. These studies further emphasize commonalities between distant AAA+ family members, including protein and DNA translocases, helicases, motor proteins, clamp loaders, and other ATP-dependent enzymes.

Nucleotide-dependent substrate handoff from the SspB adaptor to the AAA+ ClpXP protease

The SspB adaptor enhances ClpXP degradation by binding the ssrA degradation tag of substrates and the AAA+ ClpX unfoldase. To probe the mechanism of substrate delivery, we engineered a disulfide bond between the ssrA tag and SspB and demonstrated otherwise normal interactions by solving the crystal structure. Although the covalent link prevents adaptor.substrate dissociation, ClpXP degraded GFP-ssrA that was disulfide bonded to the adaptor. Thus, crosslinked substrate must be handed directly from SspB to ClpX. The ssrA tag in the covalent adaptor complex interacted with ClpX.ATPgammaS but not ClpX.ADP, suggesting that handoff occurs in the ATP bound enzyme. By contrast, SspB alone bound ClpX in both nucleotide states. Similar handoff mechanisms will undoubtedly be used by many AAA+ adaptors and enzymes, allowing assembly of delivery complexes in either nucleotide state, engagement of the recognition tag in the ATP state, and application of an unfolding force to the attached protein following hydrolysis.

Sculpting the proteome with AAA(+) proteases and disassembly machines

Machines of protein destruction-including energy-dependent proteases and disassembly chaperones of the AAA(+) ATPase family-function in all kingdoms of life to sculpt the cellular proteome, ensuring that unnecessary and dangerous proteins are eliminated and biological responses to environmental change are rapidly and properly regulated. Exciting progress has been made in understanding how AAA(+) machines recognize specific proteins as targets and then carry out ATP-dependent dismantling of the tertiary and/or quaternary structure of these molecules during the processes of protein degradation and the disassembly of macromolecular complexes.

Bivalent Tethering of SspB to ClpXP Is Required for Efficient Substrate Delivery

SspB homodimers deliver ssrA-tagged substrates to ClpXP for degradation. SspB consists of a substrate binding domain and an unstructured tail with a ClpX binding module (XB). Using computational design, we engineered an SspB heterodimer whose subunits did not form homodimers. Experiments with the designed molecule and variants lacking one or two tails demonstrate that both XB modules are required for strong binding and efficient substrate delivery to ClpXP. Assembly of stable SspB-substrate-ClpX delivery complexes requires the coupling of weak tethering interactions between ClpX and the SspB XB modules as well as interactions between ClpX and the substrate degradation tag. The ClpX hexamer contains three XB binding sites, one per N domain dimer, and thus binds strongly to just one SspB dimer at a time. Because different adaptor proteins use the same tethering sites in ClpX, those which employ bivalent tethering, like SspB, will compete more effectively for substrate delivery to ClpXP.

Flexible linkers leash the substrate binding domain of SspB to a peptide module that stabilizes delivery complexes with the AAA+ ClpXP protease

SspB dimers bind proteins bearing the ssrA-degradation tag and stimulate their degradation by the ClpXP protease. Here, E. coli SspB is shown to contain a dimeric substrate binding domain of 110-120 N-terminal residues, which binds ssrA-tagged substrates but does not stimulate their degradation. The C-terminal 40-50 residues of SspB are unstructured but are required for SspB to form substrate-delivery complexes with ClpXP. A synthetic peptide containing the 10 C-terminal residues of SspB binds ClpX, stimulates its ATPase activity, and prevents SspB-mediated delivery of GFP-ssrA for ClpXP degradation. This tripartite structure--an ssrA-tag binding and dimerization domain, a flexible linker, and a short peptide module that docks with ClpX--allows SspB to deliver tagged substrates to ClpXP without interfering with their denaturation or degradation.

Prudent modeling of core polar residues in computational protein design

Hydrogen bond interactions were surveyed in a set of protein structures. Compared to surface positions, polar side-chains at core positions form a greater number of intra-molecular hydrogen bonds. Furthermore, the majority of polar side-chains at core positions form at least one hydrogen bond to main-chain atoms that are not involved in hydrogen bonds to other main-chain atoms. Based on this structural survey, hydrogen bond rules were generated for each polar amino acid for use in protein core design. In the context of protein core design, these prudent polar rules were used to eliminate from consideration polar amino acid rotamers that do not form a minimum number of hydrogen bonds. As an initial test, the core of Escherichia coli thioredoxin was selected as a design target. For this target, the prudent polar strategy resulted in a minor increase in computational complexity compared to a strategy that did not allow polar residues. Dead-end elimination was used to identify global minimum energy conformations for the prudent polar and no polar strategies. The prudent polar strategy identified a protein sequence that was thermodynamically stabilized by 2.5 kcal/mol relative to wild-type thioredoxin and 2.2 kcal/mol relative to a thioredoxin variant whose core was designed without polar residues.

De novo design of biocatalysts

The challenging field of de novo enzyme design is beginning to produce exciting results. The application of powerful computational methods to functional protein design has recently succeeded at engineering target activities. In addition, efforts in directed evolution continue to expand the transformations that can be accomplished by existing enzymes. The engineering of completely novel catalytic activity requires traversing inactive sequence space in a fitness landscape, a feat that is better suited to computational design. Optimizing activity, which can include subtle alterations in backbone conformation and protein motion, is better suited to directed evolution, which is highly effective at scaling fitness landscapes towards maxima. Improved rational design efforts coupled with directed evolution should dramatically improve the scope of de novo enzyme design.

Enzyme-like proteins by computational design

We report the development and initial experimental validation of a computational design procedure aimed at generating enzyme-like protein catalysts called "protozymes." Our design approach utilizes a "compute and build" strategy that is based on the physical/chemical principles governing protein stability and catalytic mechanism. By using the catalytically inert 108-residue Escherichia coli thioredoxin as a scaffold, the histidine-mediated nucleophilic hydrolysis of p-nitrophenyl acetate as a model reaction, and the ORBIT protein design software to compute sequences, an active site scan identified two promising catalytic positions and surrounding active-site mutations required for substrate binding. Experimentally, both candidate protozymes demonstrated catalytic activity significantly above background. One of the proteins, PZD2, displayed "burst" phase kinetics at high substrate concentrations, consistent with the formation of a stable enzyme intermediate. The kinetic parameters of PZD2 are comparable to early catalytic Abs. But, unlike catalytic Ab design, our design procedure is independent of fold, suggesting a possible mechanism for examining the relationships between protein fold and the evolvability of protein function.

Polar residues in the protein core of Escherichia coli thioredoxin are important for fold specificity

Most globular proteins contain a core of hydrophobic residues that are inaccessible to solvent in the folded state. In general, polar residues in the core are thermodynamically unfavorable except when they are able to form intramolecular hydrogen bonds. Compared to hydrophobic interactions, polar interactions are more directional in character and may aid in fold specificity. In a survey of 263 globular protein structures, we found a strong positive correlation between the number of polar residues at core positions and protein size. To probe the importance of buried polar residues, we experimentally tested the effects of hydrophobic mutations at the five polar core residues in Escherichia coli thioredoxin. Proteins with single hydrophobic mutations (D26I, C32A, C35A, T66L, and T77V) all have cooperative unfolding transitions like the wild type (wt), as determined by chemical denaturation. Relative to wt, D26I is more stable while the other point mutants are less stable. The combined 5-fold mutant protein (IAALV) is less stable than wt and has an unfolding transition that is substantially less cooperative than that of wt. NMR spectra as well as amide deuterium exchange indicate that IAALV is likely sampling a number of low-energy structures in the folded state, suggesting that polar residues in the core are important for specifying a well-folded native structure.