Native and non-native secondary structure and dynamics in the pH 4 intermediate of apomyoglobin

The partly folded state of apomyoglobin at pH 4 represents an excellent model for an obligatory kinetic folding intermediate. The structure and dynamics of this intermediate state have been extensively examined using NMR spectroscopy. Secondary chemical shifts, (1)H-(1)H NOEs, and amide proton temperature coefficients have been used to probe residual structure in the intermediate state, and NMR relaxation parameters T(1) and T(2) and ?(1)H?-(15)N NOE have been analyzed using spectral densities to correlate motion of the polypeptide chain with these structural observations. A significant amount of helical structure remains in the pH 4 state, indicated by the secondary chemical shifts of the (13)C(alpha), (13)CO, (1)H(alpha), and (13)C(beta) nuclei, and the boundaries of this helical structure are confirmed by the locations of (1)H-(1)H NOEs. Hydrogen bonding in the structured regions is predominantly native-like according to the amide proton chemical shifts and their temperature dependence. The locations of the A, G, and H helix segments and the C-terminal part of the B helix are similar to those in native apomyoglobin, consistent with the early, complete protection of the amides of residues in these helices in quench-flow experiments. These results confirm the similarity of the equilibrium form of apoMb at pH 4 and the kinetic intermediate observed at short times in the quench-flow experiment. Flexibility in this structured core is severely curtailed compared with the remainder of the protein, as indicated by the analysis of the NMR relaxation parameters. Regions with relatively high values of J(0) and low values of J(750) correspond well with the A, B, G, and H helices, an indication that nanosecond time scale backbone fluctuations in these regions of the sequence are restricted. Other parts of the protein show much greater flexibility and much reduced secondary chemical shifts. Nevertheless, several regions show evidence of the beginnings of helical structure, including stretches encompassing the C helix-CD loop, the boundary of the D and E helices, and the C-terminal half of the E helix. These regions are clearly not well-structured in the pH 4 state, unlike the A, B, G, and H helices, which form a native-like structured core. However, the proximity of this structured core most likely influences the region between the B and F helices, inducing at least transient helical structure.     

Conformational and dynamic characterization of the molten globule state of an apomyoglobin mutant with an altered folding pathway.

Kinetic and equilibrium studies of apomyoglobin folding pathways and intermediates have provided important insights into the mechanism of protein folding. To investigate the role of intrinsic helical propensities in the apomyoglobin folding process, a mutant has been prepared in which Asn132 and Glu136 have been substituted with glycine to destabilize the H helix. The structure and dynamics of the equilibrium molten globule state formed at pH 4.1 have been examined using NMR spectroscopy. Deviations of backbone (13)C(alpha) and (13)CO chemical shifts from random coil values reveal high populations of helical structure in the A and G helix regions and in part of the B helix. However, the H helix is significantly destabilized compared to the wild-type molten globule. Heteronuclear [(1)H]-(15)N NOEs show that, although the polypeptide backbone in the H helix region is more flexible than in the wild-type protein, its motions are restricted by transient hydrophobic interactions with the molten globule core. Quench flow hydrogen exchange measurements reveal stable helical structure in the A and G helices and part of the B helix in the burst phase kinetic intermediate and confirm that the H helix is largely unstructured. Stabilization of structure in the H helix occurs during the slow folding phases, in synchrony with the C and E helices and the CD region. The kinetic and equilibrium molten globule intermediates formed by N132G/E136G are similar in structure. Although both the wild-type apomyoglobin and the mutant fold via compact helical intermediates, the structures of the intermediates and consequently the detailed folding pathways differ. Apomyoglobin is therefore capable of compensating for mutations by using alternative folding pathways within a common basic framework. Tertiary hydrophobic interactions appear to play an important role in the formation and stabilization of secondary structure in the H helix of the N132G/E136G mutant. These studies provide important insights into the interplay between secondary and tertiary structure formation in protein folding.     

