Cis peptide bonds in proteins: residues involved, their conformations, interactions and locations

Enhanced stability of cis Pro-Pro peptide bond in Pro-Pro-Phe sequence motif

Identification of sequence motifs that favor cis peptide bonds in proteins is important for understanding and designing proteins containing turns mediated by cis peptide conformations. From (1)H NMR solution studies on short peptides, we show that the Pro-Pro peptide bond in Pro-Pro-Phe almost equally populates the cis and trans isomers, with the cis isomer stabilized by a CHc...pi interaction involving the terminal Pro and Phe. We also show that Phe is over-represented at sequence positions immediately following cis Pro-Pro motifs in known protein structures. Our results demonstrate that the Pro-Pro cis conformer in Pro-Pro-Phe sequence motifs is as important as the trans conformer, both in short peptides as well as in natively folded proteins.

Puckering Transition of Proline Residue in Water

The puckering transition of the proline residue with trans and cis prolyl peptide bonds was explored by optimizations along the torsion angle chi1 of the prolyl ring using quantum-chemical methods in water. By analyzing the potential energy surfaces and local minima in water, it is observed that the puckering transition of the proline residue proceeds from a down-puckered conformation to an up-puckered one and vice versa through the transition state with an envelope form having the N atom at the top of the envelope and not a planar one, as seen in the gas phase, although the backbone conformations are different in the gas phase and in water. The barriers to the puckering transition DeltaGup-->down are estimated to be 3.12 and 3.00 kcal/mol for trans and cis conformers at the B3LYP/6-311++G(d,p) level of theory in water, respectively, which are about 1.7 kcal/mol higher than those in the gas phase. Out of 2197 prolines from the 241 high-resolution PDB chains, four transition-state-like structures with the envelope ring puckering are identified. Three of them have the trans prolyl peptide bonds and one has the cis one. The favorable or steric interactions by neighboring residues may be responsible for the stabilization of these transition-state-like ring structures in the proteins.

Geometry of nonbonded interactions involving planar groups in proteins

Although hydrophobic interaction is the main contributing factor to the stability of the protein fold, the specificity of the folding process depends on many directional interactions. An analysis has been carried out on the geometry of interaction between planar moieties of ten side chains (Phe, Tyr, Trp, His, Arg, Pro, Asp, Glu, Asn and Gln), the aromatic residues and the sulfide planes (of Met and cystine), and the aromatic residues and the peptide planes within the protein tertiary structures available in the Protein Data Bank. The occurrence of hydrogen bonds and other nonconventional interactions such as C-H...pi, C-H...O, electrophile-nucleophile interactions involving the planar moieties has been elucidated. The specific nature of the interactions constraints many of the residue pairs to occur with a fixed sequence difference, maintaining a sequential order, when located in secondary structural elements, such as alpha-helices and beta-turns. The importance of many of these interactions (for example, aromatic residues interacting with Pro or cystine sulfur atom) is revealed by the higher degree of conservation observed for them in protein structures and binding regions. The planar residues are well represented in the active sites, and the geometry of their interactions does not deviate from the general distribution. The geometrical relationship between interacting residues provides valuable insights into the process of protein folding and would be useful for the design of protein molecules and modulation of their binding properties.

Conformations of proline analogues having double bonds in the ring

The intrinsic conformational preferences of proline analogues having double bonds between carbon atoms in their rings have been investigated using quantum mechanical calculations at the B3LYP/6-31+G(d,p) level. For this purpose, the potential energy surface of the N-acety-N'-methylamide derivatives of three dehydroprolines (proline analogues unsaturated at alpha,beta; beta,gamma; and gamma,delta) and pyrrole (proline analogue with unsaturations at both alpha,beta and gamma,delta) have been explored, and the results are compared with those obtained for the derivative of the nonmodified proline. We found that the double bonds affect the ring puckering and the geometric internal parameters, even though the backbone conformation was influenced the most. Results indicate that the formation of double bonds between carbon atoms in the pyrrolidine ring should be considered as an effective procedure to restrict the conformational flexibility of prolines. Interestingly, we also found that the N-acetyl-N'-methylamide derivative of pyrrole shows a high probability of having a cis peptide bond preceding the proline analogue.

Localization of a conformational epitope common to non-native and fibrillar immunoglobulin light chains

Amyloid fibrils and partially unfolded intermediates may be distinguished serologically from native amyloidogenic precursor proteins or peptides. In this regard, we had previously reported that the IgG1 mAb 11-1F4, generated by immunizing mice with a thermally denatured variable region fragment of the human Igkappa4 Bence Jones protein Len, reacted specifically with light chain (LC) fibrils, irrespective of kappa or lambda isotype but, notably, did not with native molecules (Hrncic, R. et al. (2000) Am. J. Pathol. 157, 1239-1246). To elucidate the molecular basis of this specificity, we have used a europium-linked fluorescent immunoassay, where it was demonstrated through epitope mapping that mAb 11-1F4 recognizes a conformational determinant contained within the first (N-terminal) 18 amino acids of misfolded LCs. The nature of this epitope was evidenced in competition studies where the peptide Len (1-18), but not the intact protein or other LCs, inhibited the binding of the antibody to fibrils. This unique reactivity was dependent on the structural integrity of this portion of the molecule, particularly the presence of a highly conserved prolyl residue at position 8. On the basis of our experimental data, we posit that the mAb 11-1F4 binding site found on partially denatured and fibrillar LCs involves an inducible N-terminal main chain reversal that results in the formation of a proline anchored beta-turn. Our delineation of this LC fibril-associated epitope provides the rationale for the design of novel amyloid-reactive antibodies with diagnostic and therapeutic potential for patients with LC-associated and other forms of amyloidosis.

Interaction geometry involving planar groups in protein-protein interfaces

The geometry of interactions of planar residues is nonrandom in protein tertiary structures and gives rise to conventional, as well as nonconventional (X--H...pi, X--H...O, where X = C, N, or O) hydrogen bonds. Whether a similar geometry is maintained when the interaction is across the protein-protein interface is addressed here. The relative geometries of interactions involving planar residues, and the percentage of contacts giving rise to different types of hydrogen bonds are quite similar in protein structures and the biological interfaces formed by protein chains in homodimers and protein-protein heterocomplexes--thus pointing to the similarity of chemical interactions that occurs during protein folding and binding. However, the percentage is considerably smaller in the nonspecific and nonphysiological interfaces that are formed in crystal lattices of monomeric proteins. The C--H...O interaction linking the aromatic and the peptide groups is quite common in protein structures as well as the three types of interfaces. However, as the interfaces formed by crystal contacts are depleted in aromatic residues, the weaker hydrogen bond interactions would contribute less toward their stability.

Synthesis of a new π-deficient phenylalanine derivative from a common 1,4-diketone intermediate and study of the influence of aromatic density on prolyl amide isomer population

