Conformation biases of amino acids based on tripeptide microenvironment from PDB database

We have constructed a bank (FTTP) of tendentious factors of three states of three-peptide units from PDB database based on conformational dihedral angle library and demonstrated that amino acid biases toward protein secondary structure are present in natural protein sequences. Our research results reveal that 20 standard amino acids fall into three groups: nine residues inclined to alpha-helix with a common character (e.g. direct side chain aliphatic residues or positive/negative charged residues) arrange in three grades, viz EA, QKRLD, and MN, in turn; seven residues are apt to beta-strand with 2'-branched side chain aliphatic residues or benzyl-included residues, namely PV, IYTC, and F, in three ranks; and four residues SHWG show a double tendency to both alpha and beta. Noticeably, proline has the strongest ability to form extended conformation, especially the Re value up to 9.5298 at position 3 (Table 3). Thus, biases of codons show an evident tendency in protein folding, where GC-rich codons are mainly in charge of forming contracted conformation, especially the codon's first letter plays a dominant role in translating the genomic GC signature into protein sequences and structures. So, biases of amino acids will play an important role in protein folding, folding codons, refining domain, structure prediction, and structural genomics/proteomics.     

SMS 2.0: an updated database to study the structural plasticity of short peptide fragments in non-redundant proteins

The function of a protein molecule is greatly influenced by its three-dimensional (3D) structure and therefore structure prediction will help identify its biological function. We have updated Sequence, Motif and Structure (SMS), the database of structurally rigid peptide fragments, by combining amino acid sequences and the corre-sponding 3D atomic coordinates of non-redundant (25%) and redundant (90%) protein chains available in the Protein Data Bank (PDB). SMS 2.0 provides information pertaining to the peptide fragments of length 5-14 resi-dues. The entire dataset is divided into three categories, namely, same sequence motifs having similar, intermedi-ate or dissimilar 3D structures. Further, options are provided to facilitate structural superposition using the pro-gram structural alignment of multiple proteins (STAMP) and the popular JAVA plug-in (Jmol) is deployed for visualization. In addition, functionalities are provided to search for the occurrences of the sequence motifs in other structural and sequence databases like PDB, Genome Database (GDB), Protein Information Resource (PIR) and Swiss-Prot. The updated database along with the search engine is available over the World Wide Web through the following URL http://cluster.physics.iisc.ernet.in/sms/.     

The illusive search for the lowest free energy state of globular proteins and RNAs

As a consequence of the one-dimensional storage and transfer of genetic information, DNA → RNA → protein, the process by which globular proteins and RNAs achieve their three-dimensional structure involves folding of a linear chain. The folding process itself could create massive activation barriers that prevent the attainment of many stable protein and RNA structures. We consider several kinds of energy barriers inherent in folding that might serve as kinetic constraints to achieving the lowest energy state. Alternative approaches to forming 3D structure, where a substantial number of weak interactions would be created prior to the formation of all the peptide (or phosphodiester) bonds, might not be subjected to such high barriers. This could lead to unique 3D conformational states, potentially more stable than “native” proteins and RNAs, with new functionalities.     

Folding of peptide fragments comprising the complete sequence of proteins. Models for initiation of protein folding. I. Myohemerythrin.

In an attempt to delineate potential folding initiation sites for different protein structural motifs, we have synthesized series of peptides that span the entire length of the polypeptide chain of two proteins, and examined their conformational preferences in aqueous solution using proton nuclear magnetic resonance and circular dichroism spectroscopy. We describe here the behavior of peptides derived from a simple four-helix bundle protein, myohemerythrin. The peptides correspond to the sequences of the four long helices (the A, B, C and D helices), the N- and C-terminal loops and the connecting sequences between the helices. The peptides corresponding to the helices of the folded protein all exhibit preferences for helix-like conformations in solution. The conformational ensembles of the A- and D-helix peptides contain ordered helical forms, as shown by extensive series of medium-range nuclear Overhauser effect connectivities, while the B- and C-helix peptides exhibit conformational preferences for nascent helix. All four peptides adopt ordered helical conformations in mixtures of trifluoroethanol and water. The terminal and interconnecting loop peptides also appear to contain appreciable populations of conformers with backbone phi and psi angles in the alpha-region and include highly populated hydrophobic cluster and/or turn conformations in some cases. Trifluoroethanol is unable to drive these peptides towards helical conformations. Overall, the peptide fragments of myohemerythrin have a marked preference towards secondary structure formation in aqueous solution. In contrast, peptide fragments derived from the beta-sandwich protein plastocyanin are relatively devoid of secondary structure in aqueous solution (see accompanying paper). These results suggest that the two different protein structural motifs may require different propensities for formation of local elements of secondary structure to initiate folding, and that there is a prepartitioning of conformational space determined by the local amino acid sequence that is different for the helical and beta-sandwich structural motifs.     

Structural features in EIAV NCp11: a lentivirus nucleocapsid protein with a short linker

Lentiviral nucleocapsid proteins are a class of multifunctional proteins that play an essential role in RNA packaging and viral infectivity. They contain two CX(2)CX(4)HX(4)C zinc binding motifs connected by a basic linker of variable length. The 3D structure of a 37-aa peptide corresponding to sequence 22-58 from lentiviral EIAV nucleocapsid protein NCp11, complexed with zinc, has been determined by 2D (1)H NMR spectroscopy, simulated annealing, and molecular dynamics. The solution structure consists of two zinc binding domains held together by a five-residue basic linker Arg(38)-Ala-Pro-Lys-Val(42) that allows for spatial proximity between the two finger domains. Observed linker folding is stabilized by H bonded secondary structure elements, resulting in an Omega-shaped central region, asymmetrically centered on the linker. The conformational differences and similarities with other NC zinc binding knuckles have been systematically analyzed. The two CCHC motifs, both characterized by a peculiar Pro-Gly sequence preceding the His residue, although preserving Zn-binding geometry and chirality of other known NC proteins, exhibit local fold differences both between each other and in comparison with other previously characterized retroviral CCHC motifs.     

The HSSP database of protein structure-sequence alignments.

HSSP (homology-derived structures of proteins) is a derived database merging structural (2-D and 3-D) and sequence information (1-D). For each protein of known 3D structure from the Protein Data Bank, the database has a file with all sequence homologues, properly aligned to the PDB protein. Homologues are very likely to have the same 3D structure as the PDB protein to which they have been aligned. As a result, the database is not only a database of sequence aligned sequence families, but it is also a database of implied secondary and tertiary structures.     

PSSARD (2.0): a database server for making flexible queries relating amino acid sequences to main-chain secondary structure conformations for proteins of known three-dimensional structure and certain useful applications

We have updated the Protein Sequence-Structure Analysis Relational Database (PSSARD) first published in the Int. J. Biol. Macromol. 36 (2005) 259-262 corresponding to 1573 representative protein chains selected from the Protein Data Bank (PDB). In this, the updated and revised PSSARD (Version 2.0), we have included all proteins in the Protein Data Bank available at the time of developing this database including the NMR PDB entries. The current database corresponds to 22,752 XRAY PDB entries and 3977 NMR PDB entries and is separated accordingly in order to facilitate the appropriate database search. The representative protein chains can also be separately accessed within the current database. We have made a provision to combine more than one field to query the database and the results of any search can be used to carry out further nested searches using a combination of queries. We have provided hyperlinks to the individual PDB entries obtained as the result of any search in PSSARD in order to obtain additional details relevant to the protein structure. Certain applications useful to identify domains and structural motifs are discussed.     

Analysis of the factors that stabilize a designed two-stranded antiparallel beta-sheet

Autonomously folding beta-hairpins (two-strand antiparallel beta-sheets) have become increasingly valuable tools for probing the forces that control peptide and protein conformational preferences. We examine the effects of variations in sequence and solvent on the stability of a previously designed 12-residue peptide (1). This peptide adopts a beta-hairpin conformation containing a two-residue loop (D-Pro-Gly) and a four-residue interstrand sidechain cluster that is observed in the natural protein GB1. We show that the conformational propensity of the loop segment plays an important role in beta-hairpin stability by comparing 1 with (D)P--> N mutant 2. In addition, we show that the sidechain cluster contributes both to conformational stability and to folding cooperativity by comparing 1 with mutant 3, in which two of the four cluster residues have been changed to serine. Thermodynamic analysis suggests that the high loop-forming propensity of the (D)PG segment decreases the entropic cost of beta-hairpin formation relative to the more flexible NG segment, but that the conformational rigidity of (D)PG may prevent optimal contacts between the sidechains of the GB1-derived cluster. The enthalpic favorability of folding in these designed beta-hairpins suggests that they are excellent scaffolds for studying the fundamental mechanisms by which amino acid sidechains interact with one another in folded proteins.     

Beta-turn propensities as paradigms for the analysis of structural motifs to engineer protein stability.

