Plasmodium vivax CS peptides display conformational preferences for folded forms in solution.

Spectroscopic techniques have been used to study the conformations of several synthetic peptides with sequences corresponding to the repeat regions of the circumsporozoite proteins of Plasmodium vivax, variants vk-210 and vk-247. As has previously been shown for P. falciparum, turn-like folded conformations are observed, in rapid dynamic equilibrium with extended-chain forms. These results are consistent with the known similarity of the structural, biosynthetic and immunological properties of the circumsporozoite proteins of different plasmodial species. Additionally, the observation of folded conformers provides a rationale for the effectiveness of these peptides as immunogens and potential vaccines.

Potential bias in NMR relaxation data introduced by peak intensity analysis and curve fitting methods.

We present an evaluation of the accuracy and precision of relaxation rates calculated using a variety of methods, applied to data sets obtained for several very different protein systems. We show that common methods of data evaluation, such as the determination of peak heights and peak volumes, may be subject to bias, giving incorrect values for quantities such as R1 and R2. For example, one common method of peak-height determination, using a search routine to obtain the peak-height maximum in successive spectra, may be a source of significant systematic error in the relaxation rate. The alternative use of peak volumes or of a fixed coordinate position for the peak height in successive spectra gives more accurate results, particularly in cases where the signal/noise is low, but these methods have inherent problems of their own. For example, volumes are difficult to quantitate for overlapped peaks. We show that with any method of sampling the peak intensity, the choice of a 2- or 3-parameter equation to fit the exponential relaxation decay curves can dramatically affect both the accuracy and precision of the calculated relaxation rates. In general, a 2-parameter fit of relaxation decay curves is preferable. However, for very low intensity peaks a 3 parameter fit may be more appropriate.

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.

Backbone dynamics in dihydrofolate reductase complexes: role of loop flexibility in the catalytic mechanism.

To elucidate the influence of local motion of the polypeptide chain on the catalytic mechanism of an enzyme, we have measured (15)N relaxation data for Escherichia coli dihydrofolate reductase in three different complexes, representing different stages in the catalytic cycle of the enzyme. NMR relaxation data were analyzed by the model-free approach, corrected for rotational anisotropy, to provide insights into the backbone dynamics. There are significant differences in the backbone dynamics in the different complexes. Complexes in which the cofactor binding site is occluded by the Met20 loop display large amplitude motions on the picosecond/nanosecond time scale for residues in the Met20 loop, the adjacent betaF-betaG loop and for residues 67-69 in the adenosine binding loop. Formation of the closed Met20 loop conformation in the ternary complex with folate and NADP(+), results in attenuation of the motions in the Met20 loop and the betaF-betaG loop but leads to increased flexibility in the adenosine binding loop. New fluctuations on a microsecond/millisecond time scale are observed in the closed E:folate:NADP(+) complex in regions that form hydrogen bonds between the Met20 and the betaF-betaG loops. The data provide insights into the changes in backbone dynamics during the catalytic cycle and point to an important role of the Met20 and betaF-betaG loops in controlling access to the active site. The high flexibility of these loops in the occluded conformation is expected to promote tetrahydrofolate-assisted product release and facilitate binding of the nicotinamide ring to form the Michaelis complex. The backbone fluctuations in the Met20 loop become attenuated once it closes over the active site, thereby stabilizing the nicotinamide ring in a geometry conducive to hydride transfer. Finally, the relaxation data provide evidence for long-range motional coupling between the adenosine binding loop and distant regions of the protein.

Solution structure of Escherichia coli glutaredoxin-2 shows similarity to mammalian glutathione-S-transferases.

Glutaredoxin 2 (Grx2) from Escherichia coli is distinguished from other glutaredoxins by its larger size, low overall sequence identity and lack of electron donor activity with ribonucleotide reductase. However, catalysis of glutathione (GSH)-dependent general disulfide reduction by Grx2 is extremely efficient. The high-resolution solution structure of E. coli Grx2 shows a two-domain protein, with residues 1 to 72 forming a classical "thioredoxin-fold" glutaredoxin domain, connected by an 11 residue linker to the highly helical C-terminal domain, residues 84 to 215. The active site, Cys9-Pro10-Tyr11-Cys12, is buried in the interface between the two domains, but Cys9 is solvent-accessible, consistent with its role in catalysis. The structures reveal the hither to unknown fact that Grx2 is structurally similar to glutathione-S-transferases (GST), although there is no obvious sequence homology. The similarity of these structures gives important insights into the functional significance of a new class of mammalian GST-like proteins, the single-cysteine omega class, which have glutaredoxin oxidoreductase activity rather than GSH-S-transferase conjugating activity. E. coli Grx 2 is structurally and functionally a member of this new expanding family of large glutaredoxins. The primary function of Grx2 as a GST-like glutaredoxin is to catalyze reversible glutathionylation of proteins with GSH in cellular redox regulation including stress responses.

