Crystallographic and NMR spectroscopic protein structures: Interresidue contacts

Interresidue pair contacts were analyzed in detail for four pairs of protein structures solved using X-ray analysis (X-ray) and nuclear magnetic resonance (NMR). In the four NMR structures, at distances of ≤4.0 Å, the total number of pair contacts was 4–9% lower and, in general, the pair contacts were 0.02–0.16 Å shorter compared to the X-ray structures. Each of the four structural pairs contained 83–94% common pair contacts (CPCs), which were formed by identical residues in both structures; the other 6–17% were longer intrinsic pair contacts (IPCs) formed by different residues in NMR and X-ray structures, while the latter contained more IPC. Every NMR structure contained three types of CPC that were shorter, longer, or equal to the identical contact pairs in the X-ray structure of this protein. Methodologically different short CPCs prevailed at a known distance dependence of the interresidue contact density in 60–61 pairs of NMR/X-ray structures. Among the analyzed four structural pairs, contact shortening appeared upon the energy minimization of the crambin NMR structure and upon solving the ubiquitin, hen lysozyme, and monomeric hemoglobin NMR structures using X-PLOR software with decreased van der Waals atomic radii. The degree of contact shortening in the NMR structures diminished with an increase in the NMR data used to solve these structures. Among the 60 pairs of NMR/X-ray structures, the major difference between α-helical and β-structural proteins in the dependences on interresidue distances of average contact density appeared due to strong α/β differences in the backbone local geometry.     

The difference between protein structures that are obtained by X-ray analysis and magnetic resonance spectroscopy

We have analyzed the proteins whose structures were determined both by X-ray analysis (X-ray) and nuclear magnetic resonance (NMR) on condition that these structures do not differ greatly when spatially superimposed on each other (61 pairs of protein structures). Atom-atomic contacts (contact distances varied from 2 to 8 A) have been analyzed and it has been found that NMR structures (in comparison with X-ray ones) have more contacts in the range below 3.5 A and above 5.5 A. In the case of residue-residue contacts NMR structures have more contacts below 3 A and between 4.5 and 6.5 A. At all the other contact distances analyzed the X-ray structures have more contacts. The difference in the number of atom-atomic and residue-residue contacts is greater for internal residues, that are concealed from water, as compared to the surface residues. The other, not less important difference deals with the number of hydrogen bonds in the main chain: it is larger for the X-ray structures. The correlation between the hydrogen bonds identified in the structures obtained by both methods is no more than 32%. The consideration of a complete set of protein structures obtained by NMR results in the fact that the number of hydrogen bonds grows 1.2 times as compared to those obtained with the X-ray analysis, whereas the correlation increases only by 65%. We have also demonstrated that alpha-helices in the NMR structures are more distorted in comparison with the ideal alpha-helix, than alpha-helices in the X-ray structures.     

The difference between protein structures obtained by x-ray analysis and nuclear magnetic resonance

We have analyzed the pairs of protein structures obtained by X-ray diffraction analysis and nuclear magnetic resonance (X-ray and NMR structures) that display no major differences when superimposed on one another (61 pairs). Analyzing atom-to-atom contacts (contact distances 2–8 Å), it has been found that the NMR structures (compared to the X-ray structures) have more contacts at distances below 3.5 Å and above 5.5 Å. In the case of residue-to-residue contacts, the NMR structures have more contacts at distances below 3 Å and between 4.5 and 6.5 Å. At other distances analyzed, the X-ray structures have more contacts. The difference in the numbers of atom-to-atom and residue-to-residue contacts is greater for buried residues inaccessible to water than for surface residues. Another important difference is related to the number of hydrogen bonds in the main chain: this number is greater in the X-ray structures. The coefficient of correlation between the numbers of hydrogen bonds identified in the structures obtained by both methods is only 32%. If a complete set of NMR models of protein structure is considered, the total number of hydrogen bonds proves to be 1.2 times greater than in the X-ray structures, whereas the correlation coefficient increases to only 65%. We have also demon-strated that -helices in the NMR structures are more distorted (compared to the ideal -helix) than those in the X-ray structures.Translated from Molekulyarnaya Biologiya, Vol. 39, No. 1, 2005, pp. 129–138.Original Russian Text Copyright © 2005 by Melnik, Garbuzynskiy, Lobanov, Galzitskaya.     

