ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins

A major problem in structural biology is the recognition of errors in experimental and theoretical models of protein structures. The ProSA program (Protein Structure Analysis) is an established tool which has a large user base and is frequently employed in the refinement and validation of experimental protein structures and in structure prediction and modeling. The analysis of protein structures is generally a difficult and cumbersome exercise. The new service presented here is a straightforward and easy to use extension of the classic ProSA program which exploits the advantages of interactive web-based applications for the display of scores and energy plots that highlight potential problems spotted in protein structures. In particular, the quality scores of a protein are displayed in the context of all known protein structures and problematic parts of a structure are shown and highlighted in a 3D molecule viewer. The service specifically addresses the needs encountered in the validation of protein structures obtained from X-ray analysis, NMR spectroscopy and theoretical calculations. ProSA-web is accessible at https://prosa.services.came.sbg.ac.a

Comparative analysis of homology models of the Ah receptor ligand binding domain: Verification of structure-function predictions by site-directed mutagenesis of a nonfunctional receptor

The aryl hydrocarbon receptor (AHR) is a ligand-dependent transcription factor that mediates the biological and toxic effects of a wide variety of structurally diverse chemicals, including the toxic environmental contaminant 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). While significant interspecies differences in AHR ligand binding specificity, selectivity, and response have been observed, the structural determinants responsible for those differences have not been determined, and homology models of the AHR ligand-binding domain (LBD) are available for only a few species. Here we describe the development and comparative analysis of homology models of the LBD of 16 AHRs from 12 mammalian and nonmammalian species and identify the specific residues contained within their ligand binding cavities. The ligand-binding cavity of the fish AHR exhibits differences from those of mammalian and avian AHRs, suggesting a slightly different TCDD binding mode. Comparison of the internal cavity in the LBD model of zebrafish (zf) AHR2, which binds TCDD with high affinity, to that of zfAHR1a, which does not bind TCDD, revealed that the latter has a dramatically shortened binding cavity due to the side chains of three residues (Tyr296, Thr386, and His388) that reduce the amount of internal space available to TCDD. Mutagenesis of two of these residues in zfAHR1a to those present in zfAHR2 (Y296H and T386A) restored the ability of zfAHR1a to bind TCDD and to exhibit TCDD-dependent binding to DNA. These results demonstrate the importance of these two amino acids and highlight the predictive potential of comparative analysis of homology models from diverse species. The availability of these AHR LBD homology models will facilitate in-depth comparative studies of AHR ligand binding and ligand-dependent AHR activation and provide a novel avenue for examining species-specific differences in AHR responsiveness. © 2013 American Chemical Society

The relative orientation of the TM3 and TM4 domains varies between α1 and α3 glycine receptors

Glycine receptors (GlyRs) are anion-conducting members of the pentameric ligand-gated ion channel family. We previously showed that the dramatic difference in glycine efficacies of α1 and α3 GlyRs is largely attributable to their nonconserved TM4 domains. Because mutation of individual nonconserved TM4 residues had little effect, we concluded that the efficacy difference was a distributed effect of all nonconserved TM4 residues. We therefore hypothesized that the TM4 domains of α1 and α3 GlyRs differ in structure, membrane orientation, and/or molecular dynamic properties. Here we employed voltage-clamp fluorometry to test whether their TM4 domains interact differently with their respective TM3 domains. We found a rhodamine fluorophore covalently attached to a homologous TM4 residue in each receptor interacts differentially with a conserved TM3 residue. We conclude that the α1 and α3 GlyR TM4 domains are orientated differently relative to their TM3 domains. This may underlie their differential ability to influence glycine efficacy

PRIMO: an interactive homology modeling pipeline

The development of automated servers to predict the three-dimensional structure of proteins has seen much progress over the years. These servers make calculations simpler, but largely exclude users from the process. In this study, we present the PRotein Interactive MOdeling (PRIMO) pipeline for homology modeling of protein monomers. The pipeline eases the multi-step modeling process, and reduces the workload required by the user, while still allowing engagement from the user during every step. Default parameters are given for each step, which can either be modified or supplemented with additional external input. PRIMO has been designed for users of varying levels of experience with homology modeling. The pipeline incorporates a user-friendly interface that makes it easy to alter parameters used during modeling

Binding mode analyses and pharmacophore model development for stilbene derivatives as a novel and competitive class of α-glucosidase inhibitors

