Influence of Medium and Long Range Interactions in Different Structural Classes of Globular Proteins

An analysis of the dependence known three dimensional structure ofglobular proteins on their residue contacts and their interactions providesmuch information about their folding and stability. In this work, we analysethe residue-residue contacts and the role of medium and long rangeinteractions in globular proteins belonging to different structural classes.The results show that while medium range interactions predominate in allalpha class proteins, long range interactions predominate in all beta class.The residues Pro and Gly are found to have lowest medium range contacts,probably due to their helix breaking tendency. The hydrophobic residues Ile,Val and Tyr have higher long range contacts, and hence may serve as goodnucleation centres. Further, the role of charged residues and disulfidebridges in these interactions are also discussed.     

Selected peptide extension contacts hydrophobic patch on neighboring zinc finger and mediates dimerization on DNA.

Protein-protein interactions often play a crucial role in stabilizing protein-DNA complexes and thus facilitate site-specific DNA recognition. We have worked to incorporate such protein-protein contacts into our design and selection strategies for short peptide extensions that promote cooperative binding of zinc finger proteins to DNA. We have determined the crystal structure of one of these fusion protein-DNA complexes. The selected peptide extension was found to mediate dimerization by reaching across the dyad axis and contacting a hydrophobic patch on the surface of the zinc finger bound to the adjacent DNA site. The peptide-zinc finger protein interactions observed in this structure are similar to those of some homeodomain heterodimers. We also find that the region of the zinc finger surface contacted by the selected peptide extension corresponds to surfaces that also make key interactions in the zinc finger proteins GLI and SWI5.     

Protein contacts, inter-residue interactions and side-chain modelling

Three-dimensional structures of proteins are the support of their biological functions. Their folds are stabilized by contacts between residues. Inner protein contacts are generally described through direct atomic contacts, i.e. interactions between side-chain atoms, while contact prediction methods mainly used inter-Calpha distances. In this paper, we have analyzed the protein contacts on a recent high quality non-redundant databank using different criteria. First, we have studied the average number of contacts depending on the distance threshold to define a contact. Preferential contacts between types of amino acids have been highlighted. Detailed analyses have been done concerning the proximity of contacts in the sequence, the size of the proteins and fold classes. The strongest differences have been extracted, highlighting important residues. Then, we studied the influence of five different side-chain conformation prediction methods (SCWRL, IRECS, SCAP, SCATD and SCCOMP) on the distribution of contacts. The prediction rates of these different methods are quite similar. However, using a distance criterion between side chains, the results are quite different, e.g. SCAP predicts 50% more contacts than observed, unlike other methods that predict fewer contacts than observed. Contacts deduced are quite distinct from one method to another with at most 75% contacts in common. Moreover, distributions of amino acid preferential contacts present unexpected behaviours distinct from previously observed in the X-ray structures, especially at the surface of proteins. For instance, the interactions involving Tryptophan greatly decrease.     

Role of medium--and long-range interactions in discriminating globular and membrane proteins.

The analysis of inter-residue interactions in protein structures provides considerable insight to understand their folding and stability. We have previously analyzed the role of medium- and long-range interactions in the folding of globular proteins. In this work, we study the distinct role of such interactions in the three-dimensional structures of membrane proteins. We observed a higher number of long-range contacts in the termini of transmembrane helical (TMH) segments, implying their role in the stabilization of helix-helix interactions. The transmembrane strand (TMS) proteins are having appreciably higher long-range contacts than that in all-beta class of globular proteins, indicating closer packing of the strands in TMS proteins. The residues in membrane spanning segments of TMH proteins have 1.3 times higher medium-range contacts than long-range contacts whereas that of TMS proteins have 14 times higher long-range contacts than medium-range contacts. Residue-wise analysis indicates that in TMH proteins, the residues Cys, Glu, Gly, Pro, Gln, Ser and Tyr have higher long-range contacts than medium-range contacts in contrast with all-alpha class of globular proteins. The charged residue pairs have higher medium-range contacts in all-alpha proteins, whereas hydrophobic residue pairs are dominant in TMH proteins. The information on the preference of residue pairs to form medium-range contacts has been successfully used to discriminate the TMH proteins from all-alpha proteins. The statistical significance of the results obtained from the present study has been verified using randomized structures of TMH and TMS protein templates.     

