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

The basic differences between the 20 natural amino acid residues are due to differences in their side‐chain structures. This characteristic design of protein building blocks implies that side‐chain–side‐chain interactions play an important, even dominant role in 3D‐structural realization of amino acid codes. Here we present the results of a comparative analysis of the contributions of side‐chain–side‐chain (s‐s) and side‐chain–backbone (s‐b) interactions to the stabilization of folded protein structures within the framework of the CHARMm molecular data model. Contrary to intuition, our results suggest that side‐chain–backbone interactions play the major role in side‐chain packing, in stabilizing the folded structures, and in differentiating the folded structures from the unfolded or misfolded structures, while the interactions between side chains have a secondary effect. An additional analysis of electrostatic energies suggests that combinatorial dominance of the interactions between opposite charges makes the electrostatic interactions act as an unspecific folding force that stabilizes not only native structure, but also compact random conformations. This observation is in agreement with experimental findings that, in the denatured state, the charge–charge interactions stabilize more compact conformations. Taking advantage of the dominant role of side‐chain–backbone interactions in side‐chain packing to reduce the combinatorial problem, we developed a new algorithm, ChiRotor, for rapid prediction of side‐chain conformations. We present the results of a validation study of the method based on a set of high resolution X‐ray structures.

[1]  M. Karplus,et al.  Simulation of activation free energies in molecular systems , 1996 .

[2]  A Joshua Wand,et al.  Improved side‐chain prediction accuracy using an ab initio potential energy function and a very large rotamer library , 2004, Protein science : a publication of the Protein Society.

[3]  H. Scheraga,et al.  Medium- and long-range interaction parameters between amino acids for predicting three-dimensional structures of proteins. , 1976, Macromolecules.

[4]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[5]  B. Lee,et al.  The interpretation of protein structures: estimation of static accessibility. , 1971, Journal of molecular biology.

[6]  A. Karshikoff,et al.  Optimization of the electrostatic interactions between ionized groups and peptide dipoles in proteins , 1997, Protein science : a publication of the Protein Society.

[7]  Martin Karplus,et al.  pH-Dependence of Protein Stability: Absolute Electrostatic Free Energy Differences between Conformations† , 1997 .

[8]  M Karplus,et al.  Side-chain torsional potentials: effect of dipeptide, protein, and solvent environment. , 1979, Biochemistry.

[9]  I Lasters,et al.  All in one: a highly detailed rotamer library improves both accuracy and speed in the modelling of sidechains by dead-end elimination. , 1997, Folding & design.

[10]  B. Atanasov,et al.  Spatial optimization of electrostatic interactions between the ionized groups in globular proteins , 1994, Proteins.

[11]  A. Karshikoff,et al.  Effects of charge–charge interactions on dimensions of unfolded proteins: A Monte Carlo study , 2003 .

[12]  D. Case,et al.  Generalized born models of macromolecular solvation effects. , 2000, Annual review of physical chemistry.

[13]  Adrian A Canutescu,et al.  Access the most recent version at doi: 10.1110/ps.03154503 References , 2003 .

[14]  A. Karshikoff,et al.  Optimization of the electrostatic interactions in proteins of different functional and folding type , 1994, Protein science : a publication of the Protein Society.

[15]  F. Momany,et al.  Validation of the general purpose QUANTA ®3.2/CHARMm® force field , 1992 .

[16]  P Argos,et al.  A method to configure protein side-chains from the main-chain trace in homology modelling. , 1993, Journal of molecular biology.

[17]  Cheng-Yan Kao,et al.  GEM: A Gaussian evolutionary method for predicting protein side‐chain conformations , 2002, Protein science : a publication of the Protein Society.

[18]  Johan Desmet,et al.  The dead-end elimination theorem and its use in protein side-chain positioning , 1992, Nature.

[19]  L L Looger,et al.  Generalized dead-end elimination algorithms make large-scale protein side-chain structure prediction tractable: implications for protein design and structural genomics. , 2001, Journal of molecular biology.

[20]  C. Sander,et al.  Evaluation of protein models by atomic solvation preference. , 1992, Journal of molecular biology.

[21]  J. Ponder,et al.  Tertiary templates for proteins. Use of packing criteria in the enumeration of allowed sequences for different structural classes. , 1987, Journal of molecular biology.

[22]  A. Kidera,et al.  Determinants of protein side‐chain packing , 1994, Protein science : a publication of the Protein Society.

[23]  Kevin L. Shaw,et al.  Charge–charge interactions influence the denatured state ensemble and contribute to protein stability , 2000, Protein science : a publication of the Protein Society.

[24]  A. Godzik,et al.  Derivation and testing of pair potentials for protein folding. When is the quasichemical approximation correct? , 1997, Protein science : a publication of the Protein Society.

[25]  Roland L. Dunbrack,et al.  Backbone-dependent rotamer library for proteins. Application to side-chain prediction. , 1993, Journal of molecular biology.

[26]  R. Jernigan,et al.  Estimation of effective interresidue contact energies from protein crystal structures: quasi-chemical approximation , 1985 .

[27]  J. Straub,et al.  Orientational potentials extracted from protein structures improve native fold recognition , 2004, Protein science : a publication of the Protein Society.

[28]  B. Honig,et al.  Free energy determinants of tertiary structure and the evaluation of protein models , 2000, Protein science : a publication of the Protein Society.

[29]  N. Grishin,et al.  Side‐chain modeling with an optimized scoring function , 2002, Protein science : a publication of the Protein Society.

[30]  J. Warwicker,et al.  Simplified methods for pKa and acid pH‐dependent stability estimation in proteins: Removing dielectric and counterion boundaries , 2008, Protein science : a publication of the Protein Society.

[31]  B. Dominy,et al.  Development of a generalized Born model parameterization for proteins and nucleic acids , 1999 .

[32]  M. Gunner,et al.  Backbone dipoles generate positive potentials in all proteins: origins and implications of the effect. , 2000, Biophysical journal.

[33]  A. Elcock.  Realistic modeling of the denatured states of proteins allows accurate calculations of the pH dependence of protein stability. , 1999, Journal of molecular biology.

[34]  Guoli Wang,et al.  PISCES: a protein sequence culling server , 2003, Bioinform..

[35]  Z. Xiang,et al.  Extending the accuracy limits of prediction for side-chain conformations. , 2001, Journal of molecular biology.

[36]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .