Energy‐based reconstruction of a protein backbone from its α‐carbon trace by a Monte‐Carlo method
- 26 March 2002
- journal article
- research article
- Published by Wiley in Journal of Computational Chemistry
- Vol. 23 (7) , 715-723
- https://doi.org/10.1002/jcc.10068
Abstract
An automatic procedure is proposed for reconstruction of a protein backbone from its Cα‐trace; it is based on optimization of a simplified energy function of a peptide backbone, given its α‐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 λ(1) and λ(2) defining the rotation of peptide groups adjacent to a Cα carbon atom about the corresponding Cα … Cα virtual‐bond axes. To explore all possible rotations of peptide groups within a fixed Cα‐trace, a Monte Carlo search is carried out. The initial λ 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α‐trace geometry. The procedure was tested on model α‐helices and β‐sheets, as well as on the crystal structure of the immunoglobulin binding protein (PDB code: 1IGD, an α/β protein). In all cases, complete backbone geometry was reconstructed with a root‐mean‐square (rms) deviation of 0.5 Å from the all‐atom target structure. © 2002 Wiley Periodicals, Inc. J Comput Chem 7: 715–723, 2002Keywords
This publication has 31 references indexed in Scilit:
- Development of Physics-Based Energy Functions that Predict Medium-Resolution Structures for Proteins of the α, β, and α/β Structural ClassesThe Journal of Physical Chemistry B, 2001
- Hierarchical energy-based approach to protein-structure prediction: Blind-test evaluation with CASP3 targetsInternational Journal of Quantum Chemistry, 2000
- New developments of the electrostatically driven monte carlo method: Test on the membrane-bound portion of melittinBiopolymers, 1998
- United-residue force field for off-lattice protein-structure simulations: III. Origin of backbone hydrogen-bonding cooperativity in united-residue potentialsJournal of Computational Chemistry, 1998
- A united-residue force field for off-lattice protein-structure simulations. II. Parameterization of short-range interactions and determination of weights of energy terms by Z-score optimizationJournal of Computational Chemistry, 1997
- The Third IgG-Binding Domain from Streptococcal Protein GJournal of Molecular Biology, 1994
- Prediction of protein conformation on the basis of a search for compact structures: Test on avian pancreatic polypeptideProtein Science, 1993
- Database algorithm for generating protein backbone and side-chain co-ordinates from a Cα trace: Application to model building and detection of co-ordinate errorsJournal of Molecular Biology, 1991
- The multiple-minima problem in the conformational analysis of polypeptides. III. An Electrostatically Driven Monte Carlo Method: Tests on enkephalinProtein Journal, 1989
- The protein data bank: A computer-based archival file for macromolecular structuresJournal of Molecular Biology, 1977