LS-SNP: large-scale annotation of coding non-synonymous SNPs based on multiple information sources

MOTIVATION: The NCBI dbSNP database lists over 9 million single nucleotide polymorphisms (SNPs) in the human genome, but currently contains limited annotation information. SNPs that result in amino acid residue changes (nsSNPs) are of critical importance in variation between individuals, including disease and drug sensitivity. RESULTS: We have developed LS-SNP, a genomic scale software pipeline to annotate nsSNPs. LS-SNP comprehensively maps nsSNPs onto protein sequences, functional pathways and comparative protein structure models, and predicts positions where nsSNPs destabilize proteins, interfere with the formation of domain-domain interfaces, have an effect on protein-ligand binding or severely impact human health. It currently annotates 28,043 validated SNPs that produce amino acid residue substitutions in human proteins from the SwissProt/TrEMBL database. Annotations can be viewed via a web interface either in the context of a genomic region or by selecting sets of SNPs, genes, proteins or pathways. These results are useful for identifying candidate functional SNPs within a gene, haplotype or pathway and in probing molecular mechanisms responsible for functional impacts of nsSNPs. AVAILABILITY: http://www.salilab.org/LS-SNP CONTACT: rachelk@salilab.org SUPPLEMENTARY INFORMATION: http://salilab.org/LS-SNP/supp-info.pdf.     

A "solvated rotamer" approach to modeling water-mediated hydrogen bonds at protein-protein interfaces

Water-mediated hydrogen bonds play critical roles at protein-protein and protein-nucleic acid interfaces, and the interactions formed by discrete water molecules cannot be captured using continuum solvent models. We describe a simple model for the energetics of water-mediated hydrogen bonds, and show that, together with knowledge of the positions of buried water molecules observed in X-ray crystal structures, the model improves the prediction of free-energy changes upon mutation at protein-protein interfaces, and the recovery of native amino acid sequences in protein interface design calculations. We then describe a "solvated rotamer" approach to efficiently predict the positions of water molecules, at protein-protein interfaces and in monomeric proteins, that is compatible with widely used rotamer-based side-chain packing and protein design algorithms. Finally, we examine the extent to which the predicted water molecules can be used to improve prediction of amino acid identities and protein-protein interface stability, and discuss avenues for overcoming current limitations of the approach.     

A quantitative analysis of interfacial amino acid conservation in protein-protein hetero complexes

A long-standing question in molecular biology is whether interfaces of protein-protein complexes are more conserved than the rest of the protein surfaces. Although it has been reported that conservation can be used as an indicator for predicting interaction sites on proteins, there are recent reports stating that the interface regions are only slightly more conserved than the rest of the protein surfaces, with conservation signals not being statistically significant enough for predicting protein-protein binding sites. In order to properly address these controversial reports we have studied a set of 28 well resolved hetero complex structures of proteins that consists of transient and non-transient complexes. The surface positions were classified into four conservation classes and the conservation index of the surface positions was quantitatively analyzed. The results indicate that the surface density of highly conserved positions is significantly higher in the protein-protein interface regions compared with the other regions of the protein surface. However, the average conservation index of the patches in the interface region is not significantly higher compared with other surface regions of the protein structures. This finding demonstrates that the number of conserved residue positions is a more appropriate indicator for predicting protein-protein binding sites than the average conservation index in the interacting region. We have further validated our findings on a set of 59 benchmark complex structures. Furthermore, an analysis of 19 complexes of antigen-antibody interactions shows that there is no conservation of amino acid positions in the interacting regions of these complexes, as expected, with the variable region of the immunoglobulins interacting mostly with the antigens. Interestingly, antigen interacting regions also have a higher number of non-conserved residue positions in the interacting region than the rest of the protein surface.     

Predicting protein interaction sites: binding hot-spots in protein-protein and protein-ligand interfaces

MOTIVATION: Protein assemblies are currently poorly represented in structural databases and their structural elucidation is a key goal in biology. Here we analyse clefts in protein surfaces, likely to correspond to binding 'hot-spots', and rank them according to sequence conservation and simple measures of physical properties including hydrophobicity, desolvation, electrostatic and van der Waals potentials, to predict which are involved in binding in the native complex. RESULTS: The resulting differences between predicting binding-sites at protein-protein and protein-ligand interfaces are striking. There is a high level of prediction accuracy (< or =93%) for protein-ligand interactions, based on the following attributes: van der Waals potential, electrostatic potential, desolvation and surface conservation. Generally, the prediction accuracy for protein-protein interactions is lower, with the exception of enzymes. Our results show that the ease of cleft desolvation is strongly predictive of interfaces and strongly maintained across all classes of protein-binding interface.     