The thermodynamic stability of a protein provides an experimental metric for the relationship of protein sequence and native structure. We have investigated an approach based on an analysis of the structural database for stability engineering of an immunoglobulin variable domain. The most frequently occurring residues in specific positions of beta-turn motifs were predicted to increase the folding stability of mutants that were constructed by site-directed mutagenesis. Even in positions in which different residues are conserved in immunoglobulin sequences, the predictions were confirmed. Frequently, mutants with increased beta-turn propensities display increased folding cooperativities, suggesting pronounced effects on the unfolded state independent of the expected effect on conformational entropy. We conclude that structural motifs with predominantly local interactions can serve as templates with which patterns of sequence preferences can be extracted from the database of protein structures. Such preferences can predict the stability effects of mutations for protein engineering and design.     

COPS--cis/trans peptide bond conformation prediction of amino acids on the basis of secondary structure information

SUMMARY: COPS predicts for all 20 naturally occurring amino acids whether the peptide bond in a protein is in cis or trans conformation. The algorithm is based only on secondary structure information of amino acid triplets without considering the amino acid sequence information. Conformation parameters are derived from solved 3D structures deposited in the PDB and led to propensities based on modified Chou-Fasman parameters. COPS analyses amino acid triplets taking only their respective secondary structure into consideration and upon application of a set of rules utilizing the conformation parameters, the N-terminal peptide bond conformation of the middle residue is predicted. COPS was tested on a random selection of protein datasets. AVAILABILITY: The COPS program and further information are freely available from the FMP website at http://www.fmp-berlin.de/nmr/cops CONTACT: labudde@fmp-berlin.de.     

A rapid test for identification of autonomous folding units in proteins

The structure of a protein is dictated by a large number of weak interactions that cooperatively stabilize the native state. Usually, excised fragments smaller than a domain have little if any residual structure. When autonomous units of structure are found within domains, this challenges common assumptions about the cooperativity of protein structure. Such autonomous folding units (AFUs) are of wide interest and have applications in protein engineering and as simple model systems for studying the determinants of stability and specificity. A new method of identifying AFUs within proteins is presented here. The rapid autonomous fragment test (RAFT) identifies AFUs based on analysis of inter-residue contacts present in the three-dimensional structure of a protein. RAFT is fast enough to mine the entire PDB for AFUs and provide a library of potential small stable folds. We show that RAFT is able to predict whether a protein fragment will be structured if isolated from its parent domain.     

The chicken-egg scenario of protein folding revisited

What is the first step in protein folding - hydrophobic collapse (compaction) or secondary structure formation? It is still not clear if the major driving force in protein folding is hydrogen bonding or hydrophobic interactions or both. We analyzed data on the conformational characteristics of 41 globular proteins in native and partially folded conformational states. Our analysis shows that a good correlation exists between relative decrease in hydrodynamic volume and increase in secondary structure content. No compact equilibrium intermediates lacking secondary structure, or highly ordered non-compact species, were found. This correlation provides experimental support for the hypothesis that hydrophobic collapse occurs simultaneously with formation of secondary structure in the early stages of the protein folding.     

The SuMo server: 3D search for protein functional sites

SUMMARY: We provide the scientific community with a web server which gives access to SuMo, a bioinformatic system for finding similarities in arbitrary 3D structures or substructures of proteins. SuMo is based on a unique representation of macromolecules using selected triplets of chemical groups having their own geometry and symmetry, regardless of the restrictive notions of main chain and lateral chains of amino acids. The heuristic for extracting similar sites was used to drive two major large-scale approaches. First, searching for ligand binding sites onto a query structure has been made possible by comparing the structure against each of the ligand binding sites found in the Protein Data Bank (PDB). Second, the reciprocal process, i.e. searching for a given 3D site of interest among the structures of the PDB is also possible and helps detect cross-reacting targets in drug design projects. AVAILABILITY: The web server is freely accessible to academia through http://sumo-pbil.ibcp.fr and full support is available from MEDIT (http://www.medit.fr). CONTACT: mjambon@burnham.org.     

Het-PDB Navi.: a database for protein-small molecule interactions