Influence of medium and long range interactions in protein folding.

Protein structures are stabilized by both local and long range interactions. In this work, we analyze the residue-residue contacts and the role of medium- and long-range interactions in globular proteins belonging to different structural classes. The results show that while medium range interactions predominate in all-alpha class proteins, long-range interactions predominate in all-beta class. Based on this, we analyze the performance of several structure prediction methods in different structural classes of globular proteins and found that all the methods predict the secondary structures of all-alpha proteins more accurately than other classes. Also, we observed that the residues occurring in the range of 21-30 residues apart contributes more towards long-range contacts and about 85% of residues are involved in long-range contacts. Further, the preference of residue pairs to the folding and stability of globular proteins is discussed.     

Importance of long-range interactions in protein folding

Long-range interactions play an active role in the stability of protein molecules. In this work, we have analyzed the importance of long-range interactions in different structural classes of globular proteins in terms of residue distances. We found that 85% of residues are involved in long-range contacts. The residues occurring in the range of 4-10 residues apart contribute more towards long-range contacts in all-alpha proteins while the range is 11-20 in all-beta proteins. The hydrophobic residues Cys, Ile and Val prefer the 11-20 range and all other residues prefer the 4-10 range. The residues in all-beta proteins have an average of 3-8 long-range contacts whereas the residues in other classes have 1-4 long-range contracts. Furthermore, the preference of residue pairs to the folding and stability will be discussed.     

Clustering of non-polar contacts in proteins.

MOTIVATION: Hydrophobic or non-polar contacts in proteins are important for protein folding, protein stability and protein-protein interactions. In particular, in the interior of a protein, in the hydrophobic core, a large number of such contacts are found. The residues involved in these contacts often form a tightly packed cluster of atoms. It is useful for the understanding of protein structure to be able to identify and analyse such clusters. RESULTS: Tools for hierarchical cluster analysis of non-polar contacts in proteins are described. These tools allow for efficient identification of clusters of non-polar interactions in proteins, both internal clusters and clusters involved in protein-protein contacts. The non-polar contacts are represented by a dendrogram structure, which is a simple approach for flexible identification of clusters by visual inspection. The tools are demonstrated on the structure of crambin, the structure of the complex between human growth hormone and the human growth hormone binding protein, and a pair of lipase/esterase structures. Availability: On request from the author.     

Eukaryotic ribonucleases P/MRP: the crystal structure of the P3 domain

Ribonuclease (RNase) P is a site‐specific endoribonuclease found in all kingdoms of life. Typical RNase P consists of a catalytic RNA component and a protein moiety. In the eukaryotes, the RNase P lineage has split into two, giving rise to a closely related enzyme, RNase MRP, which has similar components but has evolved to have different specificities. The eukaryotic RNases P/MRP have acquired an essential helix‐loop‐helix protein‐binding RNA domain P3 that has an important function in eukaryotic enzymes and distinguishes them from bacterial and archaeal RNases P. Here, we present a crystal structure of the P3 RNA domain from Saccharomyces cerevisiae RNase MRP in a complex with RNase P/MRP proteins Pop6 and Pop7 solved to 2.7 Å. The structure suggests similar structural organization of the P3 RNA domains in RNases P/MRP and possible functions of the P3 domains and proteins bound to them in the stabilization of the holoenzymes' structures as well as in interactions with substrates. It provides the first insight into the structural organization of the eukaryotic enzymes of the RNase P/MRP family.     

Influence of medium- and long-range interactions in different folding types of globular proteins

Recognition of protein fold from amino acid sequence is a challenging task. The structure and stability of proteins from different fold are mainly dictated by inter-residue interactions. In our earlier work, we have successfully used the medium- and long-range contacts for predicting the protein folding rates, discriminating globular and membrane proteins and for distinguishing protein structural classes. In this work, we analyze the role of inter-residue interactions in commonly occurring folds of globular proteins in order to understand their folding mechanisms. In the medium-range contacts, the globin fold and four-helical bundle proteins have more contacts than that of DNA-RNA fold although they all belong to all-alpha class. In long-range contacts, only the ribonuclease fold prefers 4-10 range and the other folding types prefer the range 21-30 in alpha/beta class proteins. Further, the preferred residues and residue pairs influenced by these different folds are discussed. The information about the preference of medium- and long-range contacts exhibited by the 20 amino acid residues can be effectively used to predict the folding type of each protein.     

