The search for local native-like nucleation centers in the unfolded state of beta -sheet proteins

Phyloproteomic Analysis of 11780 Six-Residue-Long Motifs Occurrences

How is it possible to find good traits for phylogenetic reconstructions? Here, we present a new phyloproteomic criterion that is an occurrence of simple motifs which can be imprints of evolution history. We studied the occurrences of 11780 six-residue-long motifs consisting of two randomly located amino acids in 97 eukaryotic and 25 bacterial proteomes. For all eukaryotic proteomes, with the exception of the Amoebozoa, Stramenopiles, and Diplomonadida kingdoms, the number of proteins containing the motifs from the first group (one of the two amino acids occurs once at the terminal position) made about 20%; in the case of motifs from the second (one of two amino acids occurs one time within the pattern) and third (the two amino acids occur randomly) groups, 30% and 50%, respectively. For bacterial proteomes, this relationship was 10%, 27%, and 63%, respectively. The matrices of correlation coefficients between numbers of proteins where a motif from the set of 11780 motifs appears at least once in 9 kingdoms and 5 phyla of bacteria were calculated. Among the correlation coefficients for eukaryotic proteomes, the correlation between the animal and fungi kingdoms (0.62) is higher than between fungi and plants (0.54). Our study provides support that animals and fungi are sibling kingdoms. Comparison of the frequencies of six-residue-long motifs in different proteomes allows obtaining phylogenetic relationships based on similarities between these frequencies: the Diplomonadida kingdoms are more close to Bacteria than to Eukaryota; Stramenopiles and Amoebozoa are more close to each other than to other kingdoms of Eukaryota.

Early folding events protect aggregation-prone regions of a β-rich protein

Protein folding and aggregation inevitably compete with one another. This competition is even keener for proteins with frustrated landscapes, such as those rich in β structure. It is interesting that, despite their rugged energy landscapes and high β sheet content, intracellular lipid-binding proteins (iLBPs) appear to successfully avoid aggregation, as they are not implicated in aggregation diseases. In this study, we used a canonical iLBP, cellular retinoic acid-binding protein 1 (CRABP1), to understand better how folding is favored over aggregation. Analysis of folding kinetics of point mutants reveals that the folding pathway of CRABP1 involves early barrel closure. This folding mechanism protects sequences in CRABP1 that comprise cores of aggregates as identified by nuclear magnetic resonance. The amino acid conservation pattern in other iLBPs suggests that early barrel closure may be a general strategy for successful folding and minimization of aggregation. We suggest that folding mechanisms in general may incorporate steps that disfavor aggregation.

Temperature and urea have opposing impacts on polyproline II conformational bias

The native states of globular proteins have been accessed in atomic detail by X-ray crystallography and nuclear magnetic resonance spectroscopy, yet characterization of denatured proteins beyond global metrics has proven to be elusive. Denatured proteins have been observed to exhibit global geometric properties of a random coil polymer. However, this does not preclude the existence of nonrandom, local conformational bias that may be significant for protein folding and function. Indeed, circular dichroism (CD) spectroscopy and other methods have suggested that the denatured state contains considerable local bias to the polyproline II (PII) conformation. Here, we develop predictive models to determine the extent that temperature and the chemical denaturant urea modulate PII propensity. In agreement with our predictive model, PII propensity is observed experimentally to decrease with an increase in temperature. Conversely, urea appears to promote the PII conformation as determined by CD and isothermal titration calorimetry. Importantly, the calorimetric data are in quantitative agreement with a model that predicts the stability of the PII helix relative to other denatured state conformations based upon solvent accessible surface area and experimentally measured Gibbs transfer free energies. The ability of urea to promote the PII conformation can be attributed to the favorable interaction of urea with the peptide backbone. Thus, perturbing denatured states by temperature or cosolutes has subtle, yet opposing, impacts on local PII conformational biases. These results have implications for protein folding as well as for the function of signaling proteins that bind proline-rich targets in globular or intrinsically disordered proteins.

Disordered patterns in clustered Protein Data Bank and in eukaryotic and bacterial proteomes

We have constructed the clustered Protein Data Bank and obtained clusters of chains of different identity inside each cluster, http://bioinfo.protres.ru/st_pdb/. We have compiled the largest database of disordered patterns (141) from the clustered PDB where identity between chains inside of a cluster is larger or equal to 75% (version of 28 June 2010) by using simple rules of selection. The results of these analyses would help to further our understanding of the physicochemical and structural determinants of intrinsically disordered regions that serve as molecular recognition elements. We have analyzed the occurrence of the selected patterns in 97 eukaryotic and in 26 bacterial proteomes. The disordered patterns appear more often in eukaryotic than in bacterial proteomes. The matrix of correlation coefficients between numbers of proteins where a disordered pattern from the library of 141 disordered patterns appears at least once in 9 kingdoms of eukaryota and 5 phyla of bacteria have been calculated. As a rule, the correlation coefficients are higher inside of the considered kingdom than between them. The patterns with the frequent occurrence in proteomes have low complexity (PPPPP, GGGGG, EEEED, HHHH, KKKKK, SSTSS, QQQQQP), and the type of patterns vary across different proteomes, http://bioinfo.protres.ru/fp/search_new_pattern.html.

Residual interactions in unfolded bile acid-binding protein by 19F NMR

