Cavities and Atomic Packing in Protein Structures and Interfaces

A comparative analysis of cavities enclosed in a tertiary structure of proteins and interfaces formed by the interaction of two protein subunits in obligate and non-obligate categories (represented by homodimeric molecules and heterocomplexes, respectively) is presented. The total volume of cavities increases with the size of the protein (or the interface), though the exact relationship may vary in different cases. Likewise, for individual cavities also there is quantitative dependence of the volume on the number of atoms (or residues) lining the cavity. The larger cavities tend to be less spherical, solvated, and the interfaces are enriched in these. On average 15 Å3 of cavity volume is found to accommodate single water, with another 40–45 Å3 needed for each additional solvent molecule. Polar atoms/residues have a higher propensity to line solvated cavities. Relative to the frequency of occurrence in the whole structure (or interface), residues in β-strands are found more often lining the cavities, and those in turn and loop the least. Any depression in one chain not complemented by a protrusion in the other results in a cavity in the protein–protein interface. Through the use of the Voronoi volume, the packing of residues involved in protein–protein interaction has been compared to that in the protein interior. For a comparable number of atoms the interface has about twice the number of cavities relative to the tertiary structure

Secondary structure based analysis and classification of biological interfaces: identification of binding motifs in protein-protein interactions

Motivation: The increasing amount of data on protein-protein interaction needs to be rationalized for deriving guidelines for the alteration or design of an interface between two proteins. Results: We present a detailed structural analysis and comparison of homo- versus heterodimeric protein-protein interfaces. Regular secondary structures (helices and strands) are the main components of the former, whereas non-regular structures (turns, loops, etc.) frequently mediate interactions in the latter. Interface helices get longer with increasing interface area, but only in heterocomplexes. On average, the homodimers have longer helical segments and prominent helix-helix pairs. There is a surprising distinction in the relative orientation of interface helices, with a tendency for aligned packing in homodimers and a clear preference for packing at 90° in heterodimers. Arg and the aromatic residues have a higher preference to occur in all secondary structural elements (SSEs) in the interface. Based on the dominant SSE, the interfaces have been grouped into four classes: α, β, αβ and non-regular. Identity between protein and interface classes is the maximum for α proteins, but rather mediocre for the other protein classes. The interface classes of the two chains forming a heterodimer are often dissimilar. Eleven binding motifs can capture the prominent architectural features of most of the interfaces

Estimates of the loss of main-chain conformational entropy of different residues on protein folding

The average contribution of conformational entropy for individual amino acid residues towards the free energy of protein folding is not well understood. We have developed empirical scales for the loss of the main-chain (torsion angles, φ and ψ) conformational entropy by taking its side-chain into account. The analysis shows that the main-chain component of the total conformational entropy loss for a residue is significant and reflects intrinsic characteristics associated with individual residues. The values have direct correlation with the hydrophobicity values and this has important bearing on the folding process

Cavities in protein–DNA and protein–RNA interfaces

An analysis of cavities present in protein–DNA and protein–RNA complexes is presented. In terms of the number of cavities and their total volume, the interfaces formed in these complexes are akin to those in transient protein–protein heterocomplexes. With homodimeric proteins protein–DNA interfaces may contain cavities involving both the protein subunits and DNA, and these are more than twice as large as cavities involving a single protein subunit and DNA. A parameter, cavity index, measuring the degree of surface complementarity, indicates that the packing of atoms in protein–protein/DNA/RNA is very similar, but it is about two times less efficient in the permanent interfaces formed between subunits in homodimers. As within the tertiary structure and protein–protein interfaces, protein–DNA interfaces have a higher inclination to be lined by β-sheet residues; from the DNA side, base atoms, in particular those in minor grooves, have a higher tendency to be located in cavities. The larger cavities tend to be less spherical and solvated. A small fraction of water molecules are found to mediate hydrogen-bond interactions with both the components, suggesting their primary role is to fill in the void left due to the local non-complementary nature of the surface patches

pi-Turns: types, systematics and the context of their occurrence in protein structures

<p>Abstract</p> <p>Background</p> <p>For a proper understanding of protein structure and folding it is important to know if a polypeptide segment adopts a conformation inherent in the sequence or it depends on the context of its flanking secondary structures. Turns of various lengths have been studied and characterized starting from three-residue γ-turn to six-residue π-turn. The Schellman motif occurring at the C-terminal end of α-helices is a classical example of hydrogen bonded π-turn involving residues at (i) and (i+5) positions. Hydrogen bonded and non-hydrogen bonded β- and α-turns have been identified previously; likewise, a systematic characterization of π-turns would provide valuable insight into turn structures.</p> <p>Results</p> <p>An analysis of protein structures indicates that at least 20% of π-turns occur independent of the Schellman motif. The two categories of π-turns, designated as π-HB and SCH, have been further classified on the basis of backbone conformation and both have AAAa as the major class. They differ in the residue usage at position (i+1), the former having a large preference for Pro that is absent in the latter. As in the case of shorter length β- and α-turns, π-turns have also been identified not only on the basis of the existence of hydrogen bond, but also using the distance between terminal C^α-atoms, and this resulted in a comparable number of non-hydrogen-bonded π-turns (π-NHB). The presence of shorter β- and α-turns within all categories of π-turns, the subtle variations in backbone torsion angles along the turn residues, the location of the turns in the context of tertiary structures have been studied.</p> <p>Conclusion</p> <p>π-turns have been characterized, first using hydrogen bond and the distance between C^αatoms of the terminal residues, and then using backbone torsion angles. While the Schellman motif has a structural role in helix termination, many of the π-HB turns, being located on surface cavities, have functional role and there is also sequence conservation.</p

