Search

The Rosetta software for macromolecular modeling, docking and design is extensively used in laboratories worldwide. During two decades of development by a community of laboratories at more than 60 institutions, Rosetta has been continuously refactored and extended. Its advantages are its performance and interoperability between broad modeling capabilities. Here we review tools developed in the last 5 years, including over 80 methods. We discuss improvements to the score function, user interfaces and usability. Rosetta is available at ., This Perspective reviews tools developed over the past five years in the macromolecular modeling, docking and design software Rosetta.

Structure-based methods are powerful tools that are being exploited to unravel new functions with therapeutic advantage. Here, we report the discovery of a new class of deaminases, predominantly found in mycobacterial species that act on the commercially important s-triazine class of compounds. The enzyme Msd from Mycobacterium smegmatis was taken as a representative candidate from an evolutionarily conserved subgroup that possesses high density of Mycobacterium deaminases. Biochemical investigation reveals that Msd specifically acts on mutagenic nucleobases such as 5-azacytosine and isoguanine and does not accept natural bases as substrates. Determination of the X-ray structure of Msd to a resolution of 1.9 angstrom shows that Msd has fine-tuned its active site such that it is a hybrid of a cytosine as well as a guanine deaminase, thereby conferring Msd the ability to expand its repertoire to both purine and pyrimidine-like mutagens. Mapping of active site residues along with X-ray structures with a series of triazine analogues aids in deciphering the mechanism by which Msd proofreads the base milieu for mutagens. The genome location of the enzyme reveals that Msd is part of a conserved cluster that confers the organism with innate resistance: toward select xenobiotics by triggering their efflux.

International audience; To adapt to changing environments, bacteria have evolved numerous pathways that activate stress response genes. In Gram-positive bacteria, the stressosome, a cytoplasmic complex, relays external cues and activates the sigma B regulon. The stressosome is structurally well-characterized in Bacillus, but how it senses stress remains elusive. Here, we report a genome-wide N-terminomic approach in Listeria that strikingly led to the discovery of 19 internal translation initiation sites and 6 miniproteins, among which one, Prli42, is conserved in Firmicutes. Prli42 is membrane-anchored and interacts with orthologues of Bacillus stressosome components. We reconstituted the Listeria stressosome in vitro and visualized its supramolecular structure by electron microscopy. Analysis of a series of Prli42 mutants demonstrated that Prli42 is important for sigma B activation, bacterial growth following oxidative stress and for survival in macrophages. Taken together, our N-terminonic approach unveiled Prli42 as a long-sought link between stress and the stressosome.

Contains fulltext : 174131.pdf (Publisher’s version ) (Open Access) CorNet is a web-based tool for the analysis of co-evolving residue positions in protein super-family sequence alignments. CorNet projects external information such as mutation data extracted from literature on interactively displayed groups of co-evolving residue positions to shed light on the functions associated with these groups and the residues in them. We used CorNet to analyse six enzyme super-families and found that groups of strongly co-evolving residues tend to consist of residues involved in a same function such as activity, specificity, co-factor binding, or enantioselectivity. This finding allows to assign a function to residues for which no data is available yet in the literature. A mutant library was designed to mutate residues observed in a group of co-evolving residues predicted to be involved in enantioselectivity, but for which no literature data is available yet. The resulting set of mutations indeed showed many instances of increased enantioselectivity.

NE0047 from Nitrosomonas europaea has been annotated as a zinc-dependent deaminase; however, the substrate specificity is unknown because of the low level of structural similarity and sequence identity compared to other family members. In this study, the function of NE0047 was established as a guanine deaminase (catalytic efficiency of 1.2 X 10(5) M-1 s(-1)), exhibiting secondary activity towards ammeline. The structure of NE0047 in the presence of the substrate analogue 8-azaguanine was also determined to a resolution of 1.9 angstrom. NE0047 crystallized as a homodimer in an asymmetric unit. It was found that the extreme nine-amino acid C-terminal loop forms an active site flap; in one monomer, the flap is in the closed conformation and in the other in the open conformation with this loop region exposed to the solvent. Calorimetric data obtained using the full-length version of the enzyme fit to a sequential binding model, thus supporting a cooperative mode of ligand occupancy. In contrast, the mutant form of the enzyme (Delta C) with the deletion of the extreme nine amino acids follows an independent model of ligand occupancy. In addition, the Delta C mutant also does not exhibit any enzyme activity. Therefore, we propose that the progress of the reaction is communicated via changes in the conformation of the C-terminal flap and the closed form of the enzyme is the catalytically active form, while the open form allows for product release. The catalytic mechanism of deamination was also investigated, and we found that the mutagenesis of the highly conserved active site residues Glu79 and Glu143 resulted in a complete loss of activity and concluded that they facilitate the reaction by serving as proton shuttles.

