Distance-constraint approach to protein folding. I. Statistical analysis of protein conformations in terms of distances between residues

A statistical analysis of protein conformations in terms of the distance between residues, represented by their C atoms, is presented. We consider four factors that contribute to the determination of the distanced _i,i+k between a given pair ofith and(i+k)th residues in the native conformation of a globular protein: (1) the distancek along the chain, (2) the size of the protein, (3) the conformational states of theith to(i+k)th residues, and (4) the amino acid types of the and(i+k)th residues. In order to account for the dependence on the distancek along the chain, the statistics are taken for three ranges, viz., short, medium, and long ranges (k8; 9k20; andk21; respectively). In the statistics of short-range distances, a mean distanceD _k and its standard deviationS _k are calculated for each value ofk, with and without taking into account the conformational states of all residues fromi toi+k (factors 1 and 3). As an Appendix, the relations for converting from the distances between residues into other conformational parameters are discussed. In the statistics of long-range distances, a reduced distanced* _ij (the actual distance divided by the radius of gyration) is used to scale the data so that they become independent of protein size, and then a mean reduced distanceD _l (a, a) and its standard deviation _l (a, a) are calculated for each amino acid pair (a, a) (factors 2 and 4). The effect of the neighboring residues along the chain on the value of the distanced* _ij is explored by a linear regression analysis between the actual reduced distanced* _ij and the mean value over theD _l for all possible pairs of residues in the two segments of the (i–2)th to the (i+2)th and the (j–2)th to the (j+2)th residues. The effect is assessed in terms of the tangentA _l (a, a) of the calculated regression line for each amino acid pair (a, a). In the statistics of medium-range distances, only factors 1 and 4 are considered, to simplify the analysis. The scaled distanced _i,i+k =(d _i,i+k -D _k )/S _k is used to eliminate the dependence onk, the distance along the chain. The propertiesD _m (a, a), _m (a, a) andA _m (a, a) corresponding toD _l (a, a), _l (a, a), andA _l (a, a), and also calculated for each amino acid pair (a, a). The results are interpreted as follows: the smaller values ofD _l (a, a) andD _m (a, a) indicate a preference of the pair (a, a) for a contact (e.g., pairs between hydrophobic amino acids, and pairs of Cys with aromatic amino acids), and the larger values of these quantities indicate a preference for distant mutual location (e.g., pairs between strong hydrophilic amino acids); the smaller values of _l (a, a) and _m (a, a) indicate a strong preference for either contact or noncontact (e.g., pairs between hydrophobic amino acids, and pairs between strong hydrophobic and hydrophilic amino acids, respectively), and the larger values of these quantities indicate the ambivalent/neutral nature of the preference for contact and noncontact (e.g., pairs containing Ser or Thr); the smaller values ofA _l (a, a) andA _m (a, a) indicate that the distance of an (a, a) pair is determined independently of the amino acid character of the neighboring residues along the chain (e.g., some pairs of Cys or Met with other amino acids) and the larger values of these quantities indicare that such amino acid character contributes strongly to the determination of the distance (e.g., pairs containing Ser or Thr, and pairs between amino acids with small side chains). The difference between the statistics for the long- and medium-range distances is also discussed; the former reflect the difference between the hydrophobic and hydrophilic character of the residues, but the latter cannot be easily interpretable only in terms of hydrophobicity and hydrophilicity. The data analyzed here are used in the optimization of an object function to compute protein conformation in a subsequent paper.