Binding and cleavage of unstructured RNA by nuclear RNase P

Ribonuclease P (RNase P) is an essential endoribonuclease for which the best-characterized function is processing the 5' leader of pre-tRNAs. Compared to bacterial RNase P, which contains a single small protein subunit and a large catalytic RNA subunit, eukaryotic nuclear RNase P is more complex, containing nine proteins and an RNA subunit in Saccharomyces cerevisiae. Consistent with this, nuclear RNase P has been shown to possess unique RNA binding capabilities. To understand the unique molecular recognition of nuclear RNase P, the interaction of S. cerevisiae RNase P with single-stranded RNA was characterized. Unstructured, single-stranded RNA inhibits RNase P in a size-dependent manner, suggesting that multiple interactions are required for high affinity binding. Mixed-sequence RNAs from protein-coding regions also bind strongly to the RNase P holoenzyme. However, in contrast to poly(U) homopolymer RNA that is not cleaved, a variety of mixed-sequence RNAs have multiple preferential cleavage sites that do not correspond to identifiable consensus structures or sequences. In addition, pre-tRNA(Tyr), poly(U)(50) RNA, and mixed-sequence RNA cross-link with purified RNase P in the RNA subunit Rpr1 near the active site in "Conserved Region I," although the exact positions vary. Additional contacts between poly(U)(50) and the RNase P proteins Rpr2p and Pop4p were identified. We conclude that unstructured RNAs interact with multiple protein and RNA contacts near the RNase P RNA active site, but that cleavage depends on the nature of interaction with the active site.     

Making ends meet: the importance of the N- and C-termini for the structure, stability, and function of the third SH3 domain of CIN85

SH3 domains are common structure, interaction, and regulation modules found in more than 200 human proteins. In this report, we studied the third SH3 domain from the human CIN85 adaptor protein, which plays an important role in both receptor tyrosine kinase downregulation and phosphatidylinositol 3 kinase inhibition. The structure of this domain includes an additional 90° kink after the last canonical β-strand and features unusual interactions between the termini well outside the boundaries of the standard SH3 domain definition. The extended portions of the domain are well-structured and held together entirely by side chain-side chain interactions. Extensive expression screening showed that these additional contacts provide significantly increased stability to the domain. A similar 90° kink is found in only one other SH3 domain structure, while side chain contacts linking the termini have never been described before. As a result of the increased size of CIN85 SH3 domain C, the proximal proline rich region is positioned such that a possible intramolecular interaction is structurally inhibited. Using the key interactions of the termini as the basis for sequence analysis allowed the identification of several SH3 domains with flanking sequences that could adopt similar structures. This work illustrates the importance of careful experimental analysis of domain boundaries even for a well-characterized fold such as the SH3 domain.     

Helix-helix packing and interfacial pairwise interactions of residues in membrane proteins

Helix-helix packing plays a critical role in maintaining the tertiary structures of helical membrane proteins. By examining the overall distribution of voids and pockets in the transmembrane (TM) regions of helical membrane proteins, we found that bacteriorhodopsin and halorhodopsin are the most tightly packed, whereas mechanosensitive channel is the least tightly packed. Large residues F, W, and H have the highest propensity to be in a TM void or a pocket, whereas small residues such as S, G, A, and T are least likely to be found in a void or a pocket. The coordination number for non-bonded interactions for each of the residue types is found to correlate with the size of the residue. To assess specific interhelical interactions between residues, we have developed a new computational method to characterize nearest neighboring atoms that are in physical contact. Using an atom-based probabilistic model, we estimate the membrane helical interfacial pairwise (MHIP) propensity. We found that there are many residue pairs that have high propensity for interhelical interactions, but disulfide bonds are rarely found in the TM regions. The high propensity pairs include residue pairs between an aromatic residue and a basic residue (W-R, W-H, and Y-K). In addition, many residue pairs have high propensity to form interhelical polar-polar atomic contacts, for example, residue pairs between two ionizable residues, between one ionizable residue and one N or Q. Soluble proteins do not share this pattern of diverse polar-polar interhelical interaction. Exploratory analysis by clustering of the MHIP values suggests that residues similar in side-chain branchness, cyclic structures, and size tend to have correlated behavior in participating interhelical interactions. A chi-square test rejects the null hypothesis that membrane protein and soluble protein have the same distribution of interhelical pairwise propensity. This observation may help us to understand the folding mechanism of membrane proteins.     