Mutagenic mapping of helical structures in the transmembrane segments of the yeast alpha-factor receptor

The alpha-mating pheromone receptor encoded by the yeast STE2 gene is a G protein coupled receptor that initiates signaling via a MAP kinase pathway that prepares haploid cells for mating. To establish the range of allowed amino acid substitutions within transmembrane segments of this receptor, we conducted extensive random mutagenesis of receptors followed by screening for receptor function. A total of 157 amino acid positions in seven different mutagenic libraries corresponding to the seven predicted transmembrane segments were analyzed, yielding 390 alleles that retain at least 60 % of normal signaling function. These alleles contained a total of 576 unique amino acid substitutions, including 61 % of all the possible amino acid changes that can arise from single base substitutions. The receptor exhibits a surprising tolerance for amino acid substitutions. Every amino acid in the mutagenized regions of the transmembrane regions could be substituted by at least one other residue. Polar amino acids were tolerated in functional receptors at 115 different positions (73 % of the total). Hydrophobic amino acids were tolerated in functional receptors at all mutagenized positions. Substitutions introducing proline residues were recovered at 53 % of all positions where they could be brought about by single base changes. Residues with charged side-chains could also be tolerated at 53 % of all positions where they were accessible through single base changes. The spectrum of allowed amino acid substitutions was characterized in terms of the hydrophobicity, radius of gyration, and charge of the allowed substitutions and mapped onto alpha-helical structures. By comparing the patterns of allowed substitutions with the recently determined structure of rhodopsin, structural features indicative of helix-helix interactions can be discerned in spite of the extreme sequence divergence between these two proteins.     

Probing the roles of conserved arginine-44 of Escherichia coli dihydrofolate reductase in its function and stability by systematic sequence perturbation analysis

Sequence perturbation analysis is a powerful method to reveal roles of an amino acid residue in function and stability of a protein. By using and improving this method, we studied roles of highly conserved Arg44 of Escherichia coli dihydrofolate reductase (DHFR) in its function and stability. Here, we introduced systematic amino acid substitutions at this position and found that all 19 kinds of amino acid substitutions were tolerated, but the mutations significantly reduced the enzymatic activity and the binding affinity toward the cofactor NADPH. Moreover, the mutational effects on the cofactor binding affinity were well correlated with those on the catalytic activity, indicating that the R44X mutations affect the catalytic activity mainly by modulating the cofactor binding affinity. On the other hand, thermal denaturation measurements showed that most mutations stabilized the protein. Comparison between the mutational effects and various amino acid indices taken from the AAindex database indicated that hydrophobicity and polarity are key determinants of amino acids favorable at this position. These results suggest that through electrostatic interactions Arg44 plays a functional role in retaining the cofactor binding affinity at the cost of the protein stability.     

Comparison of epitope structures of H3HAs through protein modeling of influenza A virus hemagglutinin: mechanism for selection of antigenic variants in the presence of a monoclonal antibody

Starting with nine plaques of influenza A/Kamata/14/91(H3N2) virus, we selected mutants in the presence of monoclonal antibody 203 (mAb203). In total, amino acid substitutions were found at nine positions (77, 80, 131, 135, 141, 142, 143, 144 and 146), which localized in the antigenic site A of the hemagglutinin (HA). The escape mutants differed in the extent to which they had lost binding to mAb203. HA protein with substitutions of some amino acid residues created by site-directed mutagenesis in the escape mutants retained the ability to bind to mAb203. Changes in the amino acid character affecting charge or hydrophobicity accounted for the binding capacity to the antibody of the HA with most of the substitutions in the escape mutants and binding-positive mutants. However, the effect of some amino acid substitutions remained unexplained. A three-dimensional model of the 1991 HA was constructed and used to analyze substituted amino acids in these mutants for the accessible surface hydrophobic and hydrophilic characters. One amino acid substitution in an escape mutant and another amino acid substitution in a binding-positive mutant seemed to be explained by the changes noted on this model.     

The influence of amino acid substitutions on the conformational energy of cytochrome c.