Structural characterization of partially folded intermediates of apomyoglobin H64F

We present a detailed investigation of unfolded and partially folded states of a mutant apomyoglobin (apoMb) where the distal histidine has been replaced by phenylalanine (H64F). Previous studies have shown that substitution of His64, located in the E helix of the native protein, stabilizes the equilibrium molten globule and native states and leads to an increase in folding rate and a change in the folding pathway. Analysis of changes in chemical shift and in backbone flexibility, detected via [1H]-15N heteronuclear nuclear Overhauser effect measurements, indicates that the phenylalanine substitution has only minor effects on the conformational ensemble in the acid- and urea-unfolded states, but has a substantial effect on the structure, dynamics, and stability of the equilibrium molten globule intermediate formed near pH 4. In H64F apomyoglobin, additional regions of the polypeptide chain are recruited into the compact core of the molten globule. Since the phenylalanine substitution has negligible effect on the unfolded ensemble, its influence on folding rate and stability comes entirely from interactions within the compact folded or partly folded states. Replacement of His64 with Phe leads to favorable hydrophobic packing between the helix E region and the molten globule core and leads to stabilization of helix E secondary structure and overall thermodynamic stabilization of the molten globule. The secondary structure of the equilibrium molten globule parallels that of the burst phase kinetic intermediate; both intermediates contain significant helical structure in regions of the polypeptide that comprise the A, B, E, G, and H helices of the fully folded protein.     

Structural characterization of unfolded states of apomyoglobin using residual dipolar couplings

The conformational propensities of unfolded states of apomyoglobin have been investigated by measurement of residual dipolar couplings between (15)N and (1)H in backbone amide groups. Weak alignment of apomyoglobin in acid and urea-unfolded states was induced with both stretched and compressed polyacrylamide gels. In 8 M urea solution at pH 2.3, conditions under which apomyoglobin contains no detectable secondary or tertiary structure, significant residual dipolar couplings of uniform sign were observed for all residues. At pH 2.3 in the absence of urea, a change in the magnitude and/or sign of the residual dipolar couplings occurs in local regions of the polypeptide where there is a high propensity for helical secondary structure. These results are interpreted on the basis of the statistical properties of the unfolded polypeptide chain, viewed as a polymer of statistical segments. For a folded protein, the magnitude and sign of the residual dipolar couplings depend on the orientation of each bond vector relative to the alignment tensor of the entire molecule, which reorients as a single entity. For unfolded proteins, there is no global alignment tensor; instead, residual dipolar couplings are attributed to alignment of the statistical segments or of transient elements of secondary structure. For apomyoglobin in 8 M urea, the backbone is highly extended, with phi and psi dihedral angles favoring the beta or P(II) regions. Each statistical segment has a highly anisotropic shape, with the N-H bond vectors approximately perpendicular to the long axis, and becomes weakly aligned in the anisotropic environment of the strained acrylamide gels. Local regions of enhanced flexibility or chain compaction are characterized by a decrease in the magnitude of the residual dipolar couplings. The formation of a small population of helical structure in the acid-denatured state of apomyoglobin leads to a change in sign of the residual dipolar couplings in local regions of the polypeptide; the population of helix estimated from the residual dipolar couplings is in excellent agreement with that determined from chemical shifts. The alignment model described here for apomyoglobin can also explain the pattern of residual dipolar couplings reported previously for denatured states of staphylococcal nuclease and other proteins. In conjunction with other NMR experiments, residual dipolar couplings can provide valuable insights into the dynamic conformational propensities of unfolded and partly folded states of proteins and thereby help to chart the upper reaches of the folding landscape.     

Structural and dynamic characterization of the acid-unfolded state of hUBF HMG box 1 provides clues for the early events in protein folding

