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

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

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

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

Gravitational Waves from Gravitational Collapse

Gravitational wave emission from the gravitational collapse of massive stars has been studied for more than three decades. Current state of the art numerical investigations of collapse include those that use progenitors with realistic angular momentum profiles, properly treat microphysics issues, account for general relativity, and examine non--axisymmetric effects in three dimensions. Such simulations predict that gravitational waves from various phenomena associated with gravitational collapse could be detectable with advanced ground--based and future space--based interferometric observatories.Comment: 68 pages including 13 figures; revised version accepted for publication in Living Reviews in Relativity (http://www.livingreviews.org

Elongation of very long-chain (>C24) fatty acids in Clarias gariepinus: Cloning, functional characterization and tissue expression of elovl4 elongases

Elongation of very long-chain fatty acid 4 (Elovl4) proteins participate in the biosynthesis of very long-chain (>C24) saturated and polyunsaturated fatty acids (FA). Previous studies have shown that fish possess two different forms of Elovl4, termed Elovl4a and Elovl4b. The present study aimed to characterize both molecularly and functionally two elovl4 cDNA from the African catfish Clarias gariepinus. The results confirmed that C. gariepinus possessed two elovl4-like elongases with high homology to two previously characterized Elovl4 from Danio rerio, and thus they were termed accordingly as Elovl4a and Elovl4b. The C. gariepinus Elovl4a and Elovl4b have open reading frames (ORF) of 945 and 915 base pairs, respectively, encoding putative proteins of 314 and 304 amino acids, respectively. Functional characterization in yeast showed both Elovl4 enzymes have activity towards all the PUFA substrates assayed (18:4n-3, 18:3n-6, 20:5n-3, 20:4n-6, 22:5n-3, 22:4n-6 and 22:6n-3), producing elongated products of up to C36. Moreover, the C. gariepinus Elovl4a and Elovl4b were able to elongate very long-chain saturated FA (VLC-SFA) as denoted by increased levels of 28:0 and longer FA in yeast transformed with elovl4 ORF compared to control yeast. These results confirmed that C. gariepinus Elovl4 play important roles in the biosynthesis of very long-chain FA. Tissue distribution analysis of elovl4 mRNAs showed both genes were widely expressed in all tissues analyzed, with high expression of elovl4a in pituitary and brain, whereas female gonad and pituitary had the highest expression levels for elovl4b

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

Transient Protein-Protein Interaction of the SH3-Peptide Complex via Closely Located Multiple Binding Sites

Protein-protein interactions play an essential role in cellular processes. Certain proteins form stable complexes with their partner proteins, whereas others function by forming transient complexes. The conventional protein-protein interaction model describes an interaction between two proteins under the assumption that a protein binds to its partner protein through a single binding site. In this study, we improved the conventional interaction model by developing a Multiple-Site (MS) model in which a protein binds to its partner protein through closely located multiple binding sites on a surface of the partner protein by transiently docking at each binding site with individual binding free energies. To test this model, we used the protein-protein interaction mediated by Src homology 3 (SH3) domains. SH3 domains recognize their partners via a weak, transient interaction and are therefore promiscuous in nature. Because the MS model requires large amounts of data compared with the conventional interaction model, we used experimental data from the positionally addressable syntheses of peptides on cellulose membranes (SPOT-synthesis) technique. From the analysis of the experimental data, individual binding free energies for each binding site of peptides were extracted. A comparison of the individual binding free energies from the analysis with those from atomistic force fields gave a correlation coefficient of 0.66. Furthermore, application of the MS model to 10 SH3 domains lowers the prediction error by up to 9% compared with the conventional interaction model. This improvement in prediction originates from a more realistic description of complex formation than the conventional interaction model. The results suggested that, in many cases, SH3 domains increased the protein complex population through multiple binding sites of their partner proteins. Our study indicates that the consideration of general complex formation is important for the accurate description of protein complex formation, and especially for those of weak or transient protein complexes

Relative effectiveness and adverse effects of cervical manipulation, mobilisation and the activator instrument in patients with sub-acute non-specific neck pain: results from a stopped randomised trial

<p>Abstract</p> <p>Background</p> <p>Neck pain of a mechanical nature is a common complaint seen by practitioners of manual medicine, who use a multitude of methods to treat the condition. It is not known, however, if any of these methods are superior in treatment effectiveness. This trial was stopped due to poor recruitment. The purposes of this report are (1) to describe the trial protocol, (2) to report on the data obtained from subjects who completed the study, (3) to discuss the problems we encountered in conducting this study.</p> <p>Methods</p> <p>A pragmatic randomised clinical trial was undertaken. Patients who met eligibility criteria were randomised into three groups. One group was treated using specific segmental high velocity low amplitude manipulation (diversified), another by specific segmental mobilisation, and a third group by the Activator instrument. All three groups were also treated for any myofascial distortions and given appropriate exercises and advice. Participants were treated six times over a three-week period or until they reported being pain free. The primary outcome measure for the study was Patient Global Impression of Change (PGIC); secondary outcome measures included the Short-Form Health Survey (SF-36v₂), the neck Bournemouth Questionnaire, and the numerical rating scale for pain intensity. Participants also kept a diary of any pain medication taken and noted any perceived adverse effects of treatment. Outcomes were measured at four points: end of treatment, and 3, 6, and 12 months thereafter.</p> <p>Results</p> <p>Between January 2007 and March 2008, 123 patients were assessed for eligibility, of these 47 were considered eligible, of which 16 were allocated to manipulation, 16 to the Activator instrument and 15 to the mobilisation group. Comparison between the groups on the PGIC adjusted for baseline covariants did not show a significant difference for any of the endpoints. Within group analyses for change from baseline to the 12-month follow up for secondary outcomes were significant for all groups on the Bournemouth Questionnaire and for pain, while the mobilisation group had a significant improvement on the PCS and MCS subscales of the SF-36_v2. Finally, there were no moderate, severe, or long-lasting adverse effects reported by any participant in any group.</p> <p>Conclusions</p> <p>Although the small sample size must be taken into consideration, it appears that all three methods of treating mechanical neck pain had a long-term benefit for subacute neck pain, without moderate or serious adverse events associated with any of the treatment methods. There were difficulties in recruiting subjects to this trial. This pragmatic trial should be repeated with a larger sample size.</p