Construction of a full three-dimensional model of the transpeptidase domain of Streptococcus pneumoniae PBP2x starting from its Calpha-atom coordinates

SABBAC: online Structural Alphabet-based protein BackBone reconstruction from Alpha-Carbon trace

SABBAC is an on-line service devoted to protein backbone reconstruction from alpha-carbon trace. It is based on the assembly of fragments taken from a library of reduced size, selected from the encoding of the protein trace in a hidden Markov model-derived structural alphabet. The assembly of the fragments is achieved by a greedy algorithm, using an energy-based scoring. Alpha-carbon coordinates remain unaffected. SABBAC simply positions the missing backbone atoms, no further refinement is performed. From our tests, SABBAC performs equal or better than other similar on-line approach and is robust to deviations on the alpha-carbon coordinates. It can be accessed at http://bioserv.rpbs.jussieu.fr/SABBAC.html.

Peptide backbone reconstruction using dead-end elimination and a knowledge-based forcefield

A novel, yet simple and automated, protocol for reconstruction of complete peptide backbones from C(alpha) coordinates only is described, validated, and benchmarked. The described method collates a set of possible backbone conformations for each set of residue triads from a structural library derived from the PDB. The optimal permutation of these three residue segments of backbone conformations is determined using the dead-end elimination (DEE) algorithm. Putative conformations are evaluated using a pairwise-additive knowledge-based forcefield term and a fragment overlap term. The protocol described in this report is able to restore the full backbone coordinates to within 0.2-0.6 A of the actual crystal structure from C(alpha) coordinates only. In addition, it is insensitive to errors in the input C(alpha) coordinates with RMSDs of 3.0 A, and this is illustrated through application to deliberately distorted C(alpha) traces. The entire process, as described, is rapid, requiring of the order of a few minutes for a typical protein on a typical desktop PC. Approximations enable this to be reduced to a few seconds, although this is at the expense of prediction accuracy. This compares very favorably to previously published methods, being sufficiently fast for general use and being one of the most accurate methods. Because the method is not restricted to the reconstruction from only C(alpha) coordinates, reconstruction based on C(beta) coordinates is also demonstrated.

Energy-based reconstruction of a protein backbone from its alpha-carbon trace by a Monte-Carlo method

An automatic procedure is proposed for reconstruction of a protein backbone from its C(alpha)-trace; it is based on optimization of a simplified energy function of a peptide backbone, given its alpha-carbon trace. The energy is expressed as a sum of the energies of interaction between backbone peptide groups that are not neighbors in the sequence, the energies of local interactions within all amino acid residues, and a harmonic penalty function accounting for the conservation of standard bond angles. The energy of peptide group interactions is calculated using the assumption that each peptide group acts as a point dipole. For local interaction energy, use is made of a two-dimensional Fourier series expansion of the energies of model terminally blocked amino acid residues, calculated with the Empirical Conformational Energy Program for Peptides (ECEPP/3) force field in the angles lambda((1)) and lambda((2)) defining the rotation of peptide groups adjacent to a C(alpha) carbon atom about the corresponding C(alpha) em leader C(alpha) virtual-bond axes. To explore all possible rotations of peptide groups within a fixed C(alpha)-trace, a Monte Carlo search is carried out. The initial lambda angles are calculated by aligning the dipoles of the peptide groups that are close in space, subject to the condition of favorable local interactions. After the Monte Carlo search is accomplished with the simplified energy function, the energy of the structure is minimized with the ECEPP/3 force field, with imposition of distance constraints corresponding to the initial C(alpha)-trace geometry. The procedure was tested on model alpha-helices and beta-sheets, as well as on the crystal structure of the immunoglobulin binding protein (PDB code: 1IGD, an alpha/beta protein). In all cases, complete backbone geometry was reconstructed with a root-mean-square (rms) deviation of 0.5 A from the all-atom target structure.