Dataset Information

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.

SUBMITTER: Fogolari F

PROVIDER: S-EPMC1225047 | biostudies-other | 1996 Mar

REPOSITORIES: biostudies-other

ACCESS DATA

Publications

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

Fogolari F F Esposito G G Viglino P P Cattarinussi S S

Biophysical journal 19960301 3

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 sh ...[more]

PMID: 8785277

Similar Datasets

Project description:The implementation of molecular dynamics (MD) with our physics-based protein united-residue (UNRES) force field, described in the accompanying paper, was extended to Langevin dynamics. The equations of motion are integrated by using a simplified stochastic velocity Verlet algorithm. To compare the results to those with all-atom simulations with implicit solvent in which no explicit stochastic and friction forces are present, we alternatively introduced the Berendsen thermostat. Test simulations on the Ala(10) polypeptide demonstrated that the average kinetic energy is stable with about a 5 fs time step. To determine the correspondence between the UNRES time step and the time step of all-atom molecular dynamics, all-atom simulations with the AMBER 99 force field and explicit solvent and also with implicit solvent taken into account within the framework of the generalized Born/surface area (GBSA) model were carried out on the unblocked Ala(10) polypeptide. We found that the UNRES time scale is 4 times longer than that of all-atom MD simulations because the degrees of freedom corresponding to the fastest motions in UNRES are averaged out. When the reduction of the computational cost for evaluation of the UNRES energy function is also taken into account, UNRES (with hydration included implicitly in the side chain-side chain interaction potential) offers about at least a 4000-fold speed up of computations relative to all-atom simulations with explicit solvent and at least a 65-fold speed up relative to all-atom simulations with implicit solvent. To carry out an initial full-blown test of the UNRES/MD approach, we ran Berendsen-bath and Langevin dynamics simulations of the 46-residue B-domain of staphylococcal protein A. We were able to determine the folding temperature at which all trajectories converged to nativelike structures with both approaches. For comparison, we carried out ab initio folding simulations of this protein at the AMBER 99/GBSA level. The average CPU time for folding protein A by UNRES molecular dynamics was 30 min with a single Alpha processor, compared to about 152 h for all-atom simulations with implicit solvent. It can be concluded that the UNRES/MD approach will enable us to carry out microsecond and, possibly, millisecond simulations of protein folding and, consequently, of the folding process of proteins in real time.

Project description:Coarse-grained force fields for protein simulations are usually designed and parameterized to treat proteins composed of natural L-amino-acid residues. However, D-amino-acid residues occur in bacterial, fungal (e.g., gramicidins), as well as human-designed proteins. For this reason, we have extended the UNRES coarse-grained force field developed in our laboratory to treat systems with D-amino-acid residues. We developed the respective virtual-bond-torsional and double-torsional potentials for rotation about the C (?) · · · C (?) virtual-bond axis and two consecutive C (?) · · · C (?) virtual-bond axes, respectively, as functions of virtual-bond-dihedral angles ?. In turn, these were calculated as potentials of mean force (PMFs) from the diabatic energy surfaces of terminally-blocked model compounds for glycine, alanine, and proline. The potential-energy surfaces were calculated by using the ab initio method of molecular quantum mechanics at the Møller-Plesset (MP2) level of theory and the 6-31G(d,p) basis set, with the rotation angles of the peptide groups about [Formula: see text] and [Formula: see text] used as variables, and the energy was minimized with respect to the remaining degrees of freedom. The PMFs were calculated by numerical integration for all pairs and triplets with all possible combinations of types (glycine, alanine, and proline) and chirality (D or L); however, symmetry relations reduce the number of non-equivalent torsional potentials to 13 and the number of double-torsional potentials to 63 for a given C-terminal blocking group. Subsequently, one- (for torsional) and two-dimensional (for double-torsional potentials) Fourier series were fitted to the PMFs to obtain analytical expressions. It was found that the torsional potentials of the x-Y and X-y types, where X and Y are Ala or Pro, respectively, and a lowercase letter denotes D-chirality, have global minima for small absolute values of ?, accounting for the double-helical structure of gramicidin A, which is a dimer of two chains, each possessing an alternating D-Tyr-L-Tyr sequence, and similar peptides. The side-chain and correlation potentials for D-amino-acid residues were obtained by applying the reflection about the [Formula: see text] plane to the respective potentials for the L-amino-acid residues.

Dataset Information

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

Publications

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

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets