Search

The structure of the SH3 domain of α-spectrin (PDB code 1SHG) features Asn47 in the II' area of the Ramachandran plot, which as a rule admits only glycine residues, and this phenomenon still awaits its explanation. Here, we undertook a computational study of this particular case by means of molecular dynamics and bioinformatics approaches. We found that the region of the SH3 domain in the vicinity of Asn47 remains relatively stable during denaturing molecular dynamics simulations of the entire domain and of its parts. This increased stability may be connected with the dynamic hydrogen bonding that is susceptible to targeted in silico mutations of Arg49. Bioinformatics analysis indicated that Asn47 is in the β-turn of a distinctive structural fragment we called 'alternatingly twisted β-hairpin.' Fragments of similar conformation are quite abundant in a nonredundant set of PDB chains and are distinguished from ordinary β-hairpins by some surplus of glycine in their β-turns, lack of certain interpeptide hydrogen bonds, and an increased chirality index. Thus, the disallowed conformation of residues other than glycine is realized in the β-turns of alternatingly twisted β-hairpins.

Van der Waals interaction energy in globular proteins is presented by the interaction energies between regions of protein spatial structure with homogenous medium density distribution. We introduce a notion of the local medium permittivity as a function of absorptance of molecular groups with particular conformation. Proposed theory avoids shortcomings which are typical for the calculations on the basis of the pairwise additive approximation. The approach takes into account local peculiarities of protein spatial structure and physical-chemical characteristics of amino acid residues and molecular groups.

A complete classification of types of the protein secondary structure is developed on the basis of computer analysis of the crystallographic structural data deposited in the protein Data Bank. The majority of amino acid residues fall into five conformation types. A conclusion is drawn that the number of sequence variants of torsion angles phi, psi in globular proteins is limited and is essentially less than the number of possible amino acid sequences for this chain length. Along with alpha-helix and beta-structure, the distribution analysis assigning every maximum of distribution of amino acid conformations on Ramachandran map to a certain type of the secondary structure exposed a third type of the secondary structure that was previously neglected. This type of the structure is extended left-handed helical conformation, designated as mobile (M-) conformation. A full set of M-conformation fragments that seems to play a major role in protein globule dynamics has been obtained, a small radius of correlation for the polypeptide chain in M-conformation is demonstrated. It explains a prevalence of short segments of mobile conformation revealed in globular proteins. For secondary structure types, the frequency of occurrence of amino acid residues has been computed.