Defrosting the frozen approximation: PROSPECTOR--a new approach to threading

PROSPECTOR (PROtein Structure Predictor Employing Combined Threading to Optimize Results) is a new threading approach that uses sequence profiles to generate an initial probe-template alignment and then uses this "partly thawed" alignment in the evaluation of pair interactions. Two types of sequence profiles are used: the close set, composed of sequences in which sequence identity lies between 35% and 90%; and the distant set, composed of sequences with a FASTA E-score less than 10. Thus, a total of four scoring functions are used in a hierarchical method: the close (distant) sequence profiles screen a structural database to provide an initial alignment of the probe sequence in each of the templates. The same database is then screened with a scoring function composed of sequence plus secondary structure plus pair interaction profiles. This combined hierarchical threading method is called PROSPECTOR1. For the original Fischer database, 59 of 68 pairs are correctly identified in the top position. Next, the set of the top 20 scoring sequences (four scoring functions times the top five structures) is used to construct a protein-specific pair potential based on consensus side-chain contacts occurring in 25% of the structures. In subsequent threading iterations, this protein-specific pair potential, when combined in a composite manner, is found to be more sensitive in identifying the correct pairs than when the original statistical potential is used, and it increases the number of recognized structures for the combined scoring functions, termed PROSPECTOR2, to a total of 61 Fischer pairs identified in the top position. Application to a second, smaller Fischer database of 27 probe-template pairs places 18 (17) structures in the top position for PROSPECTOR1 (PROSPECTOR2). Overall, these studies show that the use of pair interactions as assessed by the improved Z-score enhances the specificity of probe-template matches. Thus, when the hierarchy of scoring functions is combined, the ability to identify correct probe-template pairs is significantly enhanced. Finally, a web server has been established for use by the academic community (http://bioinformatics.danforthcenter.org/services/threading.html).     

Inter-residue interactions in protein folding and stability

During the process of protein folding, the amino acid residues along the polypeptide chain interact with each other in a cooperative manner to form the stable native structure. The knowledge about inter-residue interactions in protein structures is very helpful to understand the mechanism of protein folding and stability. In this review, we introduce the classification of inter-residue interactions into short, medium and long range based on a simple geometric approach. The features of these interactions in different structural classes of globular and membrane proteins, and in various folds have been delineated. The development of contact potentials and the application of inter-residue contacts for predicting the structural class and secondary structures of globular proteins, solvent accessibility, fold recognition and ab initio tertiary structure prediction have been evaluated. Further, the relationship between inter-residue contacts and protein-folding rates has been highlighted. Moreover, the importance of inter-residue interactions in protein-folding kinetics and for understanding the stability of proteins has been discussed. In essence, the information gained from the studies on inter-residue interactions provides valuable insights for understanding protein folding and de novo protein design.     

Comparison of X‐ray and NMR structures: Is there a systematic difference in residue contacts between X‐ray‐ and NMR‐resolved protein structures?

We have compared structures of 78 proteins determined by both NMR and X-ray methods. It is shown that X-ray and NMR structures of the same protein have more differences than various X-ray structures obtained for the protein, and even more than various NMR structures of the protein. X-ray and NMR structures of 18 of these 78 proteins have obvious large-scale structural differences that seem to reflect a difference of crystal and solution structures. The other 60 pairs of structures have only small-scale differences comparable with differences between various X-ray or various NMR structures of a protein; we have analyzed these structures more attentively. One of the main differences between NMR and X-ray structures concerns the number of contacts per residue: (1) NMR structures presented in PDB have more contacts than X-ray structures at distances below 3.0 A and 4.5-6.5 A, and fewer contacts at distances of 3.0-4.5 A and 6.5-8.0 A; (2) this difference in the number of contacts is greater for internal residues than for external ones, and it is larger for beta-containing proteins than for all-alpha proteins. Another significant difference is that the main-chain hydrogen bonds identified in X-ray and NMR structures often differ. Their correlation is 69% only. However, analogous difference is found for refined and rerefined NMR structures, allowing us to suggest that the observed difference in interresidue contacts of X-ray and NMR structures of the same proteins is due mainly to a difference in mathematical treatment of experimental results.     

Identification of amino acids involved in protein structural uniqueness: implication for de novo protein design

Structural uniqueness is characteristic of native proteins and is essential to express their biological functions. The major factors that bring about the uniqueness are specific interactions between hydrophobic residues and their unique packing in the protein core. To find the origin of the uniqueness in their amino acid sequences, we analyzed the distribution of the side chain rotational isomers (rotamers) of hydrophobic amino acids in protein tertiary structures and derived deltaS(contact), the conformational-entropy changes of side chains by residue-residue contacts in each secondary structure. The deltaS(contact) values indicate distinct tendencies of the residue pairs to restrict side chain conformation by inter-residue contacts. Of the hydrophobic residues in alpha-helices, aliphatic residues (Leu, Val, Ile) strongly restrict the side chain conformations of each other. In beta-sheets, Met is most strongly restricted by contact with Ile, whereas Leu, Val and Ile are less affected by other residues in contact than those in alpha-helices. In designed and native protein variants, deltaS(contact) was found to correlate with the folding-unfolding cooperativity. Thus, it can be used as a specificity parameter for designing artificial proteins with a unique structure.     

Structure of crambin in solution, crystal and in the trajectories of molecular dynamics simulations

The mechanisms of the three-dimensional crambin structure alterations in the crystalline environments and in the trajectories of the molecular dynamics simulations in the vacuum and crystal surroundings have been analyzed. In the crystalline state and in the solution the partial regrouping of remote intramolecular packing contacts, involved in the formation and stabilization of the tertiary structure of the crambin molecule, occurs in NMR structures. In the crystalline state it is initiated by the formation of the intermolecular contacts, the conformational influence of its appearance is distributed over the structure. The changes of the conformations and positions of the residues of the loop segments, where the intermolecular contacts of the crystal surroundings are preferably concentrated, are most observable. Under the influence of these contacts the principal change of the regular secondary structure of crambin is taking place: extension of the two-strand β structure to the three-strand structure with the participation of the single last residue N46 of the C-terminal loop. In comparison with the C-terminal loop the more profound changes are observed in the conformation and the atomic positions of the backbone atoms and in the solvent accessibility of the residues of the interhelical loop. In the solution of the ensemble of the 8 NMR structures relative accessibility to the solvent differs more noticeably also in the region of the loop segments and rather markedly in the interhelical loop. In the crambin cryogenic crystal structures the positions of the atoms of the backbone and/or side chain of 14–18 of 46 residues are discretely disordered. The disorganizations of at least 8 of 14 residues occur directly in the regions of the intermolecular contacts and another 5 residues are disordered indirectly through the intramolecular contacts with the residues of the intermolecular contacts. Upon the molecular dynamics simulation in the vacuum surrounding as in the solution of the crystalline structure of crambin the essential changes of the backbone conformation are caused by the intermolecular contacts absence, but partly masked by the structure changes owing to the nonpolar H atoms absence on the simulated structure. The intermolecular contact absence is partly manifested upon the molecular dynamics simulation of the crambin crystal with one protein molecule. Compared to the crystal structure the lengths of the interpeptide hydrogen bonds and other interresidue contacts in an average solution NMR structure are somewhat shorter and accordingly the energy of the interpeptide hydrogen bonds is better. This length shortening can occur at the stage of the refinement of the NMR structures of the crambin and other proteins by its energy minimizations in the vacuum surroundings and not exist in the solution protein structures.     

The effect of backbone on the small-world properties of protein contact maps

In the last years, small-world behavior has been extensively described for proteins, when they are represented by the undirected graph defined by the inter-residue protein contacts. By adopting this representation it was possible to compute the average clustering coefficient (C) and characteristic path length (L) of protein structures, and their values were found to be similar to those of graphs characterized by small-world topology. In this comment, we analyze a large set of non-redundant protein structures (1753) and show that by randomly mimicking the protein collapse, the covalent structure of the protein chain significantly contributes to the small-world behavior of the inter-residue contact graphs. When protein graphs are generated, imposing constraints similar to those induced by the backbone connectivity, their characteristic path lengths and clustering coefficients are indistinguishable from those computed using the real contact maps showing that L and C values cannot be used for 'protein fingerprinting'. Moreover we verified that these results are independent of the selected protein representations, residue composition and protein secondary structures.     

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.     

Concordant exploration of the kinetics of RNA folding from global and local perspectives

Time-resolved small-angle X-ray scattering (SAXS) with millisecond time-resolution reveals two discrete phases of global compaction upon Mg2+-mediated folding of the Tetrahymena thermophila ribozyme. Electrostatic relaxation of the RNA occurs rapidly and dominates the first phase of compaction during which the observed radius of gyration (R(g)) decreases from 75 angstroms to 55 angstroms. A further decrease in R(g) to 45 angstroms occurs in a well-defined second phase. An analysis of mutant ribozymes shows that the latter phase depends upon the formation of long-range tertiary contacts within the P4-P6 domain of the ribozyme; disruption of the three remaining long-range contacts linking the peripheral helices has no effect on the 55-45 angstroms compaction transition. A better understanding of the role of specific tertiary contacts in compaction was obtained by concordant time-resolved hydroxyl radical (OH) analyses that report local changes in the solvent accessibility of the RNA backbone. Comparison of the global and local measures of folding shows that formation of a subset of native tertiary contacts (i.e. those defining the ribozyme core) can occur within a highly compact ensemble whose R(g) is close to that of the fully folded ribozyme. Analyses of additional ribozyme mutants and reaction conditions establish the generality of the rapid formation of a partially collapsed state with little to no detectable tertiary structure. These studies directly link global RNA compaction with formation of tertiary structure as the molecule acquires its biologically active structure, and underscore the strong dependence on salt of both local and global measures of folding kinetics.     

Prediction of protein residue contacts with a PDB-derived likelihood matrix

Proteins with similar folds often display common patterns of residue variability. A widely discussed question is how these patterns can be identified and deconvoluted to predict protein structure. In this respect, correlated mutation analysis (CMA) has shown considerable promise. CMA compares multiple members of a protein family and detects residues that remain constant or mutate in tandem. Often this behavior points to structural or functional interdependence between residues. CMA has been used to predict pairs of amino acids that are distant in the primary sequence but likely to form close contacts in the native three-dimensional structure. Until now these methods have used evolutionary or biophysical models to score the fit between residues. We wished to test whether empirical methods, derived from known protein structures, would provide useful predictive power for CMA. We analyzed 672 known protein structures, derived contact likelihood scores for all possible amino acid pairs, and used these scores to predict contacts. We then tested the method on 118 different protein families for which structures have been solved to atomic resolution. The mean performance was almost seven times better than random prediction. Used in concert with secondary structure prediction, the new CMA method could supply restraints for predicting still undetermined structures.     

Engineering a structure switching mechanism into a steroid-binding aptamer and hydrodynamic analysis of the ligand binding mechanism

The steroid binding mechanism of a DNA aptamer was studied using isothermal titration calorimetry (ITC), NMR spectroscopy, quasi-elastic light scattering (QELS), and small-angle X-ray spectroscopy (SAXS). Binding affinity determination of a series of steroid-binding aptamers derived from a parent cocaine-binding aptamer demonstrates that substituting a GA base pair with a GC base pair governs the switch in binding specificity from cocaine to the steroid deoxycholic acid (DCA). Binding of DCA to all aptamers is an enthalpically driven process with an unfavorable binding entropy. We engineered into the steroid-binding aptamer a ligand-induced folding mechanism by shortening the terminal stem by two base pairs. NMR methods were used to demonstrate that there is a transition from a state where base pairs are formed in one stem of the free aptamer, to where three stems are formed in the DCA-bound aptamer. The ability to generate a ligand-induced folding mechanism into a DNA aptamer architecture based on the three-way junction of the cocaine-binding aptamer opens the door to obtaining a series of aptamers all with ligand-induced folding mechanisms but triggered by different ligands. Hydrodynamic data from diffusion NMR spectroscopy, QELS, and SAXS show that for the aptamer with the full-length terminal stem there is a small amount of structure compaction with DCA binding. For ligand binding by the short terminal stem aptamer, we propose a binding mechanism where secondary structure forms upon DCA binding starting from a free structure where the aptamer exists in a compact form.     

The RNA i-motif

Oligodeoxynucleotides with stretches of cytidine residues associate into a four-stranded structure, the i-motif, in which two head-to-tail, intercalated, parallel-stranded duplexes are held together by hemiprotonated C.C+ pairs. We have investigated the possibility of forming an i-motif structure with C-rich ribonucleic acids. The four C-rich RNAs studied, r(UC5), r(C5), r(C5U) and r(UC3), associate into multiple intercalated structures at acidic pH. r(UC5) forms two i-motif structures that differ by their intercalation topologies. We report on a structural study of the main form and we analyze the small conformational differences found by comparison with the DNA i-motif. The stacking topology of the main structure avoids one of the six 2'-OH/2'-OH repulsive contacts expected in a fully intercalated structure. The C3'-endo pucker of the RNA sugars and the orientation of the intercalated C.C+ pairs result in a modest widening of the narrow grooves at the steps where the hydroxyl groups are in close contact. The free energy of the RNA i-motif, on average -4 kJ mol(-1) per C.C+ pair, is half of the value found in DNA i-motif structures.     

Analysis of multiple folding routes of proteins by a coarse-grained dynamics model.

Langevin dynamics of a protein molecule with Go-type potentials is developed and used to analyze long time-scale events in the folding of cytochrome c. Several trajectories are generated, starting from random coil configurations and going to the native state, that are a few angstroms root mean square deviation (RMSD) from the native structure. The dynamics is controlled, to a large scale, by the two terminal helices that are in contact in the native state. These two helices form very early during folding, and depending on the trajectory, they either stabilize rapidly or break and re-form in going over steric barriers. The extended initial chain exhibits a rapid folding transition into a relatively compact shape, after which the helices are reorganized in a highly correlated manner. The time of formation of residue pair contacts strongly points to the hierarchical nature of folding; i.e., secondary structure forms first, followed by rearrangements of larger length scales at longer times. The kinetics of formation of native contacts is also analyzed, and the onset of a stable globular configuration, referred to as the molten globule in the literature, is identified. Predictions of the model are compared with extensive experimental data on cytochrome c.     

Prediction of the location of structural domains in globular proteins

The location of structural domains in proteins is predicted from the amino acid sequence, based on the analysis of a computed contact map for the protein, the average distance map (ADM). Interactions between residues i and j in a protein are subdivided into several ranges, according to the separation |i-j| in the amino acid sequence. Within each range, average spatial distances between every pair of amino acid residues are computed from a data base of known protein structures. Infrequently occurring pairs are omitted as being statistically insignificant. The average distances are used to construct a predicted ADM. The ADM is analyzed for the occurrence of regions with high densities of contacts (compact regions). Locations of rapid changes of density between various parts of the map are determined by the use of scanning plots of contact densities. These locations serve to pinpoint the distribution of compact regions. This distribution, in turn, is used to predict boundaries of domains in the protein. The technique provides an objective method for the location of domains both on a contact map derived from a known three-dimensional protein structure, the real distance map (RDM), and on an ADM. While most other published methods for the identification of domains locate them in the known three-dimensional structure of a protein, the technique presented here also permits the prediction of domains in proteins of unknown spatial structure, as the construction of the ADM for a given protein requires knowledge of only its amino acid sequence.     

Design of an active ultrastable single-chain insulin analog: synthesis, structure, and therapeutic implications

Single-chain insulin (SCI) analogs provide insight into the inter-relation of hormone structure, function, and dynamics. Although compatible with wild-type structure, short connecting segments (<3 residues) prevent induced fit upon receptor binding and so are essentially without biological activity. Substantial but incomplete activity can be regained with increasing linker length. Here, we describe the design, structure, and function of a single-chain insulin analog (SCI-57) containing a 6-residue linker (GGGPRR). Native receptor-binding affinity (130 +/- 8% relative to the wild type) is achieved as hindrance by the linker is offset by favorable substitutions in the insulin moiety. The thermodynamic stability of SCI-57 is markedly increased (DeltaDeltaG(u) = 0.7 +/- 0.1 kcal/mol relative to the corresponding two-chain analog and 1.9 +/- 0.1 kcal/mol relative to wild-type insulin). Analysis of inter-residue nuclear Overhauser effects demonstrates that a native-like fold is maintained in solution. Surprisingly, the glycine-rich connecting segment folds against the insulin moiety: its central Pro contacts Val(A3) at the edge of the hydrophobic core, whereas the final Arg extends the A1-A8 alpha-helix. Comparison between SCI-57 and its parent two-chain analog reveals striking enhancement of multiple native-like nuclear Overhauser effects within the tethered protein. These contacts are consistent with wild-type crystal structures but are ordinarily attenuated in NMR spectra of two-chain analogs, presumably due to conformational fluctuations. Linker-specific damping of fluctuations provides evidence for the intrinsic flexibility of an insulin monomer. In addition to their biophysical interest, ultrastable SCIs may enhance the safety and efficacy of insulin replacement therapy in the developing world.     

Structure of the Wilms tumor suppressor protein zinc finger domain bound to DNA

The zinc finger domain of the Wilms tumor suppressor protein (WT1) contains four canonical Cys(2)His(2) zinc fingers. WT1 binds preferentially to DNA sequences that are closely related to the EGR-1 consensus site. We report the structure determination by both X-ray crystallography and NMR spectroscopy of the WT1 zinc finger domain in complex with DNA. The X-ray structure was determined for the complex with a cognate 14 base-pair oligonucleotide, and composite X-ray/NMR structures were determined for complexes with both the 14 base-pair and an extended 17 base-pair DNA. This combined approach allowed unambiguous determination of the position of the first zinc finger, which is influenced by lattice contacts in the crystal structure. The crystal structure shows the second, third and fourth zinc finger domains inserted deep into the major groove of the DNA where they make base-specific interactions. The DNA duplex is distorted in the vicinity of the first zinc finger, with a cytidine twisted and tilted out of the base stack to pack against finger 1 and the tip of finger 2. By contrast, the composite X-ray/NMR structures show that finger 1 continues to follow the major groove in the solution complexes. However, the orientation of the helix is non-canonical, and the fingertip and the N terminus of the helix project out of the major groove; as a consequence, the zinc finger side-chains that are commonly involved in base recognition make no contact with the DNA. We conclude that finger 1 helps to anchor WT1 to the DNA by amplifying the binding affinity although it does not contribute significantly to binding specificity. The structures provide molecular level insights into the potential consequences of mutations in zinc fingers 2 and 3 that are associated with Denys-Drash syndrome and nephritic syndrome. The mutations are of two types, and either destabilize the zinc finger structure or replace key base contact residues.     

Hairpin and duplex formation in DNA fragments CCAATTTTGG, CCAATTTTTTGG, and CCATTTTTGG: a proton NMR study

Three DNA fragments, CCAATTTTGG (1), CCAATTTTTTGG (2), and CCATTTTTGG (3), were studied by proton NMR spectroscopy in aqueous solution. All these oligodeoxyribonucleotides contain common sequences at the 5' and 3' ends (5'-CCA and TGG-3'). 2 as well as 3 forms only hairpin structures with four unpaired thymidylyl units, four and three base pair stems, respectively, in neutral solution under low and high NaCl concentrations. At high salt concentration the oligomer 1 forms a duplex structure with -TT- internal loop. On the other hand, the same oligomer forms a stable hairpin structure at low salt and low strand concentrations at pH 7. The hairpin structure of 1 has a stem containing only three base pairs (CCA.TGG) and a loop containing four nucleotides (-ATTT-) that includes a dissociated A.T base pair. The two secondary structures of 1 coexist in an aqueous solution containing 0.1 M NaCl, at pH 7. The equilibrium shifts to the hairpin side when the temperature is raised. The stabilities and base-stacking modes of all three oligonucleotides in two different structures are reported.     

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.