Conformational energies have been evaluated for each of the staggered side-chain conformations associated with the 261 amino acid substitutions known to occur among 60 eucaryotic species. At least 86% of these substitutions can be sterically accommodated (one at a time) within the structure of horse-heart cytochrome c resulting from conformational energy refinement. Simultaneous incorporation of all pertinent amino acid substitutions found in eight representative species into the refined horse-heart structure is also shown to be sterically possible, with few exceptions. In two cases (Pekin duck cytochrome with 10 substitutions and Samia cynthia cytochrome with 24 substitutions), all substitutions could be readily incorporated, and the total energies associated with their computed structures differed by less than 10 kcal/mol from that of horse-heart cytochrome c. In the cytochromes from rattlesnake (22 substitutions), tuna (18 substitutions), and Neurospora crassa (36 substitutions), tyrosine could not be substituted for phenylalanine at position 46, within the constraints of the calculations. However, when all of the remaining substitutions were incorporated into these three cytochromes, their computed conformational energies differed by less than 30 kcal/mol from that of horse-heart cytochrome c. Between two and four amino acid substitutions cause high energies in the cytochromes from human, baker's yeast, and cotton seed, but all of the remaining substitutions are consistent with a low energy conformation. These results suggest that the structures of homologous proteins may be even more similar than has previously been recognized. Substitutions of all possible amino acid types at the invariant positions (where all eucaryotic cytochromes c bear the same amino acid) have revealed some cases where different amino acids can be accommodated, thus demonstrating that the biological constraints on amino acid substitutions are often different from the purely steric constraints investigated in this work.     

Critical residues influence the affinity and selectivity of alpha-conotoxin MI for nicotinic acetylcholine receptors.

The mammalian skeletal muscle acetylcholine receptor contains two nonequivalent acetylcholine binding sites, one each at the alpha/delta and alpha/gamma subunit interfaces. Alpha-Conotoxin MI, a 14-amino acid competitive antagonist, binds at both interfaces but has approximately 10(4) higher affinity for the alpha/delta site. We performed an "alanine walk" to identify the residues in alpha-MI that contribute to this selective interaction with the alpha/delta site. Electrophysiological measurements with Xenopus oocytes expressing normal receptors or receptors lacking either the gamma or delta subunit were made to assay toxin-receptor interaction. Alanine substitutions in most amino acid positions had only modest effects on toxin potency at either binding site. However, substitutions in two positions, proline-6 and tyrosine-12, dramatically reduced toxin potency at the high-affinity alpha/delta site while having comparatively little effect on low-affinity alpha/gamma binding. When tyrosine-12 was replaced by alanine, the toxin's selectivity for the high-affinity site (relative to that for the low-affinity site) was reduced from 45,000- to 30-fold. A series of additional amino acid substitutions in this position showed that increasing side chain size/hydrophobicity increases toxin potency at the alpha/delta site without affecting alpha/gamma binding. In contrast, when tyrosine-12 is diiodinated, toxin binding is nearly irreversible at the alpha/delta site but also increases by approximately 500-fold at the alpha/gamma site. The effects of position 12 substitutions are accounted for almost entirely by changes in the rate of toxin dissociation from the high-affinity alpha/delta binding site.     

In-silico analysis of nonsynonymous genomic variants within CCM2 gene reaffirm the existence of dual cores within typical PTB domain

PurposeThe objective of this study is to validate the existence of dual cores within the typical phosphotyrosine binding (PTB) domain and to identify potentially damaging and pathogenic nonsynonymous coding single nuclear polymorphisms (nsSNPs) in the canonical PTB domain of the CCM2 gene that causes cerebral cavernous malformations (CCMs).MethodsThe nsSNPs within the coding sequence for PTB domain of human CCM2 gene, retrieved from exclusive database searches, were analyzed for their functional and structural impact using a series of bioinformatic tools. The effects of mutations on the tertiary structure of the PTB domain in human CCM2 protein were predicted to examine the effect of nsSNPs on the tertiary structure of PTB Cores.ResultsOur mutation analysis, through alignment of protein structures between wildtype CCM2 and mutant, predicted that the structural impacts of pathogenic nsSNPs is biophysically limited to only the spatially adjacent substituted amino acid site with minimal structural influence on the adjacent core of the PTB domain, suggesting both cores are independently functional and essential for proper CCM2 PTB function.ConclusionUtilizing a combination of protein conservation and structure-based analysis, we analyzed the structural effects of inherited pathogenic mutations within the CCM2 PTB domain. Our results predicted that the pathogenic amino acid substitutions lead to only subtle changes locally, confined to the surrounding tertiary structure of the PTB core within which it resides, while no structural disturbance to the neighboring PTB core was observed, reaffirming the presence of independently functional dual cores in the CCM2 typical PTB domain.     

