RosettaScripts: A Scripting Language Interface to the Rosetta Macromolecular Modeling Suite

Macromolecular modeling and design are increasingly useful in basic research, biotechnology, and teaching. However, the absence of a user-friendly modeling framework that provides access to a wide range of modeling capabilities is hampering the wider adoption of computational methods by non-experts. RosettaScripts is an XML-like language for specifying modeling tasks in the Rosetta framework. RosettaScripts provides access to protocol-level functionalities, such as rigid-body docking and sequence redesign, and allows fast testing and deployment of complex protocols without need for modifying or recompiling the underlying C++ code. We illustrate these capabilities with RosettaScripts protocols for the stabilization of proteins, the generation of computationally constrained libraries for experimental selection of higher-affinity binding proteins, loop remodeling, small-molecule ligand docking, design of ligand-binding proteins, and specificity redesign in DNA-binding proteins

Evolution of Anolis Lizard Dewlap Diversity

BACKGROUND: The dewlaps of Anolis lizards provide a classic example of a complex signaling system whose function and evolution is poorly understood. Dewlaps are flaps of skin beneath the chin that are extended and combined with head and body movements for visual signals and displays. They exhibit extensive morphological variation and are one of two cladistic features uniting anoles, yet little is known regarding their function and evolution. We quantified the diversity of anole dewlaps, investigated whether dewlap morphology was informative regarding phylogenetic relationships, and tested two separate hypotheses: (A) similar Anolis habitat specialists possess similar dewlap configurations (Ecomorph Convergence hypothesis), and (B) sympatric species differ in their dewlap morphologies to a greater extent than expected by chance (Species Recognition hypothesis). METHODOLOGY/PRINCIPAL FINDINGS: We found that dewlap configurations (sizes, patterns and colors) exhibit substantial diversity, but that most are easily categorized into six patterns that incorporate one to three of 13 recognizable colors. Dewlap morphology is not phylogenetically informative and, like other features of anoles, exhibits convergence in configurations. We found no support for the Ecomorph Convergence hypothesis; species using the same structural habitat were no more similar in dewlap configuration than expected by chance. With one exception, all sympatric species in four communities differ in dewlap configuration. However, this provides only weak support for the Species Recognition hypothesis because, due to the great diversity in dewlap configurations observed across each island, few cases of sympatric species with identical dewlaps would be expected to co-occur by chance alone. CONCLUSIONS/SIGNIFICANCE: Despite previous thought, most dewlaps exhibit easily characterizable patterns and colorations. Nevertheless, dewlap variation is extensive and explanations for the origin and evolution of this diversity are lacking. Our data do not support two hypothesized explanations for this diversity, but others such as sexual selection remain to be tested

Benchmarking and Analysis of Protein Docking Performance in Rosetta v3.2

RosettaDock has been increasingly used in protein docking and design strategies in order to predict the structure of protein-protein interfaces. Here we test capabilities of RosettaDock 3.2, part of the newly developed Rosetta v3.2 modeling suite, against Docking Benchmark 3.0, and compare it with RosettaDock v2.3, the latest version of the previous Rosetta software package. The benchmark contains a diverse set of 116 docking targets including 22 antibody-antigen complexes, 33 enzyme-inhibitor complexes, and 60 ‘other’ complexes. These targets were further classified by expected docking difficulty into 84 rigid-body targets, 17 medium targets, and 14 difficult targets. We carried out local docking perturbations for each target, using the unbound structures when available, in both RosettaDock v2.3 and v3.2. Overall the performances of RosettaDock v2.3 and v3.2 were similar. RosettaDock v3.2 achieved 56 docking funnels, compared to 49 in v2.3. A breakdown of docking performance by protein complex type shows that RosettaDock v3.2 achieved docking funnels for 63% of antibody-antigen targets, 62% of enzyme-inhibitor targets, and 35% of ‘other’ targets. In terms of docking difficulty, RosettaDock v3.2 achieved funnels for 58% of rigid-body targets, 30% of medium targets, and 14% of difficult targets. For targets that failed, we carry out additional analyses to identify the cause of failure, which showed that binding-induced backbone conformation changes account for a majority of failures. We also present a bootstrap statistical analysis that quantifies the reliability of the stochastic docking results. Finally, we demonstrate the additional functionality available in RosettaDock v3.2 by incorporating small-molecules and non-protein co-factors in docking of a smaller target set. This study marks the most extensive benchmarking of the RosettaDock module to date and establishes a baseline for future research in protein interface modeling and structure prediction

Age-Related Attenuation of Dominant Hand Superiority

