Structural compromise of disallowed conformations in peptide and protein structures

Using a data set of 454 crystal structures of peptides and 80 crystal structures of non-homologous proteins solved at ultra high resolution of 1.2 A or better we have analyzed the occurrence of disallowed Ramachandran (phi, psi) angles. Out of 1492 and 13508 non-glycyl residues in peptides and proteins respectively 12 and 76 residues in the two datasets adopt clearly disallowed combinations of Ramachandran angles. These examples include a number of conformational points which are far away from any of the allowed regions in the Ramachandran map. According to the Ramachandran map a given (phi, psi) combination is considered disallowed when two non-bonded atoms in a system of two-linked peptide units with ideal geometry are prohibitively proximal in space. However, analysis of the disallowed conformations in peptide and protein structures reveals that none of the observations of disallowed conformations in the crystal structures correspond to a short contact between non-bonded atoms. A further analysis of deviations of bond lengths and angles, from the ideal peptide geometry, at the residue positions of disallowed conformations in the crystal structures suggest that individual bond lengths and angles are all within acceptable limits. Thus, it appears that the rare tolerance of disallowed conformations is possible by gentle and acceptable deviations in a number of bond lengths and angles, from ideal geometry, over a series of bonds resulting in a net gross effect of acceptable non-bonded inter-atomic distances.     

Solution structure of native proteins with irregular folds from Raman optical activity

Raman optical activity (ROA) spectra have been measured for the proteins hen phosvitin, yeast invertase, bovine alpha-casein, soybean Bowman-Birk protease inhibitor, and rabbit Cd(7)-metallothionein, all of which have irregular folds in the native state. The results show that ROA is able to distinguish between two types of disorder. Specifically, invertase, alpha-casein, the Bowman-Birk inhibitor, and metallothionein appear to possess a "static" type of disorder similar to that in disordered states of poly(L-lysine) and poly(L-glutamic acid); whereas phosvitin appears to possess a more "dynamic" type of disorder similar to that in reduced (unfolded) lysozyme and ribonuclease A and also in molten globule protein states. In the delimiting cases, static disorder corresponds to that found in loops and turns within native proteins with well-defined tertiary folds that contain sequences of residues with fixed but nonrepetitive phi,psi angles; and dynamic disorder corresponds to that envisaged for the model random coil in which there is a distribution of Ramachandran phi,psi angles for each amino acid residue, giving rise to an ensemble of interconverting conformers. In both cases there is a propensity for the phi,psi angles to correspond to the alpha, beta and poly(L-proline) II (PPII) regions of the Ramachandran surface, as in native proteins with well-defined tertiary folds. Our results suggest that, with the exception of invertase and metallothionein, an important conformational element present in the polypeptide and protein states supporting the static type of disorder is that of the PPII helix. Long sequences of relatively unconstrained PPII helix, as in alpha-casein, may impart a plastic (rheomorphic) character to the structure.     

Polypeptide motions are dominated by peptide group oscillations resulting from dihedral angle correlations between nearest neighbors

To identify basic local backbone motions in unfolded chains, simulations are performed for a variety of peptide systems using three popular force fields and for implicit and explicit solvent models. A dominant "crankshaft-like" motion is found that involves only a localized oscillation of the plane of the peptide group. This motion results in a strong anticorrelated motion of the phi angle of the ith residue (phi(i)) and the psi angle of the residue i - 1 (psi(i-1)) on the 0.1 ps time scale. Only a slight correlation is found between the motions of the two backbone dihedral angles of the same residue. Aside from the special cases of glycine and proline, no correlations are found between backbone dihedral angles that are separated by more than one torsion angle. These short time, correlated motions are found both in equilibrium fluctuations and during the transit process between Ramachandran basins, e.g., from the beta to the alpha region. A residue's complete transit from one Ramachandran basin to another, however, occurs in a manner independent of its neighbors' conformational transitions. These properties appear to be intrinsic because they are robust across different force fields, solvent models, nonbonded interaction routines, and most amino acids.     

How different are structurally flexible and rigid binding sites? Sequence and structural features discriminating proteins that do and do not undergo conformational change upon ligand binding