Discriminating the native structure from decoys using scoring functions based on the residue packing in globular proteins

<p>Abstract</p> <p>Background</p> <p>Setting the rules for the identification of a stable conformation of a protein is of utmost importance for the efficient generation of structures in computer simulation. For structure prediction, a considerable number of possible models are generated from which the best model has to be selected.</p> <p>Results</p> <p>Two scoring functions, R_sand R_p, based on the consideration of packing of residues, which indicate if the conformation of an amino acid sequence is native-like, are presented. These are defined using the solvent accessible surface area (ASA) and the partner number (PN) (other residues that are within 4.5 Å) of a particular residue. The two functions evaluate the deviation from the average packing properties (ASA or PN) of all residues in a polypeptide chain corresponding to a model of its three-dimensional structure. While simple in concept and computationally less intensive, both the functions are at least as efficient as any other energy functions in discriminating the native structure from decoys in a large number of standard decoy sets, as well as on models submitted for the targets of CASP7. R_sappears to be slightly more effective than R_p, as determined by the number of times the native structure possesses the minimum value for the function and its separation from the average value for the decoys.</p> <p>Conclusion</p> <p>Two parameters, R_sand R_p, are discussed that can very efficiently recognize the native fold for a sequence from an ensemble of decoy structures. Unlike many other algorithms that rely on the use of composite scoring function, these are based on a single parameter, viz., the accessible surface area (or the number of residues in contact), but still able to capture the essential attribute of the native fold.</p

Molecular modeling and characterization of Vibrio cholerae transcription regulator HlyU

BACKGROUND: The SmtB/ArsR family of prokaryotic metal-regulatory transcriptional repressors represses the expression of operons linked to stress-inducing concentrations of heavy metal ions, while derepression results from direct binding of metal ions by these 'metal-sensor' proteins. The HlyU protein from Vibrio cholerae is the positive regulator of haemolysin gene, it also plays important role in the regulation of expression of the virulence genes. Despite the understanding of biochemical properties, its structure and relationship to other protein families remain unknown. RESULTS: We find that HlyU exhibits structural features common to the SmtB/ArsR family of transcriptional repressors. Analysis of the modeled structure of HlyU reveals that it does not have the key metal-sensing residues which are unique to the SmtB/ArsR family of repressors, yet the tertiary structure is very similar to the family members. HlyU is the only member that has a positive control on transcription, while all the other members in the family are repressors. An evolutionary analysis with other SmtB/ArsR family members suggests that during evolution HlyU probably occurred by gene duplication and mutational events that led to the emergence of this protein from ancestral transcriptional repressor by the loss of the metal-binding sites. CONCLUSION: The study indicates that the same protein family can contain both the positive regulator of transcription and repressors – the exact function being controlled by the absence or the presence of metal-binding sites

Disulfide bonds, their stereospecific environment and conservation in protein structures

We studied the specificity of the non-bonded interaction in the environment of 572 disulfide bonds in 247 polypeptide chains selected from the Protein Data Bank. The preferred geometry of interaction of peptide oxygen atoms is along the back of the two covalent bonds at the sulfur atom of half cystine. With aromatic residues the geometries that direct one of the sulfur lone pair of electrons into the aromatic Π-system are avoided; an orientation in which the sulfide plane is normal or inclined to the aromatic plane and on top of its edge is normally preferred. The importance of the S···aromatic interaction is manifested in the high degree of its conservation across members in homologous protein families. These interactions, while providing extra overall stability to the native fold and reducing the accessibility of the disulfide bond and thereby preventing exchange reactions, also set the orientation of the conserved aromatic rings for further interactions and binding to another molecule. The conformational features and the mode of interactions of disulfide bridges should be useful for molecular design and protein engineering experiments

Structural segments and residue propensities in protein-RNA interfaces: Comparison with protein-protein and protein-DNA complexes

The interface of a protein molecule that is involved in binding another protein, DNA or RNA has been characterized in terms of the number of unique secondary structural segments (SSSs), made up of stretches of helix, strand and non-regular (NR) regions. On average 10-11 segments define the protein interface in protein-protein (PP) and protein-DNA (PD) complexes, while the number is higher (14) for protein-RNA (PR) complexes. While the length of helical segments in PP interaction increases with the interface area, this is not the case in PD and PR complexes. The propensities of residues to occur in the three types of secondary structural elements (SSEs) in the interface relative to the corresponding elements in the protein tertiary structures have been calculated. Arg, Lys, Asn, Tyr, His and Gln are preferred residues in PR complexes; in addition, Ser and Thr are also favoured in PD interfaces

Packing of aromatic rings against tryptophan residues in proteins

The geometry of the interaction of the aromatic side chains of phenylalanine (Phe), tyrosine (Tyr), tryptophan (Trp) and histidine (His) with the indole ring of Trp has been analyzed using the structures in the Protein Data Bank in order to understand the dependence of the packing behaviour on the size and chemical nature of the aromatic rings. The Phe ring prefers to interact either perpendicularly, with its edge pointing towards the Trp face, or in an offset-stacked arrangement. The edge-to-face motif is typical of a Trp-Trp pair. While parallel stacking is the dominant feature of Trp-His interaction, Tyr packs in a more uniform manner around Trp with a higher than expected occurrence at the edge and a few cases of possible OH-Π interaction