Conserved degronome features governing quality control associated proteolysis.

The eukaryotic proteome undergoes constant surveillance by quality control systems that either sequester, refold, or eliminate aberrant proteins by ubiquitin-dependent mechanisms. Ubiquitin-conjugation necessitates the recognition of degradation determinants, termed degrons, by their cognate E3 ubiquitin-protein ligases. To learn about the distinctive properties of quality control degrons, we performed an unbiased peptidome stability screen in yeast. The search identify a large cohort of proteome-derived degrons, some of which exhibited broad E3 ligase specificity. Consequent application of a machine-learning algorithm establishes constraints governing degron potency, including the amino acid composition and secondary structure propensities. According to the set criteria, degrons with transmembrane domain-like characteristics are the most probable sequences to act as degrons. Similar quality control degrons are present in viral and human proteins, suggesting conserved degradation mechanisms. Altogether, the emerging data indicate that transmembrane domain-like degron features have been preserved in evolution as key quality control determinants of protein half-life.

Increasing protein stability by inferring substitution effects from high-throughput experiments.

We apply a computational model, global multi-mutant analysis (GMMA), to inform on effects of most amino acid substitutions from a randomly mutated gene library. Using a high mutation frequency, the method can determine mutations that increase the stability of even very stable proteins for which conventional selection systems have reached their limit. As a demonstration of this, we screened a mutant library of a highly stable and computationally redesigned model protein using an in vivo genetic sensor for folding and assigned a stability effect to 374 of 912 possible single amino acid substitutions. Combining the top 9 substitutions increased the unfolding energy 47 to 69 kJ/mol in a single engineering step. Crystal structures of stabilized variants showed small perturbations in helices 1 and 2, which rendered them closer in structure to the redesign template. This case study illustrates the capability of the method, which is applicable to any screen for protein function.

Computational and Experimental Assessment of Backbone Templates for Computational Redesign of the Thioredoxin Fold.

Computational protein design has taken big strides in recent years; however, the tools available are still not at a state where a sequence can be designed to fold into a given protein structure at will and with high probability. We have applied here a recent release of Rosetta Design to redesign a set of structurally very similar proteins belonging to the thioredoxin fold. We used a genetic screening tool to estimate solubility/folding of the designed proteins in E. coli and to select the best hits from this for further biochemical characterization. We have previously used this set of template proteins for redesign and found that success was highly dependent on template structure, a trait which was also found in this study. Nevertheless, state-of-the-art design software is now able to predict the best template, most likely due to the introduction of an energy term that reports on stress in covalent bond lengths and angles. The template that led to the greatest fraction of successful designs was the same (a thioredoxin from spinach) as that identified in our previous study. Our previously described redesign of thioredoxin, which also used the spinach protein as a template, however also performed well as a template. In the present study, both of these templates yielded proteins with compact folded structures and enforced the conclusion that any design project must carefully consider different design templates. Fortunately, selecting designs based on energies appears to correctly identify such templates.

Charge Interactions in a Highly Charge-Depleted Protein.

Electrostatic forces are important for protein folding and are favored targets of protein engineering. However, interactions between charged residues are difficult to study because of the complex network of interactions found in most proteins. We have designed a purposely simple system to investigate this problem by systematically introducing individual and pairs of charged and titratable residues in a protein otherwise free of such residues. We used constant pH molecular dynamics simulations, NMR spectroscopy, and thermodynamic double mutant cycles to probe the structure and energetics of the interaction between the charged residues. We found that the partial burial of surface charges contributes to a shift in pKa value, causing an aspartate to titrate in the neutral pH range. Additionally, the interaction between pairs of residues was found to be highly context dependent, with some pairs having no apparent preferential interaction, while other pairs would engage in coupled titration forming a highly stabilized salt bridge. We find good agreement between experiments and simulations and use the simulations to rationalize our observations and to provide a detailed mechanistic understanding of the electrostatic interactions.

Global analysis of protein stability by temperature and chemical denaturation.

The stability of a protein is a fundamental property that determines under which conditions, the protein is functional. Equilibrium unfolding with denaturants requires preparation of several samples and only provides the free energy of folding when performed at a single temperature. The typical sample requirement is around 0.5-1 mg of protein. If the stability of many proteins or protein variants needs to be determined, substantial protein production may be needed. Here we have determined the stability of acyl-coenzyme A binding protein at pH 5.3 and chymotrypsin inhibitor 2 at pH 3 and pH 6.25 by combined temperature and denaturant unfolding. We used a setup where tryptophan fluorescence is measured in quartz capillaries where only 10 µl is needed. Temperature unfolding of a series of 15 samples at increasing denaturant concentrations provided accurate and precise thermodynamic parameters. We find that the number of samples may be further reduced and less than 10 µg of protein in total are needed for reliable stability measurements. For assessment of stability of protein purified in small scale e.g. in micro plate format, our method will be highly applicable. The routine for fitting the experimental data is made available as a python notebook.

A Soluble, Folded Protein without Charged Amino Acid Residues.

Charges are considered an integral part of protein structure and function, enhancing solubility and providing specificity in molecular interactions. We wished to investigate whether charged amino acids are indeed required for protein biogenesis and whether a protein completely free of titratable side chains can maintain solubility, stability, and function. As a model, we used a cellulose-binding domain from Cellulomonas fimi, which, among proteins of more than 100 amino acids, presently is the least charged in the Protein Data Bank, with a total of only four titratable residues. We find that the protein shows a surprising resilience toward extremes of pH, demonstrating stability and function (cellulose binding) in the pH range from 2 to 11. To ask whether the four charged residues present were required for these properties of this protein, we altered them to nontitratable ones. Remarkably, this chargeless protein is produced reasonably well in Escherichia coli, retains its stable three-dimensional structure, and is still capable of strong cellulose binding. To further deprive this protein of charges, we removed the N-terminal charge by acetylation and studied the protein at pH 2, where the C-terminus is effectively protonated. Under these conditions, the protein retains its function and proved to be both soluble and have a reversible folding-unfolding transition. To the best of our knowledge, this is the first time a soluble, functional protein with no titratable side chains has been produced.

A simple probabilistic model of multibody interactions in proteins.

Protein structure prediction methods typically use statistical potentials, which rely on statistics derived from a database of know protein structures. In the vast majority of cases, these potentials involve pairwise distances or contacts between amino acids or atoms. Although some potentials beyond pairwise interactions have been described, the formulation of a general multibody potential is seen as intractable due to the perceived limited amount of data. In this article, we show that it is possible to formulate a probabilistic model of higher order interactions in proteins, without arbitrarily limiting the number of contacts. The success of this approach is based on replacing a naive table-based approach with a simple hierarchical model involving suitable probability distributions and conditional independence assumptions. The model captures the joint probability distribution of an amino acid and its neighbors, local structure and solvent exposure. We show that this model can be used to approximate the conditional probability distribution of an amino acid sequence given a structure using a pseudo-likelihood approach. We verify the model by decoy recognition and site-specific amino acid predictions. Our coarse-grained model is compared to state-of-art methods that use full atomic detail. This article illustrates how the use of simple probabilistic models can lead to new opportunities in the treatment of nonlocal interactions in knowledge-based protein structure prediction and design.