Insights into protein folding from NMR

Residue-Specific Kinetic Insights into the Transition State in Slow Polypeptide Topological Isomerization by NMR Exchange Spectroscopy

The characterization of the transition state is a central issue in biophysical studies of protein folding. NMR is a multiprobe measurement technique that provides residue-specific information. Here, we used exchange spectroscopy to characterize the transition state of the two-state slow topological isomerization of a 27-residue lantibiotic peptide. The exchange kinetic rates varied on a per-residue basis, indicating the reduced kinetic cooperativity of the two-state exchange, as well as the previously observed reduced thermodynamic cooperativity. Furthermore, temperature-dependent measurements revealed large variations in the activation enthalpy and entropy terms among residues. Interestingly, we found a linear relationship between the logarithm of the equilibrium constants and that of the exchange rates. Because the data points are derived from amino acid residues in one polypeptide chain, we refer to the linear relationship as the residue-based linear free energy relationship (rbLFER). The rbLFER offers information about the transition state of the two-state exchange.

The concept of protein folding/unfolding and its impacts on human health

Proteins have evolved in specific 3D structures and play different functions in cells and determine various reactions and pathways. The newly synthesized amino acid chains once depart ribosome must crumple into three-dimensional structures so can be biologically active. This process of protein that makes a functional molecule is called protein folding. The protein folding is both a biological and a physicochemical process that depends on the sequence of it. In fact, this process occurs more complicated and in some cases and in exposure to some molecules like glucose (glycation), mistaken folding leads to amyloid structures and fatal disorders called conformational diseases. Such conditions are detected by the quality control system of the cell and these abnormal proteins undergo renovation or degradation. This scenario takes place by the chaperones, chaperonins, and Ubiquitin-proteasome complex. Understanding of protein folding mechanisms from different views including experimental and computational approaches has revealed some intermediate ensembles such as molten globule and has been subjected to biophysical and molecular biology attempts to know more about prevalent conformational diseases.

Chemically Denatured Structures of Porcine Pepsin using Small-Angle X-ray Scattering

Porcine pepsin is a gastric aspartic proteinase that reportedly plays a pivotal role in the digestive process of many vertebrates. We have investigated the three-dimensional (3D) structure and conformational transition of porcine pepsin in solution over a wide range of denaturant urea concentrations (0–10 M) using Raman spectroscopy and small-angle X-ray scattering. Furthermore, 3D GASBOR ab initio structural models, which provide an adequate conformational description of pepsin under varying denatured conditions, were successfully constructed. It was shown that pepsin molecules retain native conformation at 0–5 M urea, undergo partial denaturation at 6 M urea, and display a strongly unfolded conformation at 7–10 M urea. According to the resulting GASBOR solution models, we identified an intermediate pepsin conformation that was dominant during the early stage of denaturation. We believe that the structural evidence presented here provides useful insights into the relationship between enzymatic activity and conformation of porcine pepsin at different states of denaturation.

Characterisation of denatured states of sensory rhodopsin II by solution-state NMR

Sensory rhodopsin II (pSRII), a retinal-binding photophobic receptor from Natronomonas pharaonis, is a novel model system for membrane protein folding studies. Recently, the SDS-denatured states and the kinetics for reversible unfolding of pSRII have been investigated, opening the door to the first detailed characterisation of denatured states of a membrane protein by solution-state nuclear magnetic resonance (NMR) using uniformly 15N-labelled pSRII. SDS denaturation and acid denaturation of pSRII both lead to fraying of helix ends but otherwise small structural changes in the transmembrane domain, consistent with little changes in secondary structure and disruption of the retinal-binding pocket and tertiary structure. Widespread changes in the backbone amide dynamics are detected in the form of line broadening, indicative of μs-to-ms timescale conformational exchange in the transmembrane region. Detailed analysis of chemical shift and intensity changes lead to high-resolution molecular insights on structural and dynamics changes in SDS- and acid-denatured pSRII, thus highlighting differences in the unfolding pathways under the two different denaturing conditions. These results will form the foundation for furthering our understanding on the folding and unfolding pathways of retinal-binding proteins and membrane proteins in general, and also for investigating the importance of ligand-binding in the folding pathways of other ligand-binding membrane proteins, such as GPCRs.

Experimental approach to study the effect of mutations on the protein folding pathway

Is it possible to compare the physicochemical properties of a wild-type protein and its mutant form under the same conditions? Provided the mutation has destabilized the protein, it may be more correct to compare the mutant protein under native conditions to the wild-type protein destabilized with a small amount of the denaturant. In general, is it appropriate to compare the properties of proteins destabilized by different treatments: mutations, pH, temperature, and denaturants like urea? These issues have compelled us to search for methods and ways of presentation of experimental results that would allow a comparison of mutant forms of proteins under different conditions and lead to conclusions on the effect of mutations on the protein folding/unfolding pathway. We have studied equilibrium unfolding of wild-type bovine carbonic anhydrase II (BCA II) and its six mutant forms using different urea concentrations. BCA II has been already studied in detail and is a good model object for validating new techniques. In this case, time-resolved fluorescence spectroscopy was chosen as the basic research method. The main features of this experimental method allowed us to compare different stages of unfolding of studied proteins and prove experimentally that a single substitution of the amino acid in three mutant forms of BCA II affected the native state of the protein but did not change its unfolding pathway. On the contrary, the inserted disulfide bridge in three other mutant forms of BCA II affected the protein unfolding pathway. An important result of this research is that we have validated the new approach allowing investigation of the effect of mutations on the folding of globular proteins, because in this way it is possible to compare proteins in the same structural states rather than under identical conditions.