Enantiopure (2S)-N-(Boc)-3-(6-methylpyridazinyl)alanine (14) has been synthesized to serve as a phenylalanine analog lacking significant pi-donor capability. Two approaches were developed to furnish the target compound from L-aspartic acid as chiral educt in respectively six and nine steps and 13% and 12% yields. In both routes, a key homoallylic ketone intermediate was synthesized by a copper-catalyzed cascade addition of vinylmagnesium bromide to a carboxylic ester. Dipeptide models Ac-Xaa-Pro-NHMe (21a-c) were prepared and the relative populations of prolyl cis- and trans-amide isomers were measured in chloroform, dimethylsulfoxide, and water by proton NMR spectroscopy in order to assess the significance of the electron density of the neighboring aromatic residue on the prolyl amide geometry.

Conformational Preferences of Non-Prolyl and Prolyl Residues

The conformational study on Ac-Ala-NHMe (the alanine dipeptide) and Ac-Pro-NHMe (the proline dipeptide) is carried out using ab initio HF and density functional methods with the self-consistent reaction field method to explore the differences in the backbone conformational preference and the cis-trans isomerization for the non-prolyl and prolyl residues in the gas phase and in the solutions (chloroform and water). For the alanine and proline dipeptides, with the increase of solvent polarity, the populations of the conformation tC with an intramolecular C(7) hydrogen bond significantly decrease, and those of the polyproline II-like conformation tF and the alpha-helical conformation tA increase, which is in good agreement with the results from circular dichroism and NMR experiments. For both the dipeptides, as the solvent polarity increases, the relative free energy of the cis conformer to the trans conformer decreases and the rotational barrier to the cis-trans isomerization increases. It is found that the cis-trans isomerization proceeds in common through only the clockwise rotation with omega' approximately +120 degrees about the non-prolyl and prolyl peptide bonds in both the gas phase and the solutions. The pertinent distance d(N...H-N(NHMe)) can successfully describe the increase in the rotational barriers for the non-prolyl and prolyl trans-cis isomerization as the solvent polarity increases and the higher barriers for the non-prolyl residue than for the prolyl residue, as seen in experimental and calculated results. By analysis of the contributions to rotational barriers, the cis-trans isomerization for the non-prolyl and prolyl peptide bonds is proven to be entirely enthalpy driven in the gas phase and in the solutions. The calculated cis populations and rotational barriers to the cis-trans isomerization for both the dipeptides in chloroform and/or water accord with the experimental values.

First-Principle Computational Study on the Full Conformational Space of l-Threonine Diamide, the Energetic Stability of Cis and Trans Isomers

First-principle computations were carried out on the conformational space of trans and cis peptide bond isomers of HCO-Thr-NH2. Using the concept of multidimensional conformational analysis (MDCA), geometry optimizations were performed at the B3LYP/6-31G(d) level of theory, and single-point energies as well as thermodynamic functions were calculated at the G3MP2B3 level of theory for the corresponding optimized structures. Two backbone Ramachandran-type potential energy surfaces (PESs) were computed, one each for the cis and trans isomers, keeping the side chain at the fully extended orientation (chi1=chi2=anti). Similarly, two side chain PESs for the cis and trans isomers were generated for the (phi=psi=anti) orientation corresponding to approximately the betaL backbone conformation. Besides correlating the relative Gibbs free energy of the various stable conformations with the number of stabilizing hydrogen bonds, the process of trans-->cis isomerization is discussed in terms of intrinsic stabilities as measured by the computed thermodynamic functions.

Conformational Preferences of Proline Oligopeptides

A conformational study on the terminally blocked proline oligopeptides, Ac-(Pro)(n)()-NMe(2) (n = 2-5), is carried out using the ab initio Hartree-Fock level of theory with the self-consistent reaction field method in the gas phase and in solutions (chloroform, 1-propanol, and water) to explore the preference and transition between polyproline II (PPII) and polyproline I (PPI) conformations depending on the chain length, the puckering, and the solvent. The mean differences in the free energy per proline of the up-puckered conformations relative to the down-puckered conformations for both diproline and triproline increases for the PPII-like conformations and decreases for the PPI-like conformations as the solvent polarity increases. These calculated results indicate that the PPII-like structures have preferentially all-down puckerings in solutions, whereas the PPI-like structures have partially mixed puckerings. The free energy difference per proline residue between the PPII- and PPI-like structures decreases as the proline chain becomes longer in the gas phase but increases as the proline chain becomes longer in solutions and the solvent polarity increases. In particular, our calculated results indicate that each of the proline oligopeptides can exist as an ensemble of conformations with the trans and cis peptide bonds in solutions, although the PPII-like structure with all-trans peptide bonds is dominantly preferred, which is reasonably consistent with the previously observed results. In diproline Ac-(Pro)(2)-NMe(2), the rotational barrier to the cis-to-trans isomerization for the first prolyl peptide bond increases as the solvent polarity increases, whereas the rotational barrier for the second prolyl peptide bond does not show the monotonic increase as the solvent polarity increases. When the rotational barriers for these two prolyl peptide bonds were compared, it could be deduced that the conformational transition from PPI with the cis peptide bond to PPII with the trans peptide bond is initiated at the C-terminus and proceeds to the N-terminus in water. This is consistent with the results from NMR experiments on polyproline in D(2)O but opposite to the results from enzymatic hydrolysis kinetics experiments on polyproline.

Impact of cis-proline analogs on peptide conformation

The beta-turn is a common motif in both proteins and peptides and often a recognition site in protein interactions. A beta-turn of four sequential residues reverses the direction of the peptide chain and is classified by the phi and psi backbone torsional angles of residues i + 1 and i + 2. The type VI turn usually contains a proline with a cis-amide bond at residue i + 2. Cis-proline analogs that constrain the peptide to adopt a type VI turn led to peptidomimetics with enhanced activity or metabolic stability. To compare the impact of different analogs on amide cis-trans isomerism and peptide conformation, the conformational preference for the cis-amide bond and the type VI turn was investigated at the MP2/6-31+G** level of theory in water (polarizable continuum water model). Analogs stabilize the cis-amide conformations through different mechanisms: (1) 5-alkylproline, with bulky hydrocarbon substituent on the C(delta) of proline, increases the cis-amide population through steric hindrance between the alkyl substituent and the N-terminal residues; (2) oxaproline or thioproline, the oxazolidine- or thiazolidine-derived proline analog, favors interactions between the dipole of the heterocyclic ring and the preceding carbonyl oxygen; and (3) azaproline, containing a nitrogen atom in place of the C(alpha) of proline, prefers the cis-amide bond by lone-pair repulsion between the alpha-nitrogen and the preceding carbonyl oxygen. Preference for the cis conformation was augmented by combining different modifications within a single proline. Azaproline and its derivatives are most effective in stabilizing cis-amide bonds without introducing additional steric bulk to compromise receptor interactions.

Effects of i and i+3 residue identity on cis-trans isomerism of the aromatic(i+1)-prolyl(i+2) amide bond: implications for type VI beta-turn formation

