Search

We have developed a protocol for identifying proteins that are predisposed to bind linear epitopes on target proteins of interest. The protocol searches through the protein database for proteins (scaffolds) that are bound to peptides with sequences similar to accessible, linear epitopes on the target protein. The sequence match is considered more significant if residues calculated to be important in the scaffold-peptide interaction are present in the target epitope. The crystal structure of the scaffold-peptide complex is then used as a template for creating a model of the scaffold bound to the target epitope. This model can then be used in conjunction with sequence optimization algorithms or directed evolution methods to search for scaffold mutations that further increase affinity for the target protein. To test the applicability of this approach we targeted three disease-causing proteins: a tuberculosis virulence factor (TVF), the apical membrane antigen (AMA) from malaria, and hemagglutinin from influenza. In each case the best scoring scaffold was tested, and binders with Kds equal to 37 μM and 50 nM for TVF and AMA, respectively, were identified. A web server (http://rosettadesign.med.unc.edu/scaffold/) has been created for performing the scaffold search process with user-defined target sequences.

Accurate energy functions are critical to macromolecular modeling and design. We describe new tools for identifying inaccuracies in energy functions and guiding their improvement, and illustrate the application of these tools to the improvement of the Rosetta energy function. The feature analysis tool identifies discrepancies between structures deposited in the PDB and low-energy structures generated by Rosetta; these likely arise from inaccuracies in the energy function. The optE tool optimizes the weights on the different components of the energy function by maximizing the recapitulation of a wide range of experimental observations. We use the tools to examine three proposed modifications to the Rosetta energy function: improving the unfolded state energy model (reference energies), using bicubic spline interpolation to generate knowledge-based torisonal potentials, and incorporating the recently developed Dunbrack 2010 rotamer library (Shapovalov & Dunbrack, 2011). © 2013 Elsevier Inc.

We computationally designed a de novo protein-protein interaction between wild-type ubiquitin and a redesigned scaffold. Our strategy was to incorporate zinc at the designed interface to promote affinity and orientation specificity. A large set of monomeric scaffold surfaces were computationally engineered with three-residue zinc coordination sites, and the ubiquitin residue H68 was docked to the open coordination site to complete a tetrahedral zinc site. This single coordination bond was intended as a hotspot and polar interaction for ubiquitin binding, and surrounding residues on the scaffold were optimized primarily as hydrophobic residues using a rotamer-based sequence design protocol in Rosetta. From thousands of independent design simulations, four sequences were selected for experimental characterization. The best performing design, called Spelter, binds tightly to zinc (Kd 10 nM) and binds ubiquitin with a Kd of 20 µM in the presence of zinc and 68 µM in the absence of zinc. Mutagenesis studies and nuclear magnetic resonance chemical shift perturbation experiments indicate that Spelter interacts with H68 and the target surface on ubiquitin; however, H68 does not form a hotspot as intended. Instead, mutation of H68 to alanine results in tighter binding. Although a 3/1 zinc coordination arrangement at an interface cannot be ruled out as a means to improve affinity, our study led us to conclude that 2/2 coordination arrangements or multiple-zinc designs are more likely to promote high-affinity protein interactions.

The accurate design of new protein–protein interactions is a longstanding goal of computational protein design. However, most computationally designed interfaces fail to form experimentally. This investigation compares five previously described successful de novo interface designs with 158 failures. Both sets of proteins were designed with the molecular modeling program Rosetta. Designs were considered a success if a high-resolution crystal structure of the complex closely matched the design model and the equilibrium dissociation constant for binding was less than 10 μM. The successes and failures represent a wide variety of interface types and design goals including heterodimers, homodimers, peptide-protein interactions, one-sided designs (i.e., where only one of the proteins was mutated) and two-sided designs. The most striking feature of the successful designs is that they have fewer polar atoms at their interfaces than many of the failed designs. Designs that attempted to create extensive sets of interface-spanning hydrogen bonds resulted in no detectable binding. In contrast, polar atoms make up more than 40% of the interface area of many natural dimers, and native interfaces often contain extensive hydrogen bonding networks. These results suggest that Rosetta may not be accurately balancing hydrogen bonding and electrostatic energies against desolvation penalties and that design processes may not include sufficient sampling to identify side chains in preordered conformations that can fully satisfy the hydrogen bonding potential of the interface.

Here we show that a recent computationally designed zinc-mediated protein interface is serendipitously capable of catalyzing carboxyester and phosphoester hydrolysis. Although the original motivation was to design a de novo zinc-mediated protein-protein interaction (called MID1-zinc), we observed in the homodimer crystal structure a small cleft and open zinc coordination site. We investigated if the cleft and zinc site at the designed interface were sufficient for formation of a primitive active site that can perform hydrolysis. MID1-zinc hydrolyzes 4-nitrophenyl acetate with a rate acceleration of 10(5) and a k(cat)/K(M) of 630 M(-1) s(-1) and 4-nitrophenyl phosphate with a rate acceleration of 10(4) and a k(cat)/K(M) of 14 M(-1) s(-1). These rate accelerations by an unoptimized active site highlight the catalytic power of zinc and suggest that the clefts formed by protein-protein interactions are well-suited for creating enzyme active sites. This discovery has implications for protein evolution and engineering: from an evolutionary perspective, three-coordinated zinc at a homodimer interface cleft represents a simple evolutionary path to nascent enzymatic activity; from a protein engineering perspective, future efforts in de novo design of enzyme active sites may benefit from exploring clefts at protein interfaces for active site placement.

SummaryPhotocontrol of functional peptides is a powerful tool for spatial and temporal control of cell signaling events. We show that the genetically encoded light-sensitive LOV2 domain of Avena Sativa phototropin 1 (AsLOV2) can be used to reversibly photomodulate the affinity of peptides for their binding partners. Sequence analysis and molecular modeling were used to embed two peptides into the Jα helix of the AsLOV2 domain while maintaining AsLOV2 structure in the dark but allowing for binding to effector proteins when the Jα helix unfolds in the light. Caged versions of the ipaA and SsrA peptides, LOV-ipaA and LOV-SsrA, bind their targets with 49- and 8-fold enhanced affinity in the light, respectively. These switches can be used as general tools for light-dependent colocalization, which we demonstrate with photo-activable gene transcription in yeast.

International audience; The CAPRI (Critical Assessment of Predicted Interactions) and CASP (Critical Assessment of protein Structure Prediction) experiments have demonstrated the power of community-wide tests of methodology in assessing the current state of the art and spurring progress in the very challenging areas of protein docking and structure prediction. We sought to bring the power of community-wide experiments to bear on a very challenging protein design problem that provides a complementary but equally fundamental test of current understanding of protein-binding thermodynamics. We have generated a number of designed protein-protein interfaces with very favorable computed binding energies but which do not appear to be formed in experiments, suggesting that there may be important physical chemistry missing in the energy calculations. A total of 28 research groups took up the challenge of determining what is missing: we provided structures of 87 designed complexes and 120 naturally occurring complexes and asked participants to identify energetic contributions and/or structural features that distinguish between the two sets. The community found that electrostatics and solvation terms partially distinguish the designs from the natural complexes, largely due to the nonpolar character of the designed interactions. Beyond this polarity difference, the community found that the designed binding surfaces were, on average, structurally less embedded in the designed monomers, suggesting that backbone conformational rigidity at the designed surface is important for realization of the designed function. These results can be used to improve computational design strategies, but there is still much to be learned; for example, one designed complex, which does form in experiments, was classified by all metrics as a nonbinder.

The inhibitory switch (IS) domain of p21-activated kinase 1 (PAK1) stabilizes full-length PAK1 in an inactive conformation by binding to the PAK1 kinase domain. Competitive binding of small guanosine triphosphatases to the IS domain disrupts the autoinhibitory interactions and exposes the IS domain binding site on the surface of the kinase domain. To build an affinity reagent that selectively binds the activated state of PAK1, we used molecular modeling to reengineer the isolated IS domain so that it was soluble and stable, did not bind to guanosine triphosphatases and bound more tightly to the PAK1 kinase domain. Three design strategies were tested: in the first and second cases, extension and redesign of the N-terminus were used to expand the hydrophobic core of the domain, and in the third case, the termini were redesigned to be adjacent in space so that the domain could be stabilized by insertion into a loop in a host cyan fluorescent protein (CFP). The best-performing design, called CFP-PAcKer, was based on the third strategy and bound the kinase domain of PAK1 with an affinity of 400 nM. CFP-PAcKer binds more tightly to a full-length variant of PAK1 that is stabilized in the “open” state (Kd = 3.3 μM) than to full-length PAK1 in the “closed” state (undetectable affinity), and binding can be monitored with fluorescence by placing an environmentally sensitive fluorescence dye on CFP-PAcKer adjacent to the binding site.

De novo protein design requires the identification of amino-acid sequences that favor the target-folded conformation and are soluble in water. One strategy for promoting solubility is to disallow hydrophobic residues on the protein surface during design. However, naturally occurring proteins often have hydrophobic amino acids on their surface that contribute to protein stability via the partial burial of hydrophobic surface area or play a key role in the formation of protein-protein interactions. A less restrictive approach for surface design that is used by the modeling program Rosetta is to parameterize the energy function so that the number of hydrophobic amino acids designed on the protein surface is similar to what is observed in naturally occurring monomeric proteins. Previous studies with Rosetta have shown that this limits surface hydrophobics to the naturally occurring frequency (∼28%), but that it does not prevent the formation of hydrophobic patches that are considerably larger than those observed in naturally occurring proteins. Here, we describe a new score term that explicitly detects and penalizes the formation of hydrophobic patches during computational protein design. With the new term, we are able to design protein surfaces that include hydrophobic amino acids at naturally occurring frequencies, but do not have large hydrophobic patches. By adjusting the strength of the new score term, the emphasis of surface redesigns can be switched between maintaining solubility and maximizing folding free energy.

B. stearothermophilus Tryptophanyl-tRNA synthetase (TrpRS) uses different conformational states to catalyze tryptophan activation. A single Mg2+ ion increases transition state stabilization by −6.5 kcal/mol for optimal catalysis. Catalytic assist occurs if, and only if, the Mg2+ interacts with the protein. We are trying to identify the metal-protein interactions that produce this catalytic effect. Physical interactions between Mg2+ and TrpRS are mediated indirectly via active-site lysines K111, K192 and K195. Mutations of these lysines showed that they all stabilize the transition state. However, their interactions with the Mg2+ significantly reduce their catalytic effects. Catalytically productive interactions between TrpRS and the Mg2+ ion must therefore arise from outside the active site. We identified a core set of residues we call the D1 switch because they move during the catalytic conformational transition. The D1 switch lies at the corner of the N-terminal s-α-s crossover distal to the active site. It is highly conserved in Rossmannoid proteins. The Rosetta Design program suggested that mutations of D1 residues could ‘‘hyperstabilize’’ the activated state observed just prior to catalysis. Multimutant thermodynamic cycles together with substitution of Mn2+ for Mg2+ and [ATP]-dependent Michaelis-Menten kinetics demonstrate significant long-range synergistic coupling between the D1 switch and the Mg2+ ion. Thus, long-range interactions to the metal likely drive catalysis indirectly, by changing an inactive Mg2+ coordination into one that can stabilize the transition state. In this way transition-state stabilization by Mg2+ occurs if, and only if, conformational changes reposition it. We suggest that other NTPase enzymes may use similar conditional activation of Mg2+ to couple catalyzed hydrolysis of their purine triphosphate substrates to conformational changes, thereby transducing chemical free energy for cellular work and signaling. Supported by NIGMS 78227, 90406.

Background Few existing protein-protein interface design methods allow for extensive backbone rearrangements during the design process. There is also a dichotomy between redesign methods, which take advantage of the native interface, and de novo methods, which produce novel binders. Methodology Here, we propose a new method for designing novel protein reagents that combines advantages of redesign and de novo methods and allows for extensive backbone motion. This method requires a bound structure of a target and one of its natural binding partners. A key interaction in this interface, the anchor, is computationally grafted out of the partner and into a surface loop on the design scaffold. The design scaffold's surface is then redesigned with backbone flexibility to create a new binding partner for the target. Careful choice of a scaffold will bring experimentally desirable characteristics into the new complex. The use of an anchor both expedites the design process and ensures that binding proceeds against a known location on the target. The use of surface loops on the scaffold allows for flexible-backbone redesign to properly search conformational space. Conclusions and Significance This protocol was implemented within the Rosetta3 software suite. To demonstrate and evaluate this protocol, we have developed a benchmarking set of structures from the PDB with loop-mediated interfaces. This protocol can recover the correct loop-mediated interface in 15 out of 16 tested structures, using only a single residue as an anchor.

The de novo design of protein-binding peptides is challenging because it requires the identification of both a sequence and a backbone conformation favorable for binding. We used a computational strategy that iterates between structure and sequence optimization to redesign the C-terminal portion of the RGS14 GoLoco motif peptide so that it adopts a new conformation when bound to Gα(i1). An X-ray crystal structure of the redesigned complex closely matches the computational model, with a backbone root-mean-square deviation of 1.1 Å.

We have recently completed a full re-architecturing of the Rosetta molecular modeling program, generalizing and expanding its existing functionality. The new architecture enables the rapid prototyping of novel protocols by providing easy to use interfaces to powerful tools for molecular modeling. The source code of this rearchitecturing has been released as Rosetta3 and is freely available for academic use. At the time of its release, it contained 470,000 lines of code. Counting currently unpublished protocols at the time of this writing, the source includes 1,285,000 lines. Its rapid growth is a testament to its ease of use. This document describes the requirements for our new architecture, justifies the design decisions, sketches out central classes, and highlights a few of the common tasks that the new software can perform.

We have recently completed a full re-architecturing of the ROSETTA molecular modeling program, generalizing and expanding its existing functionality. The new architecture enables the rapid prototyping of novel protocols by providing easy-to-use interfaces to powerful tools for molecular modeling. The source code of this rearchitecturing has been released as ROSETTA3 and is freely available for academic use. At the time of its release, it contained 470,000 lines of code. Counting currently unpublished protocols at the time of this writing, the source includes 1,285,000 lines. Its rapid growth is a testament to its ease of use. This chapter describes the requirements for our new architecture, justifies the design decisions, sketches out central classes, and highlights a few of the common tasks that the new software can perform.

During ubiquitin conjugation, the thioester bond that links ‘donor’ ubiquitin to ubiquitin-conjugating enzyme (E2) undergoes nucleophilic attack by the ε-amino group of an acceptor lysine, resulting in formation of an isopeptide bond. Models of ubiquitination have envisioned the donor ubiquitin to be a passive participant in this process. However, we show here that the I44A mutation in ubiquitin profoundly inhibits its ability to serve as a donor for ubiquitin chain initiation or elongation, but can be rescued by computationally-predicted compensatory mutations in the E2 Cdc34. The donor defect of ubiquitin-I44A can be partially suppressed either by using a low pKa amine (hydroxylamine) as the acceptor or by performing reactions at higher pH, suggesting that the discharge defect arises in part due to inefficient deprotonation of the acceptor lysine. We propose that interaction between Cdc34 and the donor ubiquitin organizes the active site to promote efficient ubiquitination of substrate.

GoLoco motif proteins bind to the inhibitory G(i) subclass of G-protein α subunits and slow the release of bound GDP; this interaction is considered critical to asymmetric cell division and neuro-epithelium and epithelial progenitor differentiation. To provide protein tools for interrogating the precise cellular role(s) of GoLoco motif/Gα(i) complexes, we have employed structure-based protein design strategies to predict gain-of-function mutations that increase GoLoco motif binding affinity. Here, we describe fluorescence polarization and isothermal titration calorimetry measurements showing three predicted Gα(i1) point mutations, E116L, Q147L, and E245L; each increases affinity for multiple GoLoco motifs. A component of this affinity enhancement results from a decreased rate of dissociation between the Gα mutants and GoLoco motifs. For Gα(i1)(Q147L), affinity enhancement was seen to be driven by favorable changes in binding enthalpy, despite reduced contributions from binding entropy. The crystal structure of Gα(i1)(Q147L) bound to the RGS14 GoLoco motif revealed disorder among three peptide residues surrounding a well defined Leu-147 side chain. Monte Carlo simulations of the peptide in this region showed a sampling of multiple backbone conformations in contrast to the wild-type complex. We conclude that mutation of Glu-147 to leucine creates a hydrophobic surface favorably buried upon GoLoco peptide binding, yet the hydrophobic Leu-147 also promotes flexibility among residues 511-513 of the RGS14 GoLoco peptide.

Metal coordination is a key structural and functional component of a large fraction of proteins. Given this dual role we considered the possibility that metal coordination may have played a templating role in the early evolution of protein folds and complexes. We describe here a rational design approach, Metal Templated Interface Redesign (MeTIR), that mimics the time course of a hypothetical evolutionary pathway for the formation of stable protein assemblies through an initial metal coordination event. Using a folded monomeric protein, cytochrome cb 562 , as a building block we show that its non-self-associating surface can be made self-associating through a minimal number of mutations that enable Zn coordination. The protein interfaces in the resulting Zn-directed, D 2 -symmetrical tetramer are subsequently redesigned, yielding unique protein architectures that self-assemble in the presence or absence of metals. Aside from its evolutionary implications, MeTIR provides a route to engineer de novo protein interfaces and metal coordination environments that can be tuned through the extensive noncovalent bonding interactions in these interfaces.

The importance of a protein-protein interaction to a signaling pathway can be established by showing that amino acid mutations that weaken the interaction disrupt signaling, and that additional mutations that rescue the interaction recover signaling. Identifying rescue mutations, often referred to as second-site suppressor mutations, controls against scenarios in which the initial deleterious mutation inactivates the protein or disrupts alternative protein-protein interactions. Here, we test a structure-based protocol for identifying second-site suppressor mutations that is based on a strategy previously described by Kortemme and Baker. The molecular modeling software Rosetta is used to scan an interface for point mutations that are predicted to weaken binding but can be rescued by mutations on the partner protein. The protocol typically identifies three types of specificity switches: knob-in-to-hole redesigns, switching hydrophobic interactions to hydrogen bond interactions, and replacing polar interactions with nonpolar interactions. Computational predictions were tested with two separate protein complexes; the G-protein Galpha(i1) bound to the RGS14 GoLoco motif, and UbcH7 bound to the ubiquitin ligase E6AP. Eight designs were experimentally tested. Swapping a buried hydrophobic residue with a polar residue dramatically weakened binding affinities. In none of these cases were we able to identify compensating mutations that returned binding to wild-type affinity, highlighting the challenges inherent in designing buried hydrogen bond networks. The strongest specificity switches were a knob-in-to-hole design (20-fold) and the replacement of a charge-charge interaction with nonpolar interactions (55-fold). In two cases, specificity was further tuned by including mutations distant from the initial design. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

Emerging evidence suggests that ubiquitination serves as a protein trafficking signal in addition to its well characterized role in promoting protein degradation. The yeast G protein α subunit Gpa1 represents a rare example of a protein that undergoes both mono- and poly-ubiquitination. Whereas mono-ubiquitinated Gpa1 is targeted to the vacuole, poly-ubiquitinated Gpa1 is directed instead to the proteasome. Here we investigate the structural requirements for mono- and poly-ubiquitination of Gpa1. We find that variants of Gpa1 engineered to be unstable are more likely to be poly-ubiquitinated and less likely to be mono-ubiquitinated. In addition, mutants that cannot be myristoylated are no longer mono-ubiquitinated but are still polyubiquitinated. Finally, we show that the ubiquitin ligase Rsp5 is necessary for Gpa1 mono-ubiquitination in vivo and that the purified enzyme is sufficient to catalyze Gpa1 mono-ubiquitination in vitro. Taken together, these data indicate that mono- and poly-ubiquitination have distinct enzyme and substrate recognition requirements; whereas poly-ubiquitination targets misfolded protein for degradation, a distinct ubiquitination apparatus targets the fully mature, fully myristoylated G protein for mono-ubiquitination and delivery to the vacuole.

The computer-based design of protein-protein interactions is a rigorous test of our understanding of molecular recognition and an attractive approach for creating novel tools for cell and molecular research. Considerable attention has been placed on redesigning the affinity and specificity of naturally occurring interactions. Several studies have shown that reducing the desolvation costs for binding while preserving shape complimentarity and hydrogen bonding is an effective strategy for improving binding affinities. In favorable cases specificity has been designed by focusing only on interactions with the target protein, while in cases with closely related off-target proteins, it has been necessary to explicitly disfavor unwanted binding partners. The rational design of protein-protein interactions from scratch is still an unsolved problem, but recent developments in flexible backbone design and energy functions hold promise for the future.

Most existing algorithms for protein design, including thosein the Rosetta molecular modeling program, precompute energies for rotamerpairs, since these energies can be examined repeatedly. Simulatedannealing algorithms, however, do not examine these energies with thesame frequency; while some are examined many times, others may not beexamined at all. This paper compares strategies for computing these energieson the fly and caching computed energy values that are likely to bereused. By avoiding the expense of computing pair energies that are notexamined by simulated annealing, we show that some caching strategiesnot only improve running time in design, but also use 90% less memory,which allows design computations to be performed on memory-limited machines.

Amino acid side chains adopt a discrete set of favorable conformations typically referred to as rotamers. The relative energies of rotamers partially determine which side chain conformations are more often observed in protein structures and accurate estimates of these energies are important for predicting protein structure and designing new proteins. Protein modelers typically calculate side chain rotamer energies by using molecular mechanics (MM) potentials or by converting rotamer probabilities from the protein database (PDB) into relative free energies. One limitation of the knowledge-based energies is that rotamer preferences observed in the PDB can reflect internal side chain energies as well as longer-range interactions with the rest of the protein. Here, we test an alternative approach for calculating rotamer energies. We use three different quantum mechanics (QM) methods (second order Moller-Plesset (MP2), density functional theory (DFT) energy calculation using the B3LYP functional, and Hartree-Fock) to calculate the energy of amino acid rotamers in a dipeptide model system, and then use these pre-calculated values in side chain placement simulations. Energies were calculated for over 36,000 different conformations of leucine, isoleucine, and valine dipeptides with backbone torsion angles from the helical and strand regions of the Ramachandran plot. In a subset of cases these energies differ significantly from those calculated with standard molecular mechanics potentials or those derived from PDB statistics. We find that in these cases the energies from the QM methods result in more accurate placement of amino acid side chains in structure prediction tests.

Despite having irregular structure, protein loops often adopt specific conformations that are critical to protein function. Most studies in de novo protein design have focused on creating proteins with regular elements of secondary structure connected by very short loops or turns. To design longer protein loops that adopt specific conformations, we have developed a protocol within the Rosetta molecular modeling program that iterates between optimizing the sequence and conformation of a loop in search of low-energy sequence–structure pairs. We have tested the procedure by designing 10-residue loops for the connection between the second and third strand in the β-sandwich protein tenascin. Three low-energy designs from 7,200 flexible backbone trajectories were selected for experimental characterization. All three designs, called LoopA, LoopB, and LoopC, adopt stable folded structures. High-resolution crystal structures of LoopA and LoopB have been solved. LoopB adopts a structure very similar to the design model (0.46 Å rmsd), and all but one of the side chains are modeled in the correct rotamers. LoopA crystallized at low pH in a structure that differs dramatically from our design model. It forms a strand-swapped dimer mediated by hydrogen bonds to protonated glutamic acids. Gel filtration indicates that the protein is not a dimer at neutral pH. These results suggest that the high-resolution design of protein loops is possible; however, they also highlight how small changes in protein energetics can dramatically perturb the low free energy structure of a protein.

The conjugation of ubiquitin to substrates requires a series of enzymatic reactions consisting of an activating enzyme (E1), conjugating enzymes (E2) and ligases (E3). Tagging the appropriate substrate with ubiquitin is achieved by specific E2-E3 and E3-substrate interactions. E6AP, a member of the HECT family of E3s, has been previously shown to bind and function with the E2s UbcH7 and UbcH8. To decipher the sequence determinants of this specificity we have developed a quantitative E2-E3 binding assay based on fluorescence polarization and used this assay to measure the affinity of wild-type and mutant E2-E6AP interactions. Alanine scanning of the E6AP-UbcH7 binding interface identified four side-chains on UbcH7 and six side-chains on E6AP that contribute more than 1 kcal/mol to the binding free energy. Two of the hot spot residues from UbcH7 (K96 and K100) are conserved in UbcH8 but vary across other E2s. To determine if these are key specificity determining residues, we attempted to induce a tighter association between the E2 UbcH5b and E6AP by mutating the corresponding positions in UbcH5b to lysine residues. Surprisingly, the mutations had little effect, but rather a mutation at UbcH7 position 4, which is not at a hot spot on the UbcH7-E6AP interface, significantly strengthened UbcH5bs affinity for E6AP. This result indicates that E2-E3 binding specificities are a function of both favorable interactions that promote binding, and unfavorable interactions that prevent binding with unwanted partners.

In its basal state, KEAP1 binds the transcription factor NRF2 (Kd = 5 nM) and promotes its degradation by ubiquitylation. Changes in the redox environment lead to modification of key cysteines within KEAP1, resulting in NRF2 protein accumulation and the transcription of genes important for restoring the cellular redox state. Using phage display and a computational loop grafting protocol, we engineered a monobody (R1) that is a potent competitive inhibitor of the KEAP1–NRF2 interaction. R1 bound to KEAP1 with a Kd of 300 pM and in human cells freed NRF2 from KEAP1 resulting in activation of the NRF2 promoter. Unlike cysteine-reactive small molecules that lack protein specificity, R1 is a genetically encoded, reversible inhibitor designed specifically for KEAP1. R1 should prove useful for studying the role of the KEAP1–NRF2 interaction in several disease states. The structure-based phage display strategy employed here is a general approach for engineering high-affinity binders that compete with naturally occurring interactions.

G-proteins cycle between an inactive GDP-bound state and an active GTP-bound state, serving as molecular switches that coordinate cellular signaling. We recently used phage display to identify a series of peptides that bind G alpha subunits in a nucleotide-dependent manner [Johnston, C. A., Willard, F. S., Jezyk, M. R., Fredericks, Z., Bodor, E. T., Jones, M. B., Blaesius, R., Watts, V. J., Harden, T. K., Sondek, J., Ramer, J. K., and Siderovski, D. P. (2005) Structure 13, 1069-1080]. Here we describe the structural features and functions of KB-1753, a peptide that binds selectively to GDP x AlF4(-)- and GTPgammaS-bound states of G alpha(i) subunits. KB-1753 blocks interaction of G alpha(transducin) with its effector, cGMP phosphodiesterase, and inhibits transducin-mediated activation of cGMP degradation. Additionally, KB-1753 interferes with RGS protein binding and resultant GAP activity. A fluorescent KB-1753 variant was found to act as a sensor for activated G alpha in vitro. The crystal structure of KB-1753 bound to G alpha(i1) x GDP x AlF4(-) reveals binding to a conserved hydrophobic groove between switch II and alpha3 helices and, along with supporting biochemical data and previous structural analyses, supports the notion that this is the site of effector interactions for G alpha(i) subunits.

The emergence of polypeptide catalysts for amino acid activation, the slowest step in protein synthesis, poses a significant puzzle associated with the origin of biology. This problem is compounded as the 20 contemporary aminoacyl-tRNA synthetases belong to two quite distinct families. We describe here the use of protein design to show experimentally that a minimal class I aminoacyl-tRNA synthetase active site might have functioned in the distant past. We deleted the anticodon binding domain from tryptophanyl-tRNA synthetase and fused the discontinuous segments comprising its active site. The resulting 130 residue minimal catalytic domain activates tryptophan. This residual catalytic activity constitutes the first experimental evidence that the conserved class I signature sequences, HIGH and KMSKS, might have arisen in-frame, opposite motifs 2 and 1 from class II, as complementary sense and antisense strands of the same ancestral gene.

The RosettaDesign server identifies low energy amino acid sequences for target protein structures (http://rosettadesign.med.unc.edu). The client provides the backbone coordinates of the target structure and specifies which residues to design. The server returns to the client the sequences, coordinates and energies of the designed proteins. The simulations are performed using the design module of the Rosetta program (RosettaDesign). RosettaDesign uses Monte Carlo optimization with simulated annealing to search for amino acids that pack well on the target structure and satisfy hydrogen bonding potential. RosettaDesign has been experimentally validated and has been used previously to stabilize naturally occurring proteins and design a novel protein structure.

Over the past 10 years there has been tremendous success in the area of computational protein design. Protein design software has been used to stabilize proteins, solubilize membrane proteins, design intermolecular interactions, and design new protein structures. A key motivation for these studies is that they test our understanding of protein energetics and structure. De novo design of novel structures is a particularly rigorous test because the protein backbone must be designed in addition to the amino acid side chains. A priori it is not guaranteed that the target backbone is even designable. To address this issue, researchers have developed a variety of methods for generating protein-like scaffolds and for optimizing the protein backbone in conjunction with the amino acid sequence. These protocols have been used to design proteins from scratch and to explore sequence space for naturally occurring protein folds.

Water-mediated hydrogen bonds play critical roles at protein-protein and protein-nucleic acid interfaces, and the interactions formed by discrete water molecules cannot be captured using continuum solvent models. We describe a simple model for the energetics of water-mediated hydrogen bonds, and show that, together with knowledge of the positions of buried water molecules observed in X-ray crystal structures, the model improves the prediction of free-energy changes upon mutation at protein-protein interfaces, and the recovery of native amino acid sequences in protein interface design calculations. We then describe a "solvated rotamer" approach to efficiently predict the positions of water molecules, at protein-protein interfaces and in monomeric proteins, that is compatible with widely used rotamer-based side-chain packing and protein design algorithms. Finally, we examine the extent to which the predicted water molecules can be used to improve prediction of amino acid identities and protein-protein interface stability, and discuss avenues for overcoming current limitations of the approach.

Protein design software places amino acid side chains by precomputing rotamer-pair energies and optimizing rotamer placement. If the software optimizes by rapid stochastic techniques, then the precomputation phase dominates run time. We present a new algorithm for rapid rotamer-pair energy computation that uses a trie data structure. The trie structure avoids redundant energy computations, and lends itself to time-saving pruning techniques based on a simple geometric criteria. With our new algorithm, we compute rotamer-pair energies nearly 4 times faster than the previous approach.

Protein structures evolved through a complex interplay of cooperative interactions and it is still very challenging to design new protein folds de novo. Here, we present a strategy to design self-assembling polypeptide nanostructured polyhedra, based on modularization using orthogonal dimerizing segments. We designed end experimentally demonstrated formation of the tetrahedron that self-assembles from a single polypeptide chain comprising 12 concatenated coiled-coil-forming segments separated by flexible peptide hinges. Path of the polypeptide chain is guided by the defined order of segments that traverse each of the 6 edges of the tetrahedron exactly twice, forming coiled-coil dimers with their corresponding partners. Coincidence of the polypeptide termini in the same vertex is demonstrated by reconstitution of the split fluorescent protein by the polypeptide with the correct tetrahedral topology, while polypeptides with a deleted or scrambled segment order fail to self-assemble correctly. This design platform provides the basis for construction of new topological polypeptide folds based on the set of orthogonal interacting polypeptide segments.

Recent advances in computational protein design have allowed exciting new insights into the sequence dependence of protein folding free energy landscapes. Whereas most previous studies have examined the sequence dependence of protein stability and folding kinetics by characterizing naturally occurring proteins and variants of these proteins that contain a small number of mutations, it is now possible to generate and characterize computationally designed proteins that differ significantly from naturally occurring proteins in sequence and/or structure. These computer-generated proteins provide insights into the determinants of protein structure, stability and folding, and make it possible to disentangle the properties of proteins that are the consequence of natural selection from those that reflect the fundamental physical chemistry of polypeptide chains.

Protein-protein interactions are critical for many biological processes, and over the past several decades, this importance has prompted researchers to investigate the physical and chemical bases of protein binding. How much do we now understand about how proteins interact? Perhaps the most rigorous way to answer this question is to endeavor to design, at an atomic level of detail, such an interaction from scratch. On page 816 of this issue, Fleishman et al. ( 1 ) take on this challenge. They used a computational method that enabled them to design two proteins that bind to a preselected surface on an influenza virus. This work demonstrates how far we have come in understanding and being able to predict protein-protein interactions, and the technique could prove useful for developing biosensors, reagents, and more effective drugs.

Background Because amino acid activation is rate-limiting for uncatalyzed protein synthesis, it is a key puzzle in understanding the origin of the genetic code. Two unrelated classes (I and II) of contemporary aminoacyl-tRNA synthetases (aaRS) now translate the code. Observing that codons for the most highly conserved, Class I catalytic peptides, when read in the reverse direction, are very nearly anticodons for Class II defining catalytic peptides, Rodin and Ohno proposed that the two superfamilies descended from opposite strands of the same ancestral gene. This unusual hypothesis languished for a decade, perhaps because it appeared to be unfalsifiable. Results The proposed sense/antisense alignment makes important predictions. Fragments that align in antiparallel orientations, and contain the respective active sites, should catalyze the same two reactions catalyzed by contemporary synthetases. Recent experiments confirmed that prediction. Invariant cores from both classes, called Urzymes after Ur = primitive, authentic, plus enzyme and representing ~20% of the contemporary structures, can be expressed and exhibit high, proportionate rate accelerations for both amino-acid activation and tRNA acylation. A major fraction (60%) of the catalytic rate acceleration by contemporary synthetases resides in segments that align sense/antisense. Bioinformatic evidence for sense/antisense ancestry extends to codons specifying the invariant secondary and tertiary structures outside the active sites of the two synthetase classes. Peptides from a designed, 46-residue gene constrained by Rosetta to encode Class I and II ATP binding sites with fully complementary sequences both accelerate amino acid activation by ATP ~400 fold. Conclusions Biochemical and bioinformatic results substantially enhance the posterior probability that ancestors of the two synthetase classes arose from opposite strands of the same ancestral gene. The remarkable acceleration by short peptides of the rate-limiting step in uncatalyzed protein synthesis, together with the synergy of synthetase Urzymes and their cognate tRNAs, introduce a new paradigm for the origin of protein catalysts, emphasize the potential relevance of an operational RNA code embedded in the tRNA acceptor stems, and challenge the RNA-World hypothesis. Reviewers This article was reviewed by Dr. Paul Schimmel (nominated by Laura Landweber), Dr. Eugene Koonin and Professor David Ardell.

Mutations that increase the dynamic range of a photoswitchable protein may be used to improve other such photoswitches and increase their potential for allosteric control of protein activity in live cells.

How large is the volume of sequence space that is compatible with a given protein structure? Starting from random sequences, low free energy sequences were generated for 108 protein backbone structures by using a Monte Carlo optimization procedure and a free energy function based primarily on Lennard–Jones packing interactions and the Lazaridis–Karplus implicit solvation model. Remarkably, in the designed sequences 51% of the core residues and 27% of all residues were identical to the amino acids in the corresponding positions in the native sequences. The lowest free energy sequences obtained for ensembles of native-like backbone structures were also similar to the native sequence. Furthermore, both the individual residue frequencies and the covariances between pairs of positions observed in the very large SH3 domain family were recapitulated in core sequences designed for SH3 domain structures. Taken together, these results suggest that the volume of sequence space optimal for a protein structure is surprisingly restricted to a region around the native sequence.

In this section we have assembled a set of reviews that focus on the design of proteins and peptides with new molecular recognition capabilities. These include proteins that promote crystallization of target proteins, peptides useful for creating nanoelectronic circuits, peptide therapeutics, and switchable enzymes that respond to novel molecular inputs. Allosteric modulation of enzyme activity is typically rapid and reversible and is used to create positive and negative feedback loops in cellular signaling networks. Marc Ostermeier reviews progress in the engineering of switchable enzymes that respond to new signaling inputs. One common strategy is to fuse a catalytic domain with a binding domain that changes conformation when bound to a target ligand. The conformational change can be used to disrupt or assemble the active site as well as relieve steric inhibition of activity. Circular permutation of the host protein allows for a variety of fusion sites, for which the best sites can be identified computationally or with high-throughput screening. There is considerable interest in the design of new switches because they can be used to probe interaction dependencies and kinetics in signaling cascades. Shohei Koide describes recent progress in engineering crystallization chaperones. These are proteins that bind to specified target molecules and promote crystal formation by reducing conformational heterogeneity of the target and providing extra surface area for crystal contact formation. They can also increase or decrease the solubility of the target protein as needed. Fragments of monoclonal antibodies were among the first proteins used as crystallization chaperones and have been used to help solve structures of several important membrane proteins. However, monoclonal antibodies are slow to create and their production at a milligram scale is expensive. Koide describes several newer strategies that do not require animal immunization, but rather use molecular display technologies to identify proteins that bind tightly to the crystallization target. One advantage of this approach is that it can be used with a variety of protein scaffolds, including single chain antibodies, ankyrin repeat proteins, and fibronectin domains. This approach has been used to crystallize the full-length form of the potassium channel KcsA as well as wild-type Polo-like kinase 1. One important drawback of the chaperone approach is that it requires the production of a new chaperone for each target, for this reason, this technology will in general be reserved for high impact targets. Protein design tests our understanding of protein structure and stability. This is particularly true of computer-based studies that use physical models of protein energetics to predict amino acid sequences that will form predetermined structures, complexes or binding sites. John Karanicolas and Brian Kuhlman review recent progress in computational design of protein–protein interactions. They focus on strategies developed for enhancing the affinity of naturally occurring interactions and redesigning protein binding specificities. Particularly exciting is progress in the use of explicit negative design to create sequences that bind to specific targets within large protein families. There is considerable interest in designing proteins that adopt specific conformations when embedded in a membrane. Critical to such efforts is a good understanding of the factors that determine the forces and sequence features that dictate the insertion of transmembrane helices into membranes. Erwin London and Khurshida Shahidullah review studies that probe the relationship between sequence and helix topography in membranes. These include the effects of overall hydrophobicity, helix length, polar residue placement, and the role of flanker sequences. They demonstrate that there is a remarkable similarity between the conclusions from experiments involving the translocon-mediated insertion of model peptides and biophysical studies of fully synthetic peptides in model membranes. These findings suggest that the interaction of potential TM segments with the lipid bilayer during their residence in the translocon, or possibly that the environment experienced by a helix in the pore of the translocon is similar to that of a phospholipid bilayer. Irrespective of the precise mechanism, this paper provides a lucid overview of the features that must be considered in designing and predicting transmembrane helices. Fred Naider and Jacob Anglister describe how structural biology studies have directed the design of peptides that bind to the HIV-1 envelope protein gp41 and inhibit fusion of the virus to target cells. One of these peptides is currently used in the clinic in cases where patients have developed resistance to more standard cocktails of protease and transcriptase inhibitors. This paper provides a particularly interesting analysis of the role of structural, thermodynamic, and kinetic studies in the design and analysis of inhibitors of gp41, and also in devising strategies to combat the growing problem of resistance to this class of pharmaceutical agents. In addition to being useful as regulators of biological systems, peptides can be engineered to form novel materials. Robert Fairman and Brian Pepe-Mooney examine recent progress toward creating nanocircuits based on peptides. They describe new strategies for depositing peptides on surfaces and methods for attaching cofactors that can be used to generate new electronic properties.

Amide H/D exchange rates have been measured for the N-terminal domain of the ribosomal protein L9, residues 1-56. The rates were measured at pD 3.91, 5.03, and 5.37. At pD 5.37, 18 amides exchange slowly enough to give reliable rate measurements. At pD 3.91, seven additional residues could be followed. The exchange is shown to occur by the EX2 mechanism for all conditions studied. The rates for the N-terminal domain are very similar to those previously measured for the corresponding region in the full-length protein (Lillemoen J et al., 1997, J Mol Biol 268:482-493). In particular, the rates for the residues that we have shown to exchange via global unfolding in the N-terminal domain agree within the experimental error with the rates measured by Hoffman and coworkers, suggesting that the structure of the domain is preserved in isolation and that the stability of the isolated domain is comparable to the stability of this domain in intact L9.

The N-terminal domain, residues 1-56, of the ribosomal protein L9 has been chemically synthesized. The isolated domain is monomeric as judged by analytical ultracentrifugation and concentration-dependent CD. Complete 1H chemical shift assignments were obtained using standard methods. 2D-NMR experiments show that the isolated domain adopts the same structure as seen in the full-length protein. It consists of a three-stranded antiparallel beta-sheet sandwiched between two helixes. Thermal and urea unfolding transitions are cooperative, and the unfolding curves generated from different experimental techniques, 1D-NMR, far-UV CD, near-UV CD, and fluorescence, are superimposable. These results suggest that the protein folds by a two-state mechanism. The thermal midpoint of folding is 77 +/- 2 degrees C at pD 8.0, and the domain has a delta G degree folding = 2.8 +/- 0.8 kcal/mol at 40 degrees C, pH 7.0. Near the thermal midpoint of the unfolding transition, the 1D-NMR peaks are significantly broadened, indicating that folding is occurring on the intermediate exchange time scale. The rate of folding was determined by fitting the NMR spectra to a two-state chemical exchange model. Similar folding rates were measured for Phe 5, located in the first beta-strand, and for Tyr 25, located in the short helix between strands two and three. The domain folds extremely rapidly with a folding rate constant of 2000 s-1 near the midpoint of the equilibrium thermal unfolding transition.

The ribosomal protein L9 has an unusual structure comprising two compact globular domains connected by a 34 residue alpha-helix. The middle 17 residues of the helix are exposed to solvent while the first seven pack against and form part of the N-terminal domain, and the last ten form part of the C-terminal domain. Here we report results which show that a peptide corresponding to the central helix of L9 is monomeric in aqueous solution and85% helical at 1 degrees C and 68(+/-7)% helical at 25 degrees C. This is considerably more helical than any other protein fragment studied to date. Another peptide corresponding to the middle 17 residues of the helix is monomeric and is 41(+/-4)% helical at 1 degrees C. Because the central helix has high intrinsic stability the globular N and C-terminal domains will likely be stabilized by their interactions with the helix. Therefore, the stability of the two terminal domains should not be completely independent because both domains gain stability from a shared structural element, the central helix. Also, the ability of the central helix to form native-like structure in isolation highlights a potential role for the helix in the early stages of the folding process.

During ubiquitin ligation, an E2 conjugating enzyme receives ubiquitin from an E1 enzyme and then interacts with an E3 ligase to modify substrates. Competitive binding experiments with three human E2-E3 protein pairs show that the binding of E1s and of E3s to E2s are mutually exclusive. These results imply that polyubiquitination requires recycling of E2 for addition of successive ubiquitins to substrate.

Reengineering protein surfaces to exhibit high net charge, referred to as "supercharging", can improve reversibility of unfolding by preventing aggregation of partially unfolded states. Incorporation of charged side chains should be optimized while considering structural and energetic consequences, as numerous mutations and accumulation of like-charges can also destabilize the native state. A previously demonstrated approach deterministically mutates flexible polar residues (amino acids DERKNQ) with the fewest average neighboring atoms per side chain atom (AvNAPSA). Our approach uses Rosetta-based energy calculations to choose the surface mutations. Both protocols are available for use through the ROSIE web server. The automated Rosetta and AvNAPSA approaches for supercharging choose dissimilar mutations, raising an interesting division in surface charging strategy. Rosetta-supercharged variants of GFP (RscG) ranging from -11 to -61 and +7 to +58 were experimentally tested, and for comparison, we re-tested the previously developed AvNAPSA-supercharged variants of GFP (AscG) with +36 and -30 net charge. Mid-charge variants demonstrated ∼3-fold improvement in refolding with retention of stability. However, as we pushed to higher net charges, expression and soluble yield decreased, indicating that net charge or mutational load may be limiting factors. Interestingly, the two different approaches resulted in GFP variants with similar refolding properties. Our results show that there are multiple sets of residues that can be mutated to successfully supercharge a protein, and combining alternative supercharge protocols with experimental testing can be an effective approach for charge-based improvement to refolding.

The Rosetta molecular modeling software package provides experimentally tested and rapidly evolving tools for the 3D structure prediction and high-resolution design of proteins, nucleic acids, and a growing number of non-natural polymers. Despite its free availability to academic users and improving documentation, use of Rosetta has largely remained confined to developers and their immediate collaborators due to the code’s difficulty of use, the requirement for large computational resources, and the unavailability of servers for most of the Rosetta applications. Here, we present a unified web framework for Rosetta applications called ROSIE (Rosetta Online Server that Includes Everyone). ROSIE provides (a) a common user interface for Rosetta protocols, (b) a stable application programming interface for developers to add additional protocols, (c) a flexible back-end to allow leveraging of computer cluster resources shared by RosettaCommons member institutions, and (d) centralized administration by the RosettaCommons to ensure continuous maintenance. This paper describes the ROSIE server infrastructure, a step-by-step ‘serverification’ protocol for use by Rosetta developers, and the deployment of the first nine ROSIE applications by six separate developer teams: Docking, RNA de novo, ERRASER, Antibody, Sequence Tolerance, Supercharge, Beta peptide design, NCBB design, and VIP redesign. As illustrated by the number and diversity of these applications, ROSIE offers a general and speedy paradigm for serverification of Rosetta applications that incurs negligible cost to developers and lowers barriers to Rosetta use for the broader biological community. ROSIE is available at http://rosie.rosettacommons.org.

Noncanonical amino acids (NCAAs) can be used in a variety of protein design contexts. For example, they can be used in place of the canonical amino acids (CAAs) to improve the biophysical properties of peptides that target protein interfaces. We describe the incorporation of 114 NCAAs into the protein-modeling suite Rosetta. We describe our methods for building backbone dependent rotamer libraries and the parameterization and construction of a scoring function that can be used to score NCAA containing peptides and proteins. We validate these additions to Rosetta and our NCAA-rotamer libraries by showing that we can improve the binding of a calpastatin derived peptides to calpain-1 by substituting NCAAs for native amino acids using Rosetta. Rosetta (executables and source), auxiliary scripts and code, and documentation can be found at (http://www.rosettacommons.org/).

This article provides information to support the database article titled “UbSRD: The Ubiquitin Structural Relational Database” (Harrison et al., 2015) [1] . The ubiquitin-like homology fold (UBL) represents a large family that encompasses both post-translational modifications, like ubiquitin (UBQ) and SUMO, and functional domains on many biologically important proteins like Parkin, UHRF1 (ubiquitin-like with PDB and RING finger domains-1), and Usp7 (ubiquitin-specific protease-7) (Zhang et al., 2015; Rothbart et al., 2013; Burroughs et al., 2012; Wauer et al., 2015) [2], [3], [4], [5]. The UBL domain can participate in several unique protein–protein interactions (PPI) since protein adducts can be attached to and removed from amino groups of lysine side chains and the N-terminus of proteins. Given the biological significance of UBL domains, many have been characterized with high-resolution techniques, and for UBQ and SUMO, many protein complexes have been characterized. We identified all the UBL domains in the PDB and created a relational database called UbSRD (Ubiquitin Structural Relational Database) by using structural analysis tools in the Rosetta (Leaver et al., 2013; O’Meara et al., 2015; Leaver-fay et al., 2011) [1], [6], [7], [8]. Querying UbSRD permitted us to report many quantitative properties of UBQ and SUMO recognition at different types interfaces (noncovalent: NC, conjugated: CJ, and deubiquitanse: DB). In this data article, we report the average number of non-UBL neighbors, secondary structure of interacting motifs, and the type of inter-molecular hydrogen bonds for each residue of UBQ and SUMO. Additionally, we used PROMALS3D to generate a multiple sequence alignment used to construct a phylogram for the entire set of UBLs (Pei and Grishin, 2014) [9]. The data described here will be generally useful to scientists studying the molecular basis for recognition of UBQ or SUMO.