Proteins are dynamic molecules and often undergo conformational change upon ligand binding. It is widely accepted that flexible loop regions have a critical functional role in enzymes. Lack of consideration of binding site flexibility has led to failures in predicting protein functions and in successfully docking ligands with protein receptors. Here we address the question: which sequence and structural features distinguish the structurally flexible and rigid binding sites? We analyze high-resolution crystal structures of ligand bound (holo) and free (apo) forms of 41 proteins where no conformational change takes place upon ligand binding, 35 examples with moderate conformational change, and 22 cases where a large conformational change has been observed. We find that the number of residue-residue contacts observed per-residue (contact density) does not distinguish flexible and rigid binding sites, suggesting a role for specific interactions and amino acids in modulating the conformational changes. Examination of hydrogen bonding and hydrophobic interactions reveals that cases that do not undergo conformational change have high polar interactions constituting the binding pockets. Intriguingly, the large, aromatic amino acid tryptophan has a high propensity to occur at the binding sites of examples where a large conformational change has been noted. Further, in large conformational change examples, hydrophobic-hydrophobic, aromatic-aromatic, and hydrophobic-polar residue pair interactions are dominant. Further analysis of the Ramachandran dihedral angles (phi, psi) reveals that the residues adopting disallowed conformations are found in both rigid and flexible cases. More importantly, the binding site residues adopting disallowed conformations clustered narrowly into two specific regions of the L-Ala Ramachandran map. Examination of the dihedral angles changes upon ligand binding shows that the magnitude of phi, psi changes are in general minimal, although some large changes particularly between right-handed alpha-helical and extended conformations are seen. Our work further provides an account of conformational changes in the dihedral angles space. The findings reported here are expected to assist in providing a framework for predicting protein-ligand complexes and for template-based prediction of protein function.     

Composites of local structure propensities: Evidence for local encoding of long-range structure

To estimate how extensively the ensemble of denatured-state conformations is constrained by local side-chain–backbone interactions, propensities of each of the 20 amino acids to occur in mono- and dipeptides mapped to discrete regions of the Ramachandran map are computed from proteins of known structure. In addition, propensities are computed for the trans, gauche−, and gauche+ rotamers, with or without consideration of the values of phi and psi. These propensities are used in scoring functions for fragment threading, which estimates the energetic favorability of fragments of protein sequence to adopt the native conformation as opposed to hundreds of thousands of incorrect conformations. As finer subdivisions of the Ramachandran plot, neighboring residue phi/psi angles, and rotamers are incorporated, scoring functions become better at ranking the native conformation as the most favorable. With the best composite propensity function, the native structure can be distinguished from 300,000 incorrect structures for 71% of the 2130 arbitrary protein segments of length 40, 48% of 2247 segments of length 30, and 20% of 2368 segments of length 20. A majority of fragments of length 30–40 are estimated to be folded into the native conformation a substantial fraction of the time. These data suggest that the variations observed in amino acid frequencies in different phi/psi/chi1 environments in folded proteins reflect energetically important local side-chain–backbone interactions, interactions that may severely restrict the ensemble of conformations populated in the denatured state to a relatively small subset with nativelike structure.     

An observation of non-superimposable stereogeometrical features in a non-chiral one-component beta-Ala model peptide

