Directionality in protein fold prediction

Background  

Ever since the ground-breaking work of Anfinsen et al. in which a denatured protein was found to refold to its native state, it has been frequently stated by the protein fold prediction community that all the information required for protein folding lies in the amino acid sequence. Recent in vitro experiments and in silico computational studies, however, have shown that cotranslation may affect the folding pathway of some proteins, especially those of ancient folds. In this paper aspects of cotranslational folding have been incorporated into a protein structure prediction algorithm by adapting the Rosetta program to fold proteins as the nascent chain elongates. This makes it possible to conduct a pairwise comparison of folding accuracy, by comparing folds created sequentially from each end of the protein.     

Implications from a network-based topological analysis of ubiquitin unfolding simulations

Background

The architectural organization of protein structures has been the focus of intense research since it can hopefully lead to an understanding of how proteins fold. In earlier works we had attempted to identify the inherent structural organization in proteins through a study of protein topology. We obtained a modular partitioning of protein structures with the modules correlating well with experimental evidence of early folding units or “foldons”. Residues that connect different modules were shown to be those that were protected during the transition phase of folding.

Methodology/Principal Findings

In this work, we follow the topological path of ubiquitin through molecular dynamics unfolding simulations. We observed that the use of recurrence quantification analysis (RQA) could lead to the identification of the transition state during unfolding. Additionally, our earlier contention that the modules uncovered through our graph partitioning approach correlated well with early folding units was vindicated through our simulations. Moreover, residues identified from native structure as connector hubs and which had been shown to be those that were protected during the transition phase of folding were indeed more stable (less flexible) well beyond the transition state. Further analysis of the topological pathway suggests that the all pairs shortest path in a protein is minimized during folding.

Conclusions

We observed that treating a protein native structure as a network by having amino acid residues as nodes and the non-covalent interactions among them as links allows for the rationalization of many aspects of the folding process. The possibility to derive this information directly from 3D structure opens the way to the prediction of important residues in proteins, while the confirmation of the minimization of APSP for folding allows for the establishment of a potentially useful proxy for kinetic optimality in the validation of sequence-structure predictions.     

Topological determinants of protein unfolding rates

For proteins that fold by two-state kinetics, the folding and unfolding processes are believed to be closely related to their native structures. In particular, folding and unfolding rates are influenced by the native structures of proteins. Thus, we focus on finding important topological quantities from a protein structure that determine its unfolding rate. After constructing graphs from protein native structures, we investigate the relationships between unfolding rates and various topological quantities of the graphs. First, we find that the correlation between the unfolding rate and the contact order is not as prominent as in the case of the folding rate and the contact order. Next, we investigate the correlation between the unfolding rate and the clustering coefficient of the graph of a protein native structure, and observe no correlation between them. Finally, we find that a newly introduced quantity, the impact of edge removal per residue, has a good overall correlation with protein unfolding rates. The impact of edge removal is defined as the ratio of the change of the average path length to the edge removal probability. From these facts, we conclude that the protein unfolding process is closely related to the protein native structure.     

Probing possible downhill folding: native contact topology likely places a significant constraint on the folding cooperativity of proteins with approximately 40 residues

Experiments point to appreciable variations in folding cooperativity among natural proteins with approximately 40 residues, indicating that the behaviors of these proteins are valuable for delineating the contributing factors to cooperative folding. To explore the role of native topology in a protein's propensity to fold cooperatively and how native topology might constrain the degree of cooperativity achievable by a given set of physical interactions, we compared folding/unfolding kinetics simulated using three classes of native-centric C^α chain models with different interaction schemes. The approach was applied to two homologous 45-residue fragments from the peripheral subunit-binding domain family and a 39-residue fragment of the N-terminal domain of ribosomal protein L9. Free-energy profiles as functions of native contact number were computed to assess the heights of thermodynamic barriers to folding. In addition, chevron plots of folding/unfolding rates were constructed as functions of native stability to facilitate comparison with available experimental data. Although common Gō-like models with pairwise Lennard-Jones-type interactions generally fold less cooperatively than real proteins, the rank ordering of cooperativity predicted by these models is consistent with experiment for the proteins investigated, showing increasing folding cooperativity with increasing nonlocality of a protein's native contacts. Models that account for water-expulsion (desolvation) barriers and models with many-body (nonadditive) interactions generally entail higher degrees of folding cooperativity indicated by more linear model chevron plots, but the rank ordering of cooperativity remains unchanged. A robust, experimentally valid rank ordering of model folding cooperativity independent of the multiple native-centric interaction schemes tested here argues that native topology places significant constraints on how cooperatively a protein can fold.     