Impact of Turn Propensity on the Folding Rates of Z34C Protein: Implications for the Folding of Helix-Turn-Helix Motif

The rate-limiting step for the folding of the helix-turn-helix (HTH) protein, Z34C, involves β-turn region ²⁰DPNL²³. This reverse turn has been observed to be part of the transition state in the folding process for Z34C, influencing its folding rates. Molecular dynamics simulations were performed on this turn peptide and its two mutants, D20A and P21A, to study turn formation using GROMOS54A7 force field. We find that this region has a turn propensity of its own, and the highest turn propensity is observed for the wild-type, which correlates well with available experimental results. We also find that a slight unfavorable change in ΔG _{turn folding} causes a drastic change in the folding rates of HTH motif and a mechanistic interpretation is given. Implications of these observations for the folding of the HTH protein Z34C are discussed.

NMR characterization of sodium carboxymethyl cellulose 2: Chemical shift assignment and conformation analysis of substituent groups

The chemical shifts of the substituent groups of sodium carboxymethyl cellulose (CMC) were assigned by examining a series of CMC samples with different degrees of substitution. Comparative analysis of the (1)H-(13)C heteronuclear single quantum coherence (HSQC) and heteronuclear multiple bond correlation (HMBC) spectra allowed the complete assignment of the substituent groups at the 2-, 3-, and 6-positions of the seven substituted monomers comprising the CMC chains, namely, 2-mono-, 3-mono-, 6-mono-, 2,3-di-, 2,6-di-, 3,6-di-, and 2,3,6-tri-substituted anhydroglucose units (AGUs). In addition, the mole fractions of the monomers were determined by lineshape analysis of the carbonyl carbon resonances. The comparison between the chemical shifts of the substituents revealed strong interactions between 2- and 3-substituents in the same AGU, and showed that the steric hindrance by a substituent at the 2- or 3-position suppresses subsequent substitution at the adjacent position.

Revisiting and re-engineering the classical zinc finger peptide: consensus peptide-1 (CP-1)

Zinc plays key structural and catalytic roles in biology. Structural zinc sites are often referred to as zinc finger (ZF) sites, and the classical ZF contains a Cys2His2 motif that is involved in coordinating Zn(II). An optimized Cys2His2 ZF, named consensus peptide 1 (CP-1), was identified more than 20 years ago using a limited set of sequenced proteins. We have reexamined the CP-1 sequence, using our current, much larger database of sequenced proteins that have been identified from high-throughput sequencing methods, and found the sequence to be largely unchanged. The CCHH ligand set of CP-1 was then altered to a CAHH motif to impart hydrolytic activity. This ligand set mimics the His2Cys ligand set of peptide deformylase (PDF), a hydrolytically active M(II)-centered (M = Zn or Fe) protein. The resultant peptide [CP-1(CAHH)] was evaluated for its ability to coordinate Zn(II) and Co(II) ions, adopt secondary structure, and promote hydrolysis. CP-1(CAHH) was found to coordinate Co(II) and Zn(II) and a pentacoordinate geometry for Co(II)-CP-1(CAHH) was implicated from UV-vis data. This suggests a His2Cys(H2O)2 environment at the metal center. The Zn(II)-bound CP-1(CAHH) was shown to adopt partial secondary structure by 1-D (1)H NMR spectroscopy. Both Zn(II)-CP-1(CAHH) and Co(II)-CP-1(CAHH) show good hydrolytic activity toward the test substrate 4-nitrophenyl acetate, exhibiting faster rates than most active synthetic Zn(II) complexes.

Real-time protein NMR spectroscopy and investigation of assisted protein folding

BACKGROUND

During protein-folding reactions toward the native structure, short-lived intermediate states can be populated. Such intermediates expose hydrophobic patches and can self-associate leading to non-productive protein misfolding. A major focus of current research is the characterization of short-lived intermediates and how molecular chaperones enable productive folding. Real-time NMR spectroscopy, together with the development of advanced methods, is reviewed here and the potential these methods have to characterize intermediate states as well as interactions with molecular chaperone proteins at single-residue resolution is highlighted.

SCOPE OF REVIEW

Various chaperone interactions can guide the protein-folding reaction and thus are important for protein structure formation, stability, and activity of their substrates. Chaperone-assisted protein folding, characterization of intermediates, and their molecular interactions using real-time NMR spectroscopy will be discussed. Additionally, recent advances in NMR methods employed for characterization of high-energy intermediates will be discussed.

MAJOR CONCLUSIONS

Real-time NMR combines high resolution with kinetic information of protein reactions, which can be employed not only for protein-folding studies and the characterization of folding intermediates but also to investigate the molecular mechanisms of assisted protein folding.

GENERAL SIGNIFICANCE

Real-time NMR spectroscopy remains an effective tool to reveal structural details about the interaction between chaperones and transient intermediates. Methodologically, it provides in-depth understanding of how kinetic intermediates and their thermodynamics contribute to the protein-folding reaction. This review summarizes the most recent advances in this field. This article is part of a Special Issue titled Proline-directed Foldases: Cell Signaling Catalysts and Drug Targets.

High-resolution characterization of intrinsic disorder in proteins: expanding the suite of (13)C-detected NMR spectroscopy experiments to determine key observables

Order in disorder: The characterization of intrinsically disordered proteins by NMR spectroscopy is a necessity on the one hand and a continuous challenge on the other. We propose two experiments that provide diagnostic parameters to monitor the degree of unfolding of a polypeptide. The test was performed on the yeast Cox17 protein, known to gain its function through maturation from an intrinsically disordered state (see figure).