The efficiency of using a variant of Hamiltonian replica-exchange molecular dynamics (Chaperone H-replica-exchange molecular dynamics [CH-REMD]) for the refinement of protein structural models generated de novo is investigated. In CH-REMD, the interaction between the protein and its environment, specifically, the electrostatic interaction between the protein and the solvating water, is varied leading to cycles of partial unfolding and refolding mimicking some aspects of folding chaperones. In 10 of the 15 cases examined, the CH-REMD approach sampled structures in which the root-mean-square deviation (RMSD) of secondary structure elements (SSE-RMSD) with respect to the experimental structure was more than 1.0 angstrom lower than the initial de novo model. In 14 of the 15 cases, the improvement was more than 0.5 angstrom. The ability of three different statistical potentials to identify near-native conformations was also examined. Little correlation between the SSE-RMSD of the sampled structures with respect to the experimental structure and any of the scoring functions tested was found. The most effective scoring function tested was the DFIRE potential. Using the DFIRE potential, the SSE-RMSD of the best scoring structures was on average 0.3 angstrom lower than the initial model. Overall the work demonstrates that targeted enhanced-sampling techniques such as CH-REMD can lead to the systematic refinement of protein structural models generated de novo but that improved potentials for the identification of near-native structures are still needed. Proteins 2012; (c) 2012 Wiley Periodicals, Inc.

A protocol is presented for the global refinement of homology models of proteins. It combines the advantages of temperature-based replica-exchange molecular dynamics (REMD) for conformational sampling and the use of statistical potentials for model selection. The protocol was tested using 21 models. Of these 14 were models of 10 small proteins for which high-resolution crystal structures were available, the remainder were targets of the recent CASPR exercise. It was found that REMD in combination with currently available force fields could sample near-native conformational states starting from high-quality homology models. Conformations in which the backbone RMSD of secondary structure elements (SSE-RMSD) was lower than the starting value by 0.5-1.0 angstrom were found for 15 out of the 21 cases (average 0.82 angstrom). Furthermore, when a simple scoring function consisting of two statistical potentials was used to rank the structures, one or more structures with SSE-RMSD of at least 0.2 angstrom lower than the starting valise was found among the five best ranked structures in 11 out of the 21 cases. The average improvement in SSE-RMSD for the best models was 0.42 angstrom. However, none of the scoring functions tested identified the structures with the lowest SSE-RMSD as the best models although all identified the native conformation as the one with lowest energy. This suggests that while the proposed protocol proved effective for the refinement of high-quality models of small proteins scoring functions remain one of the major limiting factors in structure refinement. This and other aspects by which the methodology could be further improved are discussed.

Generally, non-catalytic carbohydrate binding module (CBM) specificity has been shown to parallel the catalytic activity of the carbohydrate active enzyme (CAZyme) module it is appended to. With the rapid expansion in metagenomic sequence space for the potential discovery of new CBMs in addition to the recent emergence of several new CBM families that display diverse binding profiles and novel functions, elucidating the function of these protein modules has become a much more challenging task. This review summarizes several approaches that have been reported for using primary structure to inform CBM specificity and streamlining their biophysical characterization. In addition we discuss general trends in binding site architecture and several newly identified functions for CBMs. Streams of investigation that will facilitate the development and refinement of sequence-based prediction tools are suggested.