To understand the events that occur in the early stages of the folding of hUBF HMG box 1, we characterized its pH 2.1 unfolded state in detail with NMR. Through a triple resonance strategy, the assignments of complete backbone and some side chains were achieved. Then, significant conformational information was extracted from secondary chemical shifts, interresidual (1)H-(1)H NOEs, (3)J(HNHA) coupling constants, amide proton temperature coefficients, and (15)N relaxation data. The secondary chemical shifts for (13)CA, (13)CB, (13)CO, (1)HA, and (1)HN indicate that the residues between 64 and 78 exhibit a substantial preference for helical structure in the acid-unfolded state, which is also evidenced by the relatively more negative deviations of (3)J(HNHA) and amide proton temperature coefficients from their corresponding random-coil values and particularly confirmed by the strongest sequential d(NN)(i, i + 1) proton NOEs along the region. Following this region until residue 82 is a segment that tends to form a turn-like structure, which is unstable and exchanges between alternative states. In addition, some evidences imply that the regions 18-28 and 38-43 also possess propensities for helical structure but to a different less degree than the region 64-78. The polypeptide backbone dynamics investigated using reduced spectral density function shows apparent motional restrictions in residual structural regions and to less extent at some hydrophobic residues. On the basis of the results presented herein, we propose a potential protein-folding pathway on which these residual structures play a role of initiation site in the early folding stages.     

Native and nonnative conformational preferences in the urea-unfolded state of barstar

The refolding of barstar from its urea-unfolded state has been studied extensively using various spectroscopic probes and real-time NMR, which provide global and residue-specific information, respectively, about the folding process. Here, a preliminary structural characterization by NMR of barstar in 8 M urea has been carried out at pH 6.5 and 25 degrees C. Complete backbone resonance assignments of the urea-unfolded protein were obtained using the recently developed three-dimensional NMR techniques of HNN and HN(C)N. The conformational propensities of the polypeptide backbone in the presence of 8 M urea have been estimated by examining deviations of secondary chemical shifts from random coil values. For some residues that belong to helices in native barstar, 13C(alpha) and 13CO secondary shifts show positive deviations in the urea-unfolded state, indicating that these residues have propensities toward helical conformations. These residues are, however, juxtaposed by residues that display negative deviations indicative of propensities toward extended conformations. Thus, segments that are helical in native barstar are unlikely to preferentially populate the helical conformation in the unfolded state. Similarly, residues belonging to beta-strands 1 and 2 of native barstar do not appear to show any conformational preferences in the unfolded state. On the other hand, residues belonging to the beta-strand 3 segment show weak nonnative helical conformational preferences in the unfolded state, indicating that this segment may possess a weak preference for populating a helical conformation in the unfolded state.     

Relationship between nuclear magnetic resonance chemical shift and protein secondary structure

An analysis of the 1H nuclear magnetic resonance chemical shift assignments and secondary structure designations for over 70 proteins has revealed some very strong and unexpected relationships. Similar studies, performed on smaller databases, for 13C and 15N chemical shifts reveal equally strong correlations to protein secondary structure. Among the more interesting results to emerge from this work is the finding that all 20 naturally occurring amino acids experience a mean alpha-1H upfield shift of 0.39 parts per million (from the random coil value) when placed in a helical configuration. In a like manner, the alpha-1H chemical shift is found to move downfield by an average of 0.37 parts per million when the residue is placed in a beta-strand or extended configuration. Similar changes are also found for amide 1H, carbonyl 13C, alpha-13C and amide 15N chemical shifts. Other relationships between chemical shift and protein conformation are also uncovered; in particular, a correlation between helix dipole effects and amide proton chemical shifts as well as a relationship between alpha-proton chemical shifts and main-chain flexibility. Additionally, useful relationships between alpha-proton chemical shifts and backbone dihedral angles as well as correlations between amide proton chemical shifts and hydrogen bond effects are demonstrated.     

Conservation of folding pathways in evolutionarily distant globin sequences

