Folding rate prediction using complex network analysis for proteins with two- and three-state folding kinetics

It is a challenging task to investigate the different in- fluence of long-range and short-range interactions on two-state and three-state folding kinetics of protein. The networks of the 30 two-state proteins and 15 three-state proteins were constructed by complex networks analysis at three length scales: Protein Contact Networks, Long-range Interaction Networks and Short-range Interaction Networks. To uncover the relationship between structural properties and folding kinetics of the proteins, the correlations of protein network parameters with protein folding rate and topology parameters contact order were analyzed. The results show that Protein Contact Networks and Short-range Interaction Networks (for both two-state and three-state proteins) exhibit the “small-world” property and Long-range Interaction networks indicate “scale-free” behavior. Our results further indicate that all Protein Contact Networks and Short- range Interaction networks are assortative type. While some of Long-range Interaction Networks are of assortative type, the others are of disassortative type. For two-state proteins, the clustering coefficients of Short-range Interaction Networks show prominent correlation with folding rate and contact order. The assortativity coefficients of Short-range Interaction Networks also show remarkable correlation with folding rate and contact order. Similar correlations exist in Protein Contact Networks of three-state proteins. For two-state proteins, the correlation between contact order and folding rate is determined by the numbers of local contacts. Short- range interactions play a key role in determining the connecting trend among amino acids and they impact the folding rate of two-state proteins directly. For three-state proteins, the folding rate is determined by short-range and long-range interactions among residues together.