Important inter-residue contacts for enhancing the thermal stability of thermophilic proteins

Proteins from thermophilic organisms exhibit high thermal stability, but have structures that are very similar to their mesophilic homologues. In order to gain insight into the basis of thermostability, we have analyzed the medium- and long-range contacts in mesophilic and thermophilic proteins of 16 different families. We found that the thermophiles prefer to have contacts between residues with hydrogen-bond-forming capability. Apart from hydrophobic contacts, more contacts are observed between polar and non-polar residues in thermophiles than mesophiles. Residue-wise analysis showed that Tyr has good contacts with several other residues, and Cys has considerably higher long-range contacts in thermophiles compared with mesophiles. Furthermore, the residues occurring in the range of 31-34 residues apart in the sequence contribute significant long-range contacts to the stability of thermophilic proteins.     

Characteristic features of amino acid residues in coiled-coil protein structures

Detailed analyses of protein structures provide an opportunity to understand conformation and function in terms of amino acid sequence and composition. In this work, we have systematically analyzed the characteristic features of the amino acid residues found in alpha-helical coiled-coils and, in so doing, have developed indices for their properties, conformational parameters, surrounding hydrophobicity and flexibility. As expected, there is preference for hydrophobic (Ala, Leu), positive (Lys, Arg) and negatively (Glu) charged residues in coiled-coil domains. However, the surrounding hydrophobicity of residues in coiled-coil domains is significantly less than that for residues in other regions of coiled-coil proteins. The analysis of temperature factors in coiled-coil proteins shows that the residues in these domains are more stable than those in other regions. Further, we have delineated the medium- and long-range contacts in coiled-coil domains and compared the results with those obtained for other (non-coiled-coil) parts of the same proteins and non-coiled-coil helical segments of globular proteins. The residues in coiled-coil domains are largely influenced by medium-range contacts, whereas long-range interactions play a dominant role in other regions of these same proteins as well as in non-coiled-coil helices. We have also revealed the preference of amino acid residues to form cation-pi interactions and we found that Arg is more likely to form such interactions than Lys. The parameters developed in this work can be used to understand the folding and stability of coiled-coil proteins in general.     

Protein thermostabilization requires a fine-tuned placement of surface-charged residues

Using the information from the genome projects, recent comparative studies of thermostable proteins have revealed a certain trend of amino acid composition in which polar residues are scarce and charged residues are rich on the protein surface. To clarify experimentally the effect of the amino acid composition of surface residues on the thermostability of Escherichia coli Ribonuclease HI (RNase HI), we constructed six variants in which five to eleven polar residues were replaced by charged residues (5C, 7Ca, 7Cb, 9Ca, 9Cb and 11C). The thermal denaturation experiments indicated that all of the variant proteins are 3.2-10.1 degrees C in Tm less stable than the wild proteins. The crystal structures of resultant protein variants 7Ca, 7Cb, 9Ca and 11C closely resemble that of E. coli RNase HI in their global fold, and several different hydrogen bonding and ion-pair interactions are formed by the mutations. Comparison of the crystal structures of these variant proteins with that of E. coli RNase HI reveals that thermal destabilization is apparently related to electrostatic repulsion of the charged residues with neighbours. This result suggests that charged residues of natural thermostable proteins are strictly posted on the surface with optimal interactions and without repulsive interactions.     

Protein-RNA interactions: structural analysis and functional classes