Identification of the β-Lactamase Inhibitor Protein-II (BLIP-II) Interface Residues Essential for Binding Affinity and Specificity for Class A β-Lactamases

The interactions between β-lactamase inhibitory proteins (BLIPs) and β-lactamases have been used as model systems to understand the principles of affinity and specificity in protein-protein interactions. The most extensively studied tight binding inhibitor, BLIP, has been characterized with respect to amino acid determinants of affinity and specificity for binding β-lactamases. BLIP-II, however, shares no sequence or structural homology to BLIP and is a femtomolar to picomolar potency inhibitor, and the amino acid determinants of binding affinity and specificity are unknown. In this study, alanine scanning mutagenesis was used in combination with determinations of on and off rates for each mutant to define the contribution of residues on the BLIP-II binding surface to both affinity and specificity toward four β-lactamases of diverse sequence. The residues making the largest contribution to binding energy are heavily biased toward aromatic amino acids near the center of the binding surface. In addition, substitutions that reduce binding energy do so by increasing off rates without impacting on rates. Also, residues with large contributions to binding energy generally exhibit low temperature factors in the structures of complexes. Finally, with the exception of D206A, BLIP-II alanine substitutions exhibit a similar trend of effect for all β-lactamases, i.e., a substitution that reduces affinity for one β-lactamase usually reduces affinity for all β-lactamases tested.     

Computer-assisted mutagenesis of ecotin to engineer its secondary binding site for urokinase inhibition

Inhibitors of urokinase-type plasminogen activator (uPA) were selected in vitro from two ecotin phage-display libraries to study the effect on binding of amino acid substitutions at critical positions 108, 110, 112, and 113 within the 100s loop (RNKL, respectively, in wild type ecotin). The first, a focused library, was the result of a computation-assisted approach using the three-dimensional structure of the ecotin-trypsin complex to guide the modeling of amino acid substitutions predicted to increase affinity for uPA. The second, a complete library, allowed for all substitutions at the above identified positions. The consensus sequences selected from the focused, and complete libraries were RRWS and R(R/N)QL, respectively. Inhibition constant determinations showed ecotin variants containing these sequences to be similarly potent (K(i) = 1-2 nm). These substitutions were combined with previously identified substitutions in another critical region of ecotin. One of these combinations (D70R/M84R/RRQL) is the tightest (K(i) = 50 pm) ecotin variant inhibitor of uPA. The blending of combinatorial methods and computer algorithms designed to predict stronger binders has allowed us to obtain protein derivatives that exhibit greatly increased affinity for a predetermined target. This technology can be applied to select for enhanced binding interactions at protein-protein interfaces and accelerate the process of protease inhibitor development.     

A protein-protein docking algorithm dependent on the type of complexes

An efficient 'soft docking' algorithm is described to assist the prediction of protein-protein association using three-dimensional structures of molecules. The basic tools are the 'simplified protein' model and the docking algorithm of Wodak and Janin. The side chain flexibility of Arg, Lys, Asp, Glu and Met residues at the protein surface is taken into account. The complex type-dependent filtering technique on the basis of the geometric matching, hydrophobicity and electrostatic complementarity is used to select candidate binding modes. Subsequently, we calculate a scoring function which includes electrostatic and desolvation energy terms. In the 44 complexes tested including enzyme-inhibitor, antibody-antigen and other complexes, native-like structures were all found, of which 30 were ranked in the top 20. Thus, our soft docking algorithm has the potential to predict protein-protein recognition.     

A structural bioinformatics approach to the analysis of nonsynonymous single nucleotide polymorphisms (nsSNPs) and their relation to disease

The prediction of the effects of nonsynonymous single nucleotide polymorphisms (nsSNPs) on function depends critically on exploiting all information available on the three-dimensional structures of proteins. We describe software and databases for the analysis of nsSNPs that allow a user to move from SNP to sequence to structure to function. In both structure prediction and the analysis of the effects of nsSNPs, we exploit information about protein evolution, in particular, that derived from investigations on the relation of sequence to structure gained from the study of amino acid substitutions in divergent evolution. The techniques developed in our laboratory have allowed fast and automated sequence-structure homology recognition to identify templates and to perform comparative modeling; as well as simple, robust, and generally applicable algorithms to assess the likely impact of amino acid substitutions on structure and interactions. We describe our strategy for approaching the relationship between SNPs and disease, and the results of benchmarking our approach -- human proteins of known structure and recognized mutation.     

Computational Mapping of Anchoring Spots on Protein Surfaces