This paper describes the chemical synthesis and crystal molecular conformation of a non-chiral beta-Ala containing model peptide Boc-beta-Ala-Acc5-OCH3. The analysis revealed the existence of two crystallographically independent molecules A and B, in the asymmetric unit. Unexpectedly, while the magnitudes of the backbone torsion angles in both molecules are remarkably similar, the signs of the corresponding torsion angles are reverse therefore, inclining us to suggest the existence of non-superimposable stereogeometrical features in a non-chiral one-component beta-Ala model system. The critical mu torsion angle around CbetaH2-CalphaH2 bond of the beta-Ala residue represents a typical gauche orientation i.e., mu = 67.7 degrees in A and mu = -61.2 degrees in B, providing the molecule an overall crescent shaped topology. The observed conformation contrasts markedly to those determined for the correlated non-chiral model peptides: Boc-beta-Ala-Acc6-OCH3 and Boc-beta-Ala-Aib-OCH3 signifying the role of stereocontrolling elements since the stereochemically constrained Calpha, alpha-disubstituted glycyl residues (e.g., Acc5, Acc6, and the prototype Aib) are known to strongly restrict the peptide backbone conformations in the 3(10)/alpha-helical-regions ( phi approximately +/-60+/-20 degrees, psi approximately +/-30+/-20 degrees) of the Ramachandran map. Unpredictably, the preferred, phi, psi torsion angles of the Acc5 residue fall outside the helical regions of the Ramachandran map and exhibit opposite-handed twists for A and B. The implications of the semi-extended conformation of the Acc5 residue in the construction of backbone-modified novel scaffolds and peptides of biological relevance are highlighted. Taken together, the results indicate that in short linear beta-Ala containing peptides specific structural changes can be induced by selective substitution of non-coded linear- or cyclic symmetrically Calpha,alpha-disubstituted glycines, reinstating the hypothesis that in addition to conformational restrictions, the chemical nature of the neighboring side-chain substituents and local environments collectively influences the stabilization of folding-unfolding behavior of the two methylene units of a beta-Ala residue.     

Influence of solvent molecules on the stereochemical code of glycyl residues in proteins

The Ramachandran steric map and energy diagrams of the glycyl residue are symmetric. A plot of (phi,psi) angles of glycyl residues in 250 nonhomologous and high-resolution protein structures is also largely symmetric. However, there is a clear aberration in the symmetry. Although there is a cluster of points corresponding to the right-handed alpha-helical region, the "equivalent" cluster is clearly shifted to in and around the (phi,psi) values of (90 degrees, 0 degrees ) instead of being centered at the left-handed alpha-helical region of (60 degrees, 40 degrees ). This lack of symmetry exists even in the (phi,psi) distribution of residues from non-alpha-helical regions in proteins. Here we provide an explanation for this observation. An analysis of glycyl conformations in small peptide structures and in "coil" proteins, which are largely devoid of helical and sheet regions, shows that glycyl residues prefer to adopt conformations around (+/-90 degrees, 0 degrees ) instead of right- and left-handed alpha-helical regions. By using theoretical calculations, such conformations are shown to have highest solvent accessibility in a system of two-linked peptide units with glycyl residue at the central C(alpha) atom. This finding is consistent with the observations from 250 nonhomologous protein structures where glycyl residues with conformations close to (+/-90 degrees, 0 degrees ) are seen to have high solvent accessibility. Analysis of a subset of nonhomologous structures with very high resolution (1.5 A or better) shows that water molecules are indeed present at distances suitable for hydrogen bond interaction with glycyl residues possessing conformations close to (+/-90 degrees, 0 degrees ). It is suggested that water molecules play a key role in determining and stabilizing these conformations of glycyl residues and explain the aberration in the symmetry of glycyl conformations in proteins.     

Contribution of tertiary amides to the conformational stability of collagen triple helices

The collagen triple helix is composed of three polypeptide strands, each with a sequence of repeating (Xaa-Yaa-Gly) triplets. In these triplets, Xaa and Yaa are often tertiary amides: L-proline (Pro) and 4(R)-hydroxy-L-proline (Hyp). To determine the contribution of tertiary amides to triple-helical stability, Pro and Hyp were replaced in synthetic collagen mimics with a non-natural acyclic tertiary amide: N-methyl-L-alanine (meAla). Replacing a Pro or Hyp residue with meAla decreases triple-helical stability. Ramachandran analysis indicates that meAla residues prefer to adopt straight phi and psi angles that are dissimilar from those of the Pro and Hyp residues in the collagen triple helix. Replacement with meAla decreases triple-helical stability more than does replacement with Ala. All of the peptide bonds in triple-helical collagen are in the trans conformation. Although an Ala residue greatly prefers the trans conformation, a meAla residue exists as a nearly equimolar mixture of trans and cis conformers. These findings indicate that the favorable contribution of Pro and Hyp to the conformational stability of collagen triple helices arises from factors other than their being tertiary amides.     

CAP: conformation angles package-displaying the conformation angles of side chains in proteins

