Modeling of polypeptide chains as C alpha chains, C alpha chains with C beta, and C alpha chains with ellipsoidal lateral chains

Abstract

In an effort to reduce the number of degrees of freedom necessary to describe a polypeptide chain we analyze the statistical behavior of polypeptide chains when represented as C alpha chains, C alpha chains with C beta atoms attached, and C alpha chains with rotational ellipsoids as models of side chains. A statistical analysis on a restricted data set of 75 unrelated protein structures is performed. The comparison of the database distributions with those obtained by model calculation on very short polypeptide stretches allows the dissection of local versus nonlocal features of the distributions. The database distribution of the bend angles of polypeptide chains of pseudo bonded C alpha atoms spans a restricted range of values and shows a bimodal structure. On the other hand, the torsion angles of the C alpha chain may assume almost all possible values. The distribution is bimodal, but with a much broader probability distribution than for bend angles. The C alpha - C beta vectors may be taken as representative of the orientation of the lateral chain, as the direction of the bond is close to the direction of the vector joining C alpha to the ad hoc defined center of the "steric mass" of the side chain. Interestingly, both the bend angle defined by C alpha i-C alpha i+1-C beta i+1 and the torsional angle offset of the pseudo-dihedral C alpha i-C alpha i+1-C alpha i+2-C beta i+2 with respect to C alpha i-C alpha i+1-C alpha i+2-C alpha i+3 span a limited range of values. The latter results show that it is possible to give a more realistic representation of polypeptide chains without introducing additional degrees of freedom, i.e., by just adding to the C alpha chain a C beta with given side-chain properties. However, a more realistic description of side chains may be attained by modeling side chains as rotational ellipsoids that have roughly the same orientation and steric hindrance. To this end, we define the steric mass of an atom as proportional to its van der Waals volume and we calculate the side-chain inertia ellipsoid assuming that the steric mass of each atom is uniformly distributed within its van der Waals volume. Finally, we define the rotational ellipsoid representing the side chain as the uniform density ellipsoid possessing the same rotationally averaged inertia tensor of the side chain. The statistics of ellipsoid parameters support the possibility of representing a side chain via an ellipsoid, independently of the local conformation. To make this description useful for molecular modeling we describe ellipsoid-ellipsoid interactions via a Lennard-Jones potential that preserves the repulsive core of the interacting ellipsoids and takes into account their mutual orientation. Tests are performed for two different forms of the interaction potential on a set of high-resolution protein structures. Results are encouraging, in view of the drastic simplifications that were introduced.