Rotational restriction of nascent peptides as an essential element of co-translational protein folding: possible molecular players and structural consequences

Background

A basic tenet of protein science is that all information about the spatial structure of proteins is present in their sequences. Nonetheless, many proteins fail to attain native structure upon experimental denaturation and refolding in vitro, raising the question of the specific role of cellular machinery in protein folding in vivo. Recently, we hypothesized that energy-dependent twisting of the protein backbone is an unappreciated essential factor guiding the protein folding process in vivo. Torque force may be applied by the ribosome co-translationally, and when accompanied by simultaneous restriction of the rotational mobility of the distal part of the growing chain, the resulting tension in the protein backbone would facilitate the formation of local secondary structure and direct the folding process.

Results

Our model of the early stages of protein folding in vivo postulates that the free motion of both terminal regions of the protein during its synthesis and maturation is restricted. The long-known but unexplained phenomenon of statistical overrepresentation of protein termini on the surfaces of the protein structures may be an indication of the backbone twist-based folding mechanism; sustained maintenance of a twist requires that both ends of the protein chain are anchored in space, and if the ends are released only after the majority of folding is complete, they are much more likely to remain on the surface of the molecule. We identified the molecular components that are likely to play a role in the twisting of the nascent protein chain and in the anchoring of its N-terminus. The twist may be induced at the C-terminus of the nascent polypeptide by the peptidyltransferase center of the ribosome. Several ribosome-associated proteins, including the trigger factor in bacteria and the nascent polypeptide-associated complex in archaea and eukaryotes, may restrict the rotational mobility of the N-proximal regions of the peptides.

Conclusions

Many experimental observations are consistent with the hypothesis of co-translational twisting of the protein backbone. Several molecular players in this hypothetical mechanism of protein folding can be suggested. In addition, the new view of protein folding in vivo opens the possibility of novel potential drug targets to combat human protein folding diseases.

Reviewers

This article was reviewed by Lakshminarayan Iyer and István Simon.

     

Folding circular permutants of IL-1β: route selection driven by functional frustration

Interleukin-1β (IL-1β) is the cytokine crucial to inflammatory and immune response. Two dominant routes are populated in the folding to native structure. These distinct routes are a result of the competition between early packing of the functional loops versus closure of the β-barrel to achieve efficient folding and have been observed both experimentally and computationally. Kinetic experiments on the WT protein established that the dominant route is characterized by early packing of geometrically frustrated functional loops. However, deletion of one of the functional loops, the β-bulge, switches the dominant route to an alternative, yet, as accessible, route, where the termini necessary for barrel closure form first. Here, we explore the effect of circular permutation of the WT sequence on the observed folding landscape with a combination of kinetic and thermodynamic experiments. Our experiments show that while the rate of formation of permutant protein is always slower than that observed for the WT sequence, the region of initial nucleation for all permutants is similar to that observed for the WT protein and occurs within a similar timescale. That is, even permutants with significant sequence rearrangement in which the functional-nucleus is placed at opposing ends of the polypeptide chain, fold by the dominant WT "functional loop-packing route", despite the entropic cost of having to fold the N- and C- termini early. Taken together, our results indicate that the early packing of the functional loops dominates the folding landscape in active proteins, and, despite the entropic penalty of coalescing the termini early, these proteins will populate an entropically unfavorable route in order to conserve function. More generally, circular permutation can elucidate the influence of local energetic stabilization of functional regions within a protein, where topological complexity creates a mismatch between energetics and topology in active proteins.     

Effect of glutathione reductase knockdown in response to UVB-induced oxidative stress in human lung adenocarcinoma

Background

Glutathione reductase (GR) plays a critical role in the maintenance of physiological redox status in cells. However, the comprehensive investigations of GR-modulated oxidative stress have not been reported.

Methods

In the present study, we cultured a human lung adenocarcinoma line CL1-0 and its GR-knockdown derivative CL1-0ΔGR to evaluate their differential responses to UVB-irradiation.

Results

