The interrelationships of side-chain and main-chain conformations in proteins

Revisiting the Ramachandran plot: hard-sphere repulsion, electrostatics, and H-bonding in the alpha-helix

What determines the shape of the allowed regions in the Ramachandran plot? Although Ramachandran explained these regions in terms of 1-4 hard-sphere repulsions, there are discrepancies with the data where, in particular, the alphaR, alphaL, and beta-strand regions are diagonal. The alphaR-region also varies along the alpha-helix where it is constrained at the center and the amino terminus but diffuse at the carboxyl terminus. By analyzing a high-resolution database of protein structures, we find that certain 1-4 hard-sphere repulsions in the standard steric map of Ramachandran do not affect the statistical distributions. By ignoring these steric clashes (NH(i+1) and O(i-1)C), we identify a revised set of steric clashes (CbetaO, O(i-1)N(i+1), CbetaN(i+1), O(i-1)Cbeta, and O(i-1)O) that produce a better match with the data. We also find that the strictly forbidden region in the Ramachandran plot is excluded by multiple steric clashes, whereas the outlier region is excluded by only one significant steric clash. However, steric clashes alone do not account for the diagonal regions. Using electrostatics to analyze the conformational dependence of specific interatomic interactions, we find that the diagonal shape of the alphaR and alphaL-regions also depends on the optimization of the NH(i+1) and O(i-1)C interactions, and the diagonal beta-strand region is due to the alignment of the CO and NH dipoles. Finally, we reproduce the variation of the Ramachandran plot along the alpha-helix in a simple model that uses only H-bonding constraints. This allows us to rationalize the difference between the amino terminus and the carboxyl terminus of the alpha-helix in terms of backbone entropy.

SCit: web tools for protein side chain conformation analysis

SCit is a web server providing services for protein side chain conformation analysis and side chain positioning. Specific services use the dependence of the side chain conformations on the local backbone conformation, which is described using a structural alphabet that describes the conformation of fragments of four-residue length in a limited library of structural prototypes. Based on this concept, SCit uses sets of rotameric conformations dependent on the local backbone conformation of each protein for side chain positioning and the identification of side chains with unlikely conformations. The SCit web server is accessible at http://bioserv.rpbs.jussieu.fr/SCit.

Exploiting evolutionary relationships for predicting protein structures

In the last few years there have been many developments in computational biology, particularly with regard to novel, imaginative exploitation of genomic data. Disappointingly, there has been a lack of progress in the methodology for prediction of protein structures. In the last several years, however, promising new methods have finally begun to emerge. These methods are increasing the power and scope of the methodology, but, most importantly, they are generating new areas of investigation that we believe will accelerate progress in the field. In this review we describe recent developments and highlight the implications of their success as well as areas where efforts should be focused.

Expanded turn conformations: characterization and sequence-structure correspondence in alpha-turns with implications in helix folding

Like the beta-turns, which are characterized by a limiting distance between residues two positions apart (i, i+3), a distance criterion (involving residues at positions i and i+4) is used here to identify alpha-turns from a database of known protein structures. At least 15 classes of alpha-turns have been enumerated based on the location in the phi,psi space of the three central residues (i+1 to i+3)-one of the major being the class AAA, where the residues occupy the conventional helical backbone torsion angles. However, moving towards the C-terminal end of the turn, there is a shift in the phi,psi angles towards more negative phi, such that the electrostatic repulsion between two consecutive carbonyl oxygen atoms is reduced. Except for the last position (i+4), there is not much similarity in residue composition at different positions of hydrogen and non-hydrogen bonded AAA turns. The presence or absence of Pro at i+1 position of alpha- and beta-turns has a bearing on whether the turn is hydrogen-bonded or without a hydrogen bond. In the tertiary structure, alpha-turns are more likely to be found in beta-hairpin loops. The residue composition at the beginning of the hydrogen bonded AAA alpha-turn has similarity with type I beta-turn and N-terminal positions of helices, but the last position matches with the C-terminal capping position of helices, suggesting that the existence of a "helix cap signal" at i+4 position prevents alpha-turns from growing into helices. Our results also provide new insights into alpha-helix nucleation and folding.