Protein-protein and protein-peptide interactions are often controlled by few strong contacts that involve hot spot residues. Computational detection of such contacts, termed here anchoring spots, is important for understanding recognition processes and for predicting interactions; it is an essential step in designing interaction interfaces and therapeutic agents. We describe ANCHORSMAP, an algorithm for computational mapping of amino acid side chains on protein surfaces. The algorithm consists of two stages: A geometry based stage (LSMdet), in which sub-pockets adequate for binding single side chains are detected and amino acid probes are scattered near them, and an energy based stage in which optimal positions of the probes are determined through repeated energy minimization and clustering of nearby poses and their ΔG are calculated. ANCHORSMAP employs a new function for ΔG calculations, which is specifically designed for the context of protein-protein recognition by introducing a correction in the electrostatic energy term that compensates for the dielectric shielding exerted by a hypothetical protein bound to the probe.The algorithm successfully detects known anchoring sites and accurately positions the probes. The calculated ΔG rank high the correct anchoring spots in maps produced for unbound proteins. We find that Arg, Trp, Glu and Tyr, which are favorite hot spot residues, are also more selective of their binding environment. The usefulness of anchoring spots mapping is demonstrated by detecting the binding surfaces in the protein-protein complex barnase/barstar and the protein-peptide complex kinase/PKI, and by identifying phenylalanine anchoring sites on the surface of the nuclear transporter NTF2, C-terminus anchors on PDZ domains and phenol anchors on thermolysin. Finally, we discuss the role of anchoring spots in molecular recognition processes.     

Comprehensive and quantitative mapping of energy landscapes for protein-protein interactions by rapid combinatorial scanning

A novel, quantitative saturation (QS) scanning strategy was developed to obtain a comprehensive data base of the structural and functional effects of all possible mutations across a large protein-protein interface. The QS scan approach was applied to the high affinity site of human growth hormone (hGH) for binding to its receptor (hGHR). Although the published structure-function data base describing this system is probably the most extensive for any large protein-protein interface, it is nonetheless too sparse to accurately describe the nature of the energetics governing the interaction. Our comprehensive data base affords a complete view of the binding site and provides important new insights into the general principles underlying protein-protein interactions. The hGH binding interface is highly adaptable to mutations, but the nature of the tolerated mutations challenges generally accepted views about the evolutionary and biophysical pressures governing protein-protein interactions. Many substitutions that would be considered chemically conservative are not tolerated, while conversely, many non-conservative substitutions can be accommodated. Furthermore, conservation across species is a poor predictor of the chemical character of tolerated substitutions across the interface. Numerous deviations from generally accepted expectations indicate that mutational tolerance is highly context dependent and, furthermore, cannot be predicted by our current knowledge base. The type of data produced by the comprehensive QS scan can fill the gaps in the structure-function matrix. The compilation of analogous data bases from studies of other protein-protein interactions should greatly aid the development of computational methods for explaining and designing molecular recognition.     

Antigen-antibody interaction. Synthetic peptides define linear antigenic determinants recognized by monoclonal antibodies directed to the cytoplasmic carboxyl terminus of rhodopsin

The specificities of four monoclonal antibodies rho 1D4, 1C5, 3A6, and 3D6 prepared by immunization of rod outer segments containing rhodopsin have been defined using synthetic peptides. All of these antibodies interact within the 18 residues at the COOH terminus of rhodopsin and recognize linear antigenic determinants of 4-11 residues. Twenty-seven synthetic peptide analogs of varying lengths of native sequence or containing single amino acid substitutions at each position of the COOH-terminal 18 residues have provided some insight into the mechanism of antigen-antibody binding. Our results clearly demonstrate that antibodies can be highly specific at key positions as shown by the loss of binding on single amino acid substitutions in the binding site. In contrast single amino acid substitutions at other positions in the binding site only affect affinity for some antibodies. Ionic interactions can dominate immunogenic determinants. Immunogenic determinants are not restricted to highly charged hydrophilic regions on the surface of a protein and may be dominated by hydrophobic interactions. Although certain side chains can dominate the interaction of the antigen with antibody, our results are in agreement with the interpretation that the free energies of all the contact points are additive and a certain free energy must be present to achieve binding. Antibodies with different specificities directed to the same region of the protein antigen can be produced in an immune response. Peptide antigens representing regions of a protein antigen bind best to the anti-protein antibody when the sequence is shortened to contain only those residues binding to the specificity site in the antibody. Cross-reactivity between protein antigens can be explained by conservation of the critical residues in the combining site.