Search

In the protein sequence space, natural proteins form clusters of families which are characterized by their unique native folds whereas the great majority of random polypeptides are neither clustered nor foldable to unique structures. Since a given polypeptide can be either foldable or unfoldable, a kind of "folding transition" is expected at the boundary of a protein family in the sequence space. By Monte Carlo simulations of a statistical mechanical model of protein sequence alignment that coherently incorporates both short-range and long-range interactions as well as variable-length insertions to reproduce the statistics of the multiple sequence alignment of a given protein family, we demonstrate the existence of such transition between natural-like sequences and random sequences in the sequence subspaces for 15 domain families of various folds. The transition was found to be highly cooperative and two-state-like. Furthermore, enforcing or suppressing consensus residues on a few of the well-conserved sites enhanced or diminished, respectively, the natural-like pattern formation over the entire sequence. In most families, the key sites included ligand binding sites. These results suggest some selective pressure on the key residues, such as ligand binding activity, may cooperatively facilitate the emergence of a protein family during evolution. From a more practical aspect, the present results highlight an essential role of long-range effects in precisely defining protein families, which are absent in conventional sequence models., 13 pages, 7 figures, 2 tables (a new subsection added)

We re-examine a site-binding approach independently proposed by Schellman (Schellman, J.A. (1958) Compt. rend. Lab. Carlsberg Ser. Chim. 30, 439-449) and Aune and Tanford (Aune, K.C. and Tanford, D. (1969) Biochemistry, 8, 4586-4590) for explicitly including the denaturant concentration within the protein unfolding equilibrium. We extend and formalize the approach through development of a multi-dimensional analytical model in which the folding reaction coordinate is defined by the number of denaturant molecules bound to sites located on either the initially folded, or unfolded, states of the protein. We use the developed method to re-examine the mechanistic determinants underlying the sigmoidal shape of the unfolding transition. A natural feature of our method is that it presents a landscape picture of the denaturant induced protein unfolding reaction.

Cellular functions of an organism are maintained by protein–protein interactions. Those proteins that bind multiple partners asynchronously (date hub proteins) are important to make the interaction network coordinated. It is known that many date hub proteins bind different partners at overlapping (OV) interfaces. To understand how OV interfaces of date hub proteins can recognize multiple partners, we analyzed the difference between OV and non-overlapping (Non-OV) regions of interfaces involved in the binding of different partners. By using the structures of 16 date hub proteins with various interaction partners (ranging from 5 to 33), we compared buried surface area, compositions of amino acid residues and secondary structures, and side-chain orientations. It was found that buried interface residues are important for recognizing multiple partners, while exposed interface residues are important for determining specificity to a particular ligand. In addition, our analyses reveal that residue compositions in OV and Non-OV regions are different and that residues in OV region show diverse side-chain torsion angles to accommodate binding to multiple targets.

To elucidate the structural basis of the diversity and universality in protein-protein interactions, an exhaustive all-against-all structural comparison of all known protein interfaces in the Protein Data Bank was performed at atomic resolution. After similar interfaces were clustered, approximately 20,000 structural motifs with at least two members were identified, out of which 3678 motifs consisted of at least 10 interfaces. Except for some trivial interfaces involving single alpha helices, almost all motifs were found to be confined within single protein families. Furthermore, the interaction partners of each motif were found to be very limited, and, accordingly, the interaction networks of the motifs tend to be small and are much more restricted than the binding sites for small ligand molecules. These findings suggest that, at the level of atomic structures, protein-protein interactions are precisely designed; hence, protein interfaces with multiple interacting partners should involve incompletely overlapping multiple interfaces and/or accommodate structural changes upon binding to their targets.