Stilbene urea derivatives as a novel and competitive class of non-glycosidic α-glucosidase inhibitors are effective for the treatment of type II diabetes and obesity. The main purposes of our molecular modeling study are to explore the most suitable binding poses of stilbene derivatives with analyzing the binding affinity differences and finally to develop a pharmacophore model which would represents critical features responsible for α-glucosidase inhibitory activity. Three-dimensional structure of S. cerevisiae α-glucosidase was built by homology modeling method and the structure was used for the molecular docking study to find out the initial binding mode of compound 12, which is the most highly active one. The initial structure was subjected to molecular dynamics (MD) simulations for protein structure adjustment at compound 12-bound state. Based on the adjusted conformation, the more reasonable binding modes of the stilbene urea derivatives were obtained from molecular docking and MD simulations. The binding mode of the derivatives was validated by correlation analysis between experimental Ki value and interaction energy. Our results revealed that the binding modes of the potent inhibitors were engaged with important hydrogen bond, hydrophobic, and π-interactions. With the validated compound 12-bound structure obtained from combining approach of docking and MD simulation, a proper four featured pharmacophore model was generated. It was also validated by comparison of fit values with the Ki values. Thus, these results will be helpful for understanding the relationship between binding mode and bioactivity and for designing better inhibitors from stilbene derivatives

The Presence of the Iron-Sulfur Motif Is Important for the Conformational Stability of the Antiviral Protein, Viperin

Viperin, an antiviral protein, has been shown to contain a CX3CX2C motif, which is conserved in the radical S-adenosyl-methionine (SAM) enzyme family. A triple mutant which replaces these three cysteines with alanines has been shown to have severe deficiency in antiviral activity. Since the crystal structure of Viperin is not available, we have used a combination of computational methods including multi-template homology modeling and molecular dynamics simulation to develop a low-resolution predicted structure. The results show that Viperin is an α -β protein containing iron-sulfur cluster at the center pocket. The calculations suggest that the removal of iron-sulfur cluster would lead to collapse of the protein tertiary structure. To verify these predictions, we have prepared, expressed and purified four mutant proteins. In three mutants individual cysteine residues were replaced by alanine residues while in the fourth all the cysteines were replaced by alanines. Conformational analyses using circular dichroism and steady state fluorescence spectroscopy indicate that the mutant proteins are partially unfolded, conformationally unstable and aggregation prone. The lack of conformational stability of the mutant proteins may have direct relevance to the absence of their antiviral activity

A Score of the Ability of a Three-Dimensional Protein Model to Retrieve Its Own Sequence as a Quantitative Measure of Its Quality and Appropriateness

BACKGROUND: Despite the remarkable progress of bioinformatics, how the primary structure of a protein leads to a three-dimensional fold, and in turn determines its function remains an elusive question. Alignments of sequences with known function can be used to identify proteins with the same or similar function with high success. However, identification of function-related and structure-related amino acid positions is only possible after a detailed study of every protein. Folding pattern diversity seems to be much narrower than sequence diversity, and the amino acid sequences of natural proteins have evolved under a selective pressure comprising structural and functional requirements acting in parallel. PRINCIPAL FINDINGS: The approach described in this work begins by generating a large number of amino acid sequences using ROSETTA [Dantas G et al. (2003) J Mol Biol 332:449-460], a program with notable robustness in the assignment of amino acids to a known three-dimensional structure. The resulting sequence-sets showed no conservation of amino acids at active sites, or protein-protein interfaces. Hidden Markov models built from the resulting sequence sets were used to search sequence databases. Surprisingly, the models retrieved from the database sequences belonged to proteins with the same or a very similar function. Given an appropriate cutoff, the rate of false positives was zero. According to our results, this protocol, here referred to as Rd.HMM, detects fine structural details on the folding patterns, that seem to be tightly linked to the fitness of a structural framework for a specific biological function. CONCLUSION: Because the sequence of the native protein used to create the Rd.HMM model was always amongst the top hits, the procedure is a reliable tool to score, very accurately, the quality and appropriateness of computer-modeled 3D-structures, without the need for spectroscopy data. However, Rd.HMM is very sensitive to the conformational features of the models' backbone

Maintaining and breaking symmetry in homomeric coiled-coil assemblies

Higher order coiled coils with five or more helices can form α-helical barrels. Here the authors show that placing β-branched aliphatic residues along the lumen yields stable and open α-helical barrels, which is of interest for the rational design of functional proteins; whereas, the absence of β-branched side chains leads to unusual low-symmetry α-helical bundles