Sequence-dependent correction of random coil NMR chemical shifts.

Random coil chemical shifts are commonly used to detect secondary structure elements in proteins in chemical shift index calculations. While this technique is very reliable for folded proteins, application to unfolded proteins reveals significant deviations from measured random coil shifts for certain nuclei. While some of these deviations can be ascribed to residual structure in the unfolded protein, others are clearly caused by local sequence effects. In particular, the amide nitrogen, amide proton, and carbonyl carbon chemical shifts are highly sensitive to the local amino acid sequence. We present a detailed, quantitative analysis of the effect of the 20 naturally occurring amino acids on the random coil shifts of (15)N(H), (1)H(N), and (13)CO resonances of neighboring residues, utilizing complete resonance assignments for a set of five-residue peptides Ac-G-G-X-G-G-NH(2). The work includes a validation of the concepts used to derive sequence-dependent correction factors for random coil chemical shifts, and a comprehensive tabulation of sequence-dependent correction factors that can be applied for amino acids up to two residues from a given position. This new set of correction factors will have important applications to folded proteins as well as to short, unstructured peptides and unfolded proteins.

Solution structure of the third immunoglobulin domain of the neural cell adhesion molecule N-CAM: can solution studies define the mechanism of homophilic binding?

Homophilic binding of the neural cell adhesion molecule (N-CAM) mediates the calcium-independent cell-cell adhesion that is involved in neuronal development. Two hypotheses have been advanced for the mechanism of homophilic binding. Cell-based experiments have implicated each of the five extracellular immunoglobulin (Ig) domains of N-CAM in the homophilic adhesion interaction, and have predicted that the third domain (Ig III) self-associates. The alternative hypothesis is based on solution observations, which implicate a specific antiparallel interaction between the first two Ig domains (Ig I and Ig II). In order to test these hypotheses, we have determined a high-resolution solution structure of recombinant Ig III (sequence derived from chicken N-CAM) and examined the aggregation behavior of isolated Ig domains in solution. The structure shows that Ig III adopts a canonical Ig fold, in which the beta strands ABED and A'GFCC' form two beta sheets that are linked by a disulfide bond. In contrast to the demonstrated aggregation of Ig III on solid supports, we were unable to demonstrate self-association of Ig III under any of a variety of solution conditions. The structure shows that the surface of Ig III is dominated by two large acidic patches, which may explain our failure to observe self-association in solution. To evaluate the involvement of the Ig I-Ig II interaction in cell-cell adhesion, we designed a point mutation in Ig I (F19S) that proved sufficient to abrogate the Ig I-Ig II interaction seen in solution. However, the introduction of this mutation into full-length N-CAM expressed in COS-7 cells failed to affect N-CAM-mediated cell-cell adhesion. The inability to observe Ig III self-association in solution, combined with the failure of the F19S mutation to affect N-CAM-mediated cell-cell adhesion, suggests that, although solution studies can give important insights into the structures of individual domains, the interactions observed in solution between the domains may not be representative of the interactions that occur on the cell surface.

SANE (Structure Assisted NOE Evaluation): an automated model-based approach for NOE assignment.

A reliable automated approach for assignment of NOESY spectra would allow more rapid determination of protein structures by NMR. In this paper we describe a semi-automated procedure for complete NOESY assignment (SANE, Structure Assisted NOE Evaluation), coupled to an iterative procedure for NMR structure determination where the user is directly involved. Our method is similar to ARIA [Nilges et al. (1997) J. Mol. Biol., 269, 408-422], but is compatible with the molecular dynamics suites AMBER and DYANA. The method is ideal for systems where an initial model or crystal structure is available, but has also been used successfully for ab initio structure determination. Use of this semi-automated iterative approach assists in the identification of errors in the NOE assignments to short-cut the path to an NMR solution structure.

Genomic-scale comparison of sequence- and structure-based methods of function prediction: does structure provide additional insight?

A function annotation method using the sequence-to-structure-to-function paradigm is applied to the identification of all disulfide oxidoreductases in the Saccharomyces cerevisiae genome. The method identifies 27 sequences as potential disulfide oxidoreductases. All previously known thioredoxins, glutaredoxins, and disulfide isomerases are correctly identified. Three of the 27 predictions are probable false-positives. Three novel predictions, which subsequently have been experimentally validated, are presented. Two additional novel predictions suggest a disulfide oxidoreductase regulatory mechanism for two subunits (OST3 and OST6) of the yeast oligosaccharyltransferase complex. Based on homology, this prediction can be extended to a potential tumor suppressor gene, N33, in humans, whose biochemical function was not previously known. Attempts to obtain a folded, active N33 construct to test the prediction were unsuccessful. The results show that structure prediction coupled with biochemically relevant structural motifs is a powerful method for the function annotation of genome sequences and can provide more detailed, robust predictions than function prediction methods that rely on sequence comparison alone.

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.

NMR structural and dynamic characterization of the acid-unfolded state of apomyoglobin provides insights into the early events in protein folding.

Apomyoglobin forms a denatured state under low-salt conditions at pH 2.3. The conformational propensities and polypeptide backbone dynamics of this state have been characterized by NMR. Nearly complete backbone and some side chain resonance assignments have been obtained, using a triple-resonance assignment strategy tailored to low protein concentration (0.2 mM) and poor chemical shift dispersion. An estimate of the population and location of residual secondary structure has been made by examining deviations of (13)C(alpha), (13)CO, and (1)H(alpha) chemical shifts from random coil values, scalar (3)J(HN,H)(alpha) coupling constants and (1)H-(1)H NOEs. Chemical shifts constitute a highly reliable indicator of secondary structural preferences, provided the appropriate random coil chemical shift references are used, but in the case of acid-unfolded apomyoglobin, (3)J(HN,H)(alpha) coupling constants are poor diagnostics of secondary structure formation. Substantial populations of helical structure, in dynamic equilibrium with unfolded states, are formed in regions corresponding to the A and H helices of the folded protein. In addition, the deviation of the chemical shifts from random coil values indicates the presence of helical structure encompassing the D helix and extending into the first turn of the E helix. The polypeptide backbone dynamics of acid-unfolded apomyoglobin have been investigated using reduced spectral density function analysis of (15)N relaxation data. The spectral density J(omega(N)) is particularly sensitive to variations in backbone fluctuations on the picosecond to nanosecond time scale. The central region of the polypeptide spanning the C-terminal half of the E helix, the EF turn, and the F helix behaves as a free-flight random coil chain, but there is evidence from J(omega(N)) of restricted motions on the picosecond to nanosecond time scale in the A and H helix regions where there is a propensity to populate helical secondary structure in the acid-unfolded state. Backbone fluctuations are also restricted in parts of the B and G helices due to formation of local hydrophobic clusters. Regions of restricted backbone flexibility are generally associated with large buried surface area. A significant increase in J(0) is observed for the NH resonances of some residues located in the A and G helices of the folded protein and is associated with fluctuations on a microsecond to millisecond time scale that probably arise from transient contacts between these distant regions of the polypeptide chain. Our results indicate that the equilibrium unfolded state of apomyoglobin formed at pH 2.3 is an excellent model for the events that are expected to occur in the earliest stages of protein folding, providing insights into the regions of the polypeptide that spontaneously undergo local hydrophobic collapse and sample nativelike secondary structure.

Local structural plasticity of the prion protein. Analysis of NMR relaxation dynamics.