Cis-trans isomerization of amide bonds plays critical roles in protein molecular recognition, protein folding, protein misfolding, and disease. Aromatic-proline sequences are particularly prone to exhibit cis amide bonds. The roles of residues adjacent to a tyrosine-proline residue pair on cis-trans isomerism were examined. A short series of peptides XYPZ was synthesized and cis-trans isomerism was analyzed. Based on these initial studies, a series of peptides XYPN, X = all 20 canonical amino acids, was synthesized and analyzed by NMR for i residue effects on cis-trans isomerization. The following effects were observed: (a) aromatic residues immediately preceding Tyr-Pro disfavor cis amide bonds, with K(trans/cis)= 5.7-8.0, W > Y > F; (b) proline residues preceding Tyr-Pro lead to multiple species, exhibiting cis-trans isomerization of either or both X-Pro amide bonds; and (c) other residues exhibit similar values of K(trans/cis) (= 2.9-4.2), with Thr and protonated His exhibiting the highest fraction cis. beta-Branched and short polar residues were somewhat more favorable in stabilizing the cis conformation. Phosphorylation of serine at the i position modestly increases the stability of the cis conformer. In addition, the effect of the i+3 residue was examined in a limited series of peptides TYPZ. NMR data indicated that aromatic residues, Pro, Asn, Ala, and Val at the i+3 residue all favor cis amide bonds, with aromatic residues and Asn favoring more compact phi at Tyr(cis) and Ala and Pro favoring more extended phi at Tyr(cis). D-Alanine at the i+3 position particularly disfavors cis amide bonds.

Conservation of cis prolyl bonds in proteins during evolution

In proteins and peptides, the vast majority of peptide bonds occurs in trans conformation, but a considerable fraction (about 5%) of X-Pro bonds adopts the cis conformation. Here we study the conservation of cis prolyl residues in evolutionary related proteins. We find that overall, in contrast to local, protein sequence similarity is a clear indicator for the conformation of prolyl residues. We observe that cis prolyl residues are more often conserved than trans prolyl residues, and both are more conserved than the surrounding amino acids, which show the same extent of conservation as the whole protein. The pattern of amino acid exchanges differs between cis and trans prolyl residues. Also, the cis prolyl bond is maintained in proteins with sequence identity as low as 20%. This finding emphasizes the importance of cis peptide bonds in protein structure and function.

Ab initio conformational study of N-acetyl-L-proline-N',N'-dimethylamide: a model for polyproline

We report here the results on N-acetyl-l-proline-N',N'-dimethylamide (Ac-Pro-NMe2) as a model for polyproline at the HF/6-31+G(d) level with the conductor-like polarizable continuum model of self-consistent reaction field methods to figure out the conformational preference and cis-trans isomerization of polyproline in the gas phase, chloroform, methanol, and water. The second methyl substitution at the carboxyl amide end results in different backbone structures and their populations from those of N-acetyl-L-proline-N-methylamide (Ac-Pro-NHMe). In particular, all conformations with the C7 hydrogen bond between acetyl and amide ends, which is the most probable conformations of Ac-Pro-NHMe in the gas phase and in nonpolar solvents, disappeared for Ac-Pro-NMe2 even in the gas phase due to the lack of amide hydrogen. The dominant conformation for Ac-Pro-NMe2 is the polyproline II structure with the trans prolyl peptide bond in the gas phase and in solutions. In methanol, the population of the polyproline I structure with the cis prolyl peptide bond is calculated to be larger than that in water, which is consistent with experiments. It should be noted that Ac-Pro-NMe2 has higher rotational barriers for the cis-trans isomerization of the Ac-Pro peptide bond than Ac-Pro-NHMe in the gas phase and in solutions, which could be due to the lack of the intramolecular hydrogen bond between prolyl nitrogen and carboxyl N-H group for the transition state of Ac-Pro-NMe2. The rotational barriers for Ac-Pro-NMe2 are increased with the increase of solvent polarity, as seen for Ac-Pro-NHMe.

Main-chain conformational tendencies of amino acids

A Ramachandran plot is a visual representation of the main-chain conformational tendencies of an amino acid. Despite forty years of research, the shape of Ramachandran plots is still a matter of debate. The issue in making a Ramachandran plot based on experimental data is deciding whether sparse data represent genuine conformations. We present here a simple solution to settle the ambiguities of the sparse data, and explain how we verified the accuracies of our plots using an independent dataset. To obtain our results, we then measured the pair-wise distances of main-chain conformational tendencies among amino acids, and showed that the conformational relationships of amino acids are well preserved in a two-dimensional map, leading to the conclusion that the conformational diversity space of amino acids is largely two dimensional. We further noticed that amino acids in early and late evolutionary stages are located in different zones in the two-dimensional map. In addition to these conclusions, we here present an amino acid substitution table derived from experimental data.

Prediction of cis/trans isomerization in proteins using PSI-BLAST profiles and secondary structure information

Background

The majority of peptide bonds in proteins are found to occur in the trans conformation. However, for proline residues, a considerable fraction of Prolyl peptide bonds adopt the cis form. Proline cis/trans isomerization is known to play a critical role in protein folding, splicing, cell signaling and transmembrane active transport. Accurate prediction of proline cis/trans isomerization in proteins would have many important applications towards the understanding of protein structure and function.

Results

In this paper, we propose a new approach to predict the proline cis/trans isomerization in proteins using support vector machine (SVM). The preliminary results indicated that using Radial Basis Function (RBF) kernels could lead to better prediction performance than that of polynomial and linear kernel functions. We used single sequence information of different local window sizes, amino acid compositions of different local sequences, multiple sequence alignment obtained from PSI-BLAST and the secondary structure information predicted by PSIPRED. We explored these different sequence encoding schemes in order to investigate their effects on the prediction performance. The training and testing of this approach was performed on a newly enlarged dataset of 2424 non-homologous proteins determined by X-Ray diffraction method using 5-fold cross-validation. Selecting the window size 11 provided the best performance for determining the proline cis/trans isomerization based on the single amino acid sequence. It was found that using multiple sequence alignments in the form of PSI-BLAST profiles could significantly improve the prediction performance, the prediction accuracy increased from 62.8% with single sequence to 69.8% and Matthews Correlation Coefficient (MCC) improved from 0.26 with single local sequence to 0.40. Furthermore, if coupled with the predicted secondary structure information by PSIPRED, our method yielded a prediction accuracy of 71.5% and MCC of 0.43, 9% and 0.17 higher than the accuracy achieved based on the singe sequence information, respectively.

Conclusion

A new method has been developed to predict the proline cis/trans isomerization in proteins based on support vector machine, which used the single amino acid sequence with different local window sizes, the amino acid compositions of local sequence flanking centered proline residues, the position-specific scoring matrices (PSSMs) extracted by PSI-BLAST and the predicted secondary structures generated by PSIPRED. The successful application of SVM approach in this study reinforced that SVM is a powerful tool in predicting proline cis/trans isomerization in proteins and biological sequence analysis.

