BACKGROUND
clarifying how a protein folds into its unique 3d structure is a very significant yet unsolved problem in molecular biophysics and bioinformatics
<cit> .
in particular, some recent experimental studies have revealed that proteins sharing the same topology can take different folding pathways
<cit> .
ferredoxin-like fold proteins are well-known proteins that fold via different folding pathways.
their topology is composed of  <dig> α helices and  <dig> β strands, and the secondary structure arrangement seems similar to the β/α/β triad motif in flavodoxin
<cit>  or tim-barrel proteins
<cit> .
however, the connectivity of the secondary structures differs.
while flavodoxin or tim-barrel proteins have a parallel β sheet, ferredoxin-like proteins have an anti-parallel β sheet.
therefore, it is hard to explain the ferredoxin-like proteins’ folding behavior only with the formation of β/α/β triads and the interaction among subdomains as in the case of flavodoxin, even if it is true that most of the hydrophobic contacts are composed of ile, leu and val residues as reported in the literature
<cit> .
ferredoxin-like fold proteins are relatively small as shown in figure
<dig>  but they have interesting features in the structural transformation from denatured states to transit or native states.
for example, one ferredoxin-like protein called ribosomal protein s <dig> contains two overlapping foldons, which fold cooperatively, located at different termini, and the overlapping makes this protein fold in a two-state manner as reported by haglund et al.
<cit> .
however, other proteins such as u1a spliceosomal protein or procarboxypeptidase a <dig> fold via the n- or c-terminal foldon as reported by ternström et al.
<cit>  or villegas et al.
<cit> , respectively.
<cit> .
their names are shown above the structures.
α helices and β strands are represented with helices and arrows, respectively.
these protein models are created by visual molecular dynamics  <dig> .1
<cit> .
amino acid sequences of the study proteins.
the filled circles above them denote the location of decadal number residues.
arrows and helices below the sequences indicate the location of secondary structures.
if all the information related to the formation of 3d structures is encoded in the amino acid sequences, we should be able to decode these sequence differences to obtain their folding features.
still, this is a challenging task.
how the folding mechanism of a protein is coded in the amino acid sequence information remains an important issue to be clarified.
there are some bioinformatics approaches for predicting some aspects of folding mechanisms, like folding rate, from the amino acid sequence
<cit> .
nevertheless, it is still difficult to extract more detailed information on the folding mechanism, that is, how each protein folds.
the fact that the topologically equivalent proteins do not always fold via the same folding pathway leads us to the question of whether evolutionarily related proteins really have common folding properties.
evolutionarily related proteins have been observed to be possible to fold via different folding pathways.
for example, in spite of the fact that pdz <dig> and pdz <dig>  both of which contain mainly β sheets, are evolutionarily related and have more than 30% sequence identity, they do not share the same folding mechanism, at least in the early stage of folding
<cit> .
on the other hand, fibronectin and titin, which are evolutionarily unrelated proteins but have the same topology, share the folding mechanism involving four key residues and their peripheral residues
<cit> .
there are also some other studies focusing on the differences in the folding of topologically equivalent proteins
<cit> ; yet, these experiments were performed only for several proteins of each topology and were not applied to a whole family.
in this study, we aim to decode such information for well-studied ferredoxin-like fold
<cit>  proteins by analyzing their amino acid sequences, not only with the preiously mentioned bioinformatics approaches but also with our own analyses.
the methods we apply here are homologous sequence search, phylogenetic analysis and sequence-based analyses by means of inter-residue average distance statistics.
the methods, which are based on the inter-residue average distance statistics
<cit>  using the amino acid sequences as input, have so far provided valuable information on the initial folding segments that play crucial roles in the structural formation in the cases of lysozyme, leghemoglobin, fatty acid-binding protein, azurin, and two ancient tim-barrel proteins
<cit> .
we also apply our methods to some evolutionarily related proteins of four ferredoxin-like fold proteins to examine whether evolutionarily related proteins have common folding properties.
methods
proteins treated in this study
the proteins treated in this study are u1a spliceosomal protein
<cit> , procarboxypeptidase a <dig>
<cit> , ribosomal protein s <dig>
<cit> , and muscle-type acylphosphatase
<cit>  as shown in figure
<dig>  these were selected through the protein folding database  <dig> ,
<cit>  which provides experimental folding data on proteins
<cit> .
we call these proteins our “study proteins”.
the amino acid sequences of these proteins were obtained from the structured region in their pdb files.
their sequence identities are quite low, ranging from  <dig> to 23%.
many experimental studies on these proteins have been performed with respect to the ferredoxin-like fold, and some of these studies suggest the existence of different folding segments
<cit> .
inter-residue average distance statistics
to prepare the statistics for our methods, we calculate the average distance and its standard deviation for each inter-residue pair in  <dig> various proteins, considering the amino acid types and the sequence separation.
for the sequence separation, we simplify the sequence separations k: 1 ~  <dig>  9 ~  <dig>  21 ~  <dig>  31 ~ 40… in terms of the ranges m:  <dig>   <dig>   <dig>  4…, respectively.
the  <dig> representative proteins were carefully chosen so as not to be biased towards some specific structures and have been confirmed to extract the regions corresponding to the structural domains.
because our analysis results are strongly affected by the particular protein datasets used, we chose not to alter the datasets based on the analysis results and to instead use the same datasets as in the first paper on adm  to allow for comparability.
we present the  <dig> proteins .
average distance map analysis
the regions predicted by average distance map  analysis correspond to the regions that tend to be compact in their 3d structure.
we believe that these compact regions might be structured in the early stage of folding.
the adm analysis itself is a method for predicting the location of possible structural units in a protein by analyzing a predicted contact map
<cit>  based on inter-residue average distance statistics.
this map is referred to as an adm and is used to extract standard structural units, such as structural domains or compact regions.
in the adm, any pair of residues with smaller average distance is considered to be in contact more than the other pairs, so the segment with many such pairs is considered to form a structural unit like folding segments by mechanisms such as hydrophobic collapse.
these segments are automatically extracted by analyzing the adm .
in the current study, we use this method to extract compact regions .
for each compact region, the strength as a predicted folding segment is expressed as an η value.
the η value tends to be higher if the corresponding compact regions have many contacts within the region or with the other regions including non-local areas .
contact density has been reported to correctly identify the nucleating subdomains in t <dig> lysozyme and interleukin-1β
<cit> .
studies on flavodoxin-like proteins also suggest a relationship betrween contact density and the folding rate of the corresponding area by showing that low contact density leads to structuring late in folding
<cit> .
thus, it would be reasonable to consider the regions with many contacts  as the region structured in the early stage of folding.
finally, all high-η-value units which do not overlap with other high-η-value units are designated as predicted folding segments, except for units that cover the whole sequence.
when a predicted folding segment covers  <dig> to 100% of the whole sequence length, we conduct an additional search for folding segments overlapping in this unit.
because we could find only two folding segments for each protein, in this study, we call the segment with the higher η value the primary segment and the other, the auxiliary segment.
in additional file
1: figure s <dig>  we compare the adms of the ferredoxin-like proteins with the actual contact map constructed based on the contacts defined later.
comparison of the regions predicted by two adms
suppose that a multiple alignment of homologous sequences in a ferredoxin-like protein is obtained.
since it is convenient to define the similarity of location predicted by any two adms in the multiple alignment, we define the similarity as follows: first, two sequences are chosen from the aligned sequences, as shown in figure
<dig>  second, all sites with a gap in either one or both sequences are removed.
here, “site” refers to the common sequential number in the multiple alignment.
finally, the number of sites that are commonly included in or excluded from the regions predicted by the adms for two given sequences is calculated.
the ratio of this number to the number of all the non-gapped sites is defined as the similarity of location predicted by the two adms.
it should be noted that the similarity calculated by this procedure does not take η values and gapped sites into account: the present method is therefore not suitable in cases where the sequences in an alignment show large gaps.
having said that, the multiple alignments using the sequences obtained in this study contain few small gaps.
for this reason, we can apply this definition of the similarity to the present results.
f-value analysis
we performed additional analyses to determine the location where initial folding events, such as hydrophobic collapse, happen
<cit> .
using other statistical potentials like miyazawa-jernigan
<cit>  and skolnick
<cit>  do may return similar information, but it makes difficult to interpret the results with adm. because there is not only the average distance but also the standard deviation of distance in inter-residue average distance statistics, we expect that the potential used in our f-value analysis to better reproduce the dynamics of the denatured state ensemble compared to the potential based on the contact energy.
in f-value analysis, we use the cα bead model to represent a protein’s structure, as well as the metropolis monte carlo method with the potential energy ϵ
i,j
derived from average distance
ri,j¯ and its standard deviation σ
i,j
.
the bond and dihedral angles of the initial conformation are randomly selected.
the movement in the simulation is done as follows: the bond angle between the residue i and i +  <dig> is bent and rotated randomly from - <dig> to 10° followed by the metropolis judgment to decide if the new conformation is acceptable or not.
within a step, i = 1…n- <dig> is performed, that is, all the bond angles are altered and judged.
the probability density with the potential energy between two residues, p, is hypothesized as being equivalent to the probability density based on the standard gaussian distribution calculated with its average distance and standard deviation,
ρri,j¯,σi,j, as follows:
pϵi,j=ρri,j¯,σi,j
where this equation can be expressed by equation ;
exp-ϵi,jktz=12πσi,jexp-ri,j-ri,j¯22σi,j <dig>
finally we obtain equations  and ;
-ϵi,jkt-lnz=-ln2πσi,j-ri,j-ri,j¯22σi,j <dig>
ϵi,jkt=ri,j-ri,j¯22σi,j2-lnz2πσi,j
where kt is set so that the acceptance ratio is  <dig> .
this potential is designed to sample the ensembles which can reproduce the inter-residue average distance statistics.
from the simulation, the contact frequency, g, for each pair of residues is calculated with sampled structures generated using the potential energy function.
then we normalize the residue contact frequencies, g, in the same range m as follows:
dm=∑μ-v=mg˜μ-v∈m-gμ,vμ-v∈m2Σμ-v∈m
qi,j=gi,ji-j∈m-Σμ-v∈mgμ,vΣμ-v∈mdm
where μ or v is the residue number.
finally, the relative contact frequency, fi, is obtained by summing the normalized contact frequencies, q, from j =  <dig> to n for each residue i, where n is the protein sequence length:
fi=∑jqi,j
the peaks of the plots of the fi- or f-value peaks are thought to be located in the center of many inter-residue contacts, such as a hydrophobic cluster.
therefore, the regions around the peaks are assumed to be important for folding, especially for its initial state.
f-value analysis therefore allows us to estimate the location where a folding initiation occurs, except for the termini with their expected extreme flexibility in the simulation: due to the flexibility, the fi values at the terminal residues become unrealistically high, and this value is then considered not to be true
<cit> .
we performed this simulation with  <dig> steps  <dig> times, calculating the average fi value for residue i.
analysis of evolutionarily conserved residues
evolutionarily conserved residues maintain a protein’s function, contribute to its stability, or relate to its structural formation
<cit> .
therefore, conserved residues that have many contacts with other conserved residues in the native structure are thought to be significant indicators of potential folding segments.
based on this idea, we gather the homologous sequences for each study protein with the basic local alignment search tool
<cit>   and aligned them with the multiple sequence comparison by log-expectation tool
<cit> .
we applied the neighbor-joining method
<cit>  to construct the study protein’s phylogenetic tree for inputting into the phylogenetic analysis by maximum likelihood software package
<cit> .
with paml, for each site without any gap, we can count the number of residue substitutions by using jtt matrices
<cit>  for the substitution matrix and a poisson distribution for the substitution model.
this procedure allows us to estimate the conservation or substitution of a specific residue during evolution based on branch lengths and bifurcations in a phylogenetic tree.
because only the conservation of hydrophobic residues is taken into account in this study, the hydrophobic residues with more than 99% conservation are regarded as conserved residues, that is, we still regard a residue as conserved, when one of the hydrophobic residues a, m, w, l, f, v, i, and y has mutated to another one.
we employed the blast to identify potential homologous sequences.
only sequences that cover the whole sequence of a study protein were selected as homologous sequences based on the blast results.
the blast search identified at least  <dig> homologous sequences for each study protein were obtained .
definition of inter-residue contacts
the shrake-rupley algorithm
<cit>  was used to define a contact by the decrease in the solvent accessible surface area  upon folding.
the reduction in the surface area is calculated by the difference between a sidechain’s sasa in the presence of contacts with other residues and that in the absence of contacts.
in this study, only heavy atoms are considered, and when the decrease in the sasa reaches 27 Å <dig>  the corresponding hydrophobic residue pair is defined as being in contact.
this threshold was determined from the decrease in the sasa when two carbon atoms form a contact, namely,  <dig>  Å <dig>
summary of the experimental results from the literature
figure
<dig> provides a summary of the results reported in the literature from various Φ-value studies which provide information on structured sites in the transition state
<cit> .
we compare the regions predicted by adm with the location of secondary structures with high Φ values.
averaging the Φ values for each secondary structure is a good way to understand the differences in folding mechanisms among a set of proteins with the same topology and this method has been performed by many researchers .
we validated our predictions by comparing them to experimental results interpreted in the same manner.
.
in this study, a few secondary structures with relatively high Φ values are defined as an experimental folding segment.
even though the resolution is somewhat lower because the folding segment is defined by average Φ values, this approach is still similar in concept to the “folding nucleus” first introduced by shaknovich and his colleagues
<cit>  as a set of contacts in denatured states that are considered sufficient and necessary for transitioning or molten globule states to occur as observed by a computational technique.
the formation of these contacts is rate-limiting step and should be done by its transition state.
some later studies support the idea of a folding nucleus by means of Φ-value analysis or combining experimental Φ values with computational techniques
<cit> .
in other words, the folding segment is thought to be relatively structured compared to other regions from the denatured state to transition state.
for example, for the wild-type ribosomal protein s <dig> from thermus thermophilus, one of the ferredoxin-like folds studied by haglund et al.,
<cit>  one folding nucleus is reported to consist of β <dig>  α <dig>  and β <dig> .
however, for the circular permutants prepared by connecting n- and c-termini and disconnecting other loops between neighboring secondary structures, sometimes the folding nucleus noticeably shifts to β <dig>  α <dig>  and β4
<cit> .
therefore haglund et al.
<cit>  found that in s <dig>  there are two competing and overlapping folding nuclei, and the relative magnitude of significance for folding could be perturbed by circular permutation.
.
ternström et al.
<cit>  also performed an experimental study on u1a spliceosomal protein, which also folds through the α <dig> formation and its surrounding secondary structures in an early stage.
however, villegas et al.
<cit>  reported that procarboxypeptidase a <dig>  which has the same topology as s <dig> and u1a, folds through the α <dig> formation and part of β1/β <dig>  its folding segments seem to be in the c-terminus  like the s <dig> permutants , while there are some differences around the β strands.
specifically, there seems to be two tendencies with respect to the location of the folding segments: the n-terminus with α <dig> and the c-terminus with α <dig>  but it remains difficult to say which β strands contribute to folding.
therefore in this study, we simply chose the β strands that have higher average Φ values compared to the values of the α helix with a lower average Φ value than the other α helix.
RESULTS
folding segments predicted by adm analysis
the predicted contact maps and the location of the compact regions, namely, the predicted folding segments are shown in figure
<dig> and summarized in table
<dig>  according to these data, all four proteins have two compact regions, and each region contains one α helix with a couple of β strands.
n or c denote the n- or c-terminal borders of compact regions.
the primary compact regions are shown with n or c in the dominance column.
in the η column, compact regions extended by the 85% rule
<cit>  are identified with a superscript a; compact regions ignored due to covering more than 70% of the entire sequence are identified with a superscript b.
for u1a spliceosomal protein  and muscle-type acylphosphatase  <dig> , the n-terminal compact region has a higher η value than that of the c-terminus, which suggests that the n-terminal region is stable compared to the c-terminal region during the early stage of folding, while procarboxypeptidase a <dig>  shows the opposite trend.
as for ribosomal protein s <dig> , the two compact regions have similar η values, which is interpreted as meaning that both of these regions play equally important roles in structural formation.
it is also notable that except for the case of mtacp, β <dig> is always included in the primary folding segments  for each protein.we are interested in comparing the predicted folding segments with the secondary structures whose average Φ values are high.
figure
5a-d show the results.
according to these figures, the secondary structures within the predicted folding segments often correspond to the secondary structures with high average Φ values.
in figure
<dig>  the positions of the secondary structures with higher average Φ values than those of the α helix  are colored red in the right panel, while in the left panel, the positions of the predicted primary folding segment at the n- or c-terminus are colored yellow or green, respectively.
for s <dig>  however, we color both segments in red or yellow/green, because according to the average Φ values and η values of s <dig> , both n-terminal and c-terminal folding segments are equally important in the formation of its 3d structure.
figures
<dig>
<dig>  and
<dig> indicate that almost all the important secondary structures for folding, as defined by experimental Φ value results, are included in the folding segments predicted by the adms, although β <dig> in mtacp, which shows a relatively high Φ value, is not included in the region predicted by the adm.
evolutionary conservation analysis with f-value analysis - location of predicted hydrophobic clusters
the solid line in figure
<dig> indicates the f-value result, while the broken line shows the smoothed plots of the number of contacts with other conserved hydrophobic residues.
smoothing was performed with a gaussian kernel
<cit> .
a conserved hydrophobic contact means that a pair of conserved hydrophobic residues form a contact.
the locations of the peaks are indicated by single or double daggers for f-value or conserved hydrophobic contacts, respectively.
a number near a dagger indicates the corresponding residue number in a protein.
we follow the pdb system concerning the residue number in a protein.
the secondary structures and conserved hydrophobic residues are shown below the plot.except for several sites, most of the conserved hydrophobic residues are distributed somewhat sparsely but uniformly, which implies that it is hard to extract folding segments from only their amino acid sequences and conservation analyses.
according to figure
<dig>  most of the f-value peaks are close to those of the smoothed line within ±  <dig> residues.
this indicates that f-value peaks, which can be mainly regarded as hydrophobic clusters in the initial nucleation stage, also correspond to the region with many conserved hydrophobic contacts, which are important for the formation of a native structure.
direct comparisons of these regions with high-Φ-value sites are shown in additional file
1: figure s <dig>  for high-Φ-value sites, only the sites with a Φ value higher than the average Φ value of the protein are shown along with each site’s residue type and number.
the peaks of smoothed conserved hydrophobic contacts are marked by double daggers as in figure
<dig>  high-Φ-value sites are found to exist near the peaks of the conserved hydrophobic contacts , suggesting that some of these contacts are responsible for structural formation.
these high-Φ-value sites are also found to exist near the f-value peaks within ±  <dig> residues, as shown in additional file
1: figure s <dig>
folding segments predicted by adms in the homologous proteins of the study proteins
to confirm whether the folding segments are conserved among evolutionarily related proteins, we performed our sequence analysis on the homologues of the four study proteins.
the results of applying adm analyses to these homologues are shown in figure
<dig> this figure denotes the respective multiple alignments of the homologues.
the location of the predicted folding segments are colored dark gray: the brighter the color, the higher the region’s η value.
it can be visually confirmed in figure
7a-c that there are several bands indicating that for most of the homologues the same regions are predicted.
<cit>  is shown.
the location of the secondary structures or the ratio of being included in adms for each site is shown above  or below the results of the adm analyses , respectively.
in figure
<dig>  we ordered the sequences based on the similarity of the location of the regions predicted by the adms.
in the right column, the phylogenetic tree based on an adm similarity matrix and the neighbor-joining method is shown.
another result based on the sequence identity is shown in additional file
1: figure s <dig>  it is difficult to determine the relationship between the location of the folding segments and their evolutionary distance .
however, we can conclude that for u1a, ada2h, and s <dig>  the folding segments themselves are conserved among their homologues, while those of mtacp are not.to represent these common folding segments, we calculate the percentage of residues that are members of the predicted folding segments for each site.
the results are also shown as a histogram colored black in figure
<dig>
in the case of u1a, four secondary structures βαββ at the n-terminus form one strong folding segment for most of the homologues, while the other c-terminal region comprises a weak folding segment.
for ada2h, the c-terminal folding segment βαβ is conservative and strong, and the n-terminal folding segment is conservative but weak.
as for s <dig>  there are many homologues, and they share almost the same folding segments.
one segment consists of βαβ at the c-terminus and the other one consists of βα at the n-terminus.
the dominance of these two folding segments at the termini often differs among the homologues.
it is also notable that for some homologous proteins, the region from β <dig> to the hairpin-loop comprises the weakest folding segment, which forms a β-hairpin with β <dig> in the c-terminal folding segment.
finally, for mtacp, the folding segments are not conservative among its homologous proteins.
however, the locations of the folding segments appear similar to those of s <dig>  ada2h, and u1a.
discussion
the adm analyses of the four proteins predict two compact regions including one α helix and a couple of β strands for each protein.
these regions contain the secondary structures with high average Φ values .
therefore, we consider the predicted compact regions to correspond well to the folding initiation segments as was the case for the other proteins we treated in previous studies, including lysozyme, leghemoglobin, fatty acid binding protein, azurin, and two ancient tim-barrel proteins
<cit> .
according to the η values, mtacp and u1a have the primary predicted folding initiation segment at the n-terminus, whereas ada2h has one at the c-terminus.
on the other hand, the two folding segments of s <dig> have similar η values .
figure
<dig> shows good agreement between the adm predictions and the experimental results; however, the resolution of this analysis is too low to predict the folding mechanisms.
by means of f-value analysis, we increased the resolution of the prediction made by adm, allowing us to identify the regions that would form some hydrophobic clusters.
according to figure
<dig>  almost all the peaks are located on the secondary structures or their edges, and the highest peak is located in the primary folding segment for each protein.
for example, u1a has the primary folding segment at the n terminus from β <dig> to β <dig>  and its highest f-value peak is located in β <dig>  in addition, a peak of f values is located at a peak of the smoothed line of the conserved hydrophobic contacts within ±  <dig> residues, except for the broad peak observed at the c-terminus of ada2h, which contains several minor peaks.
these conserved contacts are thought to play important roles in the structural formation or stabilization of u1a’s native structure
<cit> .
taking these facts into account, the conserved hydrophobic residues near the f-value peaks are considered to be significant for the folding initiation.
the basis for predicting the folding mechanisms from only sequence information is the fact that the regions predicted by adm analysis contain the high-Φ-value residues measured by experiments
<cit>  and that the f-value analysis reflects the conserved hydrophobic contacts.
let us now make inferences regarding the folding mechanisms for proteins based on the results of the adms and f-value analyses.
u1a spliceosomal protein; u1a
as shown in figure
6a, the primary compact region of u1a covers β <dig>  α <dig>  β <dig>  and β <dig>  and each region of α <dig>  β <dig>  and β <dig> contains just one f-value peak.
the auxiliary compact region of u1a ranges from α <dig> to α <dig>  because the auxiliary compact region has a lower η value, it is thought to participate in the structural formation after the primary compact region has been formed.
since ternström et al.
<cit>  suggest that the region from β <dig> to β <dig> is more structured compared to α <dig> and β4
<cit> , we find that our results agree well with their experimental Φ-value analysis.
figure
8a presents the packing formed by conserved hydrophobic residues near the f-value peaks.
the residues that contribute to the hydrophobic packing are represented in the cpk model in this figure.
the regions colored yellow or green correspond to the predicted n- or c-terminal compact regions, respectively.
procarboxypeptidase a2; ada2h
the auxiliary compact region of ada2h at the n-terminus has two secondary structures, β <dig> and α <dig>  whereas the primary compact region at the c-terminus has three secondary structures, β <dig>  α <dig>  and β <dig>  the former has a high f-value peak around α <dig>  indicating that α <dig> is the center of folding within the auxiliary compact region.
on the other hand, the primary compact region has the highest and broadest peak from β <dig> to α <dig>  therefore, we can predict that after β <dig> and α <dig> form a folding segment, β <dig> and α <dig>  pack with this segment and stabilize it.
this prediction can also be validated by figure
3b.
villegas et al.
<cit>  state that the folding segment in this protein consists of α <dig> and its surrounding β strands.
this agrees with our result.
however, we could not confirm any packing between the conserved hydrophobic residues near the f-value peaks in figure
8b within the primary compact region in the native structure as observed in u1a.
this is because the region with the largest broad f-value peak in the c-terminal region seems to have only a few conserved hydrophobic residues as indicated by the smoothed plot of the conserved contacts which shows several minor peaks here.
in this case, the resolution of the f-value line is too low to detect the residues important for folding.
ribosomal protein s6; s6
the relative auxiliary folding segment of s <dig> at the n-terminus contains β <dig> and α <dig>  while the primary folding segment at the c-terminus contains β <dig>  α <dig>  and β <dig>  the η values are quite similar, so we cannot say which region folds more dominantly.
s <dig> has two significant f-value peaks within the predicted folding segments: one is around α <dig> at the n-terminus, and the other is around β <dig> at the c-terminus.
lindberg et al.
<cit>  suggest that the primary folding segment consists of β <dig>  α <dig>  and β <dig> in the early folding stage, and our results reflect this.
figure
8c shows the residue packing near the f-value peaks inside the predicted folding segments.
it is also notable that in the case of s <dig>  there is a highly frustrated region between the c-terminal unstructured coil and the β sheet based on the structure of s6
<cit> .
however, the corresponding c-terminal region does not have any specific contact with other regions in the adm. this is confirmed by its nmr structure .
at least as far as concerning the adm result and the nmr structure, the frustration between the c-terminal unstructured coil and the β sheet does not seem strong.
muscle-type acylphosphatase 2; mtacp
the primary folding segment of this protein, which has a significantly higher η value than the other folding segment, is located at the n-terminus.
this result is the same as the result from selvaraj et al.
<cit>  who suggested the existence of a hydrophobic cluster surrounding α <dig> based on the distribution of contacts.
the primary segment contains β <dig>  α <dig> and β <dig>  whereas the auxiliary segment contains α <dig>  β <dig>  and β <dig>  there are four f-value peaks near α <dig>  β <dig>  β <dig>  and α2; three of them are located in the predicted folding segments, but the peak near β <dig> is not.
therefore we propose that after α <dig> and β <dig> play a role in early structural formation, α <dig> participates in the last structural formation, followed by contributions from β <dig>  this scenario does not seem to fit the results of the average Φ-value analysis
<cit> , which indicate that β <dig> plays a more important role than β <dig> .
according to parrini et al.
<cit> , when β <dig> is forced to join the folding process by a disulfide bond between β <dig> and β <dig>  the folding rate is improved dramatically.
this result suggests that the participation of β <dig> in the folding process is rate limiting and may reflect the present findings.
for this reason, we do not consider the two analyses to conflict with each other.
the inter-segment packing between the conserved residues is represented in figure
8d.
meanwhile, it should be noted that in their analysis, chiti et al.
<cit>  ignore the highest Φ value of the 23rd residue located in α <dig> and then compare the result with that of ada2h.
in this case, the average Φ value of α <dig> is higher than that of α <dig>  indicating that the more structured secondary structure in its transition state is α <dig>  the same as for ada2h.
there are also several studies that suggest that the α <dig> in mtacp is more important than α1
<cit> .
interestingly, one of the studies refers to the large effect on α <dig> induced by point-mutation.
taddei et al.
<cit>  consider the Φ values of α <dig> to be unreliable because, according to their experiment, inducing a point-mutation on α <dig> makes mtacp form amyloid fibrils.
thus, interpreting the folding segment of mtacp is difficult.
our previous simplified go-model simulations reveal that the interactions between the folding segments in the present definition are significant in the formation of transition state ensembles
<cit> .according to the discussion above, the conserved hydrophobic residues among homologues are distributed near f-value peaks , and they seem to be involved in folding.
yet, as we mentioned in the introduction section, there are some studies that have shown that several sequences that fold into the same 3d structure have the same folding process, while other studies have shown that two evolutionarily related sequences with the same 3d structure could have different folding processes.
in this sense, whether all the homologues of our study proteins really have the same folding segments or not would be an interesting question.in the present study, we aim to analyze the conservation of the predicted folding segments in the homologues of a protein by applying adm analysis to them.
the predicted folding segments are highly conserved among the highly homologous proteins  except for mtacp and its homologous proteins, as shown in figures
<dig> and
<dig>  in figure
<dig>  the conservation of the predicted folding segments is summarized as a histogram below each result.
as seen in this figure, the histogram of mtacp is uneven compared to the other histograms.
figure
<dig> depicts this property from another aspect.
the abscissa in figure
<dig> represents the lower limit of the sequence identity for calculating the average adm similarity, that is, an average adm similarity value is calculated using homologues with a sequence identity of more than a lower limit, while the ordinate denotes the average adm similarity.
the doubled line, solid line, thick line, and dotted line correspond to u1a, ada2h, s <dig>  and mtacp, respectively.
while the other three proteins maintain their folding segment similarity of more than 75% even when the sequence identity decreases; only mtacp loses similarity down to 62%.
this result suggests a diversity of folding processes in mtacp compared to those in the other proteins, especially when sequence identity is low.
on the other hand, as we mentioned in the results section, the relationship between adm similarity and sequence identity is not parallel .
this is an unexpected result: we expected that the more similar the sequence identity is, the more similar the protein folding segments are.
yet the present results suggest that a property related to folding segments is conserved more than sequence identity.
summarizing the discussion above, there are mainly two situations.
one of them clearly comprises a main large folding segment around α <dig>  like in u1a.
the other situation comprises complex folding segments in which, one of the α helices and its surrounding β strands play a key role at first, immediately followed by the other helix and its surrounding strands, as in ada2h or s <dig>  the homologous proteins of mtacp have either property: some of them have folding segments similar to those of u1a proteins, and some others have folding segments similar to those of s <dig> or ada2h proteins.
this implies that mtacp and its homologous proteins do not seem to have any common or rigid folding segments.
CONCLUSIONS
the secondary structures that are thought to play important roles in folding as revealed by their average Φ values correspond to the folding segments predicted by adm analyses at least for the proteins treated in this study, as was the case in our previous studies
<cit> .
there are two predicted folding segments at the termini of each protein; however, which segment is primary is completely determined on a case-by-case basis.
this tendency was also in good agreement with the experimental results for the present four study proteins.
some of the conserved hydrophobic contacts considered to play important roles in structural formation
<cit>  are located near the f-value peaks.
therefore, we can predict the folding mechanisms by extracting the conserved hydrophobic residues near them.
for the four proteins we studied above, we conclude that we succeeded in predicting their folding mechanisms correctly from only their sequences.
according to the adm results of the homologues, their folding segments seem to be conserved, especially when the sequence identity is above 80%.
below this level, only mtacp represents a diversity of folding segments, whereas the other three proteins show high conservations.
our findings suggest that it should be possible to predict the folding mechanisms or properties of many other kinds of proteins from only the amino acid sequences by means of our adm analysis and f-value analysis.
abbreviations
adm: average distance map; blast: basic local alignment search tool; sasa: solvent accessible surface area; pdb: protein data bank; muscle: multiple sequence comparison by log expectation; paml: phylogenetic analysis with maximum likelihood; jtt: jones taylor thornton.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
tk conceived and designed the basis of this study.
mm performed all the calculations and data analysis.
tk and mm wrote the manuscript.
both authors read and approved the final manuscript.
supplementary material
additional file 1
details of the adm analysis and optional results are provided.
click here for file
acknowledgments
the authors wish to acknowledge asst.
prof. yosuke kawai.
this work is supported by jsps kakenhi to mm, grant-in-aid for jsps fellows, grant number  <dig>  one of the authors  expresses his gratitude to the ministry of education, culture, sports, science, and technology for the support of the present work through a program for strategic research foundations at private universities, 2010– <dig>  .

 clarifying how a protein folds into its unique 3d structure is a very significant yet unsolved problem in molecular biophysics and bioinformatics. in particular, some recent experimental studies have revealed that proteins sharing the same topology can take different folding pathways. ferredoxin-like fold proteins are well-known proteins that fold via different folding pathways. their topology is composed of α helices and β strands, and the secondary structure arrangement seems similar to the β/α/β triad motif in flavodoxin or tim-barrel proteins. however, the connectivity of the secondary structures differs. while flavodoxin or tim-barrel proteins have a parallel β sheet, ferredoxin-like proteins have an anti-parallel β sheet. therefore, it is hard to explain the ferredoxin-like proteins’ folding behavior only with the formation of β/α/β triads and the interaction among subdomains as in the case of flavodoxin, even if it is true that most of the hydrophobic contacts are composed of ile, leu and val residues as reported in the literature. ferredoxin-like fold proteins are relatively small as shown in figure  but they have interesting features in the structural transformation from denatured states to transit or native states. for example, one ferredoxin-like protein called ribosomal protein s contains two overlapping foldons, which fold cooperatively, located at different termini, and the overlapping makes this protein fold in a two-state manner as reported by haglund et al.. however, other proteins such as u1a spliceosomal protein or procarboxypeptidase a fold via the n- or c-terminal foldon as reported by ternström et al. or villegas et al., respectively.. their names are shown above the structures. α helices and β strands are represented with helices and arrows, respectively. these protein models are created by visual molecular dynamics.1. amino acid sequences of the study proteins. the filled circles above them denote the location of decadal number residues. arrows and helices below the sequences indicate the location of secondary structures. if all the information related to the formation of 3d structures is encoded in the amino acid sequences, we should be able to decode these sequence differences to obtain their folding features. still, this is a challenging task. how the folding mechanism of a protein is coded in the amino acid sequence information remains an important issue to be clarified. there are some bioinformatics approaches for predicting some aspects of folding mechanisms, like folding rate, from the amino acid sequence. nevertheless, it is still difficult to extract more detailed information on the folding mechanism, that is, how each protein folds. the fact that the topologically equivalent proteins do not always fold via the same folding pathway leads us to the question of whether evolutionarily related proteins really have common folding properties. evolutionarily related proteins have been observed to be possible to fold via different folding pathways. for example, in spite of the fact that pdz and pdz both of which contain mainly β sheets, are evolutionarily related and have more than 30% sequence identity, they do not share the same folding mechanism, at least in the early stage of folding. on the other hand, fibronectin and titin, which are evolutionarily unrelated proteins but have the same topology, share the folding mechanism involving four key residues and their peripheral residues. there are also some other studies focusing on the differences in the folding of topologically equivalent proteins; yet, these experiments were performed only for several proteins of each topology and were not applied to a whole family. in this study, we aim to decode such information for well-studied ferredoxin-like fold proteins by analyzing their amino acid sequences, not only with the preiously mentioned bioinformatics approaches but also with our own analyses. the methods we apply here are homologous sequence search, phylogenetic analysis and sequence-based analyses by means of inter-residue average distance statistics. the methods, which are based on the inter-residue average distance statistics using the amino acid sequences as input, have so far provided valuable information on the initial folding segments that play crucial roles in the structural formation in the cases of lysozyme, leghemoglobin, fatty acid-binding protein, azurin, and two ancient tim-barrel proteins. we also apply our methods to some evolutionarily related proteins of four ferredoxin-like fold proteins to examine whether evolutionarily related proteins have common folding properties. methods proteins treated in this study the proteins treated in this study are u1a spliceosomal protein, procarboxypeptidase a, ribosomal protein s, and muscle-type acylphosphatase as shown in figure  these were selected through the protein folding database, which provides experimental folding data on proteins. we call these proteins our “study proteins”. the amino acid sequences of these proteins were obtained from the structured region in their pdb files. their sequence identities are quite low, ranging from to 23%. many experimental studies on these proteins have been performed with respect to the ferredoxin-like fold, and some of these studies suggest the existence of different folding segments. inter-residue average distance statistics to prepare the statistics for our methods, we calculate the average distance and its standard deviation for each inter-residue pair in various proteins, considering the amino acid types and the sequence separation. for the sequence separation, we simplify the sequence separations k: 1 ~  9 ~  21 ~  31 ~ 40… in terms of the ranges m: 4…, respectively. the representative proteins were carefully chosen so as not to be biased towards some specific structures and have been confirmed to extract the regions corresponding to the structural domains. because our analysis results are strongly affected by the particular protein datasets used, we chose not to alter the datasets based on the analysis results and to instead use the same datasets as in the first paper on adm to allow for comparability. we present the proteins. average distance map analysis the regions predicted by average distance map analysis correspond to the regions that tend to be compact in their 3d structure. we believe that these compact regions might be structured in the early stage of folding. the adm analysis itself is a method for predicting the location of possible structural units in a protein by analyzing a predicted contact map based on inter-residue average distance statistics. this map is referred to as an adm and is used to extract standard structural units, such as structural domains or compact regions. in the adm, any pair of residues with smaller average distance is considered to be in contact more than the other pairs, so the segment with many such pairs is considered to form a structural unit like folding segments by mechanisms such as hydrophobic collapse. these segments are automatically extracted by analyzing the adm. in the current study, we use this method to extract compact regions. for each compact region, the strength as a predicted folding segment is expressed as an η value. the η value tends to be higher if the corresponding compact regions have many contacts within the region or with the other regions including non-local areas. contact density has been reported to correctly identify the nucleating subdomains in t lysozyme and interleukin-1β. studies on flavodoxin-like proteins also suggest a relationship betrween contact density and the folding rate of the corresponding area by showing that low contact density leads to structuring late in folding. thus, it would be reasonable to consider the regions with many contacts as the region structured in the early stage of folding. finally, all high-η-value units which do not overlap with other high-η-value units are designated as predicted folding segments, except for units that cover the whole sequence. when a predicted folding segment covers to 100% of the whole sequence length, we conduct an additional search for folding segments overlapping in this unit. because we could find only two folding segments for each protein, in this study, we call the segment with the higher η value the primary segment and the other, the auxiliary segment. in additional file 1: figure s we compare the adms of the ferredoxin-like proteins with the actual contact map constructed based on the contacts defined later. comparison of the regions predicted by two adms suppose that a multiple alignment of homologous sequences in a ferredoxin-like protein is obtained. since it is convenient to define the similarity of location predicted by any two adms in the multiple alignment, we define the similarity as follows: first, two sequences are chosen from the aligned sequences, as shown in figure  second, all sites with a gap in either one or both sequences are removed. here, “site” refers to the common sequential number in the multiple alignment. finally, the number of sites that are commonly included in or excluded from the regions predicted by the adms for two given sequences is calculated. the ratio of this number to the number of all the non-gapped sites is defined as the similarity of location predicted by the two adms. it should be noted that the similarity calculated by this procedure does not take η values and gapped sites into account: the present method is therefore not suitable in cases where the sequences in an alignment show large gaps. having said that, the multiple alignments using the sequences obtained in this study contain few small gaps. for this reason, we can apply this definition of the similarity to the present results. f-value analysis we performed additional analyses to determine the location where initial folding events, such as hydrophobic collapse, happen. using other statistical potentials like miyazawa-jernigan and skolnick do may return similar information, but it makes difficult to interpret the results with adm. because there is not only the average distance but also the standard deviation of distance in inter-residue average distance statistics, we expect that the potential used in our f-value analysis to better reproduce the dynamics of the denatured state ensemble compared to the potential based on the contact energy. in f-value analysis, we use the cα bead model to represent a protein’s structure, as well as the metropolis monte carlo method with the potential energy ϵ i,j derived from average distance ri,j¯ and its standard deviation σ i,j. the bond and dihedral angles of the initial conformation are randomly selected. the movement in the simulation is done as follows: the bond angle between the residue i and i +  is bent and rotated randomly from - to 10° followed by the metropolis judgment to decide if the new conformation is acceptable or not. within a step, i = 1…n- is performed, that is, all the bond angles are altered and judged. the probability density with the potential energy between two residues, p, is hypothesized as being equivalent to the probability density based on the standard gaussian distribution calculated with its average distance and standard deviation, ρri,j¯,σi,j, as follows: pϵi,j=ρri,j¯,σi,j where this equation can be expressed by equation; exp-ϵi,jktz=12πσi,jexp-ri,j-ri,j¯22σi,j finally we obtain equations and; -ϵi,jkt-lnz=-ln2πσi,j-ri,j-ri,j¯22σi,j ϵi,jkt=ri,j-ri,j¯22σi,j2-lnz2πσi,j where kt is set so that the acceptance ratio is. this potential is designed to sample the ensembles which can reproduce the inter-residue average distance statistics. from the simulation, the contact frequency, g, for each pair of residues is calculated with sampled structures generated using the potential energy function. then we normalize the residue contact frequencies, g, in the same range m as follows: dm=∑μ-v=mg˜μ-v∈m-gμ,vμ-v∈m2Σμ-v∈m qi,j=gi,ji-j∈m-Σμ-v∈mgμ,vΣμ-v∈mdm where μ or v is the residue number. finally, the relative contact frequency, fi, is obtained by summing the normalized contact frequencies, q, from j =  to n for each residue i, where n is the protein sequence length: fi=∑jqi,j the peaks of the plots of the fi- or f-value peaks are thought to be located in the center of many inter-residue contacts, such as a hydrophobic cluster. therefore, the regions around the peaks are assumed to be important for folding, especially for its initial state. f-value analysis therefore allows us to estimate the location where a folding initiation occurs, except for the termini with their expected extreme flexibility in the simulation: due to the flexibility, the fi values at the terminal residues become unrealistically high, and this value is then considered not to be true. we performed this simulation with steps times, calculating the average fi value for residue i. analysis of evolutionarily conserved residues evolutionarily conserved residues maintain a protein’s function, contribute to its stability, or relate to its structural formation. therefore, conserved residues that have many contacts with other conserved residues in the native structure are thought to be significant indicators of potential folding segments. based on this idea, we gather the homologous sequences for each study protein with the basic local alignment search tool and aligned them with the multiple sequence comparison by log-expectation tool. we applied the neighbor-joining method to construct the study protein’s phylogenetic tree for inputting into the phylogenetic analysis by maximum likelihood software package. with paml, for each site without any gap, we can count the number of residue substitutions by using jtt matrices for the substitution matrix and a poisson distribution for the substitution model. this procedure allows us to estimate the conservation or substitution of a specific residue during evolution based on branch lengths and bifurcations in a phylogenetic tree. because only the conservation of hydrophobic residues is taken into account in this study, the hydrophobic residues with more than 99% conservation are regarded as conserved residues, that is, we still regard a residue as conserved, when one of the hydrophobic residues a, m, w, l, f, v, i, and y has mutated to another one. we employed the blast to identify potential homologous sequences. only sequences that cover the whole sequence of a study protein were selected as homologous sequences based on the blast results. the blast search identified at least homologous sequences for each study protein were obtained. definition of inter-residue contacts the shrake-rupley algorithm was used to define a contact by the decrease in the solvent accessible surface area upon folding. the reduction in the surface area is calculated by the difference between a sidechain’s sasa in the presence of contacts with other residues and that in the absence of contacts. in this study, only heavy atoms are considered, and when the decrease in the sasa reaches 27 Å the corresponding hydrophobic residue pair is defined as being in contact. this threshold was determined from the decrease in the sasa when two carbon atoms form a contact, namely,  Å summary of the experimental results from the literature figure  provides a summary of the results reported in the literature from various Φ-value studies which provide information on structured sites in the transition state. we compare the regions predicted by adm with the location of secondary structures with high Φ values. averaging the Φ values for each secondary structure is a good way to understand the differences in folding mechanisms among a set of proteins with the same topology and this method has been performed by many researchers. we validated our predictions by comparing them to experimental results interpreted in the same manner.. in this study, a few secondary structures with relatively high Φ values are defined as an experimental folding segment. even though the resolution is somewhat lower because the folding segment is defined by average Φ values, this approach is still similar in concept to the “folding nucleus” first introduced by shaknovich and his colleagues as a set of contacts in denatured states that are considered sufficient and necessary for transitioning or molten globule states to occur as observed by a computational technique. the formation of these contacts is rate-limiting step and should be done by its transition state. some later studies support the idea of a folding nucleus by means of Φ-value analysis or combining experimental Φ values with computational techniques. in other words, the folding segment is thought to be relatively structured compared to other regions from the denatured state to transition state. for example, for the wild-type ribosomal protein s from thermus thermophilus, one of the ferredoxin-like folds studied by haglund et al., one folding nucleus is reported to consist of β α and β. however, for the circular permutants prepared by connecting n- and c-termini and disconnecting other loops between neighboring secondary structures, sometimes the folding nucleus noticeably shifts to β α and β4. therefore haglund et al. found that in s there are two competing and overlapping folding nuclei, and the relative magnitude of significance for folding could be perturbed by circular permutation.. ternström et al. also performed an experimental study on u1a spliceosomal protein, which also folds through the α formation and its surrounding secondary structures in an early stage. however, villegas et al. reported that procarboxypeptidase a which has the same topology as s and u1a, folds through the α formation and part of β1/β its folding segments seem to be in the c-terminus like the s permutants, while there are some differences around the β strands. specifically, there seems to be two tendencies with respect to the location of the folding segments: the n-terminus with α and the c-terminus with α but it remains difficult to say which β strands contribute to folding. therefore in this study, we simply chose the β strands that have higher average Φ values compared to the values of the α helix with a lower average Φ value than the other α helix.  folding segments predicted by adm analysis the predicted contact maps and the location of the compact regions, namely, the predicted folding segments are shown in figure  and summarized in table  according to these data, all four proteins have two compact regions, and each region contains one α helix with a couple of β strands. n or c denote the n- or c-terminal borders of compact regions. the primary compact regions are shown with n or c in the dominance column. in the η column, compact regions extended by the 85% rule are identified with a superscript a; compact regions ignored due to covering more than 70% of the entire sequence are identified with a superscript b. for u1a spliceosomal protein and muscle-type acylphosphatase, the n-terminal compact region has a higher η value than that of the c-terminus, which suggests that the n-terminal region is stable compared to the c-terminal region during the early stage of folding, while procarboxypeptidase a shows the opposite trend. as for ribosomal protein s, the two compact regions have similar η values, which is interpreted as meaning that both of these regions play equally important roles in structural formation. it is also notable that except for the case of mtacp, β is always included in the primary folding segments for each protein.we are interested in comparing the predicted folding segments with the secondary structures whose average Φ values are high. figure  5a-d show the results. according to these figures, the secondary structures within the predicted folding segments often correspond to the secondary structures with high average Φ values. in figure  the positions of the secondary structures with higher average Φ values than those of the α helix are colored red in the right panel, while in the left panel, the positions of the predicted primary folding segment at the n- or c-terminus are colored yellow or green, respectively. for s however, we color both segments in red or yellow/green, because according to the average Φ values and η values of s, both n-terminal and c-terminal folding segments are equally important in the formation of its 3d structure. figures  and indicate that almost all the important secondary structures for folding, as defined by experimental Φ value results, are included in the folding segments predicted by the adms, although β in mtacp, which shows a relatively high Φ value, is not included in the region predicted by the adm. evolutionary conservation analysis with f-value analysis - location of predicted hydrophobic clusters the solid line in figure  indicates the f-value result, while the broken line shows the smoothed plots of the number of contacts with other conserved hydrophobic residues. smoothing was performed with a gaussian kernel. a conserved hydrophobic contact means that a pair of conserved hydrophobic residues form a contact. the locations of the peaks are indicated by single or double daggers for f-value or conserved hydrophobic contacts, respectively. a number near a dagger indicates the corresponding residue number in a protein. we follow the pdb system concerning the residue number in a protein. the secondary structures and conserved hydrophobic residues are shown below the plot.except for several sites, most of the conserved hydrophobic residues are distributed somewhat sparsely but uniformly, which implies that it is hard to extract folding segments from only their amino acid sequences and conservation analyses. according to figure  most of the f-value peaks are close to those of the smoothed line within ±  residues. this indicates that f-value peaks, which can be mainly regarded as hydrophobic clusters in the initial nucleation stage, also correspond to the region with many conserved hydrophobic contacts, which are important for the formation of a native structure. direct comparisons of these regions with high-Φ-value sites are shown in additional file 1: figure s for high-Φ-value sites, only the sites with a Φ value higher than the average Φ value of the protein are shown along with each site’s residue type and number. the peaks of smoothed conserved hydrophobic contacts are marked by double daggers as in figure  high-Φ-value sites are found to exist near the peaks of the conserved hydrophobic contacts, suggesting that some of these contacts are responsible for structural formation. these high-Φ-value sites are also found to exist near the f-value peaks within ±  residues, as shown in additional file 1: figure s folding segments predicted by adms in the homologous proteins of the study proteins to confirm whether the folding segments are conserved among evolutionarily related proteins, we performed our sequence analysis on the homologues of the four study proteins. the results of applying adm analyses to these homologues are shown in figure  this figure denotes the respective multiple alignments of the homologues. the location of the predicted folding segments are colored dark gray: the brighter the color, the higher the region’s η value. it can be visually confirmed in figure  7a-c that there are several bands indicating that for most of the homologues the same regions are predicted. is shown. the location of the secondary structures or the ratio of being included in adms for each site is shown above or below the results of the adm analyses, respectively. in figure  we ordered the sequences based on the similarity of the location of the regions predicted by the adms. in the right column, the phylogenetic tree based on an adm similarity matrix and the neighbor-joining method is shown. another result based on the sequence identity is shown in additional file 1: figure s it is difficult to determine the relationship between the location of the folding segments and their evolutionary distance. however, we can conclude that for u1a, ada2h, and s the folding segments themselves are conserved among their homologues, while those of mtacp are not.to represent these common folding segments, we calculate the percentage of residues that are members of the predicted folding segments for each site. the results are also shown as a histogram colored black in figure  in the case of u1a, four secondary structures βαββ at the n-terminus form one strong folding segment for most of the homologues, while the other c-terminal region comprises a weak folding segment. for ada2h, the c-terminal folding segment βαβ is conservative and strong, and the n-terminal folding segment is conservative but weak. as for s there are many homologues, and they share almost the same folding segments. one segment consists of βαβ at the c-terminus and the other one consists of βα at the n-terminus. the dominance of these two folding segments at the termini often differs among the homologues. it is also notable that for some homologous proteins, the region from β to the hairpin-loop comprises the weakest folding segment, which forms a β-hairpin with β in the c-terminal folding segment. finally, for mtacp, the folding segments are not conservative among its homologous proteins. however, the locations of the folding segments appear similar to those of s ada2h, and u1a. discussion the adm analyses of the four proteins predict two compact regions including one α helix and a couple of β strands for each protein. these regions contain the secondary structures with high average Φ values. therefore, we consider the predicted compact regions to correspond well to the folding initiation segments as was the case for the other proteins we treated in previous studies, including lysozyme, leghemoglobin, fatty acid binding protein, azurin, and two ancient tim-barrel proteins. according to the η values, mtacp and u1a have the primary predicted folding initiation segment at the n-terminus, whereas ada2h has one at the c-terminus. on the other hand, the two folding segments of s have similar η values. figure  shows good agreement between the adm predictions and the experimental results; however, the resolution of this analysis is too low to predict the folding mechanisms. by means of f-value analysis, we increased the resolution of the prediction made by adm, allowing us to identify the regions that would form some hydrophobic clusters. according to figure  almost all the peaks are located on the secondary structures or their edges, and the highest peak is located in the primary folding segment for each protein. for example, u1a has the primary folding segment at the n terminus from β to β and its highest f-value peak is located in β in addition, a peak of f values is located at a peak of the smoothed line of the conserved hydrophobic contacts within ±  residues, except for the broad peak observed at the c-terminus of ada2h, which contains several minor peaks. these conserved contacts are thought to play important roles in the structural formation or stabilization of u1a’s native structure. taking these facts into account, the conserved hydrophobic residues near the f-value peaks are considered to be significant for the folding initiation. the basis for predicting the folding mechanisms from only sequence information is the fact that the regions predicted by adm analysis contain the high-Φ-value residues measured by experiments and that the f-value analysis reflects the conserved hydrophobic contacts. let us now make inferences regarding the folding mechanisms for proteins based on the results of the adms and f-value analyses. u1a spliceosomal protein; u1a as shown in figure  6a, the primary compact region of u1a covers β α β and β and each region of α β and β contains just one f-value peak. the auxiliary compact region of u1a ranges from α to α because the auxiliary compact region has a lower η value, it is thought to participate in the structural formation after the primary compact region has been formed. since ternström et al. suggest that the region from β to β is more structured compared to α and β4, we find that our results agree well with their experimental Φ-value analysis. figure  8a presents the packing formed by conserved hydrophobic residues near the f-value peaks. the residues that contribute to the hydrophobic packing are represented in the cpk model in this figure. the regions colored yellow or green correspond to the predicted n- or c-terminal compact regions, respectively. procarboxypeptidase a2; ada2h the auxiliary compact region of ada2h at the n-terminus has two secondary structures, β and α whereas the primary compact region at the c-terminus has three secondary structures, β α and β the former has a high f-value peak around α indicating that α is the center of folding within the auxiliary compact region. on the other hand, the primary compact region has the highest and broadest peak from β to α therefore, we can predict that after β and α form a folding segment, β and α pack with this segment and stabilize it. this prediction can also be validated by figure  3b. villegas et al. state that the folding segment in this protein consists of α and its surrounding β strands. this agrees with our result. however, we could not confirm any packing between the conserved hydrophobic residues near the f-value peaks in figure  8b within the primary compact region in the native structure as observed in u1a. this is because the region with the largest broad f-value peak in the c-terminal region seems to have only a few conserved hydrophobic residues as indicated by the smoothed plot of the conserved contacts which shows several minor peaks here. in this case, the resolution of the f-value line is too low to detect the residues important for folding. ribosomal protein s6; s6 the relative auxiliary folding segment of s at the n-terminus contains β and α while the primary folding segment at the c-terminus contains β α and β the η values are quite similar, so we cannot say which region folds more dominantly. s has two significant f-value peaks within the predicted folding segments: one is around α at the n-terminus, and the other is around β at the c-terminus. lindberg et al. suggest that the primary folding segment consists of β α and β in the early folding stage, and our results reflect this. figure  8c shows the residue packing near the f-value peaks inside the predicted folding segments. it is also notable that in the case of s there is a highly frustrated region between the c-terminal unstructured coil and the β sheet based on the structure of s6. however, the corresponding c-terminal region does not have any specific contact with other regions in the adm. this is confirmed by its nmr structure. at least as far as concerning the adm result and the nmr structure, the frustration between the c-terminal unstructured coil and the β sheet does not seem strong. muscle-type acylphosphatase 2; mtacp the primary folding segment of this protein, which has a significantly higher η value than the other folding segment, is located at the n-terminus. this result is the same as the result from selvaraj et al. who suggested the existence of a hydrophobic cluster surrounding α based on the distribution of contacts. the primary segment contains β α and β whereas the auxiliary segment contains α β and β there are four f-value peaks near α β β and α2; three of them are located in the predicted folding segments, but the peak near β is not. therefore we propose that after α and β play a role in early structural formation, α participates in the last structural formation, followed by contributions from β this scenario does not seem to fit the results of the average Φ-value analysis, which indicate that β plays a more important role than β. according to parrini et al., when β is forced to join the folding process by a disulfide bond between β and β the folding rate is improved dramatically. this result suggests that the participation of β in the folding process is rate limiting and may reflect the present findings. for this reason, we do not consider the two analyses to conflict with each other. the inter-segment packing between the conserved residues is represented in figure  8d. meanwhile, it should be noted that in their analysis, chiti et al. ignore the highest Φ value of the 23rd residue located in α and then compare the result with that of ada2h. in this case, the average Φ value of α is higher than that of α indicating that the more structured secondary structure in its transition state is α the same as for ada2h. there are also several studies that suggest that the α in mtacp is more important than α1. interestingly, one of the studies refers to the large effect on α induced by point-mutation. taddei et al. consider the Φ values of α to be unreliable because, according to their experiment, inducing a point-mutation on α makes mtacp form amyloid fibrils. thus, interpreting the folding segment of mtacp is difficult. our previous simplified go-model simulations reveal that the interactions between the folding segments in the present definition are significant in the formation of transition state ensembles.according to the discussion above, the conserved hydrophobic residues among homologues are distributed near f-value peaks, and they seem to be involved in folding. yet, as we mentioned in the introduction section, there are some studies that have shown that several sequences that fold into the same 3d structure have the same folding process, while other studies have shown that two evolutionarily related sequences with the same 3d structure could have different folding processes. in this sense, whether all the homologues of our study proteins really have the same folding segments or not would be an interesting question.in the present study, we aim to analyze the conservation of the predicted folding segments in the homologues of a protein by applying adm analysis to them. the predicted folding segments are highly conserved among the highly homologous proteins except for mtacp and its homologous proteins, as shown in figures  and in figure  the conservation of the predicted folding segments is summarized as a histogram below each result. as seen in this figure, the histogram of mtacp is uneven compared to the other histograms. figure  depicts this property from another aspect. the abscissa in figure  represents the lower limit of the sequence identity for calculating the average adm similarity, that is, an average adm similarity value is calculated using homologues with a sequence identity of more than a lower limit, while the ordinate denotes the average adm similarity. the doubled line, solid line, thick line, and dotted line correspond to u1a, ada2h, s and mtacp, respectively. while the other three proteins maintain their folding segment similarity of more than 75% even when the sequence identity decreases; only mtacp loses similarity down to 62%. this result suggests a diversity of folding processes in mtacp compared to those in the other proteins, especially when sequence identity is low. on the other hand, as we mentioned in the results section, the relationship between adm similarity and sequence identity is not parallel. this is an unexpected result: we expected that the more similar the sequence identity is, the more similar the protein folding segments are. yet the present results suggest that a property related to folding segments is conserved more than sequence identity. summarizing the discussion above, there are mainly two situations. one of them clearly comprises a main large folding segment around α like in u1a. the other situation comprises complex folding segments in which, one of the α helices and its surrounding β strands play a key role at first, immediately followed by the other helix and its surrounding strands, as in ada2h or s the homologous proteins of mtacp have either property: some of them have folding segments similar to those of u1a proteins, and some others have folding segments similar to those of s or ada2h proteins. this implies that mtacp and its homologous proteins do not seem to have any common or rigid folding segments.  the secondary structures that are thought to play important roles in folding as revealed by their average Φ values correspond to the folding segments predicted by adm analyses at least for the proteins treated in this study, as was the case in our previous studies. there are two predicted folding segments at the termini of each protein; however, which segment is primary is completely determined on a case-by-case basis. this tendency was also in good agreement with the experimental results for the present four study proteins. some of the conserved hydrophobic contacts considered to play important roles in structural formation are located near the f-value peaks. therefore, we can predict the folding mechanisms by extracting the conserved hydrophobic residues near them. for the four proteins we studied above, we conclude that we succeeded in predicting their folding mechanisms correctly from only their sequences. according to the adm results of the homologues, their folding segments seem to be conserved, especially when the sequence identity is above 80%. below this level, only mtacp represents a diversity of folding segments, whereas the other three proteins show high conservations. our findings suggest that it should be possible to predict the folding mechanisms or properties of many other kinds of proteins from only the amino acid sequences by means of our adm analysis and f-value analysis. abbreviations adm: average distance map; blast: basic local alignment search tool; sasa: solvent accessible surface area; pdb: protein data bank; muscle: multiple sequence comparison by log expectation; paml: phylogenetic analysis with maximum likelihood; jtt: jones taylor thornton. competing interests the authors declare that they have no competing interests. authors’ contributions tk conceived and designed the basis of this study. mm performed all the calculations and data analysis. tk and mm wrote the manuscript. both authors read and approved the final manuscript. supplementary material additional file 1 details of the adm analysis and optional results are provided. click here for file acknowledgments the authors wish to acknowledge asst. prof. yosuke kawai. this work is supported by jsps kakenhi to mm, grant-in-aid for jsps fellows, grant number one of the authors expresses his gratitude to the ministry of education, culture, sports, science, and technology for the support of the present work through a program for strategic research foundations at private universities, 2010–.