Mapping the conformational stability of maltose binding protein at the residue scale using nuclear magnetic resonance hydrogen exchange experiments

Being able to differentiate local fluctuations from global folding-unfolding dynamics of a protein is of major interest for improving our understanding of structure-function determinants. The maltose binding protein (MBP), a protein that belongs to the maltose transport system, has a structure composed of two globular domains separated by a rigid-body "hinge bending". Here we determined, by using hydrogen exchange (HX) nuclear magnetic resonance experiments, the apparent stabilization free energies of 101 residues of MBP bound to β-cyclodextrin (MBP-βCD) under native conditions. We observed that the last helix of MBP (helix α14) has a lower protection factor than the rest of the protein. Further, HX experiments were performed using guanidine hydrochloride under subdenaturing conditions to discriminate between local fluctuations and global unfolding events and to determine the MBP-βCD energy landscape. The results show that helix α4 and a part of helices α5 and α6 are clearly grouped into a subdenaturing folding unit and represent a partially folded intermediate under native conditions. In addition, we observed that amide protons located in the hinge between the two globular domains share similar ΔG(gu)(app) and m values and should unfold simultaneously. These observations provide new points of view for improving our understanding of the thermodynamic stability and the mechanisms that drive folding-unfolding dynamics of proteins.

Folding of an all-helical Greek-key protein monitored by quenched-flow hydrogen-deuterium exchange and NMR spectroscopy

To advance our understanding of the protein folding process, we use stopped-flow far-ultraviolet (far-UV) circular dichroism and quenched-flow hydrogen-deuterium exchange coupled with nuclear magnetic resonance (NMR) spectroscopy to monitor the formation of hydrogen-bonded secondary structure in the C-terminal domain of the Fas-associated death domain (Fadd-DD). The death domain superfamily fold consists of six α-helices arranged in a Greek-key topology, which is shared by the all-β-sheet immunoglobulin and mixed α/β-plait superfamilies. Fadd-DD is selected as our model death domain protein system because the structure of this protein has been solved by NMR spectroscopy, and both thermodynamic and kinetic analysis indicate it to be a stable, monomeric protein with a rapidly formed hydrophobic core. Stopped-flow far-UV circular dichroism spectroscopy revealed that the folding process was monophasic and the rate is 23.4 s(-1). Twenty-two amide hydrogens in the backbone of the helices and two in the backbone of the loops were monitored, and the folding of all six helices was determined to be monophasic with rate constants between 19 and 22 s(-1). These results indicate that the formation of secondary structure is largely cooperative and concomitant with the hydrophobic collapse. This study also provides unprecedented insight into the formation of secondary structure within the highly populated Greek-key fold more generally. Additional insights are gained by calculating the exchange rates of 23 residues from equilibrium hydrogen-deuterium exchange experiments. The majority of protected amide protons are found on helices 2, 4, and 5, which make up core structural elements of the Greek-key topology.

Macromolecular crowding fails to fold a globular protein in cells

Proteins perform their functions in cells where macromolecular solutes reach concentrations of >300 g/L and occupy >30% of the volume. The volume excluded by these macromolecules stabilizes globular proteins because the native state occupies less space than the denatured state. Theory predicts that crowding can increase the ratio of folded to unfolded protein by a factor of 100, amounting to 3 kcal/mol of stabilization at room temperature. We tested the idea that volume exclusion dominates the crowding effect in cells using a variant of protein L, a 7 kDa globular protein with seven lysine residues replaced by glutamic acids; 84% of the variant molecules populate the denatured state in dilute buffer at room temperature, compared with 0.1% for the wild-type protein. We then used in-cell NMR spectroscopy to show that the cytoplasm of Escherichia coli does not overcome even this modest (∼1 kcal/mol) free-energy deficit. The data are consistent with the idea that nonspecific interactions between cytoplasmic components can overcome the excluded-volume effect. Evidence for these interactions is provided by the observations that adding simple salts folds the variant in dilute solution but increasing the salt concentration inside E. coli does not fold the protein. Our data are consistent with the results of other studies of protein stability in cells and suggest that stabilizing excluded-volume effects, which must be present under crowded conditions, can be ameliorated by nonspecific interactions between cytoplasmic components.

The folding of disulfide-rich proteins

The articles in this forum issue describe various aspects of the folding of disulfide-rich proteins. They include review articles using proteins such as bovine pancreatic trypsin inhibitor as models to highlight the range of folding pathways seen in disulfide-rich proteins, along with a detailed analysis of the methods used to study them. Following two comprehensive reviews on the methods and applications of protein folding, three original articles in this issue focus on two specific classes of disulfide-rich proteins that have applications in drug design and development, namely cyclotides and conotoxins. Cyclotides are head-to-tail cyclic and disulfide-rich proteins from plants and function as a defense against insect attack. Conotoxins are the disulfide-rich components of the venom of marine cone snails that is used to capture prey. These research articles report on factors that modulate protein folding pathways in these molecules and determine the outcomes of protein folding, that is, yield and heterogeneity of products. Finally, the issue concludes with a comprehensive review on a different type of disulfide bond, namely those that have a functional rather than structural role in proteins, with a particular focus on allosteric disulfide bonds that modify protein function.

Distinct folding pathways of two homologous disulfide proteins: bovine pancreatic trypsin inhibitor and tick anticoagulant peptide