Side-chain conformation angles of amino acids: effect of temperature factor cut-off

The paper presents the analysis of the side-chain conformation angles of amino acids in 90% non-identical protein structures. The analysis has been carried out using 113,699 residues, which is higher compared to the previous studies. In the present study, one more quality check, namely, temperature factor cut-off, has been introduced in addition to resolution and R-factor cut-offs. Due to this, the present calculation reveals the approximate values for the minimum and the maximum of the three-rotameric states of chi1. In addition, the conformation angles chi2 and chi3 have been addressed with the improved data set. The results reported here could be of use in protein modeling.

The occurrence of CH…O hydrogen bonds in α-helices and helix termini in globular proteins

A comprehensive database analysis of C--H...O hydrogen bonds in 3124 alpha-helices and their corresponding helix termini has been carried out from a nonredundant data set of high-resolution globular protein structures resolved at better than 2.0 A in order to investigate their role in the helix, the important protein secondary structural element. The possible occurrence of 5 --> 1 C--H...O hydrogen bond between the ith residue CH group and (i - 4)th residue C==O with C...O < or = 3.8 A is studied, considering as potential donors the main-chain Calpha and the side-chain carbon atoms Cbeta, Cgamma, Cdelta and Cepsilon. Similar analysis has been carried out for 4 --> 1 C--H...O hydrogen bonds, since the C--H...O hydrogen bonds found in helices are predominantly of type 5 --> 1 or 4 --> 1. A total of 17,367 (9310 of type 5 --> 1 and 8057 of type 4 --> 1) C--H...O hydrogen bonds are found to satisfy the selected criteria. The average stereochemical parameters for the data set suggest that the observed C--H...O hydrogen bonds are attractive interactions. Our analysis reveals that the Cgamma and Cbeta hydrogen atom(s) are frequently involved in such hydrogen bonds. A marked preference is noticed for aliphatic beta-branched residue Ile to participate in 5 --> 1 C--H...O hydrogen bonds involving methylene Cgamma 1 atom as donor in alpha-helices. This may be an enthalpic compensation for the greater loss of side-chain conformational entropy for beta-branched amino acids due to the constraint on side-chain torsion angle, namely, chi1, when they occur in helices. The preference of amino acids for 4 --> 1 C--H...O hydrogen bonds is found to be more for Asp, Cys, and for aromatic residues Trp, Phe, and His. Interestingly, overall propensity for C--H...O hydrogen bonds shows that a majority of the helix favoring residues such as Met, Glu, Arg, Lys, Leu, and Gln, which also have large side-chains, prefer to be involved in such types of weak attractive interactions in helices. The amino acid side-chains that participate in C--H...O interactions are found to shield the acceptor carbonyl oxygen atom from the solvent. In addition, C--H...O hydrogen bonds are present along with helix stabilizing salt bridges. A novel helix terminating interaction motif, X-Gly with Gly at C(cap) position having 5 --> 1 Calpha--H...O, and a chain reversal structural motif having 1 --> 5 Calpha-H...O have been identified and discussed. Our analysis highlights that a multitude of local C--H...O hydrogen bonds formed by a variety of amino acid side-chains and Calpha hydrogen atoms occur in helices and more so at the helix termini. It may be surmised that the main-chain Calpha and the side-chain CH that participate in C--H...O hydrogen bonds collectively augment the cohesive energy and thereby contribute together with the classical N--H...O hydrogen bonds and other interactions to the overall stability of helix and therefore of proteins.

Mining a tandem mass spectrometry database to determine the trends and global factors influencing peptide fragmentation

