Four small puzzles that Rosetta doesn't solve

A complete macromolecule modeling package must be able to solve the simplest structure prediction problems. Despite recent successes in high resolution structure modeling and design, the Rosetta software suite fares poorly on deceptively small protein and RNA puzzles, some as small as four residues. To illustrate these problems, this manuscript presents extensive Rosetta results for four well-defined test cases: the 20-residue mini-protein Trp cage, an even smaller disulfide-stabilized conotoxin, the reactive loop of a serine protease inhibitor, and a UUCG RNA tetraloop. In contrast to previous Rosetta studies, several lines of evidence indicate that conformational sampling is not the major bottleneck in modeling these small systems. Instead, approximations and omissions in the Rosetta all-atom energy function currently preclude discriminating experimentally observed conformations from de novo models at atomic resolution. These molecular "puzzles" should serve as useful model systems for developers wishing to make foundational improvements to this powerful modeling suite.Comment: Published in PLoS One as a manuscript for the RosettaCon 2010 Special Collectio

Predicting the Tolerated Sequences for Proteins and Protein Interfaces Using RosettaBackrub Flexible Backbone Design

Predicting the set of sequences that are tolerated by a protein or protein interface, while maintaining a desired function, is useful for characterizing protein interaction specificity and for computationally designing sequence libraries to engineer proteins with new functions. Here we provide a general method, a detailed set of protocols, and several benchmarks and analyses for estimating tolerated sequences using flexible backbone protein design implemented in the Rosetta molecular modeling software suite. The input to the method is at least one experimentally determined three-dimensional protein structure or high-quality model. The starting structure(s) are expanded or refined into a conformational ensemble using Monte Carlo simulations consisting of backrub backbone and side chain moves in Rosetta. The method then uses a combination of simulated annealing and genetic algorithm optimization methods to enrich for low-energy sequences for the individual members of the ensemble. To emphasize certain functional requirements (e.g. forming a binding interface), interactions between and within parts of the structure (e.g. domains) can be reweighted in the scoring function. Results from each backbone structure are merged together to create a single estimate for the tolerated sequence space. We provide an extensive description of the protocol and its parameters, all source code, example analysis scripts and three tests applying this method to finding sequences predicted to stabilize proteins or protein interfaces. The generality of this method makes many other applications possible, for example stabilizing interactions with small molecules, DNA, or RNA. Through the use of within-domain reweighting and/or multistate design, it may also be possible to use this method to find sequences that stabilize particular protein conformations or binding interactions over others

Neurospora COP9 Signalosome Integrity Plays Major Roles for Hyphal Growth, Conidial Development, and Circadian Function

The COP9 signalosome (CSN) is a highly conserved multifunctional complex that has two major biochemical roles: cleaving NEDD8 from cullin proteins and maintaining the stability of CRL components. We used mutation analysis to confirm that the JAMM domain of the CSN-5 subunit is responsible for NEDD8 cleavage from cullin proteins in Neurospora crassa. Point mutations of key residues in the metal-binding motif (EXnHXHX10D) of the CSN-5 JAMM domain disrupted CSN deneddylation activity without interfering with assembly of the CSN complex or interactions between CSN and cullin proteins. Surprisingly, CSN-5 with a mutated JAMM domain partially rescued the phenotypic defects observed in a csn-5 mutant. We found that, even without its deneddylation activity, the CSN can partially maintain the stability of the SCFFWD-1 complex and partially restore the degradation of the circadian clock protein FREQUENCY (FRQ) in vivo. Furthermore, we showed that CSN containing mutant CSN-5 efficiently prevents degradation of the substrate receptors of CRLs. Finally, we found that deletion of the CAND1 ortholog in N. crassa had little effect on the conidiation circadian rhythm. Our results suggest that CSN integrity plays major roles in hyphal growth, conidial development, and circadian function in N. crassa

A Generic Program for Multistate Protein Design

Some protein design tasks cannot be modeled by the traditional single state design strategy of finding a sequence that is optimal for a single fixed backbone. Such cases require multistate design, where a single sequence is threaded onto multiple backbones (states) and evaluated for its strengths and weaknesses on each backbone. For example, to design a protein that can switch between two specific conformations, it is necessary to to find a sequence that is compatible with both backbone conformations. We present in this paper a generic implementation of multistate design that is suited for a wide range of protein design tasks and demonstrate in silico its capabilities at two design tasks: one of redesigning an obligate homodimer into an obligate heterodimer such that the new monomers would not homodimerize, and one of redesigning a promiscuous interface to bind to only a single partner and to no longer bind the rest of its partners. Both tasks contained negative design in that multistate design was asked to find sequences that would produce high energies for several of the states being modeled. Success at negative design was assessed by computationally redocking the undesired protein-pair interactions; we found that multistate design's accuracy improved as the diversity of conformations for the undesired protein-pair interactions increased. The paper concludes with a discussion of the pitfalls of negative design, which has proven considerably more challenging than positive design

Capzimin is a potent and specific inhibitor of proteasome isopeptidase Rpn11

The proteasome is a vital cellular machine that maintains protein homeostasis, which is of particular importance in multiple myeloma and possibly other cancers. Targeting of proteasome 20S peptidase activity with bortezomib and carfilzomib has been widely used to treat myeloma. However, not all patients respond to these compounds, and those who do eventually suffer relapse. Therefore, there is an urgent and unmet need to develop new drugs that target proteostasis through different mechanisms. We identified quinoline-8-thiol (8TQ) as a first-in-class inhibitor of the proteasome 19S subunit Rpn11. A derivative of 8TQ, capzimin, shows >5-fold selectivity for Rpn11 over the related JAMM proteases and >2 logs selectivity over several other metalloenzymes. Capzimin stabilized proteasome substrates, induced an unfolded protein response, and blocked proliferation of cancer cells, including those resistant to bortezomib. Proteomic analysis revealed that capzimin stabilized a subset of polyubiquitinated substrates. Identification of capzimin offers an alternative path to develop proteasome inhibitors for cancer therapy

Metal-directed, chemically tunable assembly of one-, two- and three-dimensional crystalline protein arrays.

Proteins represent the most sophisticated building blocks available to an organism and to the laboratory chemist. Yet, in contrast to nearly all other types of molecular building blocks, the designed self-assembly of proteins has largely been inaccessible because of the chemical and structural heterogeneity of protein surfaces. To circumvent the challenge of programming extensive non-covalent interactions to control protein self-assembly, we have previously exploited the directionality and strength of metal coordination interactions to guide the formation of closed, homoligomeric protein assemblies. Here, we extend this strategy to the generation of periodic protein arrays. We show that a monomeric protein with properly oriented coordination motifs on its surface can arrange, on metal binding, into one-dimensional nanotubes and two- or three-dimensional crystalline arrays with dimensions that collectively span nearly the entire nano- and micrometre scale. The assembly of these arrays is tuned predictably by external stimuli, such as metal concentration and pH