The folding pathways of disulfide proteins vary substantially (Arolas et al., Trends Biochem Sci 31: 292-301, 2006). The diversity is mainly manifested by (a) the extent of heterogeneity of folding intermediates, (b) the extent of presence of native-like intermediates, and (c) the variation of folding kinetics. Even among structurally similar proteins, the difference can be enormous. This is demonstrated in this concise review with two structurally homologous kunitz-type protease inhibitors, bovine pancreatic trypsin inhibitor and tick anticoagulant peptide, as well as a group of cystine knot proteins. The diversity of their folding mechanisms is illustrated with two different folding techniques: (a) the conventional method of disulfide oxidation (oxidative folding), and (b) the novel method of disulfide scrambling (Chang, J Biol Chem 277: 120-126, 2002). This review also highlights the convergence of folding models concluded form the conventional conformational folding and those obtained by oxidative folding.

An NMR database for simulations of membrane dynamics

Computational methods are powerful in capturing the results of experimental studies in terms of force fields that both explain and predict biological structures. Validation of molecular simulations requires comparison with experimental data to test and confirm computational predictions. Here we report a comprehensive database of NMR results for membrane phospholipids with interpretations intended to be accessible by non-NMR specialists. Experimental ¹³C-¹H and ²H NMR segmental order parameters (S(CH) or S(CD)) and spin-lattice (Zeeman) relaxation times (T(1Z)) are summarized in convenient tabular form for various saturated, unsaturated, and biological membrane phospholipids. Segmental order parameters give direct information about bilayer structural properties, including the area per lipid and volumetric hydrocarbon thickness. In addition, relaxation rates provide complementary information about molecular dynamics. Particular attention is paid to the magnetic field dependence (frequency dispersion) of the NMR relaxation rates in terms of various simplified power laws. Model-free reduction of the T(1Z) studies in terms of a power-law formalism shows that the relaxation rates for saturated phosphatidylcholines follow a single frequency-dispersive trend within the MHz regime. We show how analytical models can guide the continued development of atomistic and coarse-grained force fields. Our interpretation suggests that lipid diffusion and collective order fluctuations are implicitly governed by the viscoelastic nature of the liquid-crystalline ensemble. Collective bilayer excitations are emergent over mesoscopic length scales that fall between the molecular and bilayer dimensions, and are important for lipid organization and lipid-protein interactions. Future conceptual advances and theoretical reductions will foster understanding of biomembrane structural dynamics through a synergy of NMR measurements and molecular simulations.

Dynamics of protein folding: probing the kinetic network of folding-unfolding transitions with experiment and theory

The problem of spontaneous folding of amino acid chains into highly organized, biologically functional three-dimensional protein structures continues to challenge the modern science. Understanding how proteins fold requires characterization of the underlying energy landscapes as well as the dynamics of the polypeptide chains in all stages of the folding process. In recent years, important advances toward these goals have been achieved owing to the rapidly growing interdisciplinary interest and significant progress in both experimental techniques and theoretical methods. Improvements in the experimental time resolution led to determination of the timescales of the important elementary events in folding, such as formation of secondary structure and tertiary contacts. Sensitive single molecule methods made possible probing the distributions of the unfolded and folded states and following the folding reaction of individual protein molecules. Discovery of proteins that fold in microseconds opened the possibility of atomic-level theoretical simulations of folding and their direct comparisons with experimental data, as well as of direct experimental observation of the barrier-less folding transition. The ultra-fast folding also brought new questions, concerning the intrinsic limits of the folding rates and experimental signatures of barrier-less "downhill" folding. These problems will require novel approaches for even more detailed experimental investigations of the folding dynamics as well as for the analysis of the folding kinetic data. For theoretical simulations of folding, a main challenge is how to extract the relevant information from overwhelmingly detailed atomistic trajectories. New theoretical methods have been devised to allow a systematic approach towards a quantitative analysis of the kinetic network of folding-unfolding transitions between various configuration states of a protein, revealing the transition states and the associated folding pathways at multiple levels, from atomistic to coarse-grained representations. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Viscoelastic modeling with interfacial slip of a protein monolayer electrode-adsorbed on an acoustic wave biosensor

Transverse-shear mode acoustic wave devices have been used as real-time, label-free detectors of conformational shifts in biomolecules on surfaces. However, material changes in the biochemical monolayers and coupling between the substrate and the surrounding liquid make it difficult to isolate the desired signal, so an understanding of these phenomena is required. An important step in this understanding is knowledge of the material properties of the linker layer that attaches a biochemically selective molecule to the gold surface, in our case, neutravidin. With the goal of obtaining material properties for a neutravidin monolayer, for use in future studies, neutravidin adsorption to the gold surface of an acoustic wave biosensor is described as a viscoelastic monolayer using one-dimensional modeling. Neutravidin is described as forming hydrated, viscoelastic monolayers, and slip is allowed at all interfaces. An impedance model is numerically fit to experimental values using a two-parameter minimization algorithm and values for the shear modulus of the neutravidin monolayer, in agreement with literature values for similar proteins, are obtained. Slip is found on the electrode surface prior to neutravidin adsorption. These results will be used for future modeling studies involving this protein as a linker protein.

Biomolecular Structure and Modeling: Problem and Application Perspective

The experimental progress described in the previous chapter has been accompanied by an increasing desire to relate the complex three-dimensional (3D) shapes of biomolecules to their biological functions and interactions with other molecular systems. Structural biology, computational biology, genomics, proteomics, bioinformatics, chemoinformatics, and others are natural partner disciplines in such endeavors.

NMR investigations on residue level unfolding thermodynamics in DLC8 dimer by temperature dependent native state hydrogen exchange