We identified 18 proteins that showed significant changes under UVB-irradiation in CL1-0ΔGR cells rather than in CL1-0 cells. Several proteins involving protein folding, metabolism, protein biosynthesis and redox regulation showed significant changes in expression.

Conclusions

In summary, the current study used a comprehensive lung adenocarcinoma-based proteomic approach for the identification of GR-modulated protein expression in response to UVB-irradiation. To our knowledge, this is the first global proteomic analysis to investigate the role of GR under UVB-irradiation in mammalian cell model.     

Folding simulations and computer redesign of protein A three-helix bundle motifs

In solution, the B domain of protein A from Staphylococcus aureus (B domain) possesses a three-helix bundle structure. This simple motif has been previously reproduced by Kolinski and Skolnick (Proteins 18: 353–366, 1994) using a reduced representation lattice model of proteins with a statistical interaction scheme. In this paper, an improved version of the potential has been used, and the robustness of this result has been tested by folding from the random state a set of three-helix bundle proteins that are highly homologous to the B domain of protein A. Furthermore, an attempt to redesign the B domain native structure to its topological mirror image fold has been made by multiple mutations of the hydrophobic core and the turn region between helices I and II. A sieve method for scanning a large set of mutations to search for this desired property has been proposed. It has been shown that mutations of native B domain hydrophobic core do not introduce significant changes in the protein motif. Mutations in the turn region were also very conservative; nevertheless, a few mutants acquired the desired topological mirror image motif. A set of all atom models of the most probable mutant was reconstructed from the reduced models and refined using a molecular dynamics algorithm in the presence of water. The packing of all atom structures obtained corroborates the lattice model results. We conclude that the change in the handedness of the turn induced by the mutations, augmented by the repacking of hydrophobic core and the additional burial of the second helix N-cap side chain, are responsible for the predicted preferential adoption of the mirror image structure. © 1996 Wiley-Liss, Inc.     

A protein folding potential that places the native states of a large number of proteins near a local minimum

Background

We present a simple method to train a potential function for the protein folding problem which, even though trained using a small number of proteins, is able to place a significantly large number of native conformations near a local minimum. The training relies on generating decoys by energy minimization of the native conformations using the current potential and using a physically meaningful objective function (derivative of energy with respect to torsion angles at the native conformation) during the quadratic programming to place the native conformation near a local minimum.

Results

We also compare the performance of three different types of energy functions and find that while the pairwise energy function is trainable, a solvation energy function by itself is untrainable if decoys are generated by minimizing the current potential starting at the native conformation. The best results are obtained when a pairwise interaction energy function is used with solvation energy function.

Conclusions

We are able to train a potential function using six proteins which places a total of 42 native conformations within ~4 Å rmsd and 71 native conformations within ~6 Å rmsd of a local minimum out of a total of 91 proteins. Furthermore, the threading test using the same 91 proteins ranks 89 native conformations to be first and the other two as second.     

Experimental evolution and genome sequencing reveal variation in levels of clonal interference in large populations of bacteriophage φX174

Background

The first step of GPI anchor biosynthesis is catalyzed by PIG-A, an enzyme that transfers N -acetylglucosamine from UDP- N -acetylglucosamine to phosphatidylinositol. This protein is present in all eukaryotic organisms ranging from protozoa to higher mammals, as part of a larger complex of five to six 'accessory' proteins whose individual roles in the glycosyltransferase reaction are as yet unclear. The PIG-A gene has been shown to be an essential gene in various eukaryotes. In humans, mutations in the protein have been associated with paroxysomal noctural hemoglobuinuria. The corresponding PIG-A gene has also been recently identified in the genome of many archaeabacteria although genes of the accessory proteins have not been discovered in them. The present study explores the evolution of PIG-A and the phylogenetic relationship between this protein and other glycosyltransferases.

Results

In this paper we show that out of the twelve conserved motifs identified by us eleven are exclusively present in PIG-A and, therefore, can be used as markers to identify PIG-A from newly sequenced genomes. Three of these motifs are absent in the primitive eukaryote, G. lamblia. Sequence analyses show that seven of these conserved motifs are present in prokaryote and archaeal counterparts in rudimentary forms and can be used to differentiate PIG-A proteins from glycosyltransferases. Using partial least square regression analysis and data involving presence or absence of motifs in a range of PIG-A and glycosyltransferases we show that (i) PIG-A may have evolved from prokaryotic glycosyltransferases and lipopolysaccharide synthases, members of the GT4 family of glycosyltransferases and (ii) it is possible to uniquely classify PIG-A proteins versus glycosyltransferases.