Probing the competition between secondary structures and local preferences in gas phase isolated peptide backbones

Combining laser desorption with a supersonic expansion together with the selectivity of IR/UV double resonance spectroscopy makes it possible to isolate and characterise the gas phase of remarkable backbone conformations of short peptide chains mimicking protein segments. A systematic bottom-up approach involving a conformer-specific IR study of peptide sequences of increasing sizes has enabled us to map the spectral signatures of the intramolecular interactions, which shape the peptide backbone, in particular H-bonds. The precise data collected are directly comparable to the most sophisticated quantum chemistry calculations of these species and therefore constitute a stringent test for the theoretical methods used. One-residue chains reveal the local conformational preference of the backbone and its dependence upon the nature of the residue. The investigation of longer chains provides evidence for a competition between simple successions of local conformational preferences along the chain and more folded structures, in which a new H-bonding network, involving distant H-bonding sites along the backbone, takes place. From three residues, the issue of helical folding can also be addressed. The present review of the gas phase literature data emphasizes the observation of remarkable secondary structures of biology, including short segments of beta-strands, gamma- and beta-turns, combinations of turns, including a 3(10) helix. It also provides evidence for the flexibility of the peptide chains, i.e., a critical influence of rather minor interactions (like side-chain/backbone interactions) on the conformational stability. Finally, the paper will discuss future promising directions of the present approach.

The crystal structure of the cis-proline to glycine variant (P114G) of ribonuclease A

Replacement of a cis-proline by glycine at position 114 in ribonuclease A leads to a large decrease in thermal stability and simplifies the refolding kinetics. A crystallographic approach was used to determine whether the decrease in thermal stability results from the presence of a cis glycine peptide bond, or from a localized structural rearrangement caused by the isomerization of the mutated cis 114 peptide bond. The structure was solved at 2.0 A resolution and refined to an R-factor of 19.5% and an R(free) of 21.9%. The overall conformation of the protein was similar to that of wild-type ribonuclease A; however, there was a large localized rearrangement of the mutated loop (residues 110-117-a 9.3 A shift of the Calpha atom of residue 114). The peptide bond before Gly114 is in the trans configuration. Interestingly, a large anomalous difference density was found near residue 114, and was attributed to a bound cesium ion present in the crystallization experiment. The trans isomeric configuration of the peptide bond in the folded state of this mutant is consistent with the refolding kinetics previously reported, and the associated protein conformational change provides an explanation for the decreased thermal stability.

All-cis Cyclic Peptides

Amide bonds -NH-CO- preferentially exist in trans conformations, the cis conformation being thermodynamically unfavored with respect to the trans by about 2 kcal/mol. Yet, the main reason most proteins or peptides cannot be made from cis-peptide plaques only lies in that connecting them into open chains appears to be sterically impracticable. It is possible, however, to build all-cis cyclic peptides in which all cis-plaques are efficiently locked. The present work examines, through quantum calculations, the structural and energetic issues associated with these peculiar arrangements. Systematic exploration at DFT-B3LYP level of the potential-energy surfaces for all-cis cyclopolyglycines cG(n)(c) (n = 2-10,15), and to a lesser extent, for all-cis cyclopolyalanines and all-cis cyclopolyphenylalanines confirms that all these structures are true minima. Optimal ring size occurs around eight peptide units, resulting in planar cG7(c), cG8(c), and cG9(c). In smaller systems, the ring strain is relieved through nonplanar cup-like distortions, particularly in cG6(c). From 10 peptide units and beyond, the ring framework distorts into a saddle-edge shape. These molecules disclose some molecular flexibility, as combinatorial tilting of the plaques may give sets of minima close in energy. Indexes based on isodesmic reactions are used to estimate the energy for joining all-cis or all-trans plaques into cyclic peptides. One of them, the mean plaque-junction energy (MPJE) suggests that within sensible sizes from six peptide units and beyond, all-cis plaque association is almost equally favorable as all-trans one. The frame of radiating cis-amide bonds can be considered as defining a new kind of peptidic material, endowed with specific self-assembling properties.

Puckering Transition of 4-Substituted Proline Residues

The puckering transition of 4-substituted proline residues by electron-withdrawing groups, i.e., 4(R)-hydroxy-L-proline (Hyp) and 4(R)-fluoro-L-proline (Flp) residues, with trans and cis prolyl peptide bonds was studied by adiabatic optimizations along the torsion angle chi1 of the prolyl ring at the HF/6-31+G(d) level. By analyzing the potential energy surface and local minima, it is observed that the puckering transition of the prolyl ring for Hyp and Flp residues proceeds from a down-puckered conformation to an up-puckered one through the transition state with an envelope form having the N atom at the top of envelope and not a planar one for both trans and cis conformers, which is the same as found for the unsubstituted proline residue. At HF/6-31+G(d) and B3LYP/6-311++G(d,p) levels, the structures of the backbone and prolyl ring for local minima of Ac-Hyp-NHMe and Ac-Flp-NHMe are quite similar to those of Ac-Pro-NHMe. However, the relative stability of the up-puckered conformation to the down-puckered one is increased for Ac-Hyp-NHMe with the cis imide bond and for Ac-Flp-NHMe with the trans and cis imide bonds. In particular, the 4(R)-substitution by hydroxy and fluorine groups has brought some structural changes in the prolyl ring of the transition states and the changes in barriers for the puckering transition. The puckering transitions for Ac-Hyp-NHMe and Ac-Flp-NHMe are proven to be predominantly electronically driven by analyzing the electronic and enthalpic contributions to the barriers, as seen for Ac-Pro-NHMe.

The flexibility in the proline ring couples to the protein backbone

In proteins, the proline ring exists predominantly in two discrete states. However, there is also a small but significant amount of flexibility in the proline ring of high-resolution protein structures. We have found that this side-chain flexibility is coupled to the backbone conformation. To study this coupling, we have developed a model that is simply based on geometric and steric factors and not on energetics. We show that the coupling between phi and chi1 torsions in the proline ring can be described by an analytic equation that was developed by Bricard in 1897, and we describe a computer algorithm that implements the equation. The model predicts the observed coupling very well. The strain in the C(gamma)-C(delta)-N angle appears to be the principal barrier between the UP and DOWN pucker. This strain is relaxed to allow the proline ring to flatten in the rare PLANAR conformation.

Aromatic–aromatic and proline–aromatic interactions in endomorphin-1 and endomorphin-2

We investigated the aromatic-aromatic and proline-aromatic interactions in endomorphin-1 and endomorphin-2, and different types of these interactions were observed. For all the interacting pairs, the preferred geometric orientations were identified. We examined these interactions in the preferred secondary structures, which are different types of beta-turns and gamma-turns. These results revealed that the majority of the turn structures contained one of the interacting aromatic-aromatic or proline-aromatic pairs. On the basis of our results, it can be assumed that these interactions may be important in the determination and stabilization of the structures of endomorphins. Furthermore, our observations suggest that a conformation stabilized by an aromatic-aromatic or proline-aromatic interaction can play an important role in the association with the receptor.