To test the hypothesis that the folding pathways of evolutionarily related proteins with similar three-dimensional structures but widely different sequences should be similar, the folding pathway of apoleghemoglobin has been characterized using stopped-flow circular dichroism, heteronuclear NMR pulse labeling techniques and mass spectrometry. The pathway of folding was found to differ significantly from that of a protein of the same family, apomyoglobin, although both proteins appear to fold through helical burst phase intermediates. For leghemoglobin, the burst phase intermediate exhibits stable helical structure in the G and H helices, together with a small region in the center of the E helix. The A and B helices are not stabilized until later stages of the folding process. The structure of the burst phase folding intermediate thus differs from that of apomyoglobin, in which stable helical structure is formed in the A, B, G and H helix regions.     

Characterization of "native" apomyoglobin by molecular dynamics simulation.

We have used molecular dynamics simulation methods to study the structure and fluctuations of "native" apomyoglobin in aqueous solution for a period of greater than 0.5 nanosecond. This work was motivated by the recent attempts of Hughson et al. to characterize the structure and motion of both this molecule and the less compact, acid stabilized I stage, using methods of pulsed H/2H exchange. The study of these systems provides new insights into protein folding intermediates and our simulation has yielded a detailed model for structure and fluctuations in apomyoglobin which complements the experimental studies. We find that local (short-time) fluctuations agree well with fluctuations observed for the holoprotein in aqueous solution, as well as results from the crystallographic B-factors. In addition, the structural features we observe for native apomyoglobin are very similar to the holoprotein, in basic agreement with the findings of Hughson et al. By examining larger-scale motions, developing only over timescales in excess of a 100 picoseconds, we are able to identify conformationally "labile" and "non-labile" regions within native apomyoglobin. These regions correspond extremely well to those identified in the nuclear magnetic resonance experiments as unstable and stable "folding subdomains" in the I state of apomyoglobin. Overall we find that helices A, B, E, G and H show the least amount of motion and helices C, D and F move substantially over the timescales examined. The major motions, and the primary difference between the holo and apo structures as we have observed them, are due to the shifting motion of helices C, D and F into the vacant heme cavity. We also find that motions at the interface of helical segments can be large, with one important exception being the chain segment connecting helices G and H. This segment of chain interacts with the conformationally "non-labile" helix A to form a relatively rigid subdomain composed of helices A, G and H. We believe that these findings provide direct support for the suggestion of Hughson et al. that helices A, G and H constitute a compact subdomain that remains in a native-like conformation as the protein begins to unfold in environments of decreasing pH.     

The C terminus of apocytochrome b562 undergoes fast motions and slow exchange among ordered conformations resembling the folded state

The present work describes the dynamics of the apo form of cytochrome b(562), a small soluble protein consisting of 106 amino acid residues [Itagaki, E., and Hager, L. P. (1966) J. Biol. Chem. 241, 3687-3695]. The presence of exchange in the millisecond time scale is demonstrated for the last part of helix IV (residues 95-105 in the holo form). The chemical shift index analysis [Wishart, D. S., and Sykes, B. D. (1994) J. Biomol. NMR 4, 171-180] based on H(alpha), C(alpha), C(beta), and C' chemical shifts suggests a larger helical content than shown in the NMR structure based on NOEs. These results indicate the presence of helical-like conformations participating in the exchange process. This hypothesis is consistent with amide deuterium exchange rates and the presence of some hydrogen bonds identified from amide chemical shift temperature coefficients [Baxter, N. J., and Williamson, M. P. (1997) J. Biomol. NMR 9, 359-369]. (15)N relaxation indicates limited mobility for the amide protons of this part of the helix in the picosecond time scale. A 30 ns stochastic dynamics simulation shows small fluctuations around the helical conformation on this time scale. These fluctuations, however, do not result in a significant decrease of the calculated order parameters which are consistent with the experimental (15)N relaxation data. These results resolve an apparent discrepancy in the NMR structures between the disorder observed in helix IV due to a lack of NOEs and the secondary structure predictions based on H(alpha) chemical shifts [Feng, Y., Wand, A. J., and Sligar, S. G. (1994) Struct. Biol. 1, 30-35].     