A template-assisted conformational change of the cellular prion protein (PrP(C)) from a predominantly helical structure to an amyloid-type structure with a higher proportion of beta-sheet is thought to be the causative factor in prion diseases. Since flexibility of the polypeptide is likely to contribute to the ability of PrP(C) to undergo the conformational change that leads to the infective state, we have undertaken a comprehensive examination of the dynamics of two recombinant Syrian hamster PrP fragments, PrP(29-231) and PrP(90-231), using (15)N NMR relaxation measurements. The molecular motions of these PrP fragments have been studied in solution using (15)N longitudinal (T(1)) and transverse relaxation (T(2)) measurements as well as [(1)H]-(15)N nuclear Overhauser effects (NOE). These data have been analyzed using both reduced spectral density mapping and the Lipari-Szabo model free formalism. The relaxation properties of the common regions of PrP(29-231) and PrP(90-231) are very similar; both have a relatively inflexible globular domain (residues 128-227) with a highly flexible and largely unstructured N-terminal domain. Residues 29-89 of PrP(29-231), which include the copper-binding octarepeat sequences, are also highly flexible. Analysis of the spectral densities at each residue indicates that even within the structured core of PrP(C), a markedly diverse range of motions is observed, consistent with the inherent plasticity of the protein. The central portions of helices B and C form a relatively rigid core, which is stabilized by the presence of an interhelix disulfide bond. Of the remainder of the globular domain, the parts that are not in direct contact with the rigid region, including helix A, are more flexible. Most significantly, slow conformational fluctuations on a millisecond to microsecond time scale are observed for the small beta-sheet. These results are consistent with the hypothesis that the infectious, scrapie form of the protein PrP(Sc) could contain a helical core consisting of helices B and C, similar in structure to the cellular form PrP(C). Our results indicate that residues 90-140, which are required for prion infectivity, are relatively flexible in PrP(C), consistent with a lowered thermodynamic barrier to a template-assisted conformational change to the infectious beta-sheet-rich scrapie isoform.

Two different neurodegenerative diseases caused by proteins with similar structures.

The downstream prion-like protein (doppel, or Dpl) is a paralog of the cellular prion protein, PrP(C). The two proteins have approximately 25% sequence identity, but seem to have distinct physiologic roles. Unlike PrP(C), Dpl does not support prion replication; instead, overexpression of Dpl in the brain seems to cause a completely different neurodegenerative disease. We report the solution structure of a fragment of recombinant mouse Dpl (residues 26-157) containing a globular domain with three helices and a small amount of beta-structure. Overall, the topology of Dpl is very similar to that of PrP(C). Significant differences include a marked kink in one of the helices in Dpl, and a different orientation of the two short beta-strands. Although the two proteins most likely arose through duplication of a single ancestral gene, the relationship is now so distant that only the structures retain similarity; the functions have diversified along with the sequence.

Structure of the PHD zinc finger from human Williams-Beuren syndrome transcription factor.

The PHD (plant homeo domain) is a approximately 50-residue motif found mainly in proteins involved in eukaryotic transcription regulation. The characteristic sequence feature is a conserved Cys(4)-HisCys(3) zinc binding motif. We have determined the solution structure of the PHD motif from the human Williams-Beuren syndrome transcription factor (WSTF) protein. The domain folds into an interleaved zinc finger which binds two Zn(2+) in a similar manner to that of the RING and FYVE domains. The structure reveals a conserved zinc-binding core, together with two variable loops that are likely candidates for interactions between the various PHD domains and their specific ligands.

Solution structure and acetyl-lysine binding activity of the GCN5 bromodomain.

The solution structure of the bromodomain from the human transcriptional coactivator GCN5 has been determined using NMR methods. The structure has a left-handed four-helix bundle topology, with two short additional helices in a long connecting loop. A hydrophobic groove and deep hydrophobic cavity are formed by loops at one end of the molecule. NMR binding experiments show that the cavity forms a specific binding pocket for the acetamide moiety. Peptides containing an N(epsilon)-acetylated lysine residue bind in this pocket with modest affinity (K(D) approximately 0.9 mM); no comparable binding occurs with unacetylated peptides. The GCN5 bromodomain binds the small ligands N(omega)-acetylhistamine and N-methylacetamide, confirming specificity for the alkyl acetamide moiety and showing that the primary element of recognition is simply the sterically unhindered terminal acetamide moiety of an acetylated lysine residue. Additional experiments show that binding is enhanced if the acetyl-lysine residue occurs within the context of a basic peptide and is inhibited by the presence of nearby acidic residues and by the carboxyl group of the free acetyl-lysine amino acid. The binding of the GCN5 bromodomain to acetylated peptides appears to have little additional sequence dependence, although weak interactions with other regions of the peptide are implicated by the binding data. Discrimination between ligands of positive and negative charge is attributed to the presence of several acidic residues located on the loops that form the sides of the binding pocket. Unlike the residues forming the acetamide binding cavity, these acidic side-chains are not conserved in other bromodomain sequences, suggesting that bromodomains might display differences in substrate selectivity and specificity as well as differences in function in vivo.