A database of 5500 unique peptide tandem mass spectra acquired in an ion trap mass spectrometer was assembled for peptides derived from proteins digested with trypsin. Peptides were identified initially from their tandem mass spectra by the SEQUEST algorithm and subsequently validated manually. Two different statistical methods were used to identify sequence-dependent fragmentation patterns that could be used to improve fragmentation models incorporated into current peptide sequencing and database search algorithms. The currently accepted "mobile proton" model was expanded to derive a new classification scheme for peptide mass spectra, the "relative proton mobility" scale, which considers peptide ion charge state and amino acid composition to categorize peptide mass spectra into peptide ions containing "nonmobile", "partially mobile", or "mobile" protons. Quantitation of amide bond fragmentation, both N- and C-terminal to any given amino acid, as well as the positional effect of an amino acid in a peptide and peptide length on such fragmentation, has been determined. Peptide bond cleavage propensities, both positive (i.e., enhanced) and negative (i.e., suppressed), were determined and ranked in order of their cleavage preferences as primary, secondary, or tertiary cleavage effects. For example, primary positive cleavage effects were observed for Xaa-Pro and Asp-Xaa bond cleavage for mobile and nonmobile peptide ion categories, respectively. We also report specific pairwise interactions (e.g., Asn-Gly) that result in enhanced amide bond cleavages analogous to those observed in solution-phase chemistry. Peptides classified as nonmobile gave low or insignificant scores, below reported MS/MS score thresholds (cutoff filters), indicating that incorporation of the relative proton mobility scale classification would lead to improvements in current MS/MS scoring functions.

Structural features of transmembrane helices

A total of 160 transmembrane helices of 15 non-homologous high-resolution X-ray protein structures have been analyzed in respect of their structural features. The dihedral angles and hydrogen bonds of the helical sections that span the hydrophobic interior of the lipid bilayer have been investigated. The Ramachandran plot of protein channels and solute transporters exhibit a significant shift Delta (phi- and psi-angles) of Delta mean (+4.5 degrees and -5.4 degrees ), compared to a reference group of 151 alpha-helices of the same average length derived from water-soluble globular proteins. At the C-termini of transmembrane helices structural motifs equivalent to the Gly-caps of helices in globular proteins have been found, with two third of the transmembrane Gly-caps taking up a primary structure that is typically not found at helix termini exposed to a polar solvent. The structural particularities reported here are relevant for the three-dimensional modelling of membrane protein structures.

Beta-sheet preferences from first principles

The natural amino acids have different preferences of occurring in specific types of secondary protein structure. Simulations are performed on periodic model beta-sheets of 14 different amino acids, at the level of density functional theory, employing the generalized gradient approximation. We find that the statistically observed beta-sheet propensities correlate very well with the calculated binding energies. Analysis of the calculations shows that the beta-sheet propensities are determined by the local flexibility of the individual polypeptide strands.

Statistically Based Reduced Representation of Amino Acid Side Chains

Preferred conformations of amino acid side chains have been well established through statistically obtained rotamer libraries. Typically, these provide bond torsion angles allowing a side chain to be traced atom by atom. In cases where it is desirable to reduce the complexity of a protein representation or prediction, fixing all side-chain atoms may prove unwieldy. Therefore, we introduce a general parametrization to allow positions of representative atoms (in the present study, these are terminal atoms) to be predicted directly given backbone atom coordinates. Using a large, culled data set of amino acid residues from high-resolution protein crystal structures, anywhere from 1 to 7 preferred conformations were observed for each terminal atom of the non-glycine residues. Side-chain length from the backbone C(alpha) is one of the parameters determined for each conformation, which should itself be useful. Prediction of terminal atoms was then carried out for a second, nonredundant set of protein structures to validate the data set. Using four simple probabilistic approaches, the Monte Carlo style prediction of terminal atom locations given only backbone coordinates produced an average root mean-square deviation (RMSD) of approximately 3 A from the experimentally determined terminal atom positions. With prediction using conditional probabilities based on the side-chain chi(1) rotamer, this average RMSD was improved to 1.74 A. The observed terminal atom conformations therefore provide reasonable and potentially highly accurate representations of side-chain conformation, offering a viable alternative to existing all-atom rotamers for any case where reduction in protein model complexity, or in the amount of data to be handled, is desired. One application of this representation with strong potential is the prediction of charge density in proteins. This would likely be especially valuable on protein surfaces, where side chains are much less likely to be fixed in single rotamers. Prediction of ensembles of structures provides a method to determine the probability density of charge and atom location; such a prediction is demonstrated graphically.