The genomes of more than 100 species have been sequenced, and the biological functions of encoded proteins are now actively being researched. Protein function is based on interactions between proteins and other molecules. One approach to assuming protein function based on genomic sequence is to predict interactions between an encoded protein and other molecules. As a data source for such predictions, knowledge regarding known protein-small molecule interactions needs to be compiled. We have, therefore, surveyed interactions between proteins and other molecules in Protein Data Bank (PDB), the protein three-dimensional (3D) structure database. Among 20,685 entries in PDB (April, 2003), 4,189 types of small molecules were found to interact with proteins. Biologically relevant small molecules most often found in PDB were metal ions, such as calcium, zinc, and magnesium. Sugars and nucleotides were the next most common. These molecules are known to act as cofactors for enzymes and/or stabilizers of proteins. In each case of interactions between a protein and small molecule, we found preferred amino acid residues at the interaction sites. These preferences can be the basis for predicting protein function from genomic sequence and protein 3D structures. The data pertaining to these small molecules were collected in a database named Het-PDB Navi., which is freely available at http://daisy.nagahama-i-bio.ac.jp/golab/hetpdbnavi.html and linked to the official PDB home page.     

The nature of the free energy barriers to two-state folding

We introduce a simple procedure to analyze the temperature dependence of the folding and unfolding rates of two-state proteins. We start from the simple transition-state-like rate expression: k = D(eff)exp(-DeltaG(TS)/RT), in which upper and lower bounds for the intra-chain effective diffusion coefficient (D(eff)) are obtained empirically using the timescales of elementary processes in protein folding. From the changes in DeltaG(TS) as a function of temperature, we calculate enthalpies and heat capacities of activation, together with the more elusive entropies of activation. We then estimate the conformational entropy of the transition state by extrapolation to the temperature at which the solvation entropy vanishes by cancellation between polar and apolar terms. This approach is based on the convergence temperatures for the entropy of solvating apolar (approximately 385 K) and polar groups (approximately 335 K), the assumption that the structural properties of the transition state are somewhere in between the unfolded and folded states, and the established relationship between observed heat capacity and solvent accessibility.1 To circumvent the lack of structural information about transition states, we use the empirically determined heat capacities of activation as constraints to identify the extreme values of the transition state conformational entropy that are consistent with experiment. The application of this simple approach to six two-state folding proteins for which there is temperature-dependent data available in the literature provides important clues about protein folding. For these six proteins, we obtain an average equilibrium cost in conformational entropy of -4.3 cal x mol(-1)K(-1)per residue, which is in close agreement to previous empirical and computational estimates of the same quantity. Furthermore, we find that all these proteins have a conformationally diverse transition state, with more than half of the conformational entropy of the unfolded state. In agreement with predictions from theory and computer simulations, the transition state signals the change from a regime dominated by loss in conformational entropy to one driven by the gain in stabilization free energy (i.e., including protein interactions and solvation effects). Moreover, the height of the barrier is determined by how much stabilization free energy is realized at that point, which is related to the relative contribution of local versus non-local interactions. A remarkable observation is that the fraction of conformational entropy per residue that is present in the transition state is very similar for the six proteins in this study. Based on this commonality, we propose that the observed change in thermodynamic regime is connected to a change in the pattern of structure formation: from one driven by formation of pairwise interactions to one dominated by coupling of the networks of interactions involved in forming the protein core. In this framework, the barrier to two-state folding is crossed when the folding protein reaches a "critical native density" that allows expulsion of remaining interstitial water and consolidation of the core. The principle of critical native density should be general for all two-state proteins, but can accommodate different folding mechanisms depending on the particularities of the structure and sequence.     

The PAS fold. A redefinition of the PAS domain based upon structural prediction.

In the postgenomic era it is essential that protein sequences are annotated correctly in order to help in the assignment of their putative functions. Over 1300 proteins in current protein sequence databases are predicted to contain a PAS domain based upon amino acid sequence alignments. One of the problems with the current annotation of the PAS domain is that this domain exhibits limited similarity at the amino acid sequence level. It is therefore essential, when using proteins with low-sequence similarities, to apply profile hidden Markov model searches for the PAS domain-containing proteins, as for the PFAM database. From recent 3D X-ray and NMR structures, however, PAS domains appear to have a conserved 3D fold as shown here by structural alignment of the six representative 3D-structures from the PDB database. Large-scale modelling of the PAS sequences from the PFAM database against the 3D-structures of these six structural prototypes was performed. All 3D models generated (> 5700) were evaluated using prosaii. We conclude from our large-scale modelling studies that the PAS and PAC motifs (which are separately defined in the PFAM database) are directly linked and that these two motifs form the PAS fold. The existing subdivision in PAS and PAC motifs, as used by the PFAM and SMART databases, appears to be caused by major differences in sequences in the region connecting these two motifs. This region, as has been shown by Gardner and coworkers for human PAS kinase (Amezcua, C.A., Harper, S.M., Rutter, J. & Gardner, K.H. (2002) Structure 10, 1349-1361, [1]), is very flexible and adopts different conformations depending on the bound ligand. Some PAS sequences present in the PFAM database did not produce a good structural model, even after realignment using a structure-based alignment method, suggesting that these representatives are unlikely to have a fold resembling any of the structural prototypes of the PAS domain superfamily.     

Conformational analysis of invariant peptide sequences in bacterial genomes

The functional significance of evolutionarily conserved motifs/patterns of short regions in proteins is well documented. Although a large number of sequences are conserved, only a small fraction of these are invariant across several organisms. Here, we have examined the structural features of the functionally important peptide sequences, which have been found invariant across diverse bacterial genera. Ramachandran angles (phi,psi) have been used to analyze the conformation, folding patterns and geometrical location (buried/exposed) of these invariant peptides in different crystal structures harboring these sequences. The analysis indicates that the peptides preferred a single conformation in different protein structures, with the exception of only a few longer peptides that exhibited some conformational variability. In addition, it is noticed that the variability of conformation occurs mainly due to flipping of peptide units about the virtual C(alpha)...C(alpha) bond. However, for a given invariant peptide, the folding patterns are found to be similar in almost all the cases. Over and above, such peptides are found to be buried in the protein core. Thus, we can safely conclude that these invariant peptides are structurally important for the proteins, since they acquire unique structures across different proteins and can act as structural determinants (SD) of the proteins. The location of these SD peptides on the protein chain indicated that most of them are clustered towards the N-terminal and middle region of the protein with the C-terminal region exhibiting low preference. Another feature that emerges out of this study is that some of these SD peptides can also play the roles of "fold boundaries" or "hinge nucleus" in the protein structure. The study indicates that these SD peptides may act as chain-reversal signatures, guiding the proteins to adopt appropriate folds. In some cases the invariant signature peptides may also act as folding nuclei (FN) of the proteins.     

PSI: indexing protein structures for fast similarity search

MOTIVATION: We consider the problem of finding similarities in protein structure databases. Current techniques sequentially compare the given query protein to all of the proteins in the database to find similarities. Therefore, the cost of similarity queries increases linearly as the volume of the protein databases increase. As the sizes of experimentally determined and theoretically estimated protein structure databases grow, there is a need for scalable searching techniques. RESULTS: Our techniques extract feature vectors on triplets of SSEs (Secondary Structure Elements). Later, these feature vectors are indexed using a multidimensional index structure. For a given query protein, this index structure is used to quickly prune away unpromising proteins in the database. The remaining proteins are then aligned using a popular alignment tool such as VAST. We also develop a novel statistical model to estimate the goodness of a match using the SSEs. Experimental results show that our techniques improve the pruning time of VAST 3 to 3.5 times while maintaining similar sensitivity.     

A protein short motif search tool using amino acid sequence and their secondary structure assignment

We present the development of a web server, a protein short motif search tool that allows users to simultaneously search for a protein sequence motif and its secondary structure assignments. The web server is able to query very short motifs searches against PDB structural data from the RCSB Protein Databank, with the users defining the type of secondary structures of the amino acids in the sequence motif. The output utilises 3D visualisation ability that highlights the position of the motif in the structure and on the corresponding sequence. Researchers can easily observe the locations and conformation of multiple motifs among the results. Protein short motif search also has an application programming interface (API) for interfacing with other bioinformatics tools. AVAILABILITY: The database is available for free at http://birg3.fbb.utm.my/proteinsms.     

The dynamical contact order: Protein folding rate parameters based on quantum conformational transitions

Protein folding is regarded as a quantum transition between the torsion states of a polypeptide chain. According to the quantum theory of conformational dynamics, we propose the dynamical contact order (DCO) defined as a characteristic of the contact described by the moment of inertia and the torsion potential energy of the polypeptide chain between contact residues. Consequently, the protein folding rate can be quantitatively studied from the point of view of dynamics. By comparing theoretical calculations and experimental data on the folding rate of 80 proteins, we successfully validate the view that protein folding is a quantum conformational transition. We conclude that (i) a correlation between the protein folding rate and the contact inertial moment exists; (ii) multi-state protein folding can be regarded as a quantum conformational transition similar to that of two-state proteins but with an intermediate delay. We have estimated the order of magnitude of the time delay; (iii) folding can be classified into two types, exergonic and endergonic. Most of the two-state proteins with higher folding rate are exergonic and most of the multi-state proteins with low folding rate are endergonic. The folding speed limit is determined by exergonic folding.