Analysis of substructural variation in families of enzymatic proteins with applications to protein function prediction

Background  

Structural variations caused by a wide range of physico-chemical and biological sources directly influence the function of a protein. For enzymatic proteins, the structure and chemistry of the catalytic binding site residues can be loosely defined as a substructure of the protein. Comparative analysis of drug-receptor substructures across and within species has been used for lead evaluation. Substructure-level similarity between the binding sites of functionally similar proteins has also been used to identify instances of convergent evolution among proteins. In functionally homologous protein families, shared chemistry and geometry at catalytic sites provide a common, local point of comparison among proteins that may differ significantly at the sequence, fold, or domain topology levels.     

Exploiting residue-level and profile-level interface propensities for usage in binding sites prediction of proteins

Background  

Recognition of binding sites in proteins is a direct computational approach to the characterization of proteins in terms of biological and biochemical function. Residue preferences have been widely used in many studies but the results are often not satisfactory. Although different amino acid compositions among the interaction sites of different complexes have been observed, such differences have not been integrated into the prediction process. Furthermore, the evolution information has not been exploited to achieve a more powerful propensity.     

Sequence and structural features of binding site residues in protein-protein complexes: comparison with protein-nucleic acid complexes

Background

Protein-protein interactions are important for several cellular processes. Understanding the mechanism of protein-protein recognition and predicting the binding sites in protein-protein complexes are long standing goals in molecular and computational biology.

Methods

We have developed an energy based approach for identifying the binding site residues in protein–protein complexes. The binding site residues have been analyzed with sequence and structure based parameters such as binding propensity, neighboring residues in the vicinity of binding sites, conservation score and conformational switching.

Results

We observed that the binding propensities of amino acid residues are specific for protein-protein complexes. Further, typical dipeptides and tripeptides showed high preference for binding, which is unique to protein-protein complexes. Most of the binding site residues are highly conserved among homologous sequences. Our analysis showed that 7% of residues changed their conformations upon protein-protein complex formation and it is 9.2% and 6.6% in the binding and non-binding sites, respectively. Specifically, the residues Glu, Lys, Leu and Ser changed their conformation from coil to helix/strand and from helix to coil/strand. Leu, Ser, Thr and Val prefer to change their conformation from strand to coil/helix.

Conclusions

The results obtained in this study will be helpful for understanding and predicting the binding sites in protein-protein complexes.

     

Structural and Functional Similarity between the Bacterial Type III Secretion System Needle Protein PrgI and the Eukaryotic Apoptosis Bcl-2 Proteins

Background

Functional similarity is challenging to identify when global sequence and structure similarity is low. Active-sites or functionally relevant regions are evolutionarily more stable relative to the remainder of a protein structure and provide an alternative means to identify potential functional similarity between proteins. We recently developed the FAST-NMR methodology to discover biochemical functions or functional hypotheses of proteins of unknown function by experimentally identifying ligand binding sites. FAST-NMR utilizes our CPASS software and database to assign a function based on a similarity in the structure and sequence of ligand binding sites between proteins of known and unknown function.

Methodology/Principal Findings

The PrgI protein from Salmonella typhimurium forms the needle complex in the type III secretion system (T3SS). A FAST-NMR screen identified a similarity between the ligand binding sites of PrgI and the Bcl-2 apoptosis protein Bcl-xL. These ligand binding sites correlate with known protein-protein binding interfaces required for oligomerization. Both proteins form membrane pores through this oligomerization to release effector proteins to stimulate cell death. Structural analysis indicates an overlap between the PrgI structure and the pore forming motif of Bcl-xL. A sequence alignment indicates conservation between the PrgI and Bcl-xL ligand binding sites and pore formation regions. This active-site similarity was then used to verify that chelerythrine, a known Bcl-xL inhibitor, also binds PrgI.

Conclusions/Significance

A structural and functional relationship between the bacterial T3SS and eukaryotic apoptosis was identified using our FAST-NMR ligand affinity screen in combination with a bioinformatic analysis based on our CPASS program. A similarity between PrgI and Bcl-xL is not readily apparent using traditional global sequence and structure analysis, but was only identified because of conservation in ligand binding sites. These results demonstrate the unique opportunity that ligand-binding sites provide for the identification of functional relationships when global sequence and structural information is limited.     

PCA-MutPred: Prediction of Binding Free Energy Change Upon Missense Mutation in Protein-carbohydrate Complexes

Protein-carbohydrate interactions play an important role in several biological processes. The mutation of amino acid residues in carbohydrate-binding proteins may alter the binding affinity, affect the functions and lead to diseases. Elucidating the factors influencing the binding affinity change (ΔΔG) of protein-carbohydrate complexes upon mutation is a challenging task. In this work, we have collected the experimental data for the binding affinity change of 318 unique mutants and related with sequence and structural features of amino acid residues at the mutant sites. We found that accessible surface area, secondary structure, mutation preference, conservation score, hydrophobicity and contact energies are important to understand the binding affinity change upon mutation. We have developed multiple regression equations for predicting the binding affinity change upon mutation and our method showed an average correlation of 0.74 and a mean absolute error of 0.70 kcal/mol between experimental and predicted ΔΔG on a 10-fold cross-validation. Further, we have validated our method using an independent test data set of 124 (62 unique) mutations, which showed a correlation and MAE of 0.79 and 0.56 kcal/mol, respectively. We have developed a web server PCA-MutPred, Protein-CArbohydrate complex Mutation affinity Predictor, for predicting the change in binding affinity of protein–carbohydrate complexes and it is freely accessible at https://web.iitm.ac.in/bioinfo2/pcamutpred. We suggest that the method could be a useful resource for designing protein-carbohydrate complexes with desired affinities.     

Domain-based small molecule binding site annotation

Background  

Accurate small molecule binding site information for a protein can facilitate studies in drug docking, drug discovery and function prediction, but small molecule binding site protein sequence annotation is sparse. The Small Molecule Interaction Database (SMID), a database of protein domain-small molecule interactions, was created using structural data from the Protein Data Bank (PDB). More importantly it provides a means to predict small molecule binding sites on proteins with a known or unknown structure and unlike prior approaches, removes large numbers of false positive hits arising from transitive alignment errors, non-biologically significant small molecules and crystallographic conditions that overpredict ion binding sites.     

Sequence-based identification of interface residues by an integrative profile combining hydrophobic and evolutionary information

Background  

Protein-protein interactions play essential roles in protein function determination and drug design. Numerous methods have been proposed to recognize their interaction sites, however, only a small proportion of protein complexes have been successfully resolved due to the high cost. Therefore, it is important to improve the performance for predicting protein interaction sites based on primary sequence alone.     

Search for fucose binding domains in recently sequenced hypothetical proteins using molecular modeling techniques and structural analysis

The crystal structure of a fucose-binding lectin from the bacteria Pseudomonas aeruginosa in complex with α-L-fucose has been recently determined. It is a tetramer; each monomer displays a nine-stranded, antiparallel, β-sandwiched arrangement and contains two calcium ions that mediate the binding of fucose in a recognition mode unique among protein-carbohydrate interactions. In search of this type of unique interactions in other newly discovered protein sequences, we have used molecular modeling techniques to predict and analyze the 3-D structures of some proteins, which exhibited reasonable degree of homology with the amino acid sequence of the bacterial protein. A BLAST search with the sequence of Pseudomonas aeruginosa as query in the non-redundant sequence database identified four proteins from different species, three organisms from bacteria and one from archaea. We have modeled the structures of these proteins as well as those of the complexes with carbohydrates and studied the nature of physicochemical forces involved in the complex formation both in presence and absence of calcium. The calcium-binding loops have been found to be highly conserved both in terms of primary and tertiary structures in these proteins, although a less acidic character is observed in Photorhabdus lectin due to the absence of two aspartic acid residues on the calcium-binding loop which also resulted in lower binding affinity. All these structures exhibited highly negative electrostatic environment in the vicinity of the calcium-binding loops which was essential for neutralizing the positive charges of two closely situated Ca⁺² ions. The comparison of the binding affinities of some monosaccharides other than fucose, e.g. mannose and fructose, showed higher binding energies confirming the fucose specificity of these proteins.     