N-acetyl-N'-methylamide derivative of (2S,3S)-1-amino-2,3-diphenylcyclopropanecarboxylic acid: theoretical analysis of the conformational impact produced by the incorporation of the second phenyl group to the cyclopropane analogue of phenylalanine

The intrinsic conformational preferences of (2S,3S)-1-amino-2,3-diphenylcyclopropanecarboxylic acid, a phenylalanine cyclopropane analogue bearing two phenyl substituents, have been examined theoretically. For this purpose, its N-acetyl-N'-methylamide derivative, Ac-(2S,3S)c(3)diPhe-NHMe, has been investigated by using ab initio HF and DFT methods. Results have been compared with those previously reported for other cyclopropane analogues of phenylalanine, and with experimental data available for c(3)diPhe-containing peptides.

Stereospecific interactions of proline residues in protein structures and complexes

The constrained backbone torsion angle of a proline (Pro) residue has usually been invoked to explain its three-dimensional context in proteins. Here we show that specific interactions involving the pyrrolidine ring atoms also contribute to its location in a given secondary structure and its binding to another molecule. It is adept at participating in two rather non-conventional interactions, C-H...pi and C-H...O. The geometry of interaction between the pyrrolidine and aromatic rings, vis-à-vis the occurrence of the C-H...pi interactions has been elucidated. Some of the secondary structural elements stabilized by Pro-aromatic interactions are beta-turns, where a Pro can interact with an adjacent aromatic residue, and in antiparallel beta-sheet, where a Pro in an edge strand can interact with an aromatic residue in the adjacent strand at a non-hydrogen-bonded site. The C-H groups at the Calpha and Cdelta positions can form strong C-H...O interactions (as seen from the clustering of points) and such interactions involving a Pro residue at C' position relative to an alpha-helix can cap the hydrogen bond forming potentials of the free carbonyl groups at the helix C terminus. Functionally important Pro residues occurring at the binding site of a protein almost invariably engage aromatic residues (with one of them being held by C-H...pi interaction) from the partner molecule in the complex, and such aromatic residues are highly conserved during evolution.

Aromatic-aromatic interactions in crystal structures of helical peptide scaffolds containing projecting phenylalanine residues

Aromatic-aromatic interactions between phenylalanine side chains in peptides have been probed by the structure determination in crystals of three peptides: Boc-Val-Ala-Phe-Aib-Val-Ala-Phe-Aib-OMe, I; Boc-Val-Ala-Phe-Aib-Val-Ala-Phe-Aib-Val-Ala-Phe-Aib-OMe, II; Boc-Aib-Ala-Phe-Aib-Phe-Ala-Val-Aib-OMe, III. X-ray diffraction studies reveal that all three peptides adopt helical conformations in the solid state with the Phe side chains projecting outward. Interhelix association in the crystals is promoted by Phe-Phe interactions. A total of 15 unique aromatic pairs have been characterized in the three independent crystal structures. In peptides I and II, the aromatic side chains lie on the same face of the helix at i/i + 4 positions resulting in both intrahelix and interhelix aromatic interactions. In peptide III, the Phe side chains are placed on the opposite faces of the helix, resulting in exclusive intermolecular aromatic interactions. The distances between the centroids of aromatic pair ranges from 5.11 to 6.86 A, while the distance of closest approach of ring carbon atoms ranges from 3.27 to 4.59 A. Examples of T-shaped and parallel-displaced arrangements of aromatic pairs are observed, in addition to several examples of inclined arrangements. The results support the view that the interaction potential for a pair of aromatic rings is relatively broad and rugged with several minima of similar energies, separated by small activation barriers.