Proline Editing: A Divergent Strategy for the Synthesis of Conformationally Diverse Peptides

[reaction: see text] Strong conformational biases in peptides and proteins can be achieved with 4-substituted proline residues (cis-, trans-, or disubstituted fluoroproline or hydroxyproline). The practical, divergent synthesis of peptides containing these residues, via postsynthetic modification of a peptide containing an internal trans-hydroxyproline residue, is described. Significant differences in the conformations of the peptides Ac-TYXN-NH2 were observed, including K(trans/cis) values, which varied from 1.5 (X = cis-fluoroproline) to 7.0 (X = trans-fluoroproline).

Statistically significant dependence of the Xaa-Pro peptide bond conformation on secondary structure and amino acid sequence

Background

A reliable prediction of the Xaa-Pro peptide bond conformation would be a useful tool for many protein structure calculation methods. We have analyzed the Protein Data Bank and show that the combined use of sequential and structural information has a predictive value for the assessment of the cis versus trans peptide bond conformation of Xaa-Pro within proteins. For the analysis of the data sets different statistical methods such as the calculation of the Chou-Fasman parameters and occurrence matrices were used. Furthermore we analyzed the relationship between the relative solvent accessibility and the relative occurrence of prolines in the cis and in the trans conformation.

Results

One of the main results of the statistical investigations is the ranking of the secondary structure and sequence information with respect to the prediction of the Xaa-Pro peptide bond conformation. We observed a significant impact of secondary structure information on the occurrence of the Xaa-Pro peptide bond conformation, while the sequence information of amino acids neighboring proline is of little predictive value for the conformation of this bond.

Conclusion

In this work, we present an extensive analysis of the occurrence of the cis and trans proline conformation in proteins. Based on the data set, we derived patterns and rules for a possible prediction of the proline conformation. Upon adoption of the Chou-Fasman parameters, we are able to derive statistically relevant correlations between the secondary structure of amino acid fragments and the Xaa-Pro peptide bond conformation.

First Principle Computational Study on the Full Conformational Space of l-Proline Diamides

Ab initio molecular orbital computations were carried out at three levels of theory: RHF/3-21G, RHF/6-31G(d), and B3LYP/6-31G(d), on four model systems of the amino acid proline, HCO-Pro-NH2 [I], HCO-Pro-NH-Me [II], MeCO-Pro-NH2 [III], and MeCO-Pro-NH-Me [IV], representing a systematic variation in the protecting N- and C-terminal groups. Three previously located backbone conformations, gammaL, epsilonL, and alphaL, were characterized together with two ring-puckered forms syn (gauche+ = g+) or "DOWN" and anti (gauche- = g-) or "UP", as well as trans-trans, trans-cis, cis-trans, and cis-cis peptide bond isomers. The topologies of the conformational potential energy cross-sections (PECS) of the potential energy hypersurfaces (PEHS) for compounds [I]-[IV] were explored and analyzed in terms of potential energy curves (PEC), and HCO-Pro-NH2 [I] was also analyzed in terms of potential energy surfaces (PESs). Thermodynamic functions were also calculated for HCO-Pro-NH2 [I] at the CBS-4M and G3MP2 levels of theory. The study confirms that the use of the simplest model, compound [I] with P(N) = P(C) = H, along with the RHF/3-21G level of theory, is an acceptable practice for the analysis of peptide models because only minor differences in geometry and stability are observed.

A Plasmodium falciparum FK506-binding protein (FKBP) with peptidyl–prolyl cis–trans isomerase and chaperone activities

The immunosuppressive drugs FK506 and rapamycin have anti-malarial properties but their mechanisms of action against malaria parasites remain unknown. The pathway by which these drugs cause immunosuppression in humans is known to involve an FK506-binding protein (FKBP). Homologues of FKBPs have been identified in almost every organism in which they have been sought. Here, we describe the characterisation of the first member of the FKBP family identified in the human malarial parasite, Plasmodium falciparum. This 35-kDa protein, PfFKBP35, comprises a single, N-terminal, FKBP domain and a C-terminal tripartite tetratricopeptide repeat domain. A recombinant form of PfFKBP35, like most other FKBPs, displayed peptidyl-prolyl cis-trans isomerase activity that was inhibitable by FK506 and rapamycin. Unusually the phosphatase activity of calcineurin, the target of the FK506-FKBP complex in T-lymphocytes, was inhibited by PfFKBP35 independently of FK506 binding. PfFKBP35 also inhibited the thermal aggregation in vitro of two model substrates, suggesting that it has general chaperone properties. Analysis of the P. falciparum genome database suggested this to be the only FKBP present in the parasite. The function of this protein remains unknown but the presence of tetratricopeptide repeat motifs suggests a role in intracellular protein transport or modulation of protein function.

Intrinsic Folding of Small Peptide Chains: Spectroscopic Evidence for the Formation of β-Turns in the Gas Phase

Laser desorption of model peptides coupled to laser spectroscopic techniques enables the gas-phase observation of genuine secondary structures of biology. Spectroscopic evidence for the formation of beta-turns in gas-phase peptide chains containing glycine and phenylalanine residues establishes the intrinsic stability of these forms and their ability to compete with other stable structures. The precise characterization of local minima on the potential energy surface from IR spectroscopy constitutes an acute assessment for the state-of-the-art quantum mechanical calculations also presented. The observation of different types of beta-turns depending upon the residue order within the sequence is found to be consistent with the residue propensities in beta-turns of proteins, which suggests that the prevalence of glycine in type II and II' turns stems essentially from an energetic origin, already at play under isolated conditions.

CH/π hydrogen bonds as evidenced in the substrate specificity of acetylcholine esterase

The crystal structure of acetylcholine esterase (AchE) in complex with various inhibitors, investigated as drugs for improvement of the cognitive ability of early stage Alzheimer's disease, has been analyzed with the use of our program CHPI. A number of CH/pi hydrogen bonds have been disclosed in the binding of the inhibitors with Torpedo californica AchE. It has been demonstrated that, in order to be effective in the binding with AchE, C-H bonds in the inhibitor need not be polarized.

Context-dependent effects of proline residues on the stability and folding pathway of ubiquitin

Substitution of trans-proline at three positions in ubiquitin (residues 19, 37 and 38) produces significant context-dependent effects on protein stability (both stabilizing and destabilizing) that reflect changes to a combination of parameters including backbone flexibility, hydrophobic interactions, solvent accessibility to polar groups and intrinsic backbone conformational preferences. Kinetic analysis of the wild-type yeast protein reveals a predominant fast-folding phase which conforms to an apparent two-state folding model. Temperature-dependent studies of the refolding rate reveal thermodynamic details of the nature of the transition state for folding consistent with hydrophobic collapse providing the overall driving force. Brønsted analysis of the refolding and unfolding rates of a family of mutants with a variety of side chain substitutions for P37 and P38 reveals that the two prolines, which are located in a surface loop adjacent to the C terminus of the main alpha-helix (residues 24-33), are not significantly structured in the transition state for folding and appear to be consolidated into the native structure only late in the folding process. We draw a similar conclusion regarding position 19 in the loop connecting the N-terminal beta-hairpin to the main alpha-helix. The proline residues of ubiquitin are passive spectators in the folding process, but influence protein stability in a variety of ways.

