Methods of ab initio prediction of α helices, β sheets, and polypeptide tertiary structures

• US 6,832,162 B2
• Filed: 02/16/2001
• Issued: 12/14/2004
• Est. Priority Date: 02/16/2001
• Status: Expired due to Fees

• Current Assignee: The Trustees of Princeton University (Princeton University)
• Original Assignee: The Trustees of Princeton University (Princeton University)
• Inventors: Floudas, Christodoulos A., Klepeis, John L.
• Primary Examiner(s): BORIN, MICHAEL L
• Application Number:
  US09/788,006
• Publication Number:
  US 20030036093A1
• Time in Patent Office: 1,397 Days
• Field of Search: 702/19, 702/27
• US Class Current: 702/19
• CPC Class Codes: G16B 15/00 ICT specially adapted for analysing two-dimensional or three-dimensional molecular structures, e.g. structural or functional relations or structure alignment
• CPC Class Codes: G16B 15/20 Protein or domain folding
• CPC Class Codes: G16B 30/00 ICT specially adapted for sequence analysis involving nucleotides or amino acids
• CPC Class Codes: G16B 30/20 Sequence assembly

Abstract

The present invention provides a novel four stage ab initio approach for predicting the tertiary structure of polypeptides. The methods of the invention combine the classical and modern views of protein folding, while using free energy calculations and integer linear optimization to predict helical and β-sheet structures. Derivation of restraints, detailed atomistic modeling, and a deterministic global optimization method, αBB, coupled with torsion angle dynamics, form the basis for the final tertiary structure prediction. The performance of the methods of the invention is illustrated using several different polypeptides.

Claims

1. A method of determining the existence and location of alpha helix regions of a polypeptide having an amino acid sequence by classifying individual residues in the overall sequence as helical or non-helical comprising:

(a) defining a first segment of said amino acid sequence having a length, between 5 and about 15 residues;

(b) partitioning the amino acid sequence into overlapping segments such that consecutive segments possess common core of at least four amino acid residues;

(c) for overlapping segments comprising the common core, before “

performing atomistic modeling ”

performing atomistic modeling upon each segment to assay the free energy of each segment;

(d) generating an ensemble of low energy conformations for each segment;

(e) adjusting the ensemble for each segment to reflect a determination of entropic and free energy contributions for each segment;

(f) modifying the ensemble for each segment to reflect the contributions to free energy accorded by at least one of cavity formation, solvation, and ionization;

(g) ascertaining the equilibrium probabilities for helical clusters in each of the overlapping segments.

(h) determining probability for each residue of being involved in an alpha helix region from the average of equilibrium probabilities of those overlapping segments in which the residue constitutes a core position.

2. The method of claim 1 wherein the defining of further segments is performed iteratively and subsequent steps (c) through (h) are performed on all segments thus defined.

3. The method of claim 1, wherein each segment has a length of 5 amino acid residues.

4. The method of claim 1, wherein each segment has a length of 7 amino acid residues.

5. The method of claim 1, wherein each segment has a length of 9 amino acid residues.

6. The method of claim 2, wherein each segment has a length of 5 amino acid residues.

7. The method of claim 2, wherein each segment has a length of 7 amino acid residues.

8. The method of claim 2, wherein each segment has a length of 9 amino acid residues.

9. The method of any of claims 1-8, wherein the atomistic modeling is performed with a force field.

10. The method of any of claims 1-8, wherein the force field is an Empirical Conformational Energy Program for Peptides/3 (ECEPP/3) force field.