Conclusion

Besides identifying unique motifs and showing that PIG-A protein from G. lamblia and some putative PIG-A proteins from archaebacteria are evolutionarily closer to glycosyltransferases, these studies provide a new method for identification and classification of PIG-A proteins.     

Functional evolution of two subtly different (similar) folds

Background

The function of proteins is a direct consequence of their three-dimensional structure. The structural classification of proteins describes the ways of folding patterns all proteins could adopt. Although, the protein folds were described in many ways the functional properties of individual folds were not studied.

Results

We have analyzed two β-barrel folds generally adopted by small proteins to be looking similar but have different topology. On the basis of the topology they could be divided into two different folds named SH3-fold and OB-fold. There was no sequence homology between any of the proteins considered. The sequence diversity and loop variability was found to be important for various binding functions.

Conclusions

The function of Oligonucleotide/oligosaccharide-binding (OB) fold proteins was restricted to either DNA/RNA binding or sugar binding whereas the Src homology 3 (SH3) domain like proteins bind to a variety of ligands through loop modulations. A question was raised whether the evolution of these two folds was through DNA shuffling.     

Monte carlo simulations of protein folding. II. Application to protein A,ROP, and crambin

The hierarchy of lattice Monte Carlo models described in the accompanying paper (Kolinski, A., Skolnick, J. Monte Carlo simulations of protein folding. I. Lattice model and interaction scheme. Proteins 18:338–352, 1994) is applied to the simulation of protein folding and the prediction of 3-dimensional structure. Using sequence information alone, three proteins have been successfully folded: the B domain of staphylococcal protein A, a 120 residue, monomeric version of ROP dimer, and crambin. Starting from a random expanded conformation, the model proteins fold along relatively well-defined folding pathways. These involve a collection of early intermediates, which are followed by the final (and rate-determining) transition from compact intermediates closely resembling the molten globule state to the native-like state. The predicted structures are rather unique, with native-like packing of the side chains. The accuracy of the predicted native conformations is better than those obtained in previous folding simulations. The best (but by no means atypical) folds of protein A have a coordinate rms of 2.25 Å from the native Cα trace, and the best coordinate rms from crambin is 3.18 Å. For ROP monomer, the lowest coordinate rms from equivalent Cαs of ROP dimer is 3.65 Å. Thus, for two simple helical proteins and a small α/β protein, the ability to predict protein structure from sequence has been demonstrated. © 1994 John Wiley & Sons, Inc.     

Comparison of mcl-Poly(3-hydroxyalkanoates) synthesis by different Pseudomonas putida strains from crude glycerol: citrate accumulates at high titer under PHA-producing conditions

Background

Eukaryotic ubiquitin and SUMO are frequently used as tags to enhance the fusion protein expression in microbial host. They increase the solubility and stability, and protect the peptides from proteolytic degradation due to their stable and highly conserved structures. Few of prokaryotic ubiquitin-like proteins was used as fusion tags except ThiS, which enhances the fusion expression, however, reduces the solubility and stability of the expressed peptides in E. coli. Hence, we investigated if MoaD, a conserved small sulfur carrier in prokaryotes with the similar structure of ubiquitin, could also be used as fusion tag in heterologous expression in E. coli.

Results

Fusion of MoaD to either end of EGFP enhanced the expression yield of EGFP with a similar efficacy of ThiS. However, the major parts of the fusion proteins were expressed in the aggregated form, which was associated with the retarded folding of EGFP, similar to ThiS fusions. Fusion of MoaD to insulin chain A or B did not boost their expression as efficiently as ThiS tag did, probably due to a less efficient aggregation of products. Interestingly, fusion of MoaD to the murine ribonuclease inhibitor enhanced protein expression by completely protecting the protein from intracellular degradation in contrast to ThiS fusion, which enhanced degradation of this unstable protein when expressed in E. coli.

Conclusions

Prokaryotic ubiquitin-like protein MoaD can act as a fusion tag to promote the fusion expression with varying mechanisms, which enriches the arsenal of fusion tags in the category of insoluble expression.     

All-atom ab initio folding of a diverse set of proteins