A new protein binding pocket similarity measure based on comparison of clouds of atoms in 3D: application to ligand prediction

Background  

Predicting which molecules can bind to a given binding site of a protein with known 3D structure is important to decipher the protein function, and useful in drug design. A classical assumption in structural biology is that proteins with similar 3D structures have related molecular functions, and therefore may bind similar ligands. However, proteins that do not display any overall sequence or structure similarity may also bind similar ligands if they contain similar binding sites. Quantitatively assessing the similarity between binding sites may therefore be useful to propose new ligands for a given pocket, based on those known for similar pockets.     

Prediction of "hot spots" of aggregation in disease-linked polypeptides

Background  

The polypeptides involved in amyloidogenesis may be globular proteins with a defined 3D-structure or natively unfolded proteins. The first class includes polypeptides such as β2-microglobulin, lysozyme, transthyretin or the prion protein, whereas β-amyloid peptide, amylin or α-synuclein all belong to the second class. Recent studies suggest that specific regions in the proteins act as "hot spots" driving aggregation. This should be especially relevant for natively unfolded proteins or unfolded states of globular proteins as they lack significant secondary and tertiary structure and specific intra-chain interactions that can mask these aggregation-prone regions. Prediction of such sequence stretches is important since they are potential therapeutic targets.     

Identification of NAD interacting residues in proteins

Background  

Small molecular cofactors or ligands play a crucial role in the proper functioning of cells. Accurate annotation of their target proteins and binding sites is required for the complete understanding of reaction mechanisms. Nicotinamide adenine dinucleotide (NAD⁺ or NAD) is one of the most commonly used organic cofactors in living cells, which plays a critical role in cellular metabolism, storage and regulatory processes. In the past, several NAD binding proteins (NADBP) have been reported in the literature, which are responsible for a wide-range of activities in the cell. Attempts have been made to derive a rule for the binding of NAD⁺ to its target proteins. However, so far an efficient model could not be derived due to the time consuming process of structure determination, and limitations of similarity based approaches. Thus a sequence and non-similarity based method is needed to characterize the NAD binding sites to help in the annotation. In this study attempts have been made to predict NAD binding proteins and their interacting residues (NIRs) from amino acid sequence using bioinformatics tools.     

Markov clustering versus affinity propagation for the partitioning of protein interaction graphs

Background  

Genome scale data on protein interactions are generally represented as large networks, or graphs, where hundreds or thousands of proteins are linked to one another. Since proteins tend to function in groups, or complexes, an important goal has been to reliably identify protein complexes from these graphs. This task is commonly executed using clustering procedures, which aim at detecting densely connected regions within the interaction graphs. There exists a wealth of clustering algorithms, some of which have been applied to this problem. One of the most successful clustering procedures in this context has been the Markov Cluster algorithm (MCL), which was recently shown to outperform a number of other procedures, some of which were specifically designed for partitioning protein interactions graphs. A novel promising clustering procedure termed Affinity Propagation (AP) was recently shown to be particularly effective, and much faster than other methods for a variety of problems, but has not yet been applied to partition protein interaction graphs.     

No statistical support for correlation between the positions of protein interaction sites and alternatively spliced regions

Background  

Alternative splicing is an efficient mechanism for increasing the variety of functions fulfilled by proteins in a living cell. It has been previously demonstrated that alternatively spliced regions often comprise functionally important and conserved sequence motifs. The objective of this work was to test the hypothesis that alternative splicing is correlated with contact regions of protein-protein interactions.     

Predicting DNA-binding locations and orientation on proteins using knowledge-based learning of geometric properties

Background

DNA-binding proteins perform their functions through specific or non-specific sequence recognition. Although many sequence- or structure-based approaches have been proposed to identify DNA-binding residues on proteins or protein-binding sites on DNA sequences with satisfied performance, it remains a challenging task to unveil the exact mechanism of protein-DNA interactions without crystal complex structures. Without information from complexes, the linkages between DNA-binding proteins and their binding sites on DNA are still missing.