Understanding protein stability at residue level detail in the native state ensemble of a protein is crucial to understanding its biological function. At the same time, deriving thermodynamic parameters using conventional spectroscopic and calorimetric techniques remains a major challenge for some proteins due to protein aggregation and irreversibility of denaturation at higher temperature values. In this regard, we describe here the NMR investigations on the conformational stabilities and related thermodynamic parameters such as local unfolding enthalpies, heat capacities and transition midpoints in DLC8 dimer, by using temperature dependent native state hydrogen exchange; this protein aggregates at high (>65 degrees C) temperatures. The stability (free energy) of the native state was found to vary substantially with temperature at every residue. Significant differences were found in the thermodynamic parameters at individual residue sites indicating that the local environments in the protein structure would respond differently to external perturbations; this reflects on plasticity differences in different regions of the protein. Further, comparison of this data with similar data obtained from GdnHCl dependent native state hydrogen exchange indicated many similarities at residue level, suggesting that local unfolding transitions may be similar in both the cases. This has implications for the folding/unfolding mechanisms of the protein.

Non-natural amino acids as modulating agents of the conformational space of model glycopeptides

The synthesis and conformational analysis in aqueous solution of different alpha-methyl-alpha-amino acid diamides, derived from serine, threonine, beta-hydroxycyclobutane-alpha-amino acids, and their corresponding model beta-O-glucopeptides, are reported. The study reveals that the presence of an alpha-methyl group forces the model peptides to adopt helix-like conformations. These folded conformations are especially significant for cyclobutane derivatives. Interestingly, this feature was also observed in the corresponding model glucopeptides, thus indicating that the alpha-methyl group and not the beta-O-glucosylation process largely determines the conformational preference of the backbone in these structures. On the other hand, atypical conformations of the glycosidic linkage were experimentally determined. Therefore, when a methyl group was located at the Cbeta atom with an R configuration, the glycosidic linkage was rather rigid. Nevertheless, when the S configuration was displayed, a significant degree of flexibility was observed for the glycosidic linkage, thus showing both alternate and eclipsed conformations of the psi(s) dihedral angle. In addition, some derivatives exhibited an unusual value for the phi(s) angle, which was far from a value of -60 degrees expected for a conventional beta-O-glycosidic linkage. In this sense, the different conformations exhibited by these molecules could be a useful tool in obtaining systems with conformational preferences "à la carte".

Residue-wise conformational stability of DLC8 dimer from native-state hydrogen exchange

Dynein light chain (DLC8) is the smallest subunit of the dynein motor complex, which is known to act as a cargo adaptor in intracellular trafficking. The protein exists as a pure dimer at physiological pH and a completely folded monomer below pH 4. Here, we have determined the energy landscape of the dimeric protein using a combination of optical techniques and native-state hydrogen exchange of amide groups, the former giving the global features and the latter yielding the residue level details. The data indicated the presence of intermediates along the equilibrium unfolding transition. The hydrogen exchange data suggested that the molecule has differential stability in its various segments. We deduce from the free energy data that the antiparallel beta-sheets (beta4 and beta5) that form the hydrophobic core of the protein and the alpha2 helix, all of which are highly protected with regard to hydrogen exchange, contribute significantly to the initial step of the protein folding mechanism. Denaturant-dependent hydrogen exchange indicated further that some amides exchange via local fluctuations, whereas there are others which exchange via global unfolding events. Implications of these to cargo adaptability of the dimer are discussed.

Pockets of short-range transient order and restricted topological heterogeneity in the guanidine-denatured state ensemble of GED of dynamin

The nature and variety in the denatured state of a protein, a non-native state under a given set of conditions, has been a subject of intense debate. Here, using multidimensional NMR, we have characterized the 6 M Gdn-HCl-denatured state of GED, the assembly domain of dynamin. Even under such strongly denaturing conditions, we detected the presence of conformations in slow exchange on the NMR chemical shift time scale. Although the GED oligomer as well as the SDS-denatured monomeric GED were seen to be predominantly helical [Chugh et al. (2006) FEBS J. 273, 388-397], the 6 M Gdn-HCl-denatured GED has largely beta-structural preferences. However, against such a background, we could detect the presence of a population with a short helical stretch (Arg42-Ile47) in the ensemble. The 1H-1H NOEs suggested presence of pockets of transient short-range order along the chain. Put together these segments may lead to a rather small number of interconverting topologically distinguishable ensembles. Spectral density analysis of 15N relaxation rates and {1H}-15N NOE, measured at 600 and 800 MHz, and comparison of J(0) with hydrophobic patches calculated using AABUF approach, indicated presence of four domains of slow motions. These coincided to a large extent with those showing significant Rex. Additionally, a proline residue in the connection between two of these domains seems to cause a fast hinge motion. These observations help enhance our understanding of protein denatured states, and of folding concepts, in general.

Synthesis and structural model of an alpha(2,6)-sialyl-t glycosylated MUC1 eicosapeptide under physiological conditions

To study the effect of O-glycosylation on the conformational propensities of a peptide backbone, a 20-residue peptide (GSTAPPAHGVTSAPDTRPAP) representing the full length tandem repeat sequence of the human mucin MUC1 and its analogue glycosylated with the (2,6)-sialyl-T antigen on Thr11, were prepared and investigated by NMR and molecular modeling. The peptides contain both the GVTSAP sequence, which is an effective substrate for GalNAc transferases, and the PDTRP fragment, a known epitope recognized by several anti-MUC1 monoclonal antibodies. It has been shown that glycosylation of threonine in the GVTSAP sequence is a prerequisite for subsequent glycosylation of the serine at GVTSAP. Furthermore, carbohydrates serve as additional epitopes for MUC1 antibodies. Investigation of the solution structure of the sialyl-T glycoeicosapeptide in a H(2)O/D(2)O mixture (9:1) under physiological conditions (25 degrees C and pH 6.5) revealed that the attachment of the saccharide side-chain affects the conformational equilibrium of the peptide backbone near the glycosylated Thr11 residue. For the GVTSA region, an extended, rod-like secondary structure was found by restrained molecular dynamics simulation. The APDTR region formed a turn structure which is more flexibly organized. Taken together, the joined sequence GVTSAPDTR represents the largest structural model of MUC1 derived glycopeptides analyzed so far.