Acetoin reductase is an important enzyme for the fermentative production of 2,3-butanediol, a chemical compound with a very broad industrial use. Here, we report on the discovery and characterization of an acetoin reductase from Clostridium beijerinckii NCIMB 8052. An in silico screen of the C. beijerinckii genome revealed eight potential acetoin reductases. One of them ( CBEI_1464 ) showed substantial acetoin reductase activity after expression in Escherichia coli . The purified enzyme ( C. beijerinckii acetoin reductase [Cb-ACR]) was found to exist predominantly as a homodimer. In addition to acetoin (or 2,3-butanediol), other secondary alcohols and corresponding ketones were converted as well, provided that another electronegative group was attached to the adjacent C-3 carbon. Optimal activity was at pH 6.5 (reduction) and 9.5 (oxidation) and around 68°C. Cb-ACR accepts both NADH and NADPH as electron donors; however, unlike closely related enzymes, NADPH is preferred ( K m , 32 μM). Cb-ACR was compared to characterized close homologs, all belonging to the “threonine dehydrogenase and related Zn-dependent dehydrogenases” (COG1063). Metal analysis confirmed the presence of 2 Zn 2+ atoms. To gain insight into the substrate and cofactor specificity, a structural model was constructed. The catalytic zinc atom is likely coordinated by Cys 37 , His 70 , and Glu 71 , while the structural zinc site is probably composed of Cys 100 , Cys 103 , Cys 106 , and Cys 114 . Residues determining NADP specificity were predicted as well. The physiological role of Cb-ACR in C. beijerinckii is discussed.

The Enterococcus faecalis leucine-rich protein ElrA promotes virulence by stimulating bacterial persistence in macrophages and production of the interleukin-6 (IL-6) cytokine. The ElrA protein is encoded within an operon that is poorly expressed under laboratory conditions but induced in vivo . In this study, we identify ef2687 (renamed elrR ), which encodes a member of the Rgg ( r egulator g ene for g lucosyltransferase) family of putative regulatory proteins. Using quantitative reverse transcription-PCR, translational lacZ fusions, and electrophoretic mobility shift assays, we demonstrate that ElrR positively regulates expression of elrA . These results correlate with the attenuated virulence of the Δ elrR strain in a mouse peritonitis model. Virulence of simple and double elrR and elrA deletion mutants also suggests a remaining ElrR-independent expression of elrA in vivo and additional virulence-related genes controlled by ElrR.

The efficiency of using a variant of Hamiltonian replica-exchange molecular dynamics (Chaperone H-replica-exchange molecular dynamics [CH-REMD]) for the refinement of protein structural models generated de novo is investigated. In CH-REMD, the interaction between the protein and its environment, specifically, the electrostatic interaction between the protein and the solvating water, is varied leading to cycles of partial unfolding and refolding mimicking some aspects of folding chaperones. In 10 of the 15 cases examined, the CH-REMD approach sampled structures in which the root-mean-square deviation (RMSD) of secondary structure elements (SSE-RMSD) with respect to the experimental structure was more than 1.0 angstrom lower than the initial de novo model. In 14 of the 15 cases, the improvement was more than 0.5 angstrom. The ability of three different statistical potentials to identify near-native conformations was also examined. Little correlation between the SSE-RMSD of the sampled structures with respect to the experimental structure and any of the scoring functions tested was found. The most effective scoring function tested was the DFIRE potential. Using the DFIRE potential, the SSE-RMSD of the best scoring structures was on average 0.3 angstrom lower than the initial model. Overall the work demonstrates that targeted enhanced-sampling techniques such as CH-REMD can lead to the systematic refinement of protein structural models generated de novo but that improved potentials for the identification of near-native structures are still needed. Proteins 2012; (c) 2012 Wiley Periodicals, Inc.