Methods

While it is still difficult to acquire co-crystallized structures in an efficient way, this study proposes a knowledge-based learning method to effectively predict DNA orientation and base locations around the protein’s DNA-binding sites when given a protein structure. First, the functionally important residues of a query protein are predicted by a sequential pattern mining tool. After that, surface residues falling in the predicted functional regions are determined based on the given structure. These residues are then clustered based on their spatial coordinates and the resultant clusters are ranked by a proposed DNA-binding propensity function. Clusters with high DNA-binding propensities are treated as DNA-binding units (DBUs) and each DBU is analyzed by principal component analysis (PCA) to predict potential orientation of DNA grooves. More specifically, the proposed method is developed to predict the direction of the tangent line to the helix curve of the DNA groove where a DBU is going to bind.

Results

This paper proposes a knowledge-based learning procedure to determine the spatial location of the DNA groove with respect to the query protein structure by considering geometric propensity between protein side chains and DNA bases. The 11 test cases used in this study reveal that the location and orientation of the DNA groove around a selected DBU can be predicted with satisfied errors.

Conclusions

This study presents a method to predict the location and orientation of DNA grooves with respect to the structure of a DNA-binding protein. The test cases shown in this study reveal the possibility of imaging protein-DNA binding conformation before co-crystallized structure can be determined. How the proposed method can be incorporated with existing protein-DNA docking tools to study protein-DNA interactions deserve further studies in the near future.

     

Can Clustal-style progressive pairwise alignment of multiple sequences be used in RNA secondary structure prediction?

Background  

In ribonucleic acid (RNA) molecules whose function depends on their final, folded three-dimensional shape (such as those in ribosomes or spliceosome complexes), the secondary structure, defined by the set of internal basepair interactions, is more consistently conserved than the primary structure, defined by the sequence of nucleotides.     

A Peek into Tropomyosin Binding and Unfolding on the Actin Filament

Background

Tropomyosin is a prototypical coiled coil along its length with subtle variations in structure that allow interactions with actin and other proteins. Actin binding globally stabilizes tropomyosin. Tropomyosin-actin interaction occurs periodically along the length of tropomyosin. However, it is not well understood how tropomyosin binds actin.

Principal Findings

Tropomyosin''s periodic binding sites make differential contributions to two components of actin binding, cooperativity and affinity, and can be classified as primary or secondary sites. We show through mutagenesis and analysis of recombinant striated muscle α-tropomyosins that primary actin binding sites have a destabilizing coiled-coil interface, typically alanine-rich, embedded within a non-interface recognition sequence. Introduction of an Ala cluster in place of the native, more stable interface in period 2 and/or period 3 sites (of seven) increased the affinity or cooperativity of actin binding, analysed by cosedimentation and differential scanning calorimetry. Replacement of period 3 with period 5 sequence, an unstable region of known importance for cooperative actin binding, increased the cooperativity of binding. Introduction of the fluorescent probe, pyrene, near the mutation sites in periods 2 and 3 reported local instability, stabilization by actin binding, and local unfolding before or coincident with dissociation from actin (measured using light scattering), and chain dissociation (analyzed using circular dichroism).

Conclusions

This, and previous work, suggests that regions of tropomyosin involved in binding actin have non-interface residues specific for interaction with actin and an unstable interface that is locally stabilized upon binding. The destabilized interface allows residues on the coiled-coil surface to obtain an optimal conformation for interaction with actin by increasing the number of local substates that the side chains can sample. We suggest that local disorder is a property typical of coiled coil binding sites and proteins that have multiple binding partners, of which tropomyosin is one type.     

SUPERFICIAL – Surface mapping of proteins via structure-based peptide library design

Background  

The determination of protein surfaces and the detection of binding sites are essential to our understanding of protein-protein interactions. Such binding sites can be characterised as linear and non-linear, the non-linear sites being prevailant. Conventional mapping techniques with arrays of synthetic peptides have limitations with regard to the location of discontinuous or non-linear binding sites of proteins.