SUMMARY: A graphics package has been developed to display all side chain conformation angles of the user selected residue in a given protein structure. The proposed package is incorporated with all the protein structures (solved using X-ray diffraction and NMR spectroscopy) available in the Protein Data Bank. The package displays the multiple conformations adopted by a single amino acid residue whose structure is solved and refined at very high resolution. In addition, it shows the percentage distribution of the side chain conformation angles in different rotameric states. AVAILABILITY: http://144.16.71.146/cap/     

Chiral mutagenesis of insulin. Contribution of the B20-B23 beta-turn to activity and stability

Insulin contains a beta-turn (residues B20-B23) interposed between two receptor-binding elements, the central alpha-helix of the B chain (B9-B19) and its C-terminal beta-strand (B24-B28). The turn contains conserved glycines at B20 and B23. Although insulin exhibits marked conformational variability among crystal forms, these glycines consistently maintain positive phi dihedral angles within a classic type-I beta-turn. Because the Ramachandran conformations of GlyB20 and GlyB23 are ordinarily forbidden to L-amino acids, turn architecture may contribute to structure or function. Here, we employ "chiral mutagenesis," comparison of corresponding D- and L-Ala substitutions, to investigate this turn. Control substitutions are introduced at GluB21, a neighboring residue exhibiting a conventional (negative) phi angle. The D- and L-Ala substitutions at B23 are associated with a marked stereospecific difference in activity. Whereas the D-AlaB23 analog retains native activity, the L analog exhibits a 20-fold decrease in receptor binding. By contrast, D- and L-AlaB20 analogs each exhibit high activity. Stereospecific differences between the thermodynamic stabilities of the analogs are nonetheless more pronounced at B20 (delta deltaG(u) 2.0 kcal/mole) than at B23 (delta deltaG(u) 0.7 kcal/mole). Control substitutions at B21 are well tolerated without significant stereospecificity. Chiral mutagenesis thus defines the complementary contributions of these conserved glycines to protein stability (GlyB20) or receptor recognition (GlyB23).     

Crystallographic characterization of conformation of 1-aminocyclopropane-1-carboxylic acid residue (Ac3c) in simple derivatives and peptides

The molecular and crystal structures of the C alpha,alpha-dialkylated alpha-amino acid residue 1-aminocyclopropane-1-carboxylic acid hemihydrate (H2+-Ac3c-O-.1/2 H2O) and nine derivatives and dipeptides have been determined by X-ray diffraction. The derivatives are pBrBz-Ac3c-OH, Piv-Ac3c-OH, Z-Ac3c-OH, the alpha-and beta-forms of t-Boc-Ac3c-OH, Z-Ac3c-OMe, and the 5(4H)-oxazolone from pBrBz-Ac3c-OH; the dipeptides are H-(Ac3c)-OMe and c(Ac3c)2. The values determined for the torsion angles about the N-C alpha (phi) and C alpha-C' (psi) bonds for the single Ac3c residue of Piv-Ac3c-OH, the alpha- and beta-forms of t-Boc-Ac3-OH and Z-Ac3c-OMe, and the C-terminal Ac3c residue of H-(Ac3c)2-OMe correspond to folded conformations in the "bridge" region of the Ramachandran map. The structures of pBrBz-Ac3c-OH and Z-Ac3c-OH, however, are unusual in having a semi-extended conformation for the phi, psi angles. The N-terminal Ac3c residue of H-(Ac3c)2-OMe adopts a novel type of C5 conformation, characterized inter alia by an (amino) N. . .H-N (peptide) intramolecular hydrogen bond. While the acyl N alpha-blocking groups form trans amides (pBrBz-Ac3c-OH and Piv-Ac3c-OH), the urethane groups may adopt either the trans [Z-Ac3c-OH and t-Boc-Ac3c-OH (alpha-form)] or the cis amide conformations [t-Boc-Ac3c-OH(beta-form) and Z-Ac3c-OMe]. The five- and six-membered rings of the 5(4H)-oxazolone and the 2,5-dioxopiperazine, respectively, are planar. The four independent molecules in the asymmetric unit of the free alpha-amino acid are zwitterionic.     

Structural assessment of glycyl mutations in invariantly conserved motifs