Investigation of the Structural Stability of hUBF HMG Box 5 by Native-State Hydrogen Exchange†

HMG box 5 of human upstream binding factor (hUBF) consists of three alpha-helices arranged in an L-shape with a hydrophobic core embraced by these helices and stabilized by extensive hydrophobic interactions between nonpolar residues around the core. The GdmCl-induced equilibrium unfolding transition of HMG box 5 of hUBF was monitored by both circular dichroism (CD) and fluorescence spectra. A cooperative two-state unfolding process was observed. The unfolding free energy, DeltaGU(D2O), and the cooperativity of the unfolding reaction, m, are 4.6 +/- 0.16 kcal x mol-1 and 1.62 +/- 0.06 kcal x mol-1 x M-1, respectively. Native-state hydrogen exchange (NHX) experiments under EX2 conditions were performed. NHX results clearly show that the hydrophobic core among the three helices is a slow-exchange core. The three helices would not contribute equally to the stability of the native protein. Helix 3 appears to contribute the least to the stability. The NHX data have also allowed the local, subglobal, and global unfolding structures of hUBF HMG box 5 to be dissected, and common global and subglobal unfolding units were successfully detected.

Probing nascent structures in peptides using natural abundance 13C NMR relaxation and reduced spectral density mapping

The main chain motional properties for a series of peptides that appear to have preferred conformations in solution have been systematically studied using solution-state nuclear magnetic resonance spectroscopy. The series of peptides were derived from the N-termini of pro-inflammatory chemokine proteins and HoxB1, a transcriptional regulator. As an unstructured control, a ten residue peptide was designed, synthesized, and found to be minimally structured from solution NMR data. The dynamic properties of the main chain for the peptides were assessed through longitudinal and transverse main chain (13)Calpha relaxation rates and the heteronuclear nuclear Overhauser effect. Motional parameters were interpreted using reduced spectral density mapping and compared with those derived from an extended Lipari-Szabo model in which the rotational correlation time was calculated for each main chain site of the peptide. Comparison of spectral density and Lipari-Szabo analyses for the peptides to those of the unstructured control peptide reveals significant differences in the dynamic behavior of the peptides. The amplitude of picosecond to nanosecond timescale motions for the main chain is observed to decrease for all of the chemokine peptides and HoxB1 over the regions that show partial structure at low temperatures. Comparatively, changes in picosecond to nanosecond timescale motions for the unstructured control peptide show no correlation with sequence position. These results indicate that there are distinguishable low temperature motional differences between an intrinsically unstructured peptide and peptides that have an inherent propensity to structure.

The folding nucleus of the insulin superfamily: a flexible peptide model foreshadows the native state

Oxidative folding of insulin-like growth factor I (IGF-I) and single-chain insulin analogs proceeds via one- and two-disulfide intermediates. A predominant one-disulfide intermediate in each case contains the canonical A20-B19 disulfide bridge (cystines 18-61 in IGF-I and 19-85 in human proinsulin). Here, we describe a disulfide-linked peptide model of this on-pathway intermediate. One peptide fragment (19 amino acids) spans IGF-I residues 7-25 (canonical positions B8-B26 in the insulin superfamily); the other (18 amino acids) spans IGF-I residues 53-70 (positions A12-A21 and D1-D8). Containing only half of the IGF-I sequence, the disulfide-linked polypeptide (designated IGF-p) is not well ordered. Nascent helical elements corresponding to native alpha-helices are nonetheless observed at 4 degrees C. Furthermore, (13)C-edited nuclear Overhauser effects establish transient formation of a native-like partial core; no non-native nuclear Overhauser effects are observed. Together, these observations suggest that early events in the folding of insulin-related polypeptides are nucleated by a native-like molten subdomain containing Cys(A20) and Cys(B19). We propose that nascent interactions within this subdomain orient the A20 and B19 thiolates for disulfide bond formation and stabilize the one-disulfide intermediate once formed. Substitutions in the corresponding region of insulin are associated with inefficient chain combination and impaired biosynthetic expression. The intrinsic conformational propensities of a flexible disulfide-linked peptide thus define a folding nucleus, foreshadowing the structure of the native state.

Analysis of the kinetics of folding of proteins and peptides using circular dichroism

Circular dichroism (CD) is a useful spectroscopic technique for studying the secondary structure, folding and binding properties of proteins. This protocol covers how to use the intrinsic circular dichroic properties of proteins to follow their folding and unfolding as a function of time. Included are methods of obtaining data and for analyzing the folding and unfolding data to determine the rate constants and the order of the folding and unfolding reactions. The protocol focuses on the use of CD to follow folding when it is relatively slow, on the order of minutes to days. The methods for analyzing the data, however, can also be applied to data collected with a CD machine equipped with stopped-flow accessories in the range of milliseconds to seconds and folding analyzed by other spectroscopic methods including changes in absorption or fluorescence spectra as a function of time.

Describing partially unfolded states of proteins from sparse NMR data