A data set of 89 protein-RNA complexes has been extracted from the Protein Data Bank, and the nucleic acid recognition sites characterized through direct contacts, accessible surface area, and secondary structure motifs. The differences between RNA recognition sites that bind to RNAs in functional classes has also been analyzed. Analysis of the complete data set revealed that van der Waals interactions are more numerous than hydrogen bonds and the contacts made to the nucleic acid backbone occur more frequently than specific contacts to nucleotide bases. Of the base-specific contacts that were observed, contacts to guanine and adenine occurred most frequently. The most favored amino acid-nucleotide pairings observed were lysine-phosphate, tyrosine-uracil, arginine-phosphate, phenylalanine-adenine and tryptophan-guanine. The amino acid propensities showed that positively charged and polar residues were favored as expected, but also so were tryptophan and glycine. The propensities calculated for the functional classes showed trends similar to those observed for the complete data set. However, the analysis of hydrogen bond and van der Waal contacts showed that in general proteins complexed with messenger RNA, transfer RNA and viral RNA have more base specific contacts and less backbone contacts than expected, while proteins complexed with ribosomal RNA have less base-specific contacts than the expected. Hence, whilst the types of amino acids involved in the interfaces are similar, the distribution of specific contacts is dependent upon the functional class of the RNA bound.     

Statistical analysis of domains in interacting protein pairs

MOTIVATION: Several methods have recently been developed to analyse large-scale sets of physical interactions between proteins in terms of physical contacts between the constituent domains, often with a view to predicting new pairwise interactions. Our aim is to combine genomic interaction data, in which domain-domain contacts are not explicitly reported, with the domain-level structure of individual proteins, in order to learn about the structure of interacting protein pairs. Our approach is driven by the need to assess the evidence for physical contacts between domains in a statistically rigorous way. RESULTS: We develop a statistical approach that assigns p-values to pairs of domain superfamilies, measuring the strength of evidence within a set of protein interactions that domains from these superfamilies form contacts. A set of p-values is calculated for SCOP superfamily pairs, based on a pooled data set of interactions from yeast. These p-values can be used to predict which domains come into contact in an interacting protein pair. This predictive scheme is tested against protein complexes in the Protein Quaternary Structure (PQS) database, and is used to predict domain-domain contacts within 705 interacting protein pairs taken from our pooled data set.     

Probing the structure of complex macromolecular interactions by homolog specificity scanning: the P1 and P7 plasmid partition systems.

The P1 plasmid partition locus, P1 par, actively distributes plasmid copies to Escherichia coli daughter cells. It encodes two DNA sites and two proteins, ParA and ParB. Plasmid P7 uses a similar system, but the key macromolecular interactions are species specific. Homolog specificity scanning (HSS) exploits such specificities to map critical contact points between component macromolecules. The ParA protein contacts the par operon operator for operon autoregulation, and the ParB contacts the parS partition site during partition. Here, we refine the mapping of these contacts and extend the use of HSS to map protein-protein contacts. We found that ParB participates in autoregulation at the operator site by making a specific contact with ParA. Similarly, ParA acts in partition by making a specific contact with ParB bound at parS. Both these interactions involve contacts between a C-terminal region of ParA and the extreme N-terminus of ParB. As a single type of ParA-ParB complex appears to be involved in recognizing both DNA sites, the operator and the parS sites may both be occupied by a single protein complex during partition. The general HSS strategy may aid in solving the three-dimensional structures of large complexes of macromolecules.     

Structural analysis of residues involving cation-pi interactions in different folding types of membrane proteins

Cation-pi interactions play an important role to the stability of protein structures. In our earlier work, we have analyzed the influence and energetic contribution of cation-pi interactions in three-dimensional structures of membrane proteins. In this work, we investigate the characteristic features of residues that are involved in cation-pi interactions. We have computed several parameters, such as surrounding hydrophobicity, number of long-range contacts, conservation score and normalized B-factor for all these residues and identified their location, whether in the membrane or at surface. We found that the cation-pi interactions are mainly formed by long-range interactions. The cationic residues involved in cation-pi interactions have higher surrounding hydrophobicity than their average values in the whole dataset and an opposite trend is observed for aromatic residues. In transmembrane helical proteins, except Phe, all other residues that are responsible for cation-pi interactions are highly conserved with other related protein sequences whereas in transmembrane strand proteins, an appreciable conservation is observed only for Arg. The analysis on the flexibility of residues reveals that the cation-pi interaction forming residues are more stable than other residues. The results obtained in the present study would be helpful to understand the role of cation-pi interactions in the structure and folding of membrane proteins.