A novel method of analyzing proline synonymous codons in E. coli

Proline is a special imino acid in protein and the isomerization of the prolyl peptide bond has notable biological significance and influences the final structure of protein greatly, so the correlation between proline synonymous codon usage and local amino acid, the correlation between proline synonymous codon usage and the isomerization of the prolyl peptide bond were both investigated in the Escherichia coli genome by using a novel method based on information theory. The results show that in peptide chain, the residue at the first position C-terminal influences the usage of proline synonymous codon greatly and proline synonymous codons contain some factors influencing the isomerization of the prolyl peptide bond.

Structural mechanism governing cis and trans isomeric states and an intramolecular switch for cis/trans isomerization of a non-proline peptide bond observed in crystal structures of scorpion toxins

Non-proline cis peptide bonds have been observed in numerous protein crystal structures even though the energetic barrier to this conformation is significant and no non-prolyl-cis/trans-isomerase has been identified to date. While some external factors, such as metal binding or co-factor interaction, have been identified that appear to induce cis/trans isomerization of non-proline peptide bonds, the intrinsic structural basis for their existence and the mechanism governing cis/trans isomerization in proteins remains poorly understood. Here, we report the crystal structure of a newly isolated neurotoxin, the scorpion alpha-like toxin Buthus martensii Karsch (BmK) M7, at 1.4A resolution. BmK M7 crystallizes as a dimer in which the identical non-proline peptide bond between residues 9 and 10 exists either in the cis conformation or as a mixture of cis and trans conformations in either monomer. We also determined the crystal structures of several mutants of BmK M1, a representative scorpion alpha-like toxin that contains an identical non-proline cis peptide bond as that observed in BmK M7, in which residues within or neighboring the cis peptide bond were altered. Substitution of an aspartic acid residue for lysine at residue 8 in the BmK M1 (K8D) mutant converted the cis form of the non-proline peptide bond 9-10 into the trans form, revealing an intramolecular switch for cis-to-trans isomerization. Cis/trans interconversion of the switch residue at position 8 appears to be sequence-dependent as the peptide bond between residues 9 and 10 retains its wild-type cis conformation in the BmK M1 (K8Q) mutant structure. The structural interconversion of the isomeric states of the BmK M1 non-proline cis peptide bond may relate to the conversion of the scorpion alpha-toxins subgroups.

Cis–trans isomerization and puckering of proline residue

We report here the results on N-acetyl-L-proline-N'-methylamide (Ac-Pro-NHMe) calculated at the HF/6-31+G(d) level with the conductor-like polarizable continuum model (CPCM) of self-consistent reaction field methods to investigate the changes of backbone and prolyl ring along the cis-trans isomerization of the prolyl peptide bond. From the potential energy surface, the barrier to ring flip from the down-puckered conformation to the up-puckered one is estimated to be 2.5 and 3.2 kcal/mol for trans and cis conformers of Ac-Pro-NHMe, respectively. In particular, the ring flip seems to be inaccessible in the intermediate regions between trans and cis conformations, because of higher barriers (approximately 13-19 kcal/mol) to rotation of the prolyl peptide bond. The torsion angles for backbone and prolyl ring vary largely around the transition states at omega' approximately 120 degrees and -70 degrees for the prolyl peptide bond. Three kinds of puckering amplitudes show the same trend of puckering along the cis-trans isomerization although their absolute values are different. In particular, trans and cis conformations have the almost same degree of puckering. The cis populations and barriers to rotation of the prolyl peptide bond for Ac-Pro-NHMe are increased with the increase of solvent polarity, which is mainly ascribed to the decreases of relative free energies for cis conformations and the increase of relative free energies for transition states.

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.

NMR Solution Structures of δ-Conotoxin EVIA from Conus ermineus That Selectively Acts on Vertebrate Neuronal Na+ Channels

Delta-conotoxin EVIA, from Conus ermineus, is a 32-residue polypeptide cross-linked by three disulfide bonds forming a four-loop framework. delta-Conotoxin EVIA is the first conotoxin known to inhibit sodium channel inactivation in neuronal membranes from amphibians and mammals (subtypes rNa(v)1.2a, rNa(v)1.3, and rNa(v)1.6), without affecting rat skeletal muscle (subtype rNa(v)1.4) and human cardiac muscle (subtype hNa(v)1.5) sodium channel (Barbier, J., Lamthanh, H., Le Gall, F., Favreau, P., Benoit, E., Chen, H., Gilles, N., Ilan, N., Heinemann, S. F., Gordon, D., Ménez, A., and Molgó, J. (2004) J. Biol. Chem. 279, 4680-4685). Its structure was solved by NMR and is characterized by a 1:1 cis/trans isomerism of the Leu(12)-Pro(13) peptide bond in slow exchange on the NMR time scale. The structure of both cis and trans isomers could be calculated separately. The isomerism occurs within a specific long disordered loop 2, including residues 11-19. These contribute to an important hydrophobic patch on the surface of the toxin. The rest of the structure matches the "inhibitor cystine-knot motif" of conotoxins from the "O superfamily" with a high structural order. To probe a possible functional role of the Leu(12)-Pro(13) cis/trans isomerism, a Pro(13) --> Ala delta-conotoxin EVIA was synthesized and shown to exist only as a trans isomer. P13A delta-conotoxin EVIA was estimated only two times less active than the wild-type EVIA in binding competition to rat brain synaptosomes and when injected intracerebroventricularly into mice.

A case study of 2,2-dimethylthiazolidine as locked cis proline amide bond: synthesis, NMR and molecular modeling studies of a δ-conotoxin EVIA peptide analog

The delta-conotoxin EVIA from the Conus ermineus venom, a recently characterized toxin, exhibits cis-trans isomerism of the Leu12-Pro13 bond associated with the triggering of its biological activity. In this paper we use the pseudoproline concept to target the presumed bioactive cis conformation. We report the design and the synthesis of loop 2 analogs from residue 8 to 18 containing either the cis-inducing Cys(PsiMe,MePro)13 unit or the natural proline residue. NMR studies in water and molecular modeling allowed us to identify the amide bond "locked" in a cis conformation for as in the suggested bioactive form of the natural toxin.

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.

Role of proline residues in the folding of serine hydroxymethyltransferase