Consequences of stabilizing the natively disordered f helix for the folding pathway of apomyoglobin

The F helix region of sperm whale apomyoglobin is disordered, undergoing conformational fluctuations between a folded helical conformation and one or more locally unfolded states. To examine the effects of F helix stabilization on the folding pathway of apomyoglobin, we have introduced mutations to augment intrinsic helical structure in the F helix of the kinetic folding intermediate and to increase its propensity to fold early in the pathway, using predictions based on plots of the average area buried upon folding (AABUF) derived from the primary sequence. Two mutant proteins were prepared: a double mutant, P88K/S92K (F2), and a quadruple mutant, P88K/A90L/S92K/A94L (F4). Whereas the AABUF for F2 predicts that the F helix will not fold early in the pathway, the F helix in F4 shows a significantly increased AABUF and is therefore predicted to fold early. Protection of amide protons by formation of hydrogen-bonded helical structure during the early folding events has been analyzed by pH-pulse labeling. Consistent with the AABUF prediction, many of the F helix residues for F4 are significantly protected in the kinetic intermediate but are not protected in the F2 mutant. F4 folds via a kinetically trapped burst-phase intermediate that contains stabilized secondary structure in the A, B, F, G, and H helix regions. Rapid folding of the F helix stabilizes the central core of the misfolded intermediate and inhibits translocation of the H helix back to its native position, thereby decreasing the overall folding rate.     

An evaluation of chemical shift index-based secondary structure determination in proteins: Influence of random coil chemical shifts

Random coil chemical shifts are commonly used to detect protein secondary structural elements in chemical shift index (CSI) calculations. Though this technique is widely used and seems reliable for folded proteins, the choice of reference random coil chemical shift values can significantly alter the outcome of secondary structure estimation. In order to evaluate these effects, we present a comparison of secondary structure content calculated using CSI, based on five different reference random coil chemical shift value sets, to that derived from three-dimensional structures.Our results show that none of the reference random coil data sets chosen for evaluation fully reproduces the actual secondary structures. Among the reference values generally available to date, most tend to be good estimators only of helices. Based on our evaluation, we recommend the experimental values measured by Schwarzinger et al.(2000), and statistical values obtained by Lukin et al. (1997), as good estimators of both helical and sheet content.     

Dihydrofolate reductase: sequential resonance assignments using 2D and 3D NMR and secondary structure determination in solution

Three-dimensional (3D) heteronuclear NMR techniques have been used to make sequential 1H and 15N resonance assignments for most of the residues of Lactobacillus casei dihydrofolate reductase (DHFR), a monomeric protein of molecular mass 18,300 Da. A uniformly 15N-labeled sample of the protein was prepared and its complex with methotrexate (MTX) studied by 3D 15N/1H nuclear Overhauser-heteronuclear multiple quantum coherence (NOESY-HMQC), Hartmann-Hahn-heteronuclear multiple quantum coherence (HOHAHA-HMQC), and HMQC-NOESY-HMQC experiments. These experiments overcame most of the spectral overlap problems caused by chemical shift degeneracies in 2D spectra and allowed the 1H-1H through-space and through-bond connectivities to be identified unambiguously, leading to the resonance assignments. The novel HMQC-NOESY-HMQC experiment allows NOE cross peaks to be detected between NH protons even when their 1H chemical shifts are degenerate as long as the amide 15N chemical shifts are nondegenerate. The 3D experiments, in combination with conventional 2D NOESY, COSY, and HOHAHA experiments on unlabelled and selectively deuterated DHFR, provide backbone assignments for 146 of the 162 residues and side-chain assignments for 104 residues of the protein. Data from the NOE-based experiments and identification of the slowly exchanging amide protons provide detailed information about the secondary structure of the binary complex of the protein with methotrexate. Sequential NHi-NHi+1 NOEs define four regions with helical structure. Two of these regions, residues 44-49 and 79-89, correspond to within one amino acid to helices C and E in the crystal structure of the DHFR.methotrexate.NADPH complex [Bolin et al. (1982) J. Biol. Chem. 257, 13650-13662], while the NMR-determined helix formed by residues 26-35 is about one turn shorter at the N-terminus than helix B in the crystal structure, which spans residues 23-34. Similarly, the NMR-determined helical region comprising residues 102-110 is somewhat offset from the crystal structure's helix F, which encompasses residues 97-107. Regions of beta-sheet structure were characterized in the binary complex by strong alpha CHi-NHi+1 NOEs and by slowly exchanging amide protons. In addition, several long-range NOEs were identified linking together these stretches to form a beta-sheet. These elements align perfectly with corresponding elements in the crystal structure of the DHFR.methotrexate.NADPH complex, which contains an eight-stranded beta-sheet, indicating that the main body of the beta-sheet is preserved in the binary complex in solution.     