Natural proteins fold to a unique, thermodynamically dominant state. Modeling of the folding process and prediction of the native fold of proteins are two major unsolved problems in biophysics. Here, we show successful all-atom ab initio folding of a representative diverse set of proteins by using a minimalist transferable-energy model that consists of two-body atom-atom interactions, hydrogen bonding, and a local sequence-energy term that models sequence-specific chain stiffness. Starting from a random coil, the native-like structure was observed during replica exchange Monte Carlo (REMC) simulation for most proteins regardless of their structural classes; the lowest energy structure was close to native-in the range of 2-6 A root-mean-square deviation (rmsd). Our results demonstrate that the successful folding of a protein chain to its native state is governed by only a few crucial energetic terms.     

Thermodynamics and folding kinetics analysis of the SH3 domain form discrete molecular dynamics

We perform a detailed analysis of the thermodynamics and folding kinetics of the SH3 domain fold with discrete molecular dynamic simulations. We propose a protein model that reproduces some of the experimentally observed thermodynamic and folding kinetic properties of proteins. Specifically, we use our model to study the transition state ensemble of the SH3 fold family of proteins, a set of unstable conformations that fold to the protein native state with probability 1/2. We analyze the participation of each secondary structure element formed at the transition state ensemble. We also identify the folding nucleus of the SH3 fold and test extensively its importance for folding kinetics. We predict that a set of amino acid contacts between the RT-loop and the distal hairpin are the critical folding nucleus of the SH3 fold and propose a hypothesis that explains this result.     

Designing cooperativity into the designed protein Top7

The topology of the designed protein Top7 is not found in natural proteins. Top7 shows signatures of non‐cooperative folding in both experimental studies and computer simulations. In particular, molecular dynamics of coarse‐grained structure‐based models of Top7 show a well‐populated C‐terminal folding‐intermediate. Since most similarly sized globular proteins are cooperative folders, the non‐natural topology of Top7 has been suggested as a reason for its non‐cooperative folding. Here, we computationally examine the folding of Top7 with the intent of making it cooperative. We find that its folding cooperativity can be increased in two ways: (a) Optimization of packing interactions in the N‐terminal half of the protein enables further folding of the C‐terminal intermediate. (b) Reduction in the packing density of the C‐terminal region destabilizes the intermediate. In practice, these strategies are implemented in our Top7 model through modifications to the contact‐map. These modifications do not alter the topology of Top7 but result in cooperative folding. Amino‐acid mutations that mimic these modifications also lead to a significant increase in folding cooperativity. Finally, we devise a method to randomize the sizes of amino‐acids within the same topology, and confirm that the structure of Top7 makes its folding sensitive to packing interactions. In contrast, the ribosomal protein S6, which has secondary structure similar to Top7, has folding which is much less sensitive to packing perturbations. Thus, it should be possible to make a sequence fold cooperatively to the structure of Top7, but to do so its side‐chain packing needs to be carefully designed. Proteins 2014; 82:364–374. © 2013 Wiley Periodicals, Inc.     

Discovery of a significant,nontopological preference for antiparallel alignment of helices with parallel regions in sheets

To help elucidate the role of secondary structure packing preferences in protein folding, here we present an analysis of the packing geometry observed between alpha-helices and between alpha-helices and beta-sheets in 1316 diverse, nonredundant protein structures. Finite-length vectors were fit to the alpha-carbon atoms in each of the helices and strands, and the packing angle between the vectors, Omega, was determined at the closest point of approach within each helix-helix or helix-sheet pair. Helix-sheet interactions were found in 391 of the proteins, and the distributions of Omega values were calculated for all the helix-sheet and helix-helix interactions. The packing angle preferences for helix-helix interactions are similar to those previously observed. However, analysis of helix-strand packing preferences uncovered a remarkable tendency for helices to align antiparallel to parallel regions of beta-sheets, independent of the topological constraints or prevalence of beta-alpha-beta motifs in the proteins. This packing angle preference is significantly diminished in helix interactions involving mixed and antiparallel beta-sheets, suggesting a role for helix-sheet dipole alignment in guiding supersecondary structure formation in protein folding. This knowledge of preferred packing angles can be used to guide the engineering of stable protein modules.     

Identifying protein complexes from interaction networks based on clique percolation and distance restriction

Background