Motifs that are evolutionarily conserved in proteins are crucial to their structure and function. In one of our earlier studies, we demonstrated that the conserved motifs occurring invariantly across several organisms could act as structural determinants of the proteins. We observed the abundance of glycyl residues in these invariantly conserved motifs. The role of glycyl residues in highly conserved motifs has not been studied extensively. Thus, it would be interesting to examine the structural perturbations induced by mutation in these conserved glycyl sites. In this work, we selected a representative set of invariant signature (IS) peptides for which both the PDB structure and mutation information was available. We thoroughly analyzed the conformational features of the glycyl sites and their local interactions with the surrounding residues. Using Ramachandran angles, we showed that the glycyl residues occurring in these IS peptides, which have undergone mutation, occurred more often in the L-disallowed as compared with the L-allowed region of the Ramachandran plot. Short range contacts around the mutation site were analyzed to study the steric effects. With the results obtained from our analysis, we hypothesize that any change of activity arising because of such mutations must be attributed to the long-range interaction(s) of the new residue if the glycyl residue in the IS peptide occurred in the L-allowed region of the Ramachandran plot. However, the mutation of those conserved glycyl residues that occurred in the L-disallowed region of the Ramachandran plot might lead to an altered activity of the protein as a result of an altered conformation of the backbone in the immediate vicinity of the glycyl residue, in addition to long range effects arising from the long side chains of the new residue. Thus, the loss of activity because of mutation in the conserved glycyl site might either relate to long range interactions or to local perturbations around the site depending upon the conformational preference of the glycyl residue.     

Calibration of effective van der Waals atomic contact radii for proteins and peptides

Effective van der Waals radii were calibrated in such a way that molecular models built from standard bond lengths and bond angles reproduced the amino acid conformations observed by crystallography in proteins and peptides. The calibrations were based on the comparison of the Ramachandran plots prepared from high-resolution X-ray data of proteins and peptides with the allowed phi, psi torsional angle space for the dipeptide molecular models. The calibrated radii are useful as criteria with which to filter energetically improbable conformations in molecular modeling studies of proteins and peptides.     

Comprehensive analysis of the helix-X-helix motif in soluble proteins

Alpha-helices are the most common secondary structures found in globular proteins. In this report, we analyze the stereochemical and sequence properties of helix-X-helix (HXH) motifs in which two alpha-helices are linked by a single residue, in search of characteristic structures and sequence signals. The analysis is carried out on a database of 837 nonredundant HXH motifs. The kinks are characterized by the bend angle between the axes of the N-terminal and C-terminal helices and the wobble angle corresponding to the rotation of C-terminal helix axis on the plane perpendicular to the N-terminal one. The phi-psi dihedral angles of the linker residue are clustered in six distinct areas of the Ramachandran plot: two areas are located in the additional allowed alpha region (alpha(1) and alpha(2)), two areas are in the additional allowed beta region (beta(1) and beta(2)) and two areas have positive phi values (alpha(L) and beta(M)). Each phi/psi region corresponds to characteristic bend and wobble angles and amino acid distributions. Bend angles can vary from 0 degrees to 160 degrees. Most wobble angles correspond to a counter-clockwise rotation of the C-terminal helix. Proline residues are rigorously excluded from the linker position X but have a high propensity at position X+1 of the beta(1) and beta(2) motifs (12 and 7, respectively) and at position X+3 of the alpha(1) motifs (9). Glycine linkers are located either in the alpha(L) region (20%) or in the beta(M) region (80%). This latter conformation is characterized by a marked bend angle (124 degrees +/- 18 degrees) and a clockwise wobble. Among other amino acids, Asn is remarkable for its high propensity (>3) at the linker position of the alpha(2), beta(1), and beta(2) motifs. Stabilization of HXH motifs by H-bonds between polar side chains of the linker and polar groups of the backbone is determined. A method based on position-specific scoring matrices is developed for conformational prediction. The accuracy of the predictions reaches 80% when the method is applied to proline-induced kinks or to kinks with bend angles in the 50 degrees-100 degrees range.     

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

The accurate determination of a large number of protein structures by X-ray crystallography makes it possible to conduct a reliable statistical analysis of the distribution of the main-chain and side-chain conformational angles, how these are dependent on residue type, adjacent residue in the sequence, secondary structure, residue-residue interactions and location at the polypeptide chain termini. The interrelationship between the main-chain (phi, psi) and side-chain (chi 1) torsion angles leads to a classification of amino acid residues that simplify the folding alphabet considerably and can be a guide to the design of new proteins or mutational studies. Analyses of residues occurring with disallowed main-chain conformation or with multiple conformations shed some light on why some residues are less favoured in thermophiles.     

A Bayesian-probability-based method for assigning protein backbone dihedral angles based on chemical shifts and local sequences

Chemical shifts contain substantial information about protein local conformations. We present a method to assign individual protein backbone dihedral angles into specific regions on the Ramachandran map based on the amino acid sequences and the chemical shifts of backbone atoms of tripeptide segments. The method uses a scoring function derived from the Bayesian probability for the central residue of a query tripeptide segment to have a particular conformation. The Ramachandran map is partitioned into representative regions at two levels of resolution. The lower resolution partitioning is equivalent to the conventional definitions of different secondary structure regions on the map. At the higher resolution level, the α and β regions are further divided into subregions. Predictions are attempted at both levels of resolution. We compared our method with TALOS using the original TALOS database, and obtained comparable results. Although TALOS may produce the best results with currently available databases which are much enlarged, the Bayesian-probability-based approach can provide a quantitative measure for the reliability of predictions.     

Ramachandran plot on the web (2.0)

The Ramachandran plot displays the main chain conformation angles (Phi and Psi) of the polypeptide chain of a protein molecule. The paper reports the updated version of the Ramachandran plot web server and has several improved options for displaying the conformation angles in various regions. In addition, options are provided to display the conformation angles in various secondary structural elements and regions within the user specified Phi and Psi values in the plot. The updated version is accessible at the following URL: http://dicsoft1.physics.iisc.ernet.in/rp/.     

Implication of an unfavorable residue (Thr346) in intrinsic flexibility of firefly luciferase

In order to better understand the functional role of an unusual residue (Thr346) of firefly luciferase mutagenesis at this residue was performed. Firefly luciferase, catalyzes the bioluminescence reaction and is an excellent tool as a reporter in nano-system biology studies. Nonetheless, the enzyme rapidly loses its activity at temperatures above 30°C and this leads to reduced sensitivity and precision in analytical applications. Residue Thr346 in a connecting loop (341-348) of firefly luciferase is located in a disallowed region of Ramachandran plot. In this study, we have substituted this residue (T346) with anomalous dihedral angles with Val, Gly and Pro to clarify the role of this residue in structure and function of the enzyme using site-directed mutagenesis. Substitution of this unfavorable residue (T346) with atypical dihedral angles (ψ, φ) with other residues brought about an increase of thermostability and decrease of specific activity. Structural and functional properties of the mutants were analyzed using different spectroscopic methods. It seems that this residue is a critically conserved residue to support the functional flexibility for a fast kinetic bioluminescence reaction at the expense of lower stability.     

Thermodynamic and structural characterization of Asn and Ala residues in the disallowed II' region of the Ramachandran plot

Residue Asn47 at position L1 of a type II' beta-turn of the alpha-spectrin SH3 domain is located in a disallowed region of the Ramachandran plot (phi = 56 +/- 12, psi = -118 +/- 17). Therefore, it is expected that replacement of Asn47 by Gly should result in a considerable stabilization of the protein. Thermodynamic analysis of the N47G and N47A mutants shows that the change in free energy is small (approximately 0.7 kcal/mol; approximately 3 kJ/mol) and comparable to that found when mutating a Gly to Ala in a alpha-helix or beta-sheet. X-ray structural analysis of these mutants shows that the conformation of the beta-turn does not change upon mutation and, therefore, that there is no relaxation of the structure, nor is there any gain or loss of interactions that could explain the small energy change. Our results indicate that the energetic definition of II' region of the Ramachandran plot (phi = 60 +/- 30, psi = -115 +/- 15) should be revised for at least Ala and Asn in structure validation and protein design.