Solution structure of the osteogenic 1-31 fragment of the human parathyroid hormone

The solution conformations of a selectively osteogenic 1-31 fragment of the human parathyroid hormone (hPTH), hPTH(1-31)NH(2), have been characterized by use of very high field NMR spectroscopy at 800 MHz. The combination of the CalphaH proton and (13)Calpha chemical shifts, (3)J(NH)(alpha) coupling constants, NH proton temperature coefficients, and backbone NOEs reveals that the hPTH(1-31)NH(2) peptide has well-formed helical structures localized in two distinct segments of the polypeptide backbone. There are also many characteristic NOEs defining specific side-chain/backbone and side-chain/side-chain contacts within both helical structures. The solution structure of hPTH(1-31)NH(2) contains a short N-terminal helical segment for residues 3-11, including the helix capping residues 3 and 11 and a long C-terminal helix for residues 16-30. The two helical structures are reinforced by well-defined capping motifs and side-chain packing interactions within and at both ends of these helices. On one face of the C-terminal helix, there are side-chain pairs of Glu22-Arg25, Glu22-Lys26, and Arg25-Gln29 that can form ion-pair and/or hydrogen bonding interactions. On the opposite face of this helix, there are characteristic hydrophobic interactions involving the aromatic side chain of Trp23 packing against the aliphatic side chains of Leu15, Leu24, Lys27, and Leu28. There is also a linear array of hydrophobic residues from Val2, to Leu7, to Leu11 and continuing on to residues His14 and Leu15 in the hinge region and to Trp23 in the C-terminal helix. Capping and hydrophobic interactions at the end of the N-terminal and at the beginning of the C-terminal helix appear to consolidate the helical structures into a V-shaped overall conformation for at least the folded population of the hPTH(1-31)NH(2) peptide. Stabilization of well-folded conformations in this linear 1-31 peptide fragment and possibly other analogues of human PTH may have a significant impact on the biological activities of the PTH peptides in general and specifically for the osteogenic/anabolic activities of bone-building PTH analogues.     

Identification of initiation sites for T4 lysozyme folding using CD and NMR spectroscopy of peptide fragments

Using CD and 2D (1)H NMR spectroscopy, we have identified potential initiation sites for the folding of T4 lysozyme by examining the conformational preferences of peptide fragments corresponding to regions of secondary structure. CD spectropolarimetry showed most peptides were unstructured in water, but adopted partial helical conformations in TFE and SDS solution. This was also consistent with the (1)H NMR data which showed that the peptides were predominantly disordered in water, although in some cases, nascent or small populations of partially folded conformations could be detected. NOE patterns, coupling constants, and deviations from random coil Halpha chemical shift values complemented the CD data and confirmed that many of the peptides were helical in TFE and SDS micelles. In particular, the peptide corresponding to helix E in the native enzyme formed a well-defined helix in both TFE and SDS, indicating that helix E potentially forms an initiation site for T4 lysozyme folding. The data for the other peptides indicated that helices D, F, G, and H are dependent on tertiary interactions for their folding and/or stability. Overall, the results from this study, and those of our earlier studies, are in agreement with modeling and HD-deuterium exchange experiments, and support an hierarchical model of folding for T4 lysozyme.     