despite the simplicity of our methods, we are able to correctly determine the experimentally identified folding segments by predicting the location of the compact regions considered to play an important role in structural formation. while some studies have shown that the 3d protein structures are more conservative than their amino acid sequences, other experimental studies have shown that even if two proteins share the same topology, they may have different folding pathways. there are many studies investigating this issue with molecular dynamics or go-like model simulations, however, one should be able to obtain the same information by analyzing the proteins’ amino acid sequences, if the sequences contain all the information about the 3d structures. we focus on proteins with a ferredoxin-like fold, which has a characteristic topology.
 while some studies have shown that the 3d protein structures are more conservative than their amino acid sequences, other experimental studies have shown that even if two proteins share the same topology, they may have different folding pathways. there are many studies investigating this issue with molecular dynamics or go-like model simulations, however, one should be able to obtain the same information by analyzing the proteins’ amino acid sequences, if the sequences contain all the information about the 3d structures. in this study, we use information about protein sequences to predict the location of their folding segments. we focus on proteins with a ferredoxin-like fold, which has a characteristic topology. some of these proteins have different folding segments.  despite the simplicity of our methods, we are able to correctly determine the experimentally identified folding segments by predicting the location of the compact regions considered to play an important role in structural formation. we also apply our sequence analyses to some homologues of each protein and confirm that there are highly conserved folding segments despite the homologues’ sequence diversity. these homologues have similar folding segments even though the homology of two proteins’ sequences is not so high.  our analyses have proven useful for investigating the common or different folding features of the proteins studied. folding initiation segment predictionsequence analysisinter-residue average distance statisticsevolutionarily conserved foldingribosomal protein s6procarboxypeptidase a2u1a spliceosomal proteinmt-acylphosphatase 