Protein Science, Vol 6, Issue 3 676-688, Copyright © 1997 by Cold Spring Harbor Laboratory Press

ARTICLE

Derivation and testing of pair potentials for protein folding. When is the quasichemical approximation correct?

J. SKOLNICK, L. JAROSZEWSKI, A. KOLINSKI and A. GODZIK
Department of Molecular Biology, The Scripps Research Institute, La Jolla, California 92037

Many existing derivations of knowledge-based statistical pair potentials invoke the quasichemical approximation to estimate the expected side-chain contact frequency if there were no amino acid pair-specific interactions. At first glance, the quasichemical approximation that treats the residues in a protein as being disconnected and expresses the side-chain contact probability as being proportional to the product of the mole fractions of the pair of residues would appear to be rather severe. To investigate the validity of this approximation, we introduce two new reference states in which no specific pair interactions between amino acids are allowed, but in which the connectivity of the protein chain is retained. The first estimates the expected number of side-chain contacts by treating the protein as a Gaussian random coil polymer. The second, more realistic reference state includes the effects of chain connectivity, secondary structure, and chain compactness by estimating the expected side-chain contact probability by placing the sequence of interest in each member of a library of structures of comparable compactness to the native conformation. The side-chain contact maps are not allowed to readjust to the sequence of interest, i.e., the side chains cannot repack. This situation would hold rigorously if all amino acids were the same size. Both reference states effectively permit the factorization of the side-chain contact probability into sequence-dependent and structure-dependent terms. Then, because the sequence distribution of amino acids in proteins is random, the quasichemical approximation to each of these reference states is shown to be excellent. Thus, the range of validity of the quasichemical approximation is determined by the magnitude of the side-chain repacking term, which is, at present, unknown. Finally, the performance of these two sets of pair interaction potentials as well as side-chain contact fraction-based interaction scales is assessed by inverse folding tests both without and with allowing for gaps.
CiteULike    Connotea    Del.icio.us    Digg    Reddit    Technorati    What's this?

This article has been cited by other articles:

				J. E. Fitzgerald, A. K. Jha, A. Colubri, T. R. Sosnick, and K. F. Freed Reduced Cbeta statistical potentials can outperform all-atom potentials in decoy identification Protein Sci., October 1, 2007; 16(10): 2123 - 2139. [Abstract] [Full Text] [PDF]

				J. Cheng, J. Pei, and L. Lai A Free-Rotating and Self-Avoiding Chain Model for Deriving Statistical Potentials Based on Protein Structures Biophys. J., June 1, 2007; 92(11): 3868 - 3877. [Abstract] [Full Text] [PDF]

				V. Z. Spassov, L. Yan, and P. K. Flook The dominant role of side-chain backbone interactions in structural realization of amino acid code. ChiRotor: A side-chain prediction algorithm based on side-chain backbone interactions Protein Sci., March 1, 2007; 16(3): 494 - 506. [Abstract] [Full Text] [PDF]

				M.-y. Shen and A. Sali Statistical potential for assessment and prediction of protein structures. Protein Sci., November 1, 2006; 15(11): 2507 - 2524. [Abstract] [Full Text] [PDF]

				W. W. Chen and E. I. Shakhnovich Lessons from the design of a novel atomic potential for protein folding Protein Sci., July 1, 2005; 14(7): 1741 - 1752. [Abstract] [Full Text] [PDF]

				Y. Zhang, A. Kolinski, and J. Skolnick TOUCHSTONE II: A New Approach to Ab Initio Protein Structure Prediction Biophys. J., August 1, 2003; 85(2): 1145 - 1164. [Abstract] [Full Text] [PDF]

				B. Wallner and A. Elofsson Can correct protein models be identified? Protein Sci., May 1, 2003; 12(5): 1073 - 1086. [Abstract] [Full Text] [PDF]

				C. Zhang, J. Hou, and S.-H. Kim Fold prediction of helical proteins using torsion angle dynamics and predicted restraints PNAS, March 19, 2002; 99(6): 3581 - 3585. [Abstract] [Full Text] [PDF]

				D. K. Klimov and D. Thirumalai Native topology determines force-induced unfolding pathways in globular proteins PNAS, June 20, 2000; 97(13): 7254 - 7259. [Abstract] [Full Text] [PDF]

				U. Bastolla, M. Vendruscolo, and E.-W. Knapp A statistical mechanical method to optimize energy functions for protein folding PNAS, April 11, 2000; 97(8): 3977 - 3981. [Abstract] [Full Text] [PDF]

				F. Baud and S. Karlin Measures of residue density in protein structures PNAS, October 26, 1999; 96(22): 12494 - 12499. [Abstract] [Full Text] [PDF]

				C. Zhang and S.-H. Kim Environment-dependent residue contact energies for proteins PNAS, March 14, 2000; 97(6): 2550 - 2555. [Abstract] [Full Text] [PDF]