Background: Feruloyl esterases (FAEs) are important biomass degrading accessory enzymes due to their capability of cleaving the ester links between hemicellulose and pectin to aromatic compounds of lignin, thus enhancing the accessibility of plant tissues to cellulolytic and hemicellulolytic enzymes. FAEs have gained increased attention in the area of biocatalytic transformations for the synthesis of value added compounds with medicinal and nutritional applications. Following the increasing attention on these enzymes, a novel descriptor based classification system has been proposed for FAEs resulting into 12 distinct families and pharmacophore models for three FAE sub-families have been developed. Methodology/Principal Findings: The feruloylome of Aspergillus oryzae contains 13 predicted FAEs belonging to six sub-families based on our recently developed descriptor-based classification system. The three-dimensional structures of the 13 FAEs were modeled for structural analysis of the feruloylome. The three genes coding for three enzymes, viz., A.O.2, A.O.8 and A.O.10 from the feruloylome of A. oryzae, representing sub-families with unknown functional features, were heterologously expressed in Pichia pastoris, characterized for substrate specificity and structural characterization through CD spectroscopy. Common feature-based pharamacophore models were developed according to substrate specificity characteristics of the three enzymes. The active site residues were identified for the three expressed FAEs by determining the titration curves of amino acid residues as a function of the pH by applying molecular simulations. Conclusions/Significance: Our findings on the structure-function relationships and substrate specificity of the FAEs of A. oryzae will be instrumental for further understanding of the FAE families in the novel classification system. The developed pharmacophore models could be applied for virtual screening of compound databases for short listing the putative substrates prior to docking studies or for post-processing docking results to remove false positives. Our study exemplifies how computational predictions can complement to the information obtained through experimental methods. © 2012 Udatha et al., published_or_final_version

Contains fulltext : 186518.pdf (Publisher’s version ) (Closed access) A method is presented to refine models built by homology by the use of restricted molecular dynamics (MD) techniques. The basic idea behind this method is the use of structure validation software to determine for each residue the likelihood that it is modeled correctly. This information is used to determine constraints and restraints in an MD simulation including explicit solvent molecules, which is used for model refinement. The procedure is based on the idea that residues that the validation software identifies as correctly positioned should be strongly constrained or restrained in the MD simulations, whereas residues that are likely to be positioned wrongly should move freely. Two different protocols are compared: one (applied to CASP3 target T58) using full structural constraints with separate optimization of each short fragment and the other (applied to T47) allowing some freedom using harmonic restraining potentials, with automatic optimization of the whole molecule. Structures along the MD trajectory that scored best in structural checks were selected for the construction of models that appeared to be successful in the CASP3 competition. Model refinement with MD in general leads to a model that is less like the experimental structure (Levitt et al. Nature Struct Biol 1999;6:108-111). Actually, refined T47 was slightly improved compared to the starting model; changes in model T58 led not to further enhancement. After the X-ray structure of the modeled proteins became known, the procedure was evaluated for two targets (T47 and the CASP4 target T111) by comparing a long simulation in water with the experimental target structures. It was found that structural improvements could be obtained on a nanosecond time scale by allowing appropriate freedom in the simulation. Structural checks applied to fast fluctuations do not appear to be informative for the correctness of the structure. However, both a simple hydrogen bond count and a simple compactness measure, if averaged over times of typically 300 ps, correlate well with structural correctness and we suggest that criteria based on these properties may be used in computational folding strategies.

Background Predicting protein residue-residue contacts is an important 2D prediction task. It is useful for ab initio structure prediction and understanding protein folding. In spite of steady progress over the past decade, contact prediction remains still largely unsolved. Results Here we develop a new contact map predictor (SVMcon) that uses support vector machines to predict medium- and long-range contacts. SVMcon integrates profiles, secondary structure, relative solvent accessibility, contact potentials, and other useful features. On the same test data set, SVMcon's accuracy is 4% higher than the latest version of the CMAPpro contact map predictor. SVMcon recently participated in the seventh edition of the Critical Assessment of Techniques for Protein Structure Prediction (CASP7) experiment and was evaluated along with seven other contact map predictors. SVMcon was ranked as one of the top predictors, yielding the second best coverage and accuracy for contacts with sequence separation >= 12 on 13 de novo domains. Conclusion We describe SVMcon, a new contact map predictor that uses SVMs and a large set of informative features. SVMcon yields good performance on medium- to long-range contact predictions and can be modularly incorporated into a structure prediction pipeline.