Geometry of interaction of the histidine ring with other planar and basic residues

Among the aromatic residues in protein structures, histidine (His) is unique, as it can exist in the neutral or positively charged form at the physiological pH. As such, it can interact with other aromatic residues as well as form hydrogen bonds with polar and charged (both negative and positive) residues. We have analyzed the geometry of interaction of His residues with nine other planar side chains containing aromatic (residues Phe, Tyr, Trp, and His), carboxylate (Asp and Glu), carboxamide (Asn and Gln) and guanidinium (Arg) groups in 432 polypeptide chains. With the exception of the aspartic (Asp) and glutamic (Glu) acid side-chains, all other residues prefer to interact in a face-to-face or offset-face-stacked orientation with the His ring. Such a geometry is different from the edge-to-face relative orientation normally associated with the aromatic-aromatic interaction. His-His pair prefers to interact in a face-to-face orientation; however, when both the residues bind the same metal ion, the interplanar angle is close to 90 degrees. The occurrence of different interactions (including the nonconventional N-H...pi and C-H...pi hydrogen bonds) have been correlated with the relative orientations between the interacting residues. Several structural motifs, mostly involved in binding metal ions, have been identified by considering the cases where His residues are in contact with four other planar moieties. About 10% of His residues used here are also found in sequence patterns in PROSITE database. There are examples of the amino end of the Lys side chain interacting with His residues in such a way that it is located on an arc around a ring nitrogen atom.

Observation of residual dipolar couplings in short peptides

Residual dipolar couplings provide information on the orientation of individual bond vectors with respect to a unique set of molecular axes. We report that short peptides from 2 to 15 amino acids in length of arbitrary sequence exhibit a modest range of residual dipolar couplings when aligned in either strained polyacrylamide gels or alkyl-PEG bicelles. The absence of significant line broadening in gels suggests peptides align predominantly through steric interactions with the polyacrylamide matrix. However, broadening of NMR lines for a subset of residues aligned in bicelles indicates some peptides bind weakly to these lipid disks, yet a weak negative correlation between the couplings measured in gels and bicelles is consistent with steric hindrance playing a role in both media. The observation of dipolar couplings for peptides of length 10-15 suggests the statistical segment lengths of polypeptide chains must often be >10-15 residues, with data from denatured proteins indicating even larger values. Presumably, local side-chain backbone interactions severely restrict chain flexibility, with the cumulative effect of many such restrictions giving rise to biases in chain direction that may persist for the entire length of a protein chain. Comparison of experimental dipolar couplings for peptides with couplings calculated for ensembles of conformations generated by molecular dynamics should permit evaluation of the accuracy of molecular mechanics potentials in reproducing sequence-specific preferences for phi and psi angles.

Short peptide amyloid organization: stabilities and conformations of the islet amyloid peptide NFGAIL

Experimentally, short peptides have been shown to form amyloids similar to those of their parent proteins. Consequently, they present useful systems for studies of amyloid conformation. Here we simulate extensively the NFGAIL peptide, derived from the human islet amyloid polypeptide (residues 22-27). We simulate different possible strand/sheet organizations, from dimers to nonamers. Our simulations indicate that the most stable conformation is an antiparallel strand orientation within the sheets and parallel between sheets. Consistent with the alanine mutagenesis, we find that the driving force is the hydrophobic effect. Whereas the NFGAIL forms stable oligomers, the NAGAIL oligomer is unstable, and disintegrates very quickly after the beginning of the simulation. The simulations further identify a minimal seed size. Combined with our previous simulations of the prion-derived AGAAAAGA peptide, AAAAAAAA, and the Alzheimer Abeta fragments 16-22, 24-36, 16-35, and 10-35, and the solid-state NMR data for Abeta fragments 16-22, 10-35, and 1-40, some insight into the length and the sequence matching effects may be obtained.

Structure validation by Calpha geometry: phi,psi and Cbeta deviation