Random coil chemical shifts in acidic 8 M urea: implementation of random coil shift data in NMRView.

Studies of proteins unfolded in acid or chemical denaturant can help in unraveling events during the earliest phases of protein folding. In order for meaningful comparisons to be made of residual structure in unfolded states, it is necessary to use random coil chemical shifts that are valid for the experimental system under study. We present a set of random coil chemical shifts obtained for model peptides under experimental conditions used in studies of denatured proteins. This new set, together with previously published data sets, has been incorporated into a software interface for NMRView, allowing selection of the random coil data set that fits the experimental conditions best.

Dynamics of the metallo-beta-lactamase from Bacteroides fragilis in the presence and absence of a tight-binding inhibitor.

A significant determinant for the broad substrate specificity of the metallo-beta-lactamases from Bacteroides fragilis and other similar organisms is the presence of a plastic substrate binding site that is nevertheless capable of tight substrate binding in the Michaelis complex. To achieve these two competing ends, the molecule apparently employs a flexible flap that closes over the active site in the presence of substrate. These characteristics imply that dynamic changes are an important component of the mechanism of action of these enzymes. The backbone and tryptophan side chain dynamics of the metallo-beta-lactamase from B. fragilis have been examined using (15)N NMR relaxation measurements. Two states of the protein were examined, in the presence and absence of a tight-binding inhibitor. Relaxation measurements were analyzed by the model-free method. Overall, the metallo-beta-lactamase molecule is rigid and shows little flexibility except in loops. The flexibility of the loop that covers the active site is not unusually great as compared to the other loops of the protein. Local motion on a picosecond time scale was found to be very similar throughout the protein in the presence and absence of the inhibitor, but a significant difference was observed in the motions on a nanosecond time scale (tau(e)). Large-amplitude motions with a time constant of about 1.3 ns were observed for the flexible flap region (residues 45-55) in the absence of the inhibitor. These motions were completely damped out in the presence of the inhibitor. In addition, the motion of a tryptophan side chain at the tip of the beta-hairpin of the flap shows a very significant difference in motion on the ps time scale. These results indicate that the motions of the polypeptide chain in the flap region can be invoked to explain both the wide substrate specificity (the free form has considerable amplitude of motion in this region) and the catalytic efficiency of the metallo-beta-lactamase (the motions are damped out when the inhibitor and by implication a substrate binds in the active site).

Molecular basis for modulation of biological function by alternate splicing of the Wilms' tumor suppressor protein.

Alternate splicing, leading to the insertion of the tripeptide KTS in the linker between the third and fourth C(2)H(2) zinc fingers, changes both the DNA-binding function and the subnuclear localization of the Wilms' tumor suppressor protein (WT1). We have used NMR relaxation experiments to determine the molecular basis for the differing DNA recognition properties of the WT1-KTS and WT1+KTS isoforms. Our results show that the KTS insertion increases the flexibility of the linker between fingers 3 and 4 and abrogates binding of the fourth zinc finger to its cognate site in the DNA major groove. This represents a mechanism whereby a single zinc-finger gene can be used, through alternate splicing, to fulfill different functions in the cell.

Solution structure of the TAZ2 (CH3) domain of the transcriptional adaptor protein CBP.

The TAZ2 (CH3) domain of the transcriptional adapter protein CBP has been implicated in direct functional interactions with numerous cellular transcription factors and viral oncoproteins. The solution structure of the TAZ2 domain of murine CBP has been determined by nuclear magnetic resonance (NMR). The protein adopts a novel helical fold stabilized by three zinc ions, each of which is bound to one histidine and three cysteine ligands in HCCC-type motifs. Each zinc-binding site is formed from the carboxy terminus of an alpha-helix, a short loop, and the amino terminus of the next alpha-helix. A peptide derived from the N-terminal transactivation domain of p53 binds specifically to one face of the TAZ2 domain. The close similarities between the TAZ2 and TAZ1 (CH1 domain of CBP/p300) sequences suggest that both domains will adopt similar three-dimensional structures.