Proteins involved in signal transduction can usually be present in two states: an inactive and an active (signaling) state. In the case of photoreceptors such as PYP, it has been shown that the signaling state has a large degree of structural and dynamic disorder. Conventional structural NMR approaches present difficulties in describing such partially unfolded states. Owing to the disordered dynamical and transient nature of such states classical NOE-based information, when present, is sparse. Chemical shift changes upon partial unfolding can, however, be easily monitored from HSQC spectra. We show here that such states can be modeled by defining native-like inter-residue contacts for those residues that do not shift significantly upon partial unfolding. The feasibility of this approach is demonstrated using lysozyme as a test case and applied to model the partially unfolded signaling state (pB) of a truncated form of the photoactive yellow protein for which a "classical" NOE-based structure is available for validation. This approach should be generally applicable to systems in which part of the structure remains in a well-defined native-like conformation.

cDNA cloning and functional expression of jerdostatin, a novel RTS-disintegrin from Trimeresurus jerdonii and a specific antagonist of the alpha1beta1 integrin

Jerdostatin represents a novel RTS-containing short disintegrin cloned by reverse transcriptase-PCR from the venom gland mRNA of the Chinese Jerdons pit viper Trimeresurus jerdonii. The jerdostatins precursor cDNA contained a 333-bp open reading frame encoding a signal peptide, a pre-peptide, and a 43-amino acid disintegrin domain, whose amino acid sequence displayed 80% identity with that of the KTS-disintegrins obtustatin and viperistatin. The jerdostatin cDNA structure represents the first complete open reading frame of a short disintegrin and points to the emergence of jerdostatin from a short-coding gene. The different residues between jerdostatin and obtustatin/viperistatin are segregated within the integrin-recognition loop and the C-terminal tail. Native jerdostatin (r-jerdostatin-R21) and a R21K mutant (r-jerdostatin-K21) were produced in Escherichia coli. In each case, two conformers were isolated. One-dimensional (1)H NMR showed that conformers 1 and 2 of r-jerdostatin-R21 represent, respectively, well folded and unfolded proteins. The two conformers of the wild-type and the R21K mutant inhibited the adhesion of alpha(1)-K562 cells to collagen IV with IC(50) values of 180 and 703 nm, respectively. The IC(50) values of conformers 2 of r-jerdostatin-R21 and r-jerdostatin-K21 were, respectively, 5.95 and 12.5 microm. Neither r-jerdostatin-R21 nor r-jerdostatin-K21 showed inhibitory activity toward other integrins, including alpha(IIb)beta(3), alpha(v)beta(3), alpha(2)beta(1), alpha(5)beta(1), alpha(4)beta(1), alpha(6)beta(1), and alpha(9)beta(1) up to a concentration of 24 mum. Although the RTS motif appears to be more potent than KTS inhibiting the alpha(1)beta(1) integrin, r-jerdostatin-R21 is less active than the KTS-disintegrins, strongly suggesting that substitutions outside the integrin-binding motif and/or C-terminal proteolytic processing are responsible for the decreased inhibitory activity.

Dynamics in the Unfolded State of β2-microglobulin Studied by NMR

Many proteins form amyloid-like fibrils in vitro under conditions that favour the population of partially folded conformations or denatured state ensembles. Characterising the structural and dynamic properties of these states is crucial towards understanding the mechanisms of self-assembly in amyloidosis. The aggregation of beta2-microglobulin (beta2m) into amyloid fibrils in vivo occurs in the condition known as dialysis-related amyloidosis (DRA) and the protein has been shown to form amyloid-like fibrils under acidic conditions in vitro. We have used a number of 1H-15N nuclear magnetic resonance (NMR) experiments in conjunction with site-directed mutagenesis to study the acid-unfolded state of beta2m. 15N NMR transverse relaxation experiments reveal that the acid-denatured ensemble, although predominantly unfolded at the N and C termini, contains substantial non-native structure in the central region of the polypeptide chain, stabilised by long-range interactions between aromatic residues and by the single disulphide bond. Relaxation dispersion studies indicate that the acid-unfolded ensemble involves two or more distinct species in conformational equilibrium on the micro- to millisecond time-scale. One of these species appears to be hydrophobically collapsed, as mutations in an aromatic-rich region of the protein, including residues that are solvent-exposed in the native protein, disrupt this structure and cause a consequent decrease in the population of this conformer. Thus, acid-unfolded beta2m consists of a heterogeneous ensemble of rapidly fluctuating species, some of which contain stable, non-native hydrophobic clusters. Given that amyloid assembly of beta2m proceeds with lag kinetics under the conditions of this study, a rarely populated species such as a conformer with non-native aromatic clustering could be key to the initiation of amyloidosis.

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.

Five Hierarchical Levels of Sequence-Structure Correlation in Proteins

This article reviews recent work towards modelling protein folding pathways using a bioinformatics approach. Statistical models have been developed for sequence-structure correlations in proteins at five levels of structural complexity: (i) short motifs; (ii) extended motifs; (iii) nonlocal pairs of motifs; (iv) 3-dimensional arrangements of multiple motifs; and (v) global structural homology. We review statistical models, including sequence profiles, hidden Markov models (HMMs) and interaction potentials, for the first four levels of structural detail. The I-sites (folding Initiation sites) Library models short local structure motifs. Each succeeding level has a statistical model, as follows: HMMSTR (HMM for STRucture) is an HMM for extended motifs; HMMSTR-CM (Contact Maps) is a model for pairwise interactions between motifs; and SCALI-HMM (HMMs for Structural Core ALIgnments) is a set of HMMs for the spatial arrangements of motifs. The parallels between the statistical models and theoretical models for folding pathways are discussed in this article; however, global sequence models are not discussed because they have been extensively reviewed elsewhere. The data used and algorithms presented in this article are available at http://www.bioinfo.rpi.edu/~bystrc/ (click on "servers" or "downloads") or by request to bystrc@rpi.edu .