Identification of protein complexes in large interaction networks is crucial to understand principles of cellular organization and predict protein functions, which is one of the most important issues in the post-genomic era. Each protein might be subordinate multiple protein complexes in the real protein-protein interaction networks. Identifying overlapping protein complexes from protein-protein interaction networks is a considerable research topic.

Result

As an effective algorithm in identifying overlapping module structures, clique percolation method (CPM) has a wide range of application in social networks and biological networks. However, the recognition accuracy of algorithm CPM is lowly. Furthermore, algorithm CPM is unfit to identifying protein complexes with meso-scale when it applied in protein-protein interaction networks. In this paper, we propose a new topological model by extending the definition of k-clique community of algorithm CPM and introduced distance restriction, and develop a novel algorithm called CP-DR based on the new topological model for identifying protein complexes. In this new algorithm, the protein complex size is restricted by distance constraint to conquer the shortcomings of algorithm CPM. The algorithm CP-DR is applied to the protein interaction network of Sacchromyces cerevisiae and identifies many well known complexes.

Conclusion

The proposed algorithm CP-DR based on clique percolation and distance restriction makes it possible to identify dense subgraphs in protein interaction networks, a large number of which correspond to known protein complexes. Compared to algorithm CPM, algorithm CP-DR has more outstanding performance.

     

A kingdom-specific protein domain HMM library for improved annotation of fungal genomes

Background

Mimivirus isolated from A. polyphaga is the largest virus discovered so far. It is unique among all the viruses in having genes related to translation, DNA repair and replication which bear close homology to eukaryotic genes. Nevertheless, only a small fraction of the proteins (33%) encoded in this genome has been assigned a function. Furthermore, a large fraction of the unassigned protein sequences bear no sequence similarity to proteins from other genomes. These sequences are referred to as ORFans. Because of their lack of sequence similarity to other proteins, they can not be assigned putative functions using standard sequence comparison methods. As part of our genome-wide computational efforts aimed at characterizing Mimivirus ORFans, we have applied fold-recognition methods to predict the structure of these ORFans and further functions were derived based on conservation of functionally important residues in sequence-template alignments.

Results

Using fold recognition, we have identified highly confident computational 3D structural assignments for 21 Mimivirus ORFans. In addition, highly confident functional predictions for 6 of these ORFans were derived by analyzing the conservation of functional motifs between the predicted structures and proteins of known function. This analysis allowed us to classify these 6 previously unannotated ORFans into their specific protein families: carboxylesterase/thioesterase, metal-dependent deacetylase, P-loop kinases, 3-methyladenine DNA glycosylase, BTB domain and eukaryotic translation initiation factor eIF4E.

Conclusion

Using stringent fold recognition criteria we have assigned three-dimensional structures for 21 of the ORFans encoded in the Mimivirus genome. Further, based on the 3D models and an analysis of the conservation of functionally important residues and motifs, we were able to derive functional attributes for 6 of the ORFans. Our computational identification of important functional sites in these ORFans can be the basis for a subsequent experimental verification of our predictions. Further computational and experimental studies are required to elucidate the 3D structures and functions of the remaining Mimivirus ORFans.     

Decomposing the space of protein quaternary structures with the interface fragment pair library

Background

The physical interactions between proteins constitute the basis of protein quaternary structures. They dominate many biological processes in living cells. Deciphering the structural features of interacting proteins is essential to understand their cellular functions. Similar to the space of protein tertiary structures in which discrete patterns are clearly observed on fold or sub-fold motif levels, it has been found that the space of protein quaternary structures is highly degenerate due to the packing of compact secondary structure elements at interfaces. Therefore, it is necessary to further decompose the protein quaternary structural space into a more local representation.

Results

Here we constructed an interface fragment pair library from the current structure database of protein complexes. After structural-based clustering, we found that more than 90% of these interface fragment pairs can be represented by a limited number of highly abundant motifs. These motifs were further used to guide complex assembly. A large-scale benchmark test shows that the native-like binding is highly likely in the structural ensemble of modeled protein complexes that were built through the library.

Conclusions

Our study therefore presents supportive evidences that the space of protein quaternary structures can be represented by the combination of a small set of secondary-structure-based packing at binding interfaces. Finally, after future improvements such as adding sequence profiles, we expect this new library will be useful to predict structures of unknown protein-protein interactions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0437-4) contains supplementary material, which is available to authorized users.