Geometrical validation around the Calpha is described, with a new Cbeta measure and updated Ramachandran plot. Deviation of the observed Cbeta atom from ideal position provides a single measure encapsulating the major structure-validation information contained in bond angle distortions. Cbeta deviation is sensitive to incompatibilities between sidechain and backbone caused by misfit conformations or inappropriate refinement restraints. A new phi,psi plot using density-dependent smoothing for 81,234 non-Gly, non-Pro, and non-prePro residues with B < 30 from 500 high-resolution proteins shows sharp boundaries at critical edges and clear delineation between large empty areas and regions that are allowed but disfavored. One such region is the gamma-turn conformation near +75 degrees,-60 degrees, counted as forbidden by common structure-validation programs; however, it occurs in well-ordered parts of good structures, it is overrepresented near functional sites, and strain is partly compensated by the gamma-turn H-bond. Favored and allowed phi,psi regions are also defined for Pro, pre-Pro, and Gly (important because Gly phi,psi angles are more permissive but less accurately determined). Details of these accurate empirical distributions are poorly predicted by previous theoretical calculations, including a region left of alpha-helix, which rates as favorable in energy yet rarely occurs. A proposed factor explaining this discrepancy is that crowding of the two-peptide NHs permits donating only a single H-bond. New calculations by Hu et al. [Proteins 2002 (this issue)] for Ala and Gly dipeptides, using mixed quantum mechanics and molecular mechanics, fit our nonrepetitive data in excellent detail. To run our geometrical evaluations on a user-uploaded file, see MOLPROBITY (http://kinemage.biochem.duke.edu) or RAMPAGE (http://www-cryst.bioc.cam.ac.uk/rampage).

Sequence and structure patterns in proteins from an analysis of the shortest helices: implications for helix nucleation

The shortest helices (three-length 3(10) and four-length alpha), most abundant among helices of different lengths, have been analyzed from a database of protein structures. A characteristic feature of three-length 3(10)-helices is the shifted backbone conformation for the C-terminal residue (phi,psi angles: -95 degrees,0 degrees ), compared to the rest of the helix (-62 degrees,-24 degrees ). The deviation can be attributed to the release of electrostatic repulsion between the carbonyl oxygen atoms at the two C-terminal residues and further stabilization (due to a more linear geometry) of an intrahelical hydrogen bond. A consequence of this non-canonical C-terminal backbone conformation can be a potential origin of helix kinks when a 3(10)-helix is sequence-contiguous at the alpha-helix N-terminal. An analysis of hydrogen bonding, as well as hydrophobic interactions in the shortest helices shows that capping interactions, some of them not observed for longer helices, dominate at the N termini. Further, consideration of the distribution of amino acid residues indicates that the shortest helices resemble the N-terminal end of alpha-helices rather than the C terminus, implying that the folding of helices may be initiated at the N-terminal end, which does not get propagated in the case of the shortest helices. Finally, pairwise comparison of beta-turns and the shortest helices, based on correlation matrices of site-specific amino acid composition, and the relative abundance of these short secondary structural elements, leads to a helix nucleation scheme that considers the formation of an isolated beta-turn (and not an alpha-turn) as the helix nucleation step, with shortest 3(10)-helices as intermediates between the shortest alpha-helix and the beta-turn. Our results ascribe an important role played by shortest 3(10)-helices in proteins with important structural and folding implications.

CADB: Conformation Angles DataBase of proteins

Conformation Angles DataBase (CADB) provides an online resource to access data on conformation angles (both main-chain and side-chain) of protein structures in two data sets corresponding to 25% and 90% sequence identity between any two proteins, available in the Protein Data Bank. In addition, the database contains the necessary crystallographic parameters. The package has several flexible options and display facilities to visualize the main-chain and side-chain conformation angles for a particular amino acid residue. The package can also be used to study the interrelationship between the main-chain and side-chain conformation angles. A web based JAVA graphics interface has been deployed to display the user interested information on the client machine. The database is being updated at regular intervals and can be accessed over the World Wide Web interface at the following URL: http://144.16.71.148/cadb/.

Variants of 3(10)-helices in proteins

An analysis of the shortest 3(10)-helices, containing three helical residues and two flanking capping residues that participate in two consecutive i + 3 --> i hydrogen bonds, shows that not all helices belong to the classic 3(10)-helix, where the three central residues adopt the right-handed helical conformation (alpha(R)). Three variants identified are: 3L10-helix with all residues in the left-handed helical region (alpha(L)), 3EL10-helix where the first residue is in the extended region followed by two residues in the alpha(L) conformation, and its mirror-image, the 3E'R10-helix. In the context of these helices, as well as the equivalent variants of alpha-helices, the length dependence of the handedness of secondary structures in protein structure is discussed. There are considerable differences in the amino acid preferences at different positions in the various types of 3(10)-helices. Each type of 3(10)-helix can be thought to be made up of an extension of a particular type of beta-turn (made up of residues i to i + 3) such that the (i + 3)th residue assumes the same conformation as the preceding residue. Distinct residue preferences at i and i + 3 positions seem to decide whether a particular stretch of four residues will be a beta-turn or a 3(10)-helix in the folded structure.

Rotamer libraries in the 21st century

Rotamer libraries are widely used in protein structure prediction, protein design, and structure refinement. As the size of the structure data base has increased rapidly in recent years, it has become possible to derive well-refined rotamer libraries using strict criteria for data inclusion and for studying dependence of rotamer populations and dihedral angles on local structural features.

On residues in the disallowed region of the Ramachandran map

An analysis of the occurrence of nonglycyl residues in conformations disallowed in the Ramachandran plot is presented. Ser, Asn, Thr, and Cys have the highest propensities to exhibit such conformations, and the branched aliphatic residues the lowest. Residues cluster in five regions and there are some trends in the types of residues and their side-chain conformations (chi(1)) occupying these. Majority of the residues are found at the edge of helices and strands and in short loops, and are involved in different types of weak, stabilizing interactions. A structural motif has been identified where a residue in disallowed conformation occurs as the first residue of a short 3(10)-helix. On the basis of the types of neighboring residues, the location in the three-dimensional structure and accessibility, there are similarities with the occurrence of cis peptide bonds in protein structures.

Aromatic-aromatic interactions in and around alpha-helices

To understand the role of aromatic-aromatic interactions in imparting specificity to the folding process, the geometries of four aromatic residues with different sequence spacing, located in alpha-helices or five residues from helical ends, interacting with each other have been elucidated. The geometry is found to depend on the sequence difference. Specific interactions (C-H...pi and N-H...pi) which result from this geometry may cause a given pair of residues (such as Phe-His) with a particular sequence difference to occur more than expected. The most conspicuous residue in an aromatic pair in the context of helix stability is His, which is found at the last (C1) position or the two positions (Ncap and Ccap) immediately flanking the helix. An alpha-helix and a contiguous 3(10)-helix or two helices separated by a non-helical residue can have interacting aromatic pairs, the geometry of interaction and the relative orientation between the helices being rather fixed. Short helices can also have interacting residues from either side.

The conformations of polypeptide chains where the main-chain parts of successive residues are enantiomeric. Their occurrence in cation and anion-binding regions of proteins

We have investigated the shapes of polypeptides where successive residues have main-chain phi,psi conformations of opposite hand. A graph not unlike a Ramachandran plot is presented illustrating the various possible conformations. All are ring-shaped or extended. Some of these conformations occur in native proteins, the commonest approximating to a feature we propose calling a nest, described in the accompanying paper, where the main-chain NH groups point inwards relative to the ring and give rise to an anion-binding site. Another conformation is related but more extended and is found uniquely in the four stretches of polypeptide that line the tetrameric K(+) channel; their CO groups bind the K ions in the channel. In a different ring-shaped conformation that we propose calling a catgrip, the main-chain CO groups point into the ring; this is employed for specific Ca ion binding in the annexin, phospholipase A2 and subtilisin loops, and the regularly arranged beta-roll loops of the serralysin protease family.