Previous studies on the folding mechanism of Escherichia coli serine hydroxymethyltransferase (SHMT) showed that the final rate determining folding step was from an intermediate that contained two fully folded domains with N-terminal segments of approximately 55 residues and interdomain segments of approximately 50 residues that were still solvent exposed and subject to proteolysis. The interdomain segment contains 3 Pro residues near its N terminus and 2 Pro residues near its C terminus. The 5 Pro residues were each mutated to both a Gly and Ala residue, and each mutant SHMT was purified and characterized with respect to kinetic properties, stability, secondary structure, and folding mechanism. The results showed that Pro214 and Pro218 near the N terminus of the interdomain segment are not critical for folding, stability, or activity. The P216A mutant also retained most of the characteristics of the native enzyme, but its folding rate was altered. However, the P216G mutant was severely compromised in folding into a catalytically competent enzyme. Mutation of both Pro258 and Pro264 had altered folding kinetics and resulted in enzymes that expressed little catalytic activity. The Phe257-Pro258 bond is cis in its configuration, and the P258A mutant SHMT showed reduced thermal stability. Pro216, Pro258, and Pro264 are conserved in all 53 known sequences of this enzyme. The results are discussed in terms of the role of each Pro residue in maintaining the structure and function of SHMT and a possible role in pyridoxal 5'-phosphate addition to the apo-enzyme.

Deciphering a novel thioredoxin-like fold family

Sequence--and structure-based searching strategies have proven useful in the identification of remote homologs and have facilitated both structural and functional predictions of many uncharacterized protein families. We implement these strategies to predict the structure of and to classify a previously uncharacterized cluster of orthologs (COG3019) in the thioredoxin-like fold superfamily. The results of each searching method indicate that thioltransferases are the closest structural family to COG3019. We substantiate this conclusion using the ab initio structure prediction method rosetta, which generates a thioredoxin-like fold similar to that of the glutaredoxin-like thioltransferase (NrdH) for a COG3019 target sequence. This structural model contains the thiol-redox functional motif CYS-X-X-CYS in close proximity to other absolutely conserved COG3019 residues, defining a novel thioredoxin-like active site that potentially binds metal ions. Finally, the rosetta-derived model structure assists us in assembling a global multiple-sequence alignment of COG3019 with two other thioredoxin-like fold families, the thioltransferases and the bacterial arsenate reductases (ArsC).

The nature of the turn in omega loops of proteins

An analysis of Omega loops in a nonredundant set of protein structures from the Protein Data Bank has been carried out to determine the nature of the "turn elements" present. Because Omega loops essentially reverse their direction in three-dimensional space, this analysis was made with respect to four turn elements identified as (1) Gly; (2) Pro; (3) a residue with alpha-helical phi,psi angles, termed a helical residue; and (4) a cis peptide. A set of 1079 Omega loops from a set of 680 proteins were used for the analysis. Apart from other criteria that define Omega loops, the selection of an Omega loop from a cluster of loops is based on an exposure index. In this study, analyses have been made with two sets of data: (1) Omega loops arising from a minimum exposure index indicative of a less exposed loop (xmin set) and (2) Omega loops with a maximum exposure index indicative of a relatively exposed loop (xmax set). Overall residue preferences and positional preferences have been examined. Positions of the turn elements for Omega loops of varying length have also been studied. Specific positional preferences are observed for particular turn elements with regard to the length of Omega loops. Analysis in terms of the turn elements can provide guidelines for modeling of loops in proteins. Apart from Pro, which has the natural tendency to form cis peptide bonds, a higher occurrence of non-Pro cis peptide bonds is observed. Torsion angles in Omega loops also indicate the occurrence of a large number of residues with helical phi,psi angles, necessary for the turn in the loop structures.

Parallel channels and rate-limiting steps in complex protein folding reactions: prolyl isomerization and the alpha subunit of Trp synthase, a TIM barrel protein

A kinetic folding mechanism for the alpha subunit of tryptophan synthase (alphaTS) from Escherichia coli, involving four parallel channels with multiple native, intermediate and unfolded forms, has recently been proposed. The hypothesis that cis/trans isomerization of several Xaa-Pro peptide bonds is the source of the multiple folding channels was tested by measuring the sensitivity of the three rate-limiting phases (tau(1), tau(2), tau(3)) to catalysis by cyclophilin, a peptidyl-prolyl isomerase. Although the absence of catalysis for the tau(1) (fast) phase leaves its assignment ambiguous, our previous mutational analysis demonstrated its connection to the unique cis peptide bond preceding proline 28. The acceleration of the tau(2) (medium) and tau(3) (slow) refolding phases by cyclophilin demonstrated that cis/trans prolyl isomerization is also the source of these phases. A collection of proline mutants, which covered all of the remaining 18 trans proline residues of alphaTS, was constructed to obtain specific assignments for these phases. Almost all of the mutant proteins retained the complex equilibrium and kinetic folding properties of wild-type alphaTS; only the P217A, P217G and P261A mutations caused significant changes in the equilibrium free energy surface. Both the P78A and P96A mutations selectively eliminated the tau(1) folding phase, while the P217M and P261A mutations eliminated the tau(2) and tau(3) folding phases, respectively. The redundant assignment of the tau(1) phase to Pro28, Pro78 and Pro96 may reflect their mutual interactions in non-random structure in the unfolded state. The non-native cis isomers for Pro217 and Pro261 may destabilize an autonomous C-terminal folding unit, thereby giving rise to kinetically distinct unfolded forms. The nature of the preceding amino acid, the solvent exposure, or the participation in specific elements of secondary structure in the native state, in general, are not determinative of the proline residues whose isomerization reactions can limit folding.

A cis-prolyl peptide bond isomerization dominates the folding of the alpha subunit of Trp synthase, a TIM barrel protein

The cis/trans isomerization of prolyl peptide bonds has been suggested to dominate the folding of the alpha subunit of tryptophan synthase from Escherichia coli (alphaTS). To test the role of the unique cis isomer between Asp27 and Pro28, the folding properties of P28A, P28G and G(3)P28G, a three-glycine insertion mutant between Asp27 and Gly28, were investigated using urea as a denaturant. Circular dichroism analysis demonstrated that none of the mutations perturb the secondary structure significantly, although the aromatic side-chain packing is altered for P28A and P28G. All three mutant proteins inherited the three-state thermodynamic behavior observed in wild-type alphaTS, ensuring that the fundamental features of the energy surface are intact. Kinetic studies showed that neither alanine nor glycine substitutions at Pro28 results in the elimination of any slow-refolding phases. By contrast, the G(3)P28G mutant eliminates the fastest of the slow-refolding phases and one of the two unfolding phases. Double-jump experiments on G(3)P28G confirm the assignment of the missing refolding phase to the isomerization of the Asp27-Pro28 peptide bond. These results imply that the local stability conveyed by the tight, overlapping turns containing the cis peptide bond is sufficient to favor the cis isomer for several non-prolyl residues. The free energy required to drive the isomerization reaction is provided by the formation of the stable intermediate, demonstrating that the acquisition of structure and stability is required to induce subsequent rate-limiting steps in the folding of alphaTS.

Crystal structure of the jacalin-T-antigen complex and a comparative study of lectin-T-antigen complexes

