Comprehensive in silico mutagenesis highlights functionally important residues in proteins.
Ontology highlight
ABSTRACT: Motivation
Mutating residues into alanine (alanine scanning) is one of the fastest experimental means of probing hypotheses about protein function. Alanine scans can reveal functional hot spots, i.e. residues that alter function upon mutation. In vitro mutagenesis is cumbersome and costly: probing all residues in a protein is typically as impossible as substituting by all non-native amino acids. In contrast, such exhaustive mutagenesis is feasible in silico.Results
Previously, we developed SNAP to predict functional changes due to non-synonymous single nucleotide polymorphisms. Here, we applied SNAP to all experimental mutations in the ASEdb database of alanine scans; we identi.ed 70% of the hot spots (>or=1 kCal/mol change in binding energy); more severe changes were predicted more accurately. Encouraged, we carried out a complete all-against-all in silico mutagenesis for human glucokinase. Many of the residues predicted as functionally important have indeed been con.rmed in the literature, others await experimental veri.cation, and our method is ready to aid in the design of in vitro mutagenesis.Availability
ASEdb and glucokinase scores are available at http://www.rostlab.org/services/SNAP. For submissions of large/whole proteins for processing please contact the author.
SUBMITTER: Bromberg Y
PROVIDER: S-EPMC2597370 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The FireDB database is a databank for functional information relating to proteins with known structures. It contains the most comprehensive and detailed repository of known functionally important residues, bringing together both ligand binding and catalytic residues in one site. The platform integrates biologically relevant data filtered from the close atomic contacts in Protein Data Bank crystal structures and reliably annotated catalytic residues from the Catalytic Site Atlas. The interface allows users to make queries by protein, ligand or keyword. Relevant biologically important residues are displayed in a simple and easy to read manner that allows users to assess binding site similarity across homologous proteins. Binding site residue variations can also be viewed with molecular visualization tools. The database is available at http://firedb.bioinfo.cnio.es.
| S-EPMC1716728 | biostudies-literature
Project description:We analyzed functionality and relative distribution of genetic variants across the complete Oryza sativa genome, using the 40 million single nucleotide polymorphisms (SNPs) dataset from the 3,000 Rice Genomes Project (http://snp-seek.irri.org), the largest and highest density SNP collection for any higher plant. We have shown that the DNA-binding transcription factors (TFs) are the most conserved group of genes, whereas kinases and membrane-localized transporters are the most variable ones. TFs may be conserved because they belong to some of the most connected regulatory hubs that modulate transcription of vast downstream gene networks, whereas signaling kinases and transporters need to adapt rapidly to changing environmental conditions. In general, the observed profound patterns of nucleotide variability reveal functionally important genomic regions. As expected, nucleotide diversity is much higher in intergenic regions than within gene bodies (regions spanning gene models), and protein-coding sequences are more conserved than untranslated gene regions. We have observed a sharp decline in nucleotide diversity that begins at about 250 nucleotides upstream of the transcription start and reaches minimal diversity exactly at the transcription start. We found the transcription termination sites to have remarkably symmetrical patterns of SNP density, implying presence of functional sites near transcription termination. Also, nucleotide diversity was significantly lower near 3' UTRs, the area rich with regulatory regions.
| S-EPMC5075931 | biostudies-literature
Project description:BackgroundWell-performing automated protein function recognition approaches usually comprise several complementary techniques. Beside constructing better consensus, their predictive power can be improved by either adding or refining independent modules that explore orthogonal features of proteins. In this work, we demonstrated how the exploration of global atomic distributions can be used to indicate functionally important residues.ResultsUsing a set of carefully selected globular proteins, we parametrized continuous probability density functions describing preferred central distances of individual protein atoms. Relative preferred burials were estimated using mixture models of radial density functions dependent on the amino acid composition of a protein under consideration. The unexpectedness of extraordinary locations of atoms was evaluated in the information-theoretic manner and used directly for the identification of key amino acids. In the validation study, we tested capabilities of a tool built upon our approach, called SurpResi, by searching for binding sites interacting with ligands. The tool indicated multiple candidate sites achieving success rates comparable to several geometric methods. We also showed that the unexpectedness is a property of regions involved in protein-protein interactions, and thus can be used for the ranking of protein docking predictions. The computational approach implemented in this work is freely available via a Web interface at http://www.bioinformatics.org/surpresi.ConclusionsProbabilistic analysis of atomic central distances in globular proteins is capable of capturing distinct orientational preferences of amino acids as resulting from different sizes, charges and hydrophobic characters of their side chains. When idealized spatial preferences can be inferred from the sole amino acid composition of a protein, residues located in hydrophobically unfavorable environments can be easily detected. Such residues turn out to be often directly involved in binding ligands or interfacing with other proteins.
| S-EPMC3188475 | biostudies-literature
Project description:Lack of crystal structure data of folate binding proteins has left so many questions unanswered (for example, important residues in active site, binding domain, important amino acid residues involved in interactions between ligand and receptor). With sequence alignment and PROSITE motif identification, we attempted to answer evolutionarily significant residues that are of functional importance for ligand binding and that form catalytic sites. We have analyzed 46 different FRs and FBP sequences of various organisms obtained from Genbank. Multiple sequence alignment identified 44 highly conserved identical amino acid residues with 10 cysteine residues and 12 motifs including ECSPNLGPW (which might help in the structural stability of FR).
| S-EPMC2255074 | biostudies-literature
Project description:Crystallization of the zebrafish P2X4 receptor in both open and closed states revealed conformational differences in the ectodomain structures, including the dorsal fin and left flipper domains. Here, we focused on the role of these domains in receptor activation, responsiveness to orthosteric ATP analogue agonists, and desensitization. Alanine scanning mutagenesis of the R203-L214 (dorsal fin) and the D280-N293 (left flipper) sequences of the rat P2X4 receptor showed that ATP potency/efficacy was reduced in 15 out of 26 alanine mutants. The R203A, N204A, and N293A mutants were essentially non-functional, but receptor function was restored by ivermectin, an allosteric modulator. The I205A, T210A, L214A, P290A, G291A, and Y292A mutants exhibited significant changes in the responsiveness to orthosteric analog agonists 2-(methylthio)adenosine 5'-triphosphate, adenosine 5'-(?-thio)triphosphate, 2'(3'-O-(4-benzoylbenzoyl)adenosine 5'-triphosphate, and ?,?-methyleneadenosine 5'-triphosphate. In contrast, the responsiveness of L206A, N208A, D280A, T281A, R282A, and H286A mutants to analog agonists was comparable to that of the wild type receptor. Among these mutants, D280A, T281A, R282A, H286A, G291A, and Y292A also exhibited increased time-constant of the desensitizing current response. These experiments, together with homology modeling, indicate that residues located in the upper part of the dorsal fin and left flipper domains, relative to distance from the channel pore, contribute to the organization of the ATP binding pocket and to the initiation of signal transmission towards residues in the lower part of both domains. The R203 and N204 residues, deeply buried in the protein, may integrate the output signal from these two domains towards the gate. In addition, the left flipper residues predominantly account for the control of transition of channels from an open to a desensitized state.
| S-EPMC4232510 | biostudies-literature
Project description:Thrombopoietin (TPO) is a haematopoietic growth factor responsible for megakaryocyte progenitor proliferation and differentiation. It belongs to the four-helix-bundle cytokine family and exerts its biological effects through binding to a specific receptor, c-Mpl. With the use of site-directed mutagenesis we have generated 20 TPO mutants. Each of the TPO mutants was produced in a eukaryotic expression system and the mutants' ability to induce the proliferation of factor-dependent c-Mpl-expressing megakaryoblastic M-O7e cells was compared with that of wild-type TPO. Among the mutations studied, 10 lead to a significant decrease in TPO bioactivity. Of these ten residues, three are located in helix A of the protein (Arg10, Lys14 and Arg17) and four in helix D (His133, Gln132, Lys138 and Phe141), indicating that in TPO, as in other cytokines, these two helices are important for functional cytokine/receptor interactions. Surprisingly, mutant Arg10-->Ala (R10A) lacked any proliferative activity, despite the fact that this mutation was recently reported to have no effect on TPO/c-Mpl binding in a TPO phage ELISA [Pearce, Potts, Presta, Bald, Fendly and Wells (1997) J. Biol. Chem. 272, 20595-20602]. The lack of M-O7e proliferation is probably due to an inability of R10A mutant to promote receptor dimerization and thus receptor activation. Moreover we found that the Arg10 and Arg17 residues of TPO seem to be specific determinants for TPO/c-Mpl recognition. We also demonstrate that the O-glycosylation site located at position 110 of TPO is not necessary for the bioactivity of the cytokine.
| S-EPMC1219638 | biostudies-other
Project description:Here we present firestar, an expert system for predicting ligand-binding residues in protein structures. The server provides a method for extrapolating from the large inventory of functionally important residues organized in the FireDB database and adds information about the local conservation of potential-binding residues. The interface allows users to make queries by protein sequence or structure. The user can access pairwise and multiple alignments with structures that have relevant functionally important binding sites. The results are presented in a series of easy to read displays that allow users to compare binding residue conservation across homologous proteins. The binding site residues can also be viewed with molecular visualization tools. One feature of firestar is that it can be used to evaluate the biological relevance of small molecule ligands present in PDB structures. With the server it is easy to discern whether small molecule binding is conserved in homologous structures. We found this facility particularly useful during the recent assessment of CASP7 function prediction.http://firedb.bioinfo.cnio.es/Php/FireStar.php.
| S-EPMC1933227 | biostudies-literature
Project description:The rapidly increasing volume of sequence and structure information available for proteins poses the daunting task of determining their functional importance. Computational methods can prove to be very useful in understanding and characterizing the biochemical and evolutionary information contained in this wealth of data, particularly at functionally important sites. Therefore, we perform a detailed survey of compositional and evolutionary constraints at the molecular and biological function level for a large set of known functionally important sites extracted from a wide range of protein families. We compare the degree of conservation across different functional categories and provide detailed statistical insight to decipher the varying evolutionary constraints at functionally important sites. The compositional and evolutionary information at functionally important sites has been compiled into a library of functional templates. We developed a module that predicts functionally important columns (FIC) of an alignment based on the detection of a significant "template match score" to a library template. Our template match score measures an alignment column's similarity to a library template and combines a term explicitly representing a column's residue composition with various evolutionary conservation scores (information content and position-specific scoring matrix-derived statistics). Our benchmarking studies show good sensitivity/specificity for the prediction of functional sites and high accuracy in attributing correct molecular function type to the predicted sites. This prediction method is based on information derived from homologous sequences and no structural information is required. Therefore, this method could be extremely useful for large-scale functional annotation.
| S-EPMC2222836 | biostudies-literature
Project description:Further improvement of the thermostability of inherently thermostable proteins is an attractive challenge because more thermostable proteins are industrially more useful and serve as better scaffolds for protein engineering. To establish guidelines that can be applied for the rational design of hyperthermostable proteins, we compared the amino acid sequences of two ancestral nucleoside diphosphate kinases, Arc1 and Bac1, reconstructed in our previous study. Although Bac1 is a thermostable protein whose unfolding temperature is around 100°C, Arc1 is much more thermostable with an unfolding temperature of 114°C. However, only 12 out of 139 amino acids are different between the two sequences. In this study, one or a combination of amino acid(s) in Bac1 was/were substituted by a residue(s) found in Arc1 at the same position(s). The best mutant, which contained three amino acid substitutions (S108D, G116A and L120P substitutions), showed an unfolding temperature more than 10°C higher than that of Bac1. Furthermore, a combination of the other nine amino acid substitutions also led to improved thermostability of Bac1, although the effects of individual substitutions were small. Therefore, not only the sum of the contributions of individual amino acids, but also the synergistic effects of multiple amino acids are deeply involved in the stability of a hyperthermostable protein. Such insights will be helpful for future rational design of hyperthermostable proteins.
| S-EPMC8530338 | biostudies-literature
Project description:Permeases belonging to the equilibrative nucleoside transporter family promote uptake of nucleosides and/or nucleobases into a wide range of eukaryotes and mediate the uptake of a variety of drugs used in the treatment of cancer, heart disease, AIDS, and parasitic infections. No experimental three-dimensional structure exists for any of these permeases, and they are not present in prokaryotes, the source of many membrane proteins used in crystal structure determination. To generate a structural model for such a transporter, the LdNT1.1 nucleoside permease from the parasitic protozoan Leishmania donovani was modeled using ab initio computation. Site-directed mutations that strongly impair transport or that alter substrate specificity map to the central pore of the ab initio model, whereas mutations that have less pronounced phenotypes map to peripheral positions. The model suggests that aromatic residues present in transmembrane helices 1, 2, and 7 may interact to form an extracellular gate that closes the permeation pathway in the inward oriented conformation. Mutation of two of these three residues abrogated transport activity, consistent with the prediction of the model. The ab initio model is similar to one derived previously using threading analysis, a distinct computational approach, supporting the overall accuracy of both models. However, significant differences in helix orientation and residue position between the two models are apparent, and the mutagenesis data suggest that the ab initio model represents an improvement regarding structural details over the threading model. The putative gating interaction may also help explain differences in substrate specificity between members of this family.
| S-EPMC2707223 | biostudies-literature