1H, 13C, and 15N NMR backbone assignments and secondary structure of human interferon-gamma.

1H, 13C, and 15N NMR assignments of the protein backbone of human interferon-gamma, a homodimer of 31.4 kDa, have been made using the recently introduced three-dimensional (3D) triple-resonance NMR techniques. It is shown that, despite the approximately 40-50-Hz 13C alpha and 1H alpha line widths of this high molecular weight dimer and the extensive overlap in the 1H alpha and 13C alpha spectral regions, unique sequential assignments can be made on the basis of combined use of the 3D HNCO, HNCA, HN(CO)CA, and HCACO constant-time experiments, the 15N-separated 3D NOESY-HMQC, and the 3D HOHAHA-HMQC experiments. Analysis of the 15N-separated 3D NOESY-HMQC and 13C/15N-separated four-dimensional (4D) NOESY-HMQC spectra together with the secondary C alpha and C beta chemical shifts yielded extensive secondary structure information. The NMR-derived secondary structure essentially confirms results of a recently published low-resolution crystal structure [Ealick et al. (1991) Science 252, 698-702], i.e., six helices in the monomer which are mostly alpha-helical in nature, no beta-sheets, a long flexible loop between helices A and B, and a very hydrophobic helix C. The functionally important carboxy terminus, which was not observed in the X-ray study, does not adopt a rigid conformation in solution. A high degree of internal mobility, starting at Pro-123, gives rise to significantly narrower resonance line widths for these carboxy-terminal residues compared to the rest of the protein.     

Structural and dynamic characterization of an unfolded state of poplar apo-plastocyanin formed under nondenaturing conditions

Plastocyanin is a predominantly beta-sheet protein containing a type I copper center. The conformational ensemble of a denatured state of apo-plastocyanin formed in solution under conditions of low salt and neutral pH has been investigated by multidimensional heteronuclear NMR spectroscopy. Chemical shift assignments were obtained by using three-dimensional triple-resonance NMR experiments to trace through-bond heteronuclear connectivities along the backbone and side chains. The (3)J(HN,Halpha) coupling constants, (15)N-edited proton-proton nuclear Overhauser effects (NOEs), and (15)N relaxation parameters were also measured for the purpose of structural and dynamic characterization. Most of the residues corresponding to beta-strands in the folded protein exhibit small upfield shifts of the (13)C(alpha) and (13)CO resonances relative to random coil values, suggesting a slight preference for backbone dihedral angles in the beta region of (phi,psi) space. This is further supported by the presence of strong sequential d(alphaN)(i, i + 1) NOEs throughout the sequence. The few d(NN)(i, i + 1) proton NOEs that are observed are mostly in regions that form loops in the native plastocyanin structure. No medium or long-range NOEs were observed. A short sequence, between residues 59 and 63, was found to populate a nonnative helical conformation in the unfolded state, as indicated by the shift of the (13)C(alpha), (13)CO, and (1)H(alpha) resonances relative to random coil values and by the decreased values of the (3)J(HN,Halpha) coupling constants. The (15)N relaxation parameters indicate restriction of motions on a nanosecond timescale in this region. Intriguingly, this helical conformation is present in a sequence that is close to but not in the same location as the single short helix in the native folded protein. The results are consistent with earlier NMR studies of peptide fragments of plastocyanin and confirm that the regions of the sequence that form beta-strands in the native protein spontaneously populate the beta-region of (phi,psi) space under folding conditions, even in the absence of stabilizing tertiary interactions. We conclude that the state of apo-plastocyanin present under nondenaturing conditions is a noncompact unfolded state with some evidence of nativelike and nonnative local structuring that may be initiation sites for folding of the protein.     