Comprehensive knowledge of protein-ligand interactions should provide a useful basis for annotating protein functions, studying protein evolution, engineering enzymatic activity, and designing drugs. To investigate the diversity and universality of ligand binding sites in protein structures, we conducted the all-against-all atomic-level structural comparison of over 180,000 ligand binding sites found in all the known structures in the Protein Data Bank by using a recently developed database search and alignment algorithm. By applying a hybrid top-down-bottom-up clustering analysis to the comparison results, we determined approximately 3000 well-defined structural motifs of ligand binding sites. Apart from a handful of exceptions, most structural motifs were found to be confined within single families or superfamilies, and to be associated with particular ligands. Furthermore, we analyzed the components of the similarity network and enumerated more than 4000 pairs of ligand binding sites that were shared across different protein folds., Comment: 13 pages, 8 figures

Human transcriptional regulation factors, such as activators, repressors, and enhancer-binding factors are quite different from their prokaryotic counterparts in two respects: the average sequence in human is more than twice as long as that in prokaryotes, while the fraction of sequence aligned to domains of known structure is 31% in human transcription factors (TFs), less than half of that in bacterial TFs (72%). Intrinsically disordered (ID) regions were identified by a disorder-prediction program, and were found to be in good agreement with available experimental data. Analysis of 401 human TFs with experimental evidence from the Swiss-Prot database showed that as high as 49% of the entire sequence of human TFs is occupied by ID regions. More than half of the human TFs consist of a small DNA binding domain (DBD) and long ID regions frequently sandwiching unassigned regions. The remaining TFs have structural domains in addition to DBDs and ID regions. Experimental studies, particularly those with NMR, revealed that the transactivation domains in unbound TFs are usually unstructured, but become structured upon binding to their partners. The sequences of human and mouse TF orthologues are 90.5% identical despite a high incidence of ID regions, probably reflecting important functional roles played by ID regions. In general ID regions occupy a high fraction in TFs of eukaryotes, but not in prokaryotes. Implications of this dichotomy are discussed in connection with their functional roles in transcriptional regulation and evolution.

The living cell is inherently crowded with proteins and macromolecules. To avoid aggregation of denatured proteins in the living cell, molecular chaperones play important roles. Here we introduce a simple model to describe crowded protein solutions with chaperone-like species based on a dynamic density functional theory. As predicted by others, our simulations show that macromolecular crowding enhances the association of proteins and chaperones. However, when the intrinsic folding rate of the protein is slow, it is possible that crowding also enhances aggregation of proteins. The results of simulation suggest that, when the concentration of the crowding agent is as high as that in the cell, the association of the protein and unbound chaperone becomes correlated with the aggregation process, and that the protein-bound chaperones efficiently destroy the potential nuclei of aggregates and thus prevent the aggregation.

3D complex structures of proteins and other molecules provide a clue for mechanism of interaction. However, few servers for searching complex structures are available. We are developing an updated server HOMCOS (HOmology Modeling of Complex Structure) with more services than the previous version. It separates component molecules of PDB files of complexes, such as proteins, nucleic acids, small chemical compounds, stores their binding relationships. It searches these molecules by BLAST and our chemical structure comparison program dkcombu. Based on found similar complexes, simple template-based models of the complex can be generated by replacing the template with the query amino acid sequence. The service ‘Modeling Homo Protein Multimer’ searches homologous protein complexes with one query protein sequence, whereas the service ‘Modeling Hetero Protein Multimer’ searches with two sequences. The service ‘Modeling Compound-Protein Complexes’ searches with one sequence and one chemical structure. For modeling a bound conformation of a target compound, our flexible overlay program fkcombu is performed to superimpose the target compound onto the template compound bound to a protein. The program transforms a pose and torsion angles of the target molecule to superimpose atom pairs of the target and template molecule determined by 2D-MCS. The service ‘Searching Contacting Molecules for Query Protein’ searches contacting molecules with homologous proteins to a given query sequence, and provide putative interacting molecules with putative binding sites. It is useful to estimate mutational effects on molecular interactions. The service ‘Searching Contacting Proteins for Query Compound’ searching contacting proteins with similar compounds to a given query chemical compound. It may be useful to predict unintended interacting proteins, which cause side effects of the compound. We hope this server may contribute to the platform for drug discoveries and life sciences.