The decline of motor performance of the human hand-arm system with age is well-documented. While dominant hand performance is superior to that of the non-dominant hand in young individuals, little is known of possible age-related changes in hand dominance. We investigated age-related alterations of hand dominance in 20 to 90 year old subjects. All subjects were unambiguously right-handed according to the Edinburgh Handedness Inventory. In Experiment 1, motor performance for aiming, postural tremor, precision of arm-hand movement, speed of arm-hand movement, and wrist-finger speed tasks were tested. In Experiment 2, accelerometer-sensors were used to obtain objective records of hand use in everyday activities

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

Anchored Design of Protein-Protein Interfaces

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.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.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

RosettaRemodel: A Generalized Framework for Flexible Backbone Protein Design

We describe RosettaRemodel, a generalized framework for flexible protein design that provides a versatile and convenient interface to the Rosetta modeling suite. RosettaRemodel employs a unified interface, called a blueprint, which allows detailed control over many aspects of flexible backbone protein design calculations. RosettaRemodel allows the construction and elaboration of customized protocols for a wide range of design problems ranging from loop insertion and deletion, disulfide engineering, domain assembly, loop remodeling, motif grafting, symmetrical units, to de novo structure modeling

Structure-Based Design of Non-Natural Amino Acid Inhibitors of Amyloid Fibrillation

Many globular and natively disordered proteins can convert into amyloid fibers. These fibers are associated with numerous pathologies1 as well as with normal cellular functions2,3, and frequently form during protein denaturation4,5. Inhibitors of pathological amyloid fibers could serve as leads for therapeutics, provided the inhibitors were specific enough to avoid interfering with normal processes. Here we show that computer-aided, structure-based design can yield highly specific peptide inhibitors of amyloid formation. Using known atomic structures of segments of amyloid fibers as templates, we have designed and characterized an all D-amino acid inhibitor of fibrillation of the tau protein found in Alzheimer’s disease, and a non-natural L-amino acid inhibitor of an amyloid fiber that enhances sexual transmission of HIV. Our results indicate that peptides from structure-based designs can disrupt the fibrillation of full-length proteins, including those like tau that lack fully ordered native structures.We thank M.I. Ivanova, J. Corn, T. Kortemme, D. Anderson, M.R. Sawaya, M. Phillips, S. Sambashivan, J. Park, M. Landau, Q. Zhang, R. Clubb, F. Guo, T. Yeates, J. Nowick, J. Zheng, and M.J. Thompson for discussions, HHMI, NIH, NSF, the GATES foundation, and the Joint Center for Translational Medicine for support, R. Peterson for help with NMR experiments, E. Mandelkow for providing tau constructs, R. Riek for providing amyloid beta, J. Stroud for amyloid beta preparation. Support for JK was from the Damon Runyon Cancer Research Foundation, for HWC by the Ruth L. Kirschstein National Research Service Award, for JM from the programme for junior-professors by the ministry of science, Baden-Württemberg, and for SAS by a UCLA-IGERT bioinformatics traineeship

Atomic-accuracy prediction of protein loop structures through an RNA-inspired ansatz

Consistently predicting biopolymer structure at atomic resolution from sequence alone remains a difficult problem, even for small sub-segments of large proteins. Such loop prediction challenges, which arise frequently in comparative modeling and protein design, can become intractable as loop lengths exceed 10 residues and if surrounding side-chain conformations are erased. This article introduces a modeling strategy based on a 'stepwise ansatz', recently developed for RNA modeling, which posits that any realistic all-atom molecular conformation can be built up by residue-by-residue stepwise enumeration. When harnessed to a dynamic-programming-like recursion in the Rosetta framework, the resulting stepwise assembly (SWA) protocol enables enumerative sampling of a 12 residue loop at a significant but achievable cost of thousands of CPU-hours. In a previously established benchmark, SWA recovers crystallographic conformations with sub-Angstrom accuracy for 19 of 20 loops, compared to 14 of 20 by KIC modeling with a comparable expenditure of computational power. Furthermore, SWA gives high accuracy results on an additional set of 15 loops highlighted in the biological literature for their irregularity or unusual length. Successes include cis-Pro touch turns, loops that pass through tunnels of other side-chains, and loops of lengths up to 24 residues. Remaining problem cases are traced to inaccuracies in the Rosetta all-atom energy function. In five additional blind tests, SWA achieves sub-Angstrom accuracy models, including the first such success in a protein/RNA binding interface, the YbxF/kink-turn interaction in the fourth RNA-puzzle competition. These results establish all-atom enumeration as a systematic approach to protein structure that can leverage high performance computing and physically realistic energy functions to more consistently achieve atomic resolution.Comment: Identity of four-loop blind test protein and parts of figures 5 have been omitted in this preprint to ensure confidentiality of the protein structure prior to its public releas