Cell cycle regulatory protein p27KIP1 is a substrate and interacts with the protein kinase CK2

The protein kinase CK2 is constituted by two catalytic (alpha and/or alpha') and two regulatory (beta) subunits. CK2 phosphorylates more than 300 proteins with important functions in the cell cycle. This study has looked at the relation between CK2 and p27(KIP1), which is a regulator of the cell cycle and a known inhibitor of cyclin-dependent kinases (Cdk). We demonstrated that in vitro recombinant Xenopus laevis CK2 can phosphorylate recombinant human p27(KIP1), but this phosphorylation occurs only in the presence of the regulatory beta subunit. The principal site of phosphorylation is serine-83. Analysis using pull down and surface plasmon resonance (SPR) techniques showed that p27(KIP1) interacts with the beta subunit through two domains present in the amino and carboxyl ends, while CD spectra showed that p27(KIP1) phosphorylation by CK2 affects its secondary structure. Altogether, these results suggest that p27(KIP1) phosphorylation by CK2 probably involves a docking event mediated by the CK2beta subunit. The phosphorylation of p27(KIP1) by CK2 may affect its biological activity.

Intimate View of a Kinetic Protein Folding Intermediate: Residue-resolved Structure, Interactions, Stability, Folding and Unfolding Rates, Homogeneity

A cytochrome c kinetic folding intermediate was studied by hydrogen exchange (HX) pulse labeling. Advances in the technique and analysis made it possible to define the structured and unstructured regions, equilibrium stability, and kinetic opening and closing rates, all at an amino acid-resolved level. The entire N-terminal and C-terminal helices are formed and docked together at their normal native positions. They fray in both directions from the interaction region, due to a progression in both unfolding and refolding rates, leading to the surprising suggestion that helix propagation may proceed very slowly in the condensed milieu. Several native-like beta turns are formed. Some residues in the segment that will form the native 60s helix are protected but others are not, suggesting energy minimization to some locally non-native conformation in the transient intermediate. All other regions are unprotected, presumably dynamically disordered. The intermediate resembles a partially constructed native state. It is early, on-pathway, and all of the refolding molecules pass through it. These and related results consistently point to distinct, homogeneous, native-like intermediates in a stepwise sequential pathway, guided by the same factors that determine the native structure. Previous pulse labeling efforts have always assumed EX2 exchange during the labeling pulse, often leading to the suggestion of heterogeneous intermediates in alternative parallel pathways. The present work reveals a dominant role for EX1 exchange in the high pH labeling pulse, which will mimic heterogeneous behavior when EX2 exchange is assumed. The general problem of homogeneous versus heterogeneous intermediates and pathways is discussed.

Effects of domain dissection on the folding and stability of the 43 kDa protein PGK probed by NMR

The characterization of early folding intermediates is key to understanding the protein folding process. Previous studies of the N-domain of phosphoglycerate kinase (PGK) from Bacillus stearothermophilus combined equilibrium amide exchange data with a kinetic model derived from stopped-flow kinetics. Together, these implied the rapid formation of an intermediate with extensive native-like hydrogen bonding. However, there was an absence of protection in the region proximal to the C-domain in the intact protein. We now report data for the intact PGK molecule, which at 394 residues constitutes a major extension to the protein size for which such data can be acquired. The methods utilised to achieve the backbone assignment are described in detail, including a semi-automated protocol based on a simulated annealing Monte Carlo technique. A substantial increase in the stability of the contact region is observed, allowing protection to be inferred on both faces of the beta-sheet in the intermediate. Thus, the entire N-domain acts concertedly in the formation of the kinetic refolding intermediate rather than there existing a distinct local folding nucleus.

Tyrosine phosphorylation of the well packed ephrinB cytoplasmic beta-hairpin for reverse signaling. Structural consequences and binding properties

Tyrosine phosphorylation of the 22-residue cytoplasmic region of ephrinB induces its binding to the SH2 domain of Grb4, thus initiating reverse signaling pathways controlling cytoskeleton assembly and remodeling. Recently, the region corresponding to this 22-residue motif was demonstrated to adopt a well packed beta-hairpin structure with a high conformational stability in the unphosphorylated cytoplasmic subdomain. However, because the binding to Grb4 is phosphorylation-dependent and the hairpin contains three conserved tyrosine residues that may be phosphorylated, the key events remain unknown as to how tyrosine phosphorylation affects the structure of this well packed beta-hairpin and which phosphorylation site is relevant to SH2 domain binding. By characterizing the structural and binding properties of six 22-residue SH2 domain-binding motifs with different phosphorylated sites, the present study reveals that, as shown by circular dichroism and NMR, the unphosphorylated 22-residue motif adopts a well formed beta-hairpin structure in isolation from the ephrinB cytoplasmic subdomain. However, this beta-hairpin is radically abolished by tyrosine phosphorylation, regardless of the relative location and number of Tyr residues. Unexpectedly, the peptides with either Tyr304 or Tyr316 phosphorylated show high affinity binding to SH2 domain, whereas the peptide with Tyr311 phosphorylated has no detectable binding. This implies that ephrinB with Tyr311 phosphorylated might have a currently unidentified binding partner distinct from the Grb4 protein, because Tyr311 is known to be phosphorylated in vivo. Based on the results above, it is thus proposed that the disruption of the tight side-chain packing by tyrosine phosphorylation in the well structured region of a signaling protein may represent a general activation mechanism by which a cryptic binding site is disclosed for new protein-protein interactions.