Display options
Share it on

J Chem Theory Comput. 2010 Jul 13;6(7):2191-203. doi: 10.1021/ct100007y.

Energy Matrix of Structurally Important Side-Chain/Side-Chain Interactions in Proteins.

Journal of chemical theory and computation

Karel Berka, Roman A Laskowski, Pavel Hobza, Jiří Vondrášek

Affiliations

  1. Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic and Center for Biomolecules and Complex Molecular Systems, Flemingovo nám. 2, Prague, Czech Republic, Institute of Biotechnology, Academy of Sciences of the Czech Republic, Videnska 1083, 142 00 Prague, Czech Republic, Palacký University, Department of Physical Chemistry, Faculty of Science, tø. 17. listopadu 12, 771 46, Olomouc, Czech Republic, and EMBL Outstation - Hinxton, European Bioinformatics Institute, Welcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SD, United Kingdom.

PMID: 26615945 DOI: 10.1021/ct100007y

Abstract

The interactions between amino acid side chains in proteins are generally considered to be the most important stabilizing factor controlling the precise arrangement of the polypeptide chain into a well-defined spatial structure. We used the RI-DFT-D method to calculate the full 20 × 20 matrix of interaction energies between all pairs of amino acid side chains. For each pair, we used a representative 3D conformation extracted from an analysis of known protein structures from Protein Data Bank (PDB). The representative comes from the largest cluster of relative orientations of the two side chains. We find that all of the calculated interaction energies between selected pairs of amino acids are attractive in the gas phase with the exception of side chain pairs having the same total charge. We compared these data with those calculated by the parm03 and OPLS-AA/L force fields to investigate the reliability of simple methods in modeling biomolecules and their behavior. The force fields yield good overall interaction energies for our set but have problems in evaluation of some particular interactions which could be of principal importance for protein stability. We then looked in detail at the 20 side chain interactions involving tryptophan. The histograms of interaction energies showed that the distributions of the interaction energies are neither normal nor Boltzmann-like and that our representative geometries correspond mostly to the minimum energy geometry which is rather poorly populated in the whole pairwise energy distribution. We concluded that cluster representatives obtained by the clusterization algorithm based on geometry criteria cannot be considered as a typical interaction for the whole side chain/side chain interaction distribution. They seem to epitomize the strongest interactions in a protein and are often functionally or structurally important.

Publication Types