Thomsen-Friedenreich antigen (Galbeta1-3GalNAc), generally known as T-antigen, is expressed in more than 85% of human carcinomas. Therefore, proteins which specifically bind T-antigen have potential diagnostic value. Jacalin, a lectin from jack fruit (Artocarpus integrifolia) seeds, is a tetramer of molecular mass 66kDa. It is one of the very few proteins which are known to bind T-antigen. The crystal structure of the jacalin-T-antigen complex has been determined at 1.62A resolution. The interactions of the disaccharide at the binding site are predominantly through the GalNAc moiety, with Gal interacting only through water molecules. They include a hydrogen bond between the anomeric oxygen of GalNAc and the pi electrons of an aromatic side-chain. Several intermolecular interactions involving the bound carbohydrate contribute to the stability of the crystal structure. The present structure, along with that of the Me-alpha-Gal complex, provides a reasonable qualitative explanation for the known affinities of jacalin to different carbohydrate ligands and a plausible model of the binding of the lectin to T-antigen O-linked to seryl or threonyl residues. Including the present one, the structures of five lectin-T-antigen complexes are available. GalNAc occupies the primary binding site in three of them, while Gal occupies the site in two. The choice appears to be related to the ability of the lectin to bind sialylated sugars. In either case, most of the lectin-disaccharide interactions are at the primary binding site. The conformation of T-antigen in the five complexes is nearly the same.

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.

Leishmania donovani phosphofructokinase. Gene characterization, biochemical properties and structure-modeling studies

The characterization of the gene encoding Leishmania donovani phosphofructokinase (PFK) and the biochemical properties of the expressed enzyme are reported. L. donovani has a single PFK gene copy per haploid genome that encodes a polypeptide with a deduced molecular mass of 53 988 and a pI of 9.26. The predicted amino acid sequence contains a C-terminal tripeptide that conforms to an established signal for glycosome targeting. L. donovani PFK showed most sequence similarity to inorganic pyrophosphate (PPi)-dependent PFKs, despite being ATP-dependent. It thereby resembles PFKs from other Kinetoplastida such as Trypanosoma brucei, Trypanoplasma borreli (characterized in this study), and a PFK found in Entamoeba histolytica. It exhibited hyperbolic kinetics with respect to ATP whereas the binding of the other substrate, fructose 6-phosphate, showed slight positive cooperativity. PPi, even at high concentrations, did not have any effect. AMP acted as an activator of PFK, shifting its kinetics for fructose 6-phosphate from slightly sigmoid to hyperbolic, and increasing considerably the affinity for this substrate, whereas GDP did not have any effect. Modelling studies and site-directed mutagenesis were employed to shed light on the structural basis for the AMP effector specificity and on ATP/PPi specificity among PFKs.

Crystal structure of the autocatalytic initiator of glycogen biosynthesis, glycogenin

Glycogen is an important storage reserve of glucose present in many organisms, from bacteria to humans. Its biosynthesis is initiated by a specialized protein, glycogenin, which has the unusual property of transferring glucose from UDP-glucose to form an oligosaccharide covalently attached to itself at Tyr194. Glycogen synthase and the branching enzyme complete the synthesis of the polysaccharide. The structure of glycogenin was solved in two different crystal forms. Tetragonal crystals contained a pentamer of dimers in the asymmetric unit arranged in an improper non-crystallographic 10-fold relationship, and orthorhombic crystals contained a monomer in the asymmetric unit that is arranged about a 2-fold crystallographic axis to form a dimer. The structure was first solved to 3.4 A using the tetragonal crystal form and a three-wavelength Se-Met multi-wavelength anomalous diffraction (MAD) experiment. Subsequently, an apo-enzyme structure and a complex between glycogenin and UDP-glucose/Mn2+ were solved by molecular replacement to 1.9 A using the orthorhombic crystal form. Glycogenin contains a conserved DxD motif and an N-terminal beta-alpha-beta Rossmann-like fold that are common to the nucleotide-binding domains of most glycosyltransferases. Although sequence identity amongst glycosyltransferases is minimal, the overall folds are similar. In all of these enzymes, the DxD motif is essential for coordination of the catalytic divalent cation, most commonly Mn2+. We propose a mechanism in which the Mn2+ that associates with the UDP-glucose molecule functions as a Lewis acid to stabilize the leaving group UDP and to facilitate the transfer of the glucose moiety to an intermediate nucleophilic acceptor in the enzyme active site, most likely Asp162. Following transient transfer to Asp162, the glucose moiety is then delivered to the final acceptor, either directly to Tyr194 or to glucose residues already attached to Tyr194. The positioning of the bound UDP-glucose far from Tyr194 in the glycogenin structure raises questions as to the mechanism for the attachment of the first glucose residues. Possibly the initial glucosylation is via inter-dimeric catalysis with an intra-molecular mechanism employed later in oligosaccharide synthesis.

Local structural changes caused by peptidyl-prolyl cis/trans isomerization in the native state of proteins

Peptidyl-prolyl cis/trans isomerization, observed in the native state of an increasing number of proteins, is of considerable biological significance. The first evidence for an asymmetric transmission along the polypeptide chain of the structural effects of prolyl isomerization is now derived from the statistics of the C(alpha)/C(alpha)-atom distance distributions in the crystal structures of 848 non-homologous proteins. More detailed information on how isomerization affects segments adjacent to proline is obtained from crystal structures of proteins, that are more than 95% homologous, and that exhibit two different states of isomerization at a particular prolyl bond. The resulting 64 cases, which represent 3.8% of the database used, form pairs of coordinates which were analyzed for the existence of isomer-specific intramolecular nonbonded C(alpha)/C(alpha)-atom distances around the critical proline, and for the positional preferences for particular amino acids in the isomeric sequence segment. The probability that a native protein exhibits both prolyl isomers in the crystalline state increases in particular with a Pro at the third position N-terminal to the isomeric bond (-3 position), and with Ser, Gly and Asp at the position preceding the isomeric bond (-1 position). Structural alignment of matched pairs of isomeric proteins generates three classes with respect to position-specific distribution of C(alpha)-atom displacements around an isomeric proline imide bond. In the majority of cases the distribution of these intermolecular isomer-specific C(alpha)-atom distances shows a symmetric behavior for the N-terminal and C-terminal segment flanking the proline residue, and the magnitude did not exceed 1.3+/-0.6 A including the C(alpha) atoms in proximity to the prolyl bond. However, in the remaining 12 protein pairs the structural changes are unidirectional relative to the isomerizing bond whereby the magnitude of the isomer-specific effect exceeds 3.0+/-2.0 A even at positions remote to proline. Interestingly, the magnitude of the intramolecular isomer-specific C(alpha) atom displacements reveals a lever-arm amplification of the isomerization-mediated structural changes in a protein backbone. The observed backbone effects provide a structural basis for isomer-specific reactions of proline-containing polypeptides, and thus may play a role in biological recognition and regulation.