Probability-based protein secondary structure identification using combined NMR chemical-shift data

For a long time, NMR chemical shifts have been used to identify protein secondary structures. Currently, this is accomplished through comparing the observed (1)H(alpha), (13)C(alpha), (13)C(beta), or (13)C' chemical shifts with the random coil values. Here, we present a new protocol, which is based on the joint probability of each of the three secondary structural types (beta-strand, alpha-helix, and random coil) derived from chemical-shift data, to identify the secondary structure. In combination with empirical smooth filters/functions, this protocol shows significant improvements in the accuracy and the confidence of identification. Updated chemical-shift statistics are reported, on the basis of which the reliability of using chemical shift to identify protein secondary structure is evaluated for each nucleus. The reliability varies greatly among the 20 amino acids, but, on average, is in the order of: (13)C(alpha)>(13)C'>(1)H(alpha)>(13)C(beta)>(15)N>(1)H(N) to distinguish an alpha-helix from a random coil; and (1)H(alpha)>(13)C(beta) >(1)H(N) approximately (13)C(alpha) approximately (13)C' approximately (15)N for a beta-strand from a random coil. Amide (15)N and (1)H(N) chemical shifts, which are generally excluded from the application, in fact, were found to be helpful in distinguishing a beta-strand from a random coil. In addition, the chemical-shift statistical data are compared with those reported previously, and the results are discussed. A JAVA User Interface program has been developed to make the entire procedure fully automated and is available via http://ccsr3150-p3.stanford.edu.     

Random coil structures in bacterial proteins. Relationships of their amino acid compositions to flanking structures and corresponding genic base compositions

In this study we classified regions of random coil into four types: coil between alpha helix and beta strand, coil between beta strand and alpha helix, coil between two alpha helices and coil between two beta strands. This classification may be considered as natural. We used 610 3D structures of proteins collected from the Protein Data Bank from bacteria with low, average and high genomic GC-content. Relatively short regions of coil are not random: certain amino acid residues are more or less frequent in each of the types of coil. Namely, hydrophobic amino acids with branched side chains (Ile, Val and Leu) are rare in coil between two beta strands, unlike some acrophilic amino acids (Asp, Asn and Gly). In contrast, coil between two alpha helices is enriched by Leu. Regions of coil between alpha helix and beta strand are enriched by positively charged amino acids (Arg and Lys), while the usage of residues with side chains possessing hydroxyl group (Ser and Thr) is low in them, in contrast to the regions of coil between beta strand and alpha helix. Regions of coil between beta strand and alpha helix are significantly enriched by Cys residues. The response to the symmetric mutational pressure (AT-pressure or GC-pressure) is also quite different for four types of coil. The most conserved regions of coil are “connecting bridges” between beta strand and alpha helix, since their amino acid content shows less strong dependence on GC-content of genes than amino acid contents of other three types of coil. Possible causes and consequences of the described differences in amino acid content distribution between different types of random coil have been discussed.     

Support vector machines for prediction of dihedral angle regions

MOTIVATION: Most secondary structure prediction programs target only alpha helix and beta sheet structures and summarize all other structures in the random coil pseudo class. However, such an assignment often ignores existing local ordering in so-called random coil regions. Signatures for such ordering are distinct dihedral angle pattern. For this reason, we propose as an alternative approach to predict directly dihedral regions for each residue as this leads to a higher amount of structural information. RESULTS: We propose a multi-step support vector machine (SVM) procedure, dihedral prediction (DHPRED), to predict the dihedral angle state of residues from sequence. Trained on 20,000 residues our approach leads to dihedral region predictions, that in regions without alpha helices or beta sheets is higher than those from secondary structure prediction programs. AVAILABILITY: DHPRED has been implemented as a web service, which academic researchers can access from our webpage http://www.fz-juelich.de/nic/cbb