The folding initiation mechanism of human bile acid-binding protein (BABP) has been examined by (19) F NMR. Equilibrium unfolding studies of BABP labeled with fluorine at all eight of its phenylalanine residues showed that at least two sites experience changes in solvent exposure at high denaturant concentrations. Peak assignments were made by site-specific 4FPhe incorporation. The resonances for proteins specifically labeled at Phe17, Phe47, and Phe63 showed changes in chemical shift at denaturant concentrations at which the remaining five phenylalanine residues appear to be fully solvent-exposed. Phe17 is a helical residue that was not expected to participate in a folding initiation site. Phe47 and Phe63 form part of a hydrophobic core region that may be conserved as a site for folding initiation in the intracellular lipid-binding protein family.

Study and prediction of secondary structure for membrane proteins

In this paper we present a novel approach to membrane protein secondary structure prediction based on the statistical stepwise discriminant analysis method. A new aspect of our approach is the possibility to derive physical-chemical properties that may affect the formation of membrane protein secondary structure. The certain physical-chemical properties of protein chains can be used to clarify the formation of the secondary structure types under consideration. Another aspect of our approach is that the results of multiple sequence alignment, or the other kinds of sequence alignment, are not used in the frame of the method. Using our approach, we predicted the formation of three main secondary structure types (alpha-helix, beta-structure and coil) with high accuracy, that is Q(3) = 76%. Predicting the formation of alpha-helix and non-alpha-helix states we reached the accuracy which was measured as Q(2) = 86%. Also we have identified certain protein chain properties that affect the formation of membrane protein secondary structure. These protein properties include hydrophobic properties of amino acid residues, presence of Gly, Ala and Val amino acids, and the location of protein chain end.

Water and urea interactions with the native and unfolded forms of a beta-barrel protein

A fundamental understanding of protein stability and the mechanism of denaturant action must ultimately rest on detailed knowledge about the structure, solvation, and energetics of the denatured state. Here, we use (17)O and (2)H magnetic relaxation dispersion (MRD) to study urea-induced denaturation of intestinal fatty acid-binding protein (I-FABP). MRD is among the few methods that can provide molecular-level information about protein solvation in native as well as denatured states, and it is used here to simultaneously monitor the interactions of urea and water with the unfolding protein. Whereas CD shows an apparently two-state transition, MRD reveals a more complex process involving at least two intermediates. At least one water molecule binds persistently (with residence time >10 nsec) to the protein even in 7.5 M urea, where the large internal binding cavity is disrupted and CD indicates a fully denatured protein. This may be the water molecule buried near the small hydrophobic folding core at the D-E turn in the native protein. The MRD data also provide insights about transient (residence time <1 nsec) interactions of urea and water with the native and denatured protein. In the denatured state, both water and urea rotation is much more retarded than for a fully solvated polypeptide. The MRD results support a picture of the denatured state where solvent penetrates relatively compact clusters of polypeptide segments.

Sequence and structural analysis of cellular retinoic acid-binding proteins reveals a network of conserved hydrophobic interactions

Proteins in the intracellular lipid-binding protein (iLBP) family show remarkably high structural conservation despite their low-sequence identity. A multiple-sequence alignment using 52 sequences of iLBP family members revealed 15 fully conserved positions, with a disproportionately high number of these (n=7) located in the relatively small helical region. The conserved positions displayed high structural conservation based on comparisons of known iLBP crystal structures. It is striking that the beta-sheet domain had few conserved positions, despite its high structural conservation. This observation prompted us to analyze pair-wise interactions within the beta-sheet region to ask whether structural information was encoded in interacting amino acid pairs. We conducted this analysis on the iLBP family member, cellular retinoic acid-binding protein I (CRABP I), whose folding mechanism is under study in our laboratory. Indeed, an analysis based on a simple classification of hydrophobic and polar amino acids revealed a network of conserved interactions in CRABP I that cluster spatially, suggesting a possible nucleation site for folding. Significantly, a small number of residues participated in multiple conserved interactions, suggesting a key role for these sites in the structure and folding of CRABP I. The results presented here correlate well with available experimental evidence on folding of CRABPs and their family members and suggest future experiments. The analysis also shows the usefulness of considering pair-wise conservation based on a simple classification of amino acids, in analyzing sequences and structures to find common core regions among homologues.

Possible locally driven folding pathways of TC5b, a 20-residue protein

A novel computational procedure for modeling possible locally driven folding pathways by stepwise elongations of the peptide chain was successfully applied to TC5b, a 20-residue miniprotein. Systematic exploration of the possible locally driven pathways showed that the Trp-cage structure of TC5b could be obtained by stepwise elongation starting from the noncentral local nucleation centers preexisting in the unfolded state of TC5b. The probable locally driven folding pathway starts with folding of alpha-helical fragment 4-9, followed by formation of the proper three-dimensional structure of fragment 4-12, and then 4-18. Accordingly, the Trp-cage-forming interactions emerge successively, first Trp(6)-Pro(12), then Trp(6)-Pro(18), and then Trp(6)-Tyr(3). The Trp-cage-like structures of TC5b found in this study by independent energy calculations are in excellent agreement with the NMR experimental data. The same procedure rationalizes the incomplete Trp-cage formation observed for two analogs of TC5b. Generally, the success of this novel approach is encouraging and provides some justification for the use of computational simulations of locally driven protein folding.