|
|
||||||||
T.C. Jenkins Department of Biophysics, Johns Hopkins University, Baltimore, Maryland 21218, USA
Reprint requests to: George D. Rose, Jenkins Department of Biophysics, Johns Hopkins University, 3400 N. Charles Street, Baltimore, MD 21218, USA; e-mail: grose{at}jhu.edu; fax: (410) 516-4118.
(RECEIVED June 2, 2005; FINAL REVISION June 2, 2005; ACCEPTED August 14, 2005)
| Abstract |
|---|
|
|
|---|
-Turns are sites at which proteins change their overall chain direction, and they occur with high frequency in globular proteins. The Protein Data Bank has many instances of conformations that resemble
-turns but lack the characteristic NH(i)
O=C(i 3) hydrogen bond of an authentic
-turn. Here, we identify potential hydrogen-bonded
-turns in the coil library, a Web-accessible database utility comprised of all residues not in repetitive secondary structure, neither
-helix nor
-sheet (http://www.roselab.jhu.edu/coil). In particular, candidate turns were identified as four-residue segments satisfying highly relaxed geometric criteria but lacking a strictly defined hydrogen bond. Such candidates were then subjected to a minimization protocol to determine whether slight changes in torsion angles are sufficient to shift the conformation into reference-quality geometry without deviating significantly from the original structure. This approach of applying constrained minimization to known structures reveals a substantial population of previously unidentified, stringently defined, hydrogen-bonded
-turns. In particular, 33% of coil library residues were classified as
-turns prior to minimization. After minimization, 45% of such residues could be classified as
-turns, with another 8% in 310 helixes (which closely resemble type III
-turns). Of the remaining coil library residues, 37% have backbone dihedral angles in left-handed polyproline II structure. Keywords: protein structure/folding; structure; computational analysis of protein structure; hydrogen bonds; thermodynamics; hydrodynamics
Article and publication are at http://www.proteinscience.org/cgi/doi/10.1110/ps.051625305.
| Introduction |
|---|
|
|
|---|
-Turns (Venkatachalam 1968) represent the largest category of nonrepetitive protein secondary structure (Rose et al. 1985). Several
-turn subcategories have been identified (Lewis et al. 1971; Kuntz 1972; Chou and Fasman 1977; Richardson 1981; Ramakrishnan and Soman 1982; Kabsch and Sander 1983; Wilmot and Thornton 1988, 1990; Efimov 1993); all are four residues in length with an NH(i)
O=C(i 3) hydrogen bond. Upon adopting this motif, the polypeptide chain reverses its overall direction, a frequent occurrence in globular proteins.
Typically, turns are identified from X-ray coordinates based on investigator-defined geometric criteria (Ceccarelli et al. 1981; Jeffrey and Maluszynska 1982; Taylor and Kennard 1983, 1984; Rose et al. 1985; Dasgupta et al. 2004), a highly parameter-sensitive approach near the defined threshold values used to discriminate turns from non-turns (see Table V and Fig. 26 in Rose et al. 1985). In one widely used definition, a four-residue segment with a C
(i)
C
(i 3) distance of 7.0 Å might be classified as a turn while a virtually indistinguishable structure with a corresponding distance of 7.01 Å might not.
We seek to take an accurate census of the hydrogen-bonded turns in proteins by testing whether likely turn conformers can be successfully minimized into tighter geometric criteria without deviating significantly from the original structure. In pursuit of this goal, we take advantage of the recently released coil library (Fitzkee et al. 2005b), a Web-accessible database utility comprised of all residues that do not participate in repetitive secondary structure, neither
-helix nor
-sheet. In a Pisces list of proteins (Wang 2003), this coil library accounts for 44% of total protein structure.
After minimization, we find that 53% of all residues in the coil library participate in
-turns, including 310 helices, which closely resemble type III turns. The largest remaining category is comprised of residues in polyproline II (PII) conformation (Stapley and Creamer 1999). In all, ~80% of protein residues are involved in four hydrogen-bonded backbone structures:
-helices, 310 helices,
-sheet, and
-turns.
| Results |
|---|
|
|
|---|
-turns were identified, comprising 142,069 (33.4%) of the coil library residues (Table 1
-turns that satisfy relaxed criteria and survive the minimization protocol (see Materials and Methods), bringing our total census to 192,241 (45.2%) coil library residues (Table 1
-turns or 310 helices.
|
|
-turns surviving our protocol, differences in backbone conformation upon minimization were small (Table 3
|
|
,
angles [
(
1,
1) = 14.55°, 9.29°;
(
2,
2) = 7.74°, 6.48°;
(
3 ,
3) = 24.88°, 22.06°;
(
4,
4) = 27.35°, 3.04°]. The minimized conformer is a type III turn with an NH(213)
O=C(210) H-bond donor to acceptor distance of 3.45 Å. The RMSD of all backbone atoms between the initial and the minimized four-residue structure is 0.47 Å, with the largest difference resulting from a 1.12 Å shift in the carbonyl oxygen of residue 211.
|
| Discussion |
|---|
|
|
|---|
-turns.
Prior to minimization, we find that 25% of geometrically identified
-turns do not participate in an i
(i 3) hydrogen bond, as assessed by stringent hydrogen-bonding criteria. Following minimization, 87% of these can be reclassified as having a hydrogen bond. The increase in identified turns is accompanied by a corresponding reduction in the remaining fraction of residues that cannot be classified as participating in a hydrogen-bonded back-bone structure (
-helix,
-sheet, 310 helix, and
-turn).
Our study did not include type VIII turns (Wilmot and Thornton 1988, 1990), which are not hydrogen-bonded. When the hydrogen-bonded
-turns identified in our study are supplemented by type VIII turns, only 33% of the coil library (14% of total protein structure) remains unclassified.
There is an extensive literature on hydrogen-bond energy, which is better represented by a potential than a geometric threshold (see, for example, Fig. 2
.1 in Jeffrey and Saenger 1991). Indeed, the energy of interaction between donor and acceptor does not drop to zero beyond some cutoff-value that was selected for computational convenience. Depending on how the figure of merit is defined, slight changes in backbone torsion angles that enhance hydrogen-bond energy may actually result in an improved experimental structure.
After subtracting all hydrogen-bonded
-turns from the coil library, 37% of the remaining residues are within the polyproline II (PII) region of the dipeptide map (Fig. 3
), a conformation thought to optimize backbone exposure to hydrogen bonds with solvent water (Pappu and Rose 2002). This finding underscores the tendency of residues in folded proteins to participate either in hydrogen-bonded backbone structures or in hydrogen bonds with solvent. The results presented here lend further support to the hydrogen-bonding hypothesis, our proposal that essentially all potential hydrogen-bond donors and acceptors in proteins are satisfied a significant fraction of the time, either by intramolecular hydrogen bonds or by hydrogen bonds to solvent water (Fitzkee et al. 2005a; Fleming and Rose 2005).
|
-helix nor
-sheet and are included in the coil library (Fitzkee 2005b). In our analysis, the coil library was filtered initially to count and remove all 310 helices, after which
-turns were identified and tallied using stringent geometric criteria. Remaining residues were then screened using relaxed turn criteria and candidate structures were minimized; successfully minimized candidates were added to the turn census. These steps are now described in greater detail.
Identification of 310 helices
Residues participating in 310 helices were identified prior to assessing turn candidates. A 310 helix was defined as three or more consecutive residues with backbone dihedral angles of
,
= 60 ± 30°, 30 ± 30°.
Identification of hydrogen-bonded turns using stringent criteria
Our turn definition is based on distance and angular constraints taken from well-resolved, small-molecule X-ray structures (Taylor and Kennard 1983, 1984). Specifically, a
-turn was identified as a four-residue segment in which backbone dihedral angles for the middle two residues (i + 1, i + 2) were within 30° of a canonical turn type (Rose et al. 1985), with an NH (i + 3)
O=C(i) hydrogen bond (Stickle et al. 1992). Detailed hydrogen-bonding criteria are shown in Figure 4
. In those rare cases where a candidate turn satisfied both type I and type III definitions, the turn was classified as type III. In the case of a type VI turn, the
angle was required to lie within ±30° of planarity.
|
Counting turns
The two residues preceding and following a candidate structure were used in the identification of
-turns. For two- or three-residue coil segments, these flanking residues will adjoin an
-helix or
-strand external to the coil library. In such cases, only residues within the coil library proper were added to the count. When two adjacent turn segments had residues in common, the overlapping residues were counted only once.
Minimization protocol
Every potential turn candidate identified using relaxed criteria was subjected to a minimization protocol to determine whether it could be "tweaked" into a similar structure satisfying stringent criteria, including hydrogen-bond geometry derived from small molecule X-ray studies, canonical
,
backbone angles, and satisfaction of hydrogen-bonding requirements for all polar atoms.
Up to 500 attempts of the minimization protocol were performed on an all-atom, nine-residue segment that included the four candidate turn residues. This nine-residue segment was excised from the original X-ray crystal structure, beginning with the residue preceding the potential turn (labeled residue i) and continuing through the four residues following the potential turn.
To determine whether a nearby backbone conformation satisfied strictly defined hydrogen-bonded turn criteria, a Monte Carlo search was performed (Srinivasan et al. 2004), with coordinates of residue i held constant and
,
angles of residues (i + 1)(i + 4) allowed to vary at random within ±30° of their experimentally determined X-ray values;
was also allowed to vary by ±5°. Trial conformations were identified as those satisfying stringent hydrogen-bond criteria and within 30° of a canonical turn type (Rose et al. 1985). A total of 1000 trial conformations were generated for further testing.
To maintain downstream conformational integrity, trial conformations were eliminated if the RMSD between the last four
-carbons of the experimental structure and the trial structure exceeded 0.7 Å. Remaining trial conformations were then evaluated for steric clash using published hard-sphere radii (Srinivasan and Rose 1995). Clashes were relieved, if possible, by adjusting the positions of side chain atoms, using a combination of steepest descent and conjugate gradient minimization. Finally, surviving clash-free structures were evaluated using the CHASA procedure (Fleming et al. 2005) to assure that the hydrogen-bonding requirements of all polar atoms were satisfied (Fleming and Rose 2005), either by another H-bond partner within the protein or by a water molecule.
| Acknowledgments |
|---|
| References |
|---|
|
|
|---|
Ceccarelli, C., Jeffrey, G.A., and Taylor, R. 1981. A survey of O-H- - -O hydrogen-bond geometries determined by neutron diffraction. J. Mol. Struct. 70: 255271.[CrossRef]
Chou, P.Y. and Fasman, G.D. 1977.
etaturns in proteins. J. Mol. Biol. 115: 135175.[CrossRef][Medline]
Dasgupta, B., Pal, L., Basu, G., and Chakrabarti, P. 2004. Expanded turn conformations: Characterization and sequence-structure correspondence in
-turns with implications in helix folding. Proteins 55: 305315.[CrossRef][Medline]
Efimov, A.V. 1993. Standard structures in proteins. Prog. Biophys. Mol. Biol. 60: 201239.[CrossRef][Medline]
Fabiola, F., Bertram, R., Korostelev, A., and Chapman, M.S. 2002. An improved hydrogen bond potential: Impact on medium resolution protein structures. Protein Sci. 11: 14151423.
Fitzkee, N.C., Fleming, P.J., Gong, H., Panasik Jr., N., Street, T.O., and Rose, G.D. 2005a. Are proteins made from a limited parts list? Trends Biochem. Sci. 30: 7380.[CrossRef][Medline]
Fitzkee, N.C., Fleming, P.J., and Rose, G.D. 2005b. The Protein Coil Library: A structural database of nonhelix, nonstrand fragments derived from the PDB. Proteins 58: 852854.[CrossRef][Medline]
Fleming, P.J. and Rose, G.D. 2005. Do all backbone polar groups in proteins form hydrogen bonds? Protein Sci. 14: 19111917.
Fleming, P.J., Fitzkee, N.C., Mezei, M., Srinivasan, R., and Rose, G.D. 2005. A novel method reveals that solvent water favors polyproline II over
-strand conformation in peptides and unfolded proteins: Conditional hydrophobic accessible surface area (CHASA). Protein Sci. 14: 111118.
Hutchinson, E.G. and Thornton, J.M. 1994. A revised set of potentials for
-turn formation in proteins. Protein Sci. 3: 22072216.[Abstract]
Jeffrey, G.A. and Maluszynska, H. 1982. A survey of hydrogen-bond geometries in the crystal structures of amino acids. Int. J. Biol. Macromol. 4: 173185.
Jeffrey, G.A. and Saenger, W. 1991. Hydrogen bonding in biological structures. Springer-Verlag, Berlin.
Kabsch, W. and Sander, C. 1983. Dictionary of protein secondary structure: Pattern recognition of hydrogen-bonded and geometrical features. Biopolymers 22: 25772637.[CrossRef][Medline]
Kortemme, T., Morozov, A.V., and Baker, D. 2003. An orientation-dependent hydrogen bonding potential improves prediction of specificity and structure for proteins and proteinprotein complexes. J. Mol. Biol. 326: 12391259.[CrossRef][Medline]
Kuntz, I.D. 1972. Protein folding. J. Am. Chem. Soc. 94: 40094012.[CrossRef][Medline]
Lewis, P.N., Momany, F.A., and Scheraga, H.A. 1971. Folding of polypeptide chains in proteins: A proposed mechanism for folding. Proc. Natl. Acad. Sci. 68: 22932297.
Pappu, R.V. and Rose, G.D. 2002. A simple model for polyproline II structurein unfolded states of alanine-based peptides. Protein Sci. 11: 24372455.
Ramakrishnan, C. and Soman, K.V. 1982. Identification of secondary structures in globular proteinsA new algorithm. Int. J. Peptide Protein Res. 20: 218237.[Medline]
Richardson, J.S. 1981. The anatomy and taxonomy of protein structure. Adv. Protein Chem. 34: 168340.
Rose, G.D., Gierasch, L.M., and Smith, J.A. 1985. Turns in peptides and proteins. Adv. Protein Chem. 37: 1109.[Medline]
Srinivasan, R. and Rose, G.D. 1995. LINUSA simple algorithm to predict the fold of a protein. Proteins 22: 8199.[CrossRef][Medline]
Srinivasan, R., Fleming, P.J., and Rose, G.D. 2004. Ab initio protein folding using LINUS. Methods Enzymol. 383: 4866.[Medline]
Stapley, B.J. and Creamer, T.P. 1999. A survey of left-handed polyproline II helices. Protein Sci. 8: 587595.[Abstract]
Stickle, D.F., Presta, L.G., Dill, K.A., and Rose, G.D. 1992. Hydrogen bonding in globular proteins. J. Mol. Biol. 226: 11431159.[CrossRef][Medline]
Taylor, R. and Kennard, O. 1983. Comparison of X-ray and neutron diffraction results for the NH- - -O=C hydrogen-bond. Acta Crystallogr. 39: 133138.
. 1984. Hydrogen-bond geometry in organic crystals. Acc. Chem. Res. 17: 320326.[CrossRef]
Venkatachalam, C.M. 1968. Stereochemical criteria for polypeptides and proteins. V. Conformation of a system of three linked peptide units. Biopolymers 6: 14251436.[CrossRef][Medline]
Wang, G. and Dunbrack Jr., R.L. 2003. PISCES: A protein sequence culling server. Bioinformatics 19: 15891591.
Wilmot, C.M. and Thornton, J.M. 1988. Analysis and prediction of the different types of
-turns in proteins. J. Mol. Biol. 203: 221232.[CrossRef][Medline]
. 1990.
-turns and their distortions: A proposed new nomenclature. Protein Eng. 3: 479493.
![]()
CiteULike
Connotea
Del.icio.us
Digg
Reddit
Technorati What's this?
This article has been cited by other articles:
![]() |
L. L. Perskie, T. O. Street, and G. D. Rose Structures, basins, and energies: A deconstruction of the Protein Coil Library Protein Sci., July 1, 2008; 17(7): 1151 - 1161. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. O. Street, N. C. Fitzkee, L. L. Perskie, and G. D. Rose Physical-chemical determinants of turn conformations in globular proteins Protein Sci., August 1, 2007; 16(8): 1720 - 1727. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. D. Rose, P. J. Fleming, J. R. Banavar, and A. Maritan A backbone-based theory of protein folding PNAS, November 7, 2006; 103(45): 16623 - 16633. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. W. Buchko, S. Ni, H. Robinson, E. A. Welsh, H. B. Pakrasi, and M. A. Kennedy Characterization of two potentially universal turn motifs that shape the repeated five-residues fold--Crystal structure of a lumenal pentapeptide repeat protein from Cyanothece 51142. Protein Sci., November 1, 2006; 15(11): 2579 - 2595. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Avbelj, S. G. Grdadolnik, J. Grdadolnik, and R. L. Baldwin Intrinsic backbone preferences are fully present in blocked amino acids PNAS, January 31, 2006; 103(5): 1272 - 1277. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |