Aromatic and methyl NOEs highlight hydrophobic clustering in the unfolded state of an SH3 domain

Denatured State Conformational Biases in Three-Helix Bundles Containing Divergent Sequences Localize near Turns and Helix Capping Residues

Rhodopseudomonas palustris cytochrome c', a four-helix bundle, and the second ubiquitin-associated domain, UBA(2), a three-helix bundle from the human homologue of yeast Rad23, HHR23A, deviate from random coil behavior under denaturing conditions in a fold-specific manner. The random coil deviations in each of these folds occur near interhelical turns and loops in their tertiary structures. Here, we examine an additional three-helix bundle with an identical fold to UBA(2), but a highly divergent sequence, the first ubiquitin-associated domain, UBA(1), of HHR23A. We use histidine-heme loop formation methods, employing eight single histidine variants, to probe for denatured state conformational bias of a UBA(1) domain fused to the N-terminus of iso-1-cytochrome c (iso-1-Cytc). Guanidine hydrochloride (GuHCl) denaturation shows that the iso-1-Cytc domain unfolds first, followed by the UBA(1) domain. Denatured state (4 and 6 M GuHCl) histidine-heme loop formation studies show that as the size of the histidine-heme loop increases, loop stability decreases, as expected for the Jacobson-Stockmayer relationship. However, loops formed with His35, His31, and His15, of UBA(1), are 0.6-1.1 kcal/mol more stable than expected from the Jacobson-Stockmayer relationship, confirming the importance of deviations of the denatured state from random coil behavior near interhelical turns of helical domains for facilitating folding to the correct topology. For UBA(1) and UBA(2), hydrophobic clusters on either side of the turns partially explain deviations from random coil behavior; however, helix capping also appears to be important.

Extensively sparse 13C labeling to simplify solid-state NMR 13C spectra of membrane proteins

Solid-state Nuclear Magnetic Resonance (ssNMR) is an emerging technique to investigate the structures and dynamics of membrane proteins in an artificial or native membrane environment. However, the structural studies of proteins by ssNMR are usually prolonged or impeded by signal assignments, especially the assignments of signals for collection of distance restraints, because of serious overlapping of signals in 2D ¹³C-¹³C spectra. Sparse labeling of ¹³C spins is an effective approach to simplify the ¹³C spectra and facilitate the extractions of distance restraints. Here, we propose a new reverse labeling combination of six types of amino acid residues (Ile, Leu, Phe, Trp, Tyr and Lys), and show a clean reverse labeling effect on a model membrane protein E. coli aquaporin Z (AqpZ). We further combine this reverse labeling combination and alternate ¹³C-¹²C labeling, and demonstrate an enhanced dilution effect in ¹³C-¹³C spectra. In addition, the influences of reverse labeling on the labeling of the other types of residues are quantitatively analyzed in the two strategies (1, reverse labeling and 2, reverse labeling combining alternate ¹³C-¹²C labeling). The signal intensities of some other types of residues in 2D ¹³C-¹³C spectra are observed to be 20-50% weaker because of the unwanted reverse labeling. The extensively sparse ¹³C labeling proposed in this study is expected to be useful in the collection of distance restraints using 2D ¹³C-¹³C spectra of membrane proteins.

Extended Experimental Inferential Structure Determination Method in Determining the Structural Ensembles of Disordered Protein States

Proteins with intrinsic or unfolded state disorder comprise a new frontier in structural biology, requiring the characterization of diverse and dynamic structural ensembles. We introduce a comprehensive Bayesian framework, the Extended Experimental Inferential Structure Determination (X-EISD) method, that calculates the maximum log-likelihood of a disordered protein ensemble. X-EISD accounts for the uncertainties of a range of experimental data and back-calculation models from structures, including NMR chemical shifts, J-couplings, Nuclear Overhauser Effects (NOEs), paramagnetic relaxation enhancements (PREs), residual dipolar couplings (RDCs), hydrodynamic radii (R h ), single molecule fluorescence Förster resonance energy transfer (smFRET) and small angle X-ray scattering (SAXS). We apply X-EISD to the joint optimization against experimental data for the unfolded drkN SH3 domain and find that combining a local data type, such as chemical shifts or J-couplings, paired with long-ranged restraints such as NOEs, PREs or smFRET, yields structural ensembles in good agreement with all other data types if combined with representative IDP conformers.

Collapse Transitions of Proteins and the Interplay Among Backbone, Sidechain, and Solvent Interactions

Proteins can collapse into compact globules or form expanded, solvent-accessible, coil-like conformations. Additionally, they can fold into well-defined three-dimensional structures or remain partially or entirely disordered. Recent discoveries have shown that the tendency for proteins to collapse or remain expanded is not intrinsically coupled to their ability to fold. These observations suggest that proteins do not have to form compact globules in aqueous solutions. They can be intrinsically disordered, collapsed, or expanded, and even form well-folded, elongated structures. This ability to decouple collapse from folding is determined by the sequence details of proteins. In this review, we highlight insights gleaned from studies over the past decade. Using a polymer physics framework, we explain how the interplay among sidechains, backbone units, and solvent determines the driving forces for collapsed versus expanded states in aqueous solvents.

Probing Denatured State Conformational Bias in a Three-Helix Bundle, UBA(2), Using a Cytochrome c Fusion Protein

Previous work with the four-helix-bundle protein cytochrome c' from Rhodopseudomonas palustris using histidine-heme loop formation methods revealed fold-specific deviations from random coil behavior in its denatured state ensemble. To examine the generality of this finding, we extend this work to a three-helix-bundle polypeptide, the second ubiquitin-associated domain, UBA(2), of the human DNA excision repair protein. We use yeast iso-1-cytochrome c as a scaffold, fusing the UBA(2) domain at the N-terminus of iso-1-cytochrome c. We have engineered histidine into highly solvent accessible positions of UBA(2), creating six single histidine variants. Guanidine hydrochloride denaturation studies show that the UBA(2)-cytochrome c fusion protein unfolds in a three-state process with iso-1-cytochrome c unfolding first. Furthermore, engineered histidine residues in UBA(2) strongly destabilize the iso-1-cytochrome c domain. Equilibrium and kinetic histidine-heme loop formation measurements in the denatured state at 4 and 6 M guanidine hydrochloride show that loop stability decreases as the size of the histidine-heme loop increases, in accord with the Jacobson-Stockmayer equation. However, we observe that the His27-heme loop is both more stable than expected from the Jacobson-Stockmayer relationship and breaks more slowly than expected. These results show that the sequence near His27, which is in the reverse turn between helices 2 and 3 of UBA(2), is prone to persistent interactions in the denatured state. Therefore, consistent with our results for cytochrome c', this reverse turn sequence may help to establish the topology of this fold by biasing the conformational distribution of the denatured state.

Helical Propensity Affects the Conformational Properties of the Denatured State of Cytochrome c'

Changing the helical propensity of a polypeptide sequence might be expected to affect the conformational properties of the denatured state of a protein. To test this hypothesis, alanines at positions 83 and 87 near the center of helix 3 of cytochrome c' from Rhodopseudomonas palustris were mutated to serine to decrease the stability of this helix. A set of 13 single histidine variants in the A83S/A87S background were prepared to permit assessment of the conformational properties of the denatured state using histidine-loop formation in 3 M guanidine hydrochloride. The data are compared with previous histidine-heme loop formation data for wild-type cytochrome c'. As expected, destabilization of helix 3 decreases the global stabilities of the histidine variants in the A83S/A87S background relative to the wild-type background. Loop stability versus loop size data yields a scaling exponent of 2.1 ± 0.2, similar to the value of 2.3 ± 0.2 obtained for wild-type cytochrome c'. However, the stabilities of all histidine-heme loops, which contain the helix 3 sequence segment, are increased in the A83S/A87S background compared to the wild-type background. Rate constants for histidine-heme loop breakage are similar for the wild-type and A83S/A87S variants. However, for histidine-heme loops that contain the helix 3 sequence segment, the rate constants for loop formation increase in the A83S/A87S background compared to the wild-type background. Thus, residual helical structure appears to stiffen the polypeptide chain slowing loop formation in the denatured state. The implications of these results for protein folding mechanisms are discussed.

Effect of an Imposed Contact on Secondary Structure in the Denatured State of Yeast Iso-1-cytochrome c

There is considerable evidence that long-range interactions stabilize residual protein structure under denaturing conditions. However, evaluation of the effect of a specific contact on structure in the denatured state has been difficult. Iso-1-cytochrome c variants with a Lys54 → His mutation form a particularly stable His-heme loop in the denatured state, suggestive of loop-induced residual structure. We have used multidimensional nuclear magnetic resonance methods to assign ¹H and ¹⁵N backbone amide and ¹³C backbone and side chain chemical shifts in the denatured state of iso-1-cytochrome c carrying the Lys54 → His mutation in 3 and 6 M guanidine hydrochloride and at both pH 6.4, where the His54-heme loop is formed, and pH 3.6, where the His54-heme loop is broken. Using the secondary structure propensity score, with the 6 M guanidine hydrochloride chemical shift data as a random coil reference state for data collected in 3 M guanidine hydrochloride, we found residual helical structure in the denatured state for the 60s helix and the C-terminal helix, but not in the N-terminal helix in the presence or absence of the His54-heme loop. Non-native helical structure is observed in two regions that form Ω-loops in the native state. There is more residual helical structure in the C-terminal helix at pH 6.4 when the loop is formed. Loop formation also appears to stabilize helical structure near His54, consistent with induction of helical structure observed when His-heme bonds form in heme-peptide model systems. The results are discussed in the context of the folding mechanism of cytochrome c.

Finding Our Way in the Dark Proteome

The traditional structure-function paradigm has provided significant insights for well-folded proteins in which structures can be easily and rapidly revealed by X-ray crystallography beamlines. However, approximately one-third of the human proteome is comprised of intrinsically disordered proteins and regions (IDPs/IDRs) that do not adopt a dominant well-folded structure, and therefore remain "unseen" by traditional structural biology methods. This Perspective considers the challenges raised by the "Dark Proteome", in which determining the diverse conformational substates of IDPs in their free states, in encounter complexes of bound states, and in complexes retaining significant disorder requires an unprecedented level of integration of multiple and complementary solution-based experiments that are analyzed with state-of-the art molecular simulation, Bayesian probabilistic models, and high-throughput computation. We envision how these diverse experimental and computational tools can work together through formation of a "computational beamline" that will allow key functional features to be identified in IDP structural ensembles.

Structural mapping of a chaperone-substrate interaction surface

NMR spectroscopy is used to detect site-specific intermolecular short-range contacts in a membrane-protein-chaperone complex. This is achieved by an "orthogonal" isotope-labeling scheme that permits the unambiguous detection of intermolecular NOEs between the well-folded chaperone and the unfolded substrate ensemble. The residues involved in these contacts are part of the chaperone-substrate contact interface. The approach is demonstrated for the 70 kDa bacterial Skp-tOmpA complex.

The RING domain of the scaffold protein Ste5 adopts a molten globular character with high thermal and chemical stability

Ste5 is a scaffold protein that controls the pheromone response of the MAP-kinase cascade in yeast cells. Upon pheromone stimulation, Ste5 (through its RING-H2 domain) interacts with the β and γ subunits of an activated heterodimeric G protein and promotes activation of the MAP-kinase cascade. With structural and biophysical studies, we show that the Ste5 RING-H2 domain exists as a molten globule under native buffer conditions, in yeast extracts, and even in denaturing conditions containing urea (7 M). Furthermore, it exhibits high thermal stability in native conditions. Binding of the Ste5 RING-H2 domain to the physiological Gβ/γ (Ste4/Ste18) ligand is accompanied by a conformational transition into a better folded, more globular structure. This study reveals novel insights into the folding mechanism and recruitment of binding partners by the Ste5 RING-H2 domain. We speculate that many RING domains may share a similar mechanism of substrate recognition and molten-globule-like character.

Structural characterization of intrinsically disordered proteins by NMR spectroscopy

Recent advances in NMR methodology and techniques allow the structural investigation of biomolecules of increasing size with atomic resolution. NMR spectroscopy is especially well-suited for the study of intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) which are in general highly flexible and do not have a well-defined secondary or tertiary structure under functional conditions. In the last decade, the important role of IDPs in many essential cellular processes has become more evident as the lack of a stable tertiary structure of many protagonists in signal transduction, transcription regulation and cell-cycle regulation has been discovered. The growing demand for structural data of IDPs required the development and adaption of methods such as ¹³C-direct detected experiments, paramagnetic relaxation enhancements (PREs) or residual dipolar couplings (RDCs) for the study of ‘unstructured’ molecules in vitro and in-cell. The information obtained by NMR can be processed with novel computational tools to generate conformational ensembles that visualize the conformations IDPs sample under functional conditions. Here, we address NMR experiments and strategies that enable the generation of detailed structural models of IDPs.

Observation of persistent α-helical content and discrete types of backbone disorder during a molten globule to ordered peptide transition via deep-UV resonance Raman spectroscopy

The molten globule state can aide in the folding of a protein to a functional structure and is loosely defined as an increase in structural disorder with conservation of the ensemble secondary structure content. Simultaneous observation of persistent secondary structure content with increased disorder has remained experimentally problematic. As a consequence, modeling how the molten globule state remains stable and how it facilitates proper folding remains difficult due to a lack of amenable spectroscopic techniques to characterize this class of partially unfolded proteins. Previously, deep-UV resonance Raman (dUVRR) spectroscopy has proven useful in the resolution of global and local structural fluctuations in the secondary structure of proteins. In this work, dUVRR was employed to study the molten globule to ordered transition of a model four-helix bundle protein, HP7. Both the average ensemble secondary structure and types of local disorder were monitored, without perturbation of the solvent, pH, or temperature. The molten globule to ordered transition is induced by stepwise coordination of two heme molecules. Persistent dUVRR spectral features in the amide III region at 1295-1301 and 1335-1338 cm-1 confirm previous observations that HP7 remains predominantly helical in the molten globule versus the fully ordered state. Additionally, these spectra represent the first demonstration of conserved helical content in a molten globule protein. With successive heme binding significant losses are observed in the spectral intensity of the amide III3 and S regions (1230-1260 and 1390 cm-1, respectively), which are known to be sensitive to local disorder. These observations indicate that there is a decrease in the structural populations able to explore various extended conformations, with successive heme binding events. DUVRR spectra indicate that the first heme coordination between two helical segments diminishes exploration of more elongated backbone structural conformations in the inter-helical regions. A second heme coordination by the remaining two helices further restricts protein motion.

Tryptophan stabilizes His-heme loops in the denatured state only when it is near a loop end

We use a host-guest approach to evaluate the effect of Trp guest residues relative to Ala on the kinetics and thermodynamics of formation of His-heme loops in the denatured state of iso-1-cytochrome c at 1.5, 3.0, and 6.0 M guanidine hydrochloride (GdnHCl). Trp guest residues are inserted into an alanine-rich segment placed after a unique His near the N-terminus of iso-1-cytochrome c. Trp guest residues are either 4 or 10 residues from the His end of the 28-residue loop. We find the guest Trp stabilizes the His-heme loop at all GdnHCl concentrations when it is the 4th, but not the 10th, residue from the His end of the loop. Thus, residues near loop ends are most important in developing topological constraints in the denatured state that affect protein folding. In 1.5 M GdnHCl, the loop stabilization is ~0.7 kcal/mol, providing a thermodynamic rationale for the observation that Trp often mediates residual structure in the denatured state. Measurement of loop breakage rate constants, k(b,His), indicates that loop stabilization by the Trp guest residues occurs completely after the transition state for loop formation in 6.0 M GdnHCl. Under poorer solvent conditions, approximately half of the stabilization of the loop develops in the transition state, consistent with contacts in the denatured state being energetically downhill and providing evidence for funneling even near the rim of the folding funnel.

Ensemble modeling of protein disordered states: experimental restraint contributions and validation

Disordered states of proteins include the biologically functional intrinsically disordered proteins and the unfolded states of normally folded proteins. In recent years, ensemble-modeling strategies using various experimental measurements as restraints have emerged as powerful means for structurally characterizing disordered states. However, these methods are still in their infancy compared with the structural determination of folded proteins. Here, we have addressed several issues important to ensemble modeling using our ENSEMBLE methodology. First, we assessed how calculating ensembles containing different numbers of conformers affects their structural properties. We find that larger ensembles have very similar properties to smaller ensembles fit to the same experimental restraints, thus allowing a considerable speed improvement in our calculations. In addition, we analyzed the contributions of different experimental restraints to the structural properties of calculated ensembles, enabling us to make recommendations about the experimental measurements that should be made for optimal ensemble modeling. The effects of different restraints, most significantly from chemical shifts, paramagnetic relaxation enhancements and small-angle X-ray scattering, but also from other data, underscore the importance of utilizing multiple sources of experimental data. Finally, we validate our ENSEMBLE methodology using both cross-validation and synthetic experimental restraints calculated from simulated ensembles. Our results suggest that secondary structure and molecular size distribution can generally be modeled very accurately, whereas the accuracy of calculated tertiary structure is dependent on the number of distance restraints used.

Measuring 1HN temperature coefficients in invisible protein states by relaxation dispersion NMR spectroscopy

A method based on the Carr-Purcell-Meiboom-Gill relaxation dispersion experiment is presented for measuring the temperature coefficients of amide proton chemical shifts of low populated 'invisible' protein states that exchange with a 'visible' ground state on the millisecond time-scale. The utility of the approach is demonstrated with an application to an I58D mutant of the Pfl6 Cro protein that undergoes exchange between the native, folded state and a cold denatured, unfolded conformational ensemble that is populated at a level of 6% at 2.5°C. A wide distribution of amide temperature coefficients is measured for the unfolded state. The distribution is centered about -5.6 ppb/K, consistent with an absence of intra-molecular hydrogen bonds, on average. However, the large range of values (standard deviation of 2.1 ppb/K) strongly supports the notion that the unfolded state of the protein is not a true random coil polypeptide chain.

Propensities of aromatic amino acids versus leucine and proline to induce residual structure in the denatured-state ensemble of iso-1-cytochrome c

Histidine-heme loop formation in the denatured state of a protein is a sensitive means for probing residual structure under unfolding conditions. In this study, we use a host-guest approach to investigate the relative tendencies of different amino acids to promote residual structure under denaturing conditions. The host for this work is a 6-amino-acid insert of five alanines, followed by a lysine engineered immediately following a unique histidine near the N-terminus of yeast iso-1-cytochrome c. We substitute the fourth alanine in this sequence HAAAXAK (with X=Trp, Phe, Tyr, and Leu). The effects of proline are tested with substitutions at positions 1 and 5 in the insert (HPAAAAK and HAAAAPK, respectively). Thermodynamic studies on His-heme loop formation in 3 M guanidine hydrochloride reveal significant stabilization of residual structure by aromatic amino acids, particularly Trp and Phe, and minimal stabilization of residual structure by Leu. Prolines slightly disfavor His-heme loop formation, presumably due to enhanced chain stiffness. Kinetic studies reveal that much of the change in His-heme loop stability for the aromatic amino acids is caused by a slowdown in the rate of His-heme loop breakage, indicating that residual structure is preferentially stabilized in the closed-loop form of the denatured state.

Sequence determinants of compaction in intrinsically disordered proteins

Intrinsically disordered proteins (IDPs), which lack folded structure and are disordered under nondenaturing conditions, have been shown to perform important functions in a large number of cellular processes. These proteins have interesting structural properties that deviate from the random-coil-like behavior exhibited by chemically denatured proteins. In particular, IDPs are often observed to exhibit significant compaction. In this study, we have analyzed the hydrodynamic radii of a number of IDPs to investigate the sequence determinants of this compaction. Net charge and proline content are observed to be strongly correlated with increased hydrodynamic radii, suggesting that these are the dominant contributors to compaction. Hydrophobicity and secondary structure, on the other hand, appear to have negligible effects on compaction, which implies that the determinants of structure in folded and intrinsically disordered proteins are profoundly different. Finally, we observe that polyhistidine tags seem to increase IDP compaction, which suggests that these tags have significant perturbing effects and thus should be removed before any structural characterizations of IDPs. Using the relationships observed in this analysis, we have developed a sequence-based predictor of hydrodynamic radius for IDPs that shows substantial improvement over a simple model based upon chain length alone.

Cooperative formation of native-like tertiary contacts in the ensemble of unfolded states of a four-helix protein

In studies of the ensembles of unfolded structures of a four-helix bundle protein, we have detected the presence of potential precursors of native tertiary structures. These observations were based on the perturbation of NMR chemical shifts of the protein backbone atoms by single site mutations. Some mutations change the chemical shifts of residues remote from the site of mutation indicating the presence of an interaction between the mutated and the remote residues, suggesting that the formation of helix segments and helix-helix interactions is cooperative. We can begin to track down the folding mechanism of this protein using only experimental data by combining the information available for the rate limiting structure formation during the folding process with measurements of the site specific hydrogen bond formation in the burst phase, and with the existence prior to the folding reaction of tertiary structures in the ensemble of otherwise unfolded structures observed in the present study.

15N H/D-SOLEXSY experiment for accurate measurement of amide solvent exchange rates: application to denatured drkN SH3

Amide solvent exchange rates are regarded as a valuable source of information on structure/dynamics of unfolded (disordered) proteins. Proton-based saturation transfer experiments, normally used to measure solvent exchange, are known to meet some serious difficulties. The problems mainly arise from the need to (1) manipulate water magnetization and (2) discriminate between multiple magnetization transfer pathways that occur within the proton pool. Some of these issues are specific to unfolded proteins. For example, the compensation scheme used to cancel the Overhauser effect in the popular CLEANEX experiment is not designed for use with unfolded proteins. In this report we describe an alternative experimental strategy, where amide (15)N is used as a probe of solvent exchange. The experiment is performed in 50% H(2)O-50% D(2)O solvent and is based on the (HACACO)NH pulse sequence. The resulting spectral map is fully equivalent to the conventional HSQC. To fulfill its purpose, the experiment monitors the conversion of deuterated species, (15)N(D), into protonated species, (15)N(H), as effected by the solvent exchange. Conceptually, this experiment is similar to EXSY which prompted the name of (15)N(H/D)-SOLEXSY (SOLvent EXchange SpectroscopY). Of note, our experimental scheme, which relies on nitrogen rather than proton to monitor solvent exchange, is free of the complications described above. The developed pulse sequence was used to measure solvent exchange rates in the chemically denatured state of the drkN SH3 domain. The results were found to correlate well with the CLEANEX-PM data, r = 0.97, thus providing a measure of validation for both techniques. When the experimentally measured exchange rates are converted into protection factors, most of the values fall in the range 0.5-2, consistent with random-coil behavior. However, elevated values, ca. 5, are obtained for residues R38 and A39, as well as the side-chain indole of W36. This is surprising, given that high protection factors imply hydrogen bonding or hydrophobic burial not expected to occur in a chemically denatured state of a protein. We, therefore, hypothesized that elevated protection factors are an artefact arising from the calculation of the reference (random-coil) exchange rates. To confirm this hypothesis, we prepared samples of several short peptides derived from the sequence of the drkN SH3 domain; these samples were used to directly measure the reference exchange rates. The revised protection factors obtained in this manner proved to be close to 1.0. These results also have implications for the more compact unfolded state of drkN SH3, which appears to be fully permeable to water as well, with no manifestations of hydrophobic burial.

Insights into protein aggregation by NMR characterization of insoluble SH3 mutants solubilized in salt-free water

Protein aggregation in vivo has been extensively associated with a large spectrum of human diseases. On the other hand, mechanistic insights into protein aggregation in vitro were incomplete due to the inability in solubilizing insoluble proteins for high-resolution biophysical investigations. However, a new avenue may be opened up by our recent discovery that previously-thought insoluble proteins can in fact be solubilized in salt-free water. Here we use this approach to study the NMR structural and dynamic properties of an insoluble SH3 mutant with a naturally-occurring insertion of Val22 at the tip of the diverging turn. The obtained results reveal: 1) regardless of whether the residue is Val, Ala, Asp or Arg, the insertion will render the first hNck2 SH3 domain to be insoluble in buffers. Nevertheless, all four mutants could be solubilized in salt-free water and appear to be largely unfolded as evident from their CD and NMR HSQC spectra. 2) Comparison of the chemical shift deviations reveals that while in V22-SH3 the second helical region is similarly populated as in the wild-type SH3 at pH 2.0, the first helical region is largely unformed. 3) In V22-SH3, many non-native medium-range NOEs manifest to define non-native helical conformations. In the meanwhile a small group of native-like long-range NOEs still persists, indicating the existence of a rudimentary native-like tertiary topology. 4) Although overall, V22-SH3 has significantly increased backbone motions on the ps-ns time scale, some regions still own restricted backbone motions as revealed by analyzing ¹⁵N relaxation data. Our study not only leads to the establishment of the first high-resolution structural and dynamic picture for an insoluble protein, but also shed more light on the molecular events for the nonhierarchical folding mechanism. Furthermore, a general mechanism is also proposed for in vivo protein aggregation triggered by the genetic mutation and posttranslational modification.

Structure-microbicidal activity relationship of synthetic fragments derived from the antibacterial alpha-helix of human lactoferrin

There is a need for new microbicidal agents with therapeutic potential due to antibiotic resistance in bacteria and fungi. In this study, the structure-microbicidal activity relationship of amino acid residues 14 to 31 (sequence 14-31) from the N-terminal end, corresponding to the antibacterial alpha-helix of human lactoferrin (LF), was investigated by downsizing, alanine scanning, and substitution of amino acids. Microbicidal analysis (99% killing) was performed by a microplate assay using Escherichia coli, Staphylococcus aureus, and Candida albicans as test organisms. Starting from the N-terminal end, downsizing of peptide sequence 14-31 showed that the peptide sequence 19-31 (KCFQWQRNMRKVR, HL9) was the optimal length for antimicrobial activity. Furthermore, HL9 bound to lipid A/lipopolysaccharide, as shown by neutralizing endotoxic activity in a Limulus assay. Alanine scanning of peptide sequence 20-31 showed that Cys20, Trp23, Arg28, Lys29, or Arg31 was important for expressing full killing activity, particularly against C. albicans. Substituting the neutral hydrophilic amino acids Gln24 and Asn26 for Lys and Ala (HLopt2), respectively, enhanced microbicidal activity significantly against all test organisms compared to the amino acids natural counterpart, also, in comparison with HL9, HLopt2 had more than 10-fold-stronger fungicidal activity. Furthermore, HLopt2 was less affected by metallic salts than HL9. The microbicidal activity of HLopt2 was slightly reduced only at pH 7.0, as tested in the pH range of 4.5 to 7.5. The results showed that the microbicidal activity of synthetic peptide sequences, based on the antimicrobial alpha-helix region of LF, can be significantly enhanced by optimizing the length and substitution of neutral amino acids at specific positions, thus suggesting a sequence lead with therapeutic potential.

Paramagnetic relaxation enhancements in unfolded proteins: theory and application to drkN SH3 domain

Site-directed spin labeling in combination with paramagnetic relaxation enhancement (PRE) measurements is one of the most promising techniques for studying unfolded proteins. Since the pioneering work of Gillespie and Shortle (J Mol Biol 1997;268:158), PRE data from unfolded proteins have been interpreted using the theory that was originally developed for rotational spin relaxation. At the same time, it can be readily recognized that the relative motion of the paramagnetic tag attached to the peptide chain and the reporter spin such as (1)H(N) is best described as a translation. With this notion in mind, we developed a number of models for the PRE effect in unfolded proteins: (i) mutual diffusion of the two tethered spheres, (ii) mutual diffusion of the two tethered spheres subject to a harmonic potential, (iii) mutual diffusion of the two tethered spheres subject to a simulated mean-force potential (Smoluchowski equation); (iv) explicit-atom molecular dynamics simulation. The new models were used to predict the dependences of the PRE rates on the (1)H(N) residue number and static magnetic field strength; the results are appreciably different from the Gillespie-Shortle model. At the same time, the Gillespie-Shortle approach is expected to be generally adequate if the goal is to reconstruct the distance distributions between (1)H(N) spins and the paramagnetic center (provided that the characteristic correlation time is known with a reasonable accuracy). The theory has been tested by measuring the PRE rates in three spin-labeled mutants of the drkN SH3 domain in 2M guanidinium chloride. Two modifications introduced into the measurement scheme-using a reference compound to calibrate the signals from the two samples (oxidized and reduced) and using peak volumes instead of intensities to determine the PRE rates-lead to a substantial improvement in the quality of data. The PRE data from the denatured drkN SH3 are mostly consistent with the model of moderately expanded random-coil protein, although part of the data point toward a more compact structure (local hydrophobic cluster). At the same time, the radius of gyration reported by Choy et al. (J Mol Biol 2002;316:101) suggests that the protein is highly expanded. This seemingly contradictory evidence can be reconciled if one assumes that denatured drkN SH3 forms a conformational ensemble that is dominated by extended conformations, yet also contains compact (collapsed) species. Such behavior is apparently more complex than predicted by the model of a random-coil protein in good solvent/poor solvent.

Non-sequence-specific interactions can account for the compaction of proteins unfolded under "native" conditions

Proteins unfolded by high concentrations of chemical denaturants adopt expanded, largely structure-free ensembles of conformations that are well approximated as random coils. In contrast, globular proteins unfolded under less denaturing conditions (via mutations, or transiently unfolded after a rapid jump to native conditions) and molten globules (arising due to mutations or cosolvents) are often compact. Here we explore the origins of this compaction using a truncated equilibrium-unfolded variant of the 57-residue FynSH3 domain. As monitored by far-UV circular dichroism, NMR spectroscopy, and hydrogen-exchange kinetics, CDelta4 (a 4-residue carboxy-terminal deletion variant of FynSH3) appears to be largely unfolded even in the absence of denaturant. Nevertheless, CDelta4 is quite compact under these conditions, with a hydrodynamic radius only slightly larger than that of the native protein. In order to understand the origins of this molten-globule-like compaction, we have characterized a random sequence polypeptide of identical amino acid composition to CDelta4. Notably, we find that the hydrodynamic radius of this random sequence polypeptide also approaches that of the native protein. Thus, while native-like interactions may contribute to the formation of compact "unfolded" states, it appears that non-sequence-specific monomer-monomer interactions can also account for the dramatic compaction observed for molten globules and the "physiological" unfolded state.

Structure and disorder in an unfolded state under nondenaturing conditions from ensemble models consistent with a large number of experimental restraints

Obtaining detailed structural models of disordered states of proteins under nondenaturing conditions is important for a better understanding of both functional intrinsically disordered proteins and unfolded states of folded proteins. Extensive experimental characterization of the drk N-terminal SH3 domain unfolded state has shown that, although it appears to be highly disordered, it possesses significant nonrandom secondary and tertiary structure. In our previous attempts to generate structural models of the unfolded state using the program ENSEMBLE, we were limited by insufficient experimental restraints and conformational sampling. In this study, we have vastly expanded our experimental restraint set to include (1)H-(15)N residual dipolar couplings, small-angle X-ray scattering measurements, nitroxide paramagnetic relaxation enhancements, O(2)-induced (13)C paramagnetic shifts, hydrogen-exchange protection factors, and (15)N R(2) data, in addition to the previously used nuclear Overhauser effects, amino terminal Cu(2+)-Ni(2+) binding paramagnetic relaxation enhancements, J-couplings, chemical shifts, hydrodynamic radius, and solvent accessibility restraints. We have also implemented a new ensemble calculation methodology that uses iterative conformational sampling and seeks to calculate the simplest possible ensemble models. As a result, we can now generate ensembles that are consistent with much larger experimental data sets than was previously possible. Although highly heterogeneous and having broad molecular size distributions, the calculated drk N-terminal SH3 domain unfolded-state ensembles have very different properties than expected for random or statistical coils and possess significant nonnative alpha-helical structure and both native-like and nonnative tertiary structure.

Biophysical characterization of intrinsically disordered proteins

The challenges associated with the structural characterization of disordered proteins have resulted in the application of a host of biophysical methods to such systems. NMR spectroscopy is perhaps the most readily suited technique for providing high-resolution structural information on disordered protein states in solution. Optical methods, solid state NMR, ESR and X-ray scattering can also provide valuable information regarding the ensemble of conformations sampled by disordered states. Finally, computational studies have begun to assume an increasingly important role in interpreting and extending the impact of experimental data obtained for such systems. This article discusses recent advances in the applications of these methods to intrinsically disordered proteins.

Sequence-Specific Resonance Assignment of Soluble Nonglobular Proteins by 7D APSY-NMR Spectroscopy

Based on sequence-specific resonance assignments, NMR is the method of choice for obtaining atomic-resolution experimental data on soluble nonglobular proteins. So far, however, NMR assignment of unfolded polypeptides in solution has been a time-consuming task, mainly due to the small chemical shift dispersion, which has limited practical applications of the NMR approach. This paper presents an efficient, fully automated method for sequence-specific backbone and beta-carbon NMR assignment of soluble nonglobular proteins with sizes up to at least 150 residues. The procedure is based on new APSY (automated projection spectroscopy) experiments which benefit from the short effective rotational correlation times in soluble nonglobular polypeptides to record five- to seven-dimensional NMR data sets, which reliably resolves chemical shift degeneracies. Fully automated sequence-specific resonance assignments of the backbone nuclei and C(beta) are described for the uniformly (13)C,(15)N-labeled urea-denatured 148-residue outer membrane protein X (OmpX) from E. coli. The method is generally applicable to systems with similar spectroscopic properties as unfolded OmpX, and we anticipate that this paper may open the door for extensive atomic-resolution studies of chemical denaturant-unfolded proteins, as well as some classes of functional nonglobular polypeptides in solution.

On the NMR analysis of pKa values in the unfolded state of proteins by extrapolation to zero denaturant

Detailed knowledge of the pH-dependence in both folded and unfolded states of proteins is essential to understand the role of electrostatics in protein stability. The increasing number of natively disordered proteins constitutes an excellent source for the NMR analysis of pKa values in the unfolded state of proteins. However, the tendency of many natively disordered proteins to aggregate via intermolecular hydrophobic clusters limits their NMR analysis over a wide pH range. To assess whether the pKa values in natively disordered polypeptides can be extrapolated from NMR measurements in the presence of denaturants, the natively disordered backbone of the C-terminal fragment 75 to 105 of Human Thioredoxin was studied. First, assignments using triple resonance experiments were performed to confirm lack of secondary structure. Then the pH-dependence of the amides and carboxylate side chains of Glu residues (Glu88, Glu95, Glu98, and Glu103) in the pH range from 2.0 to 7.0 was monitored using 2D 1H15N HSQC and 3D C(CO)NH experiments, and the behavior of their amides and corresponding carboxyl groups was compared to confirm the absence of nonlocal interactions. Lastly, the effect of increasing dimethyl urea concentration on the pKa values of these Glu residues was monitored. The results indicate that: (i) the dispersion in the pKa of carboxyl groups and the pH midpoints of amides in Glu residues is about 0.5 pH units and 0.6 pH units, respectively; (ii) the backbone amides of the Glu residues exhibit pH midpoints which are within 0.2 pH units from those of their carboxylates; (iii) the addition of denaturant produces upshifts in the pKa values of Glu residues that are nearly independent of their position in the sequence; and (iv) these upshifts show a nonlinear behavior in denaturant concentration, complicating the extrapolation to zero denaturant. Nevertheless, the relative ordering of the pKa values of Glu residues is preserved over the whole range of denaturant concentrations indicating that measurements at high denaturant concentration (e.g. 4 M dimethyl urea) can yield a qualitatively correct ranking of the pKa of these residues in natively disordered proteins whose pH-dependence cannot be monitored directly by NMR.

A pre-existing hydrophobic collapse in the unfolded state of an ultrafast folding protein

Insights into the conformational passage of a polypeptide chain across its free energy landscape have come from the judicious combination of experimental studies and computer simulations. Even though some unfolded and partially folded proteins are now known to possess biological function or to be involved in aggregation phenomena associated with disease states, experimentally derived atomic-level information on these structures remains sparse as a result of conformational heterogeneity and dynamics. Here we present a technique that can provide such information. Using a 'Trp-cage' miniprotein known as TC5b (ref. 5), we report photochemically induced dynamic nuclear polarization NMR pulse-labelling experiments that involve rapid in situ protein refolding. These experiments allow dipolar cross-relaxation with hyperpolarized aromatic side chain nuclei in the unfolded state to be identified and quantified in the resulting folded-state spectrum. We find that there is residual structure due to hydrophobic collapse in the unfolded state of this small protein, with strong inter-residue contacts between side chains that are relatively distant from one another in the native state. Prior structuring, even with the formation of non-native rather than native contacts, may be a feature associated with fast folding events in proteins.

Improved Structural Characterizations of the drkN SH3 Domain Unfolded State Suggest a Compact Ensemble with Native-like and Non-native Structure

Due to their dynamic ensemble nature and a deficiency of experimental restraints, disordered states of proteins are difficult to characterize structurally. Here, we have expanded upon our previous work on the unfolded state of the Drosophila drk N-terminal (drkN) SH3 domain with our program ENSEMBLE, which assigns population weights to pregenerated conformers in order to calculate ensembles of structures whose properties are collectively consistent with experimental measurements. The experimental restraint set has been enlarged with newly measured paramagnetic relaxation enhancements from Cu(2+) bound to an amino terminal Cu(2+)-Ni(2+) binding (ATCUN) motif as well as nuclear Overhauser effect (NOE) and hydrogen exchange data from recent studies. In addition, two new pseudo-energy minimization algorithms have been implemented that have dramatically improved the speed of ENSEMBLE population weight assignment. Finally, we have greatly improved our conformational sampling by utilizing a variety of techniques to generate both random structures and structures that are biased to contain elements of native-like or non-native structure. Although it is not possible to uniquely define a representative structural ensemble, we have been able to assess various properties of the drkN SH3 domain unfolded state by performing ENSEMBLE minimizations of different conformer pools. Specifically, we have found that the experimental restraint set enforces a compact structural distribution that is not consistent with an overall native-like topology but shows preference for local non-native structure in the regions corresponding to the diverging turn and the beta5 strand of the folded state and for local native-like structure in the region corresponding to the beta6 and beta7 strands. We suggest that this approach could be generally useful for the structural characterization of disordered states.

Role of the Charge–Charge Interactions in Defining Stability and Halophilicity of the CspB Proteins

Charge-charge interactions on the surface of native proteins are important for protein stability and can be computationally redesigned in a rational way to modulate protein stability. Such computational effort led to an engineered protein, CspB-TB that has the same core as the mesophilic cold shock protein CspB-Bs from Bacillus subtilis, but optimized distribution of charge-charge interactions on the surface. The CspB-TB protein shows an increase in the transition temperature by 20 degrees C relative to the unfolding temperature of CspB-Bs. The CspB-TB and CspB-Bs protein pair offers a unique opportunity to further explore the energetics of charge-charge interactions as the substitutions at the same sequence positions are done in largely similar structural but different electrostatic environments. In particular we addressed two questions. What is the contribution of charge-charge interactions in the unfolded state to the protein stability and how amino acid substitutions modulate the effect of increase in ionic strength on protein stability (i.e. protein halophilicity). To this end, we experimentally measured the stabilities of over 100 variants of CspB-TB and CspB-Bs proteins with substitutions at charged residues. We also performed computational modeling of these protein variants. Analysis of the experimental and computational data allowed us to conclude that the charge-charge interactions in the unfolded state of two model proteins CspB-Bs and CspB-TB are not very significant and computational models that are based only on the native state structure can adequately, i.e. qualitatively (stabilizing versus destabilizing) and semi-quantitatively (relative rank order), predict the effects of surface charge neutralization or reversal on protein stability. We also show that the effect of ionic strength on protein stability (protein halophilicity) appears to be mainly due to the screening of the long-range charge-charge interactions.

Oxygen as a Paramagnetic Probe of Clustering and Solvent Exposure in Folded and Unfolded States of an SH3 Domain

The N-terminal SH3 domain of the Drosophila modular protein Drk undergoes slow exchange between a folded (Fexch) and highly populated unfolded (Uexch) state under nondenaturing buffer conditions, enabling both Fexch and Uexch states to be simultaneously monitored. The addition of dissolved oxygen, equilibrated to a partial pressure of either 30 atm or 60 atm, provides the means to study solvent exposure with atomic resolution via 13C NMR paramagnetic shifts in 1H,13C HSQC (heteronuclear single quantum coherence) spectra. Absolute differences in these paramagnetic shifts between the Fexch and Uexch states allow the discrimination of regions of the protein which undergo change in solvent exposure upon unfolding. Contact with dissolved oxygen for both the Fexch and Uexch states could also be assessed through 13C paramagnetic shifts which were normalized based on the corresponding paramagnetic shifts seen in the free amino acids. In the Fexch state, the 13C nuclei belonging to the hydrophobic core of the protein exhibited very weak normalized paramagnetic shifts while those with greater solvent accessible surface area exhibited significantly larger normalized shifts. The Uexch state displayed less varied 13C paramagnetic shifts although distinct regions of protection from solvent exposure could be identified by a lack of such shifts. These regions, which included Phe9, Thr12, Ala13, Lys21, Thr22, Ile24, Ile27, and Arg38, overlapped with those found to have residual nativelike and non-native structures in previous studies and in some cases provided novel information. Thus, the paramagnetic shifts from dissolved oxygen are highly useful in the study of a transient structure or clustering in disordered systems, where conventional NMR measurements (couplings, chemical shift deviations from random coil values, and NOEs) may give little information.

Atomic-level characterization of disordered protein ensembles

The roles of unfolded states of proteins in normal folding and in diseases involving aggregation, as well as the prevalence and regulatory functions of intrinsically disordered proteins, have become increasingly recognized. The structural representation of these disordered states as ensembles of interconverting conformers can therefore provide critical insights. Experimental methods can be used to probe ensemble-averaged structural properties of disordered states and computational approaches generate representative ensembles of conformers using experimental restraints. In particular, NMR and small-angle X-ray scattering provide quantitative data that can readily be incorporated into calculations. These techniques have gleaned structural information about denatured, unfolded and intrinsically disordered proteins. The use of experimental data in different computational approaches, including ensemble molecular dynamics simulations and algorithms that assign populations to pregenerated conformers, has highlighted the presence of both local and long-range structure, and the occurrence of native-like and non-native interactions in unfolded and denatured states. Analysis of the resulting ensembles has suggested important implications of this fluctuating structure for folding, aggregation and binding.

Thermodynamics of protein denatured states

Recent work on the thermodynamics of protein denatured states is providing insight into the stability of residual structure and the conformational constraints that affect the disordered states of proteins. Current data from native state hydrogen exchange and the pH dependence of protein stability indicate that residual structure can modulate the stability of the denatured state by up to 4 kcal mol(-1). NMR structural data have emphasized the role of hydrophobic clusters in stabilizing denatured state residual structures, however recent results indicate that electrostatic interactions, both favorable and unfavorable, are also important modulators of the stability of the denatured state. Thermodynamics methods that take advantage of histidine-heme ligation chemistry have also been developed to probe the conformational constraints that act on denatured states. These methods have provided insights into the role of excluded volume, chain stiffness, and loop persistence in modulating the conformational preferences of highly disordered proteins. New insights into protein folding and novel methods to manipulate protein stability are emerging from this work.

Characterisation of the Conformational Properties of Urea-unfolded Im7: Implications for the Early Stages of Protein Folding

The colicin immunity protein Im7 folds from its unfolded state in 6 M urea to its native four-helix structure through an on-pathway intermediate that lacks one of the helices of the native structure (helix III). In order to further characterize the folding mechanism of Im7, we have studied the conformational properties of the protein unfolded in 6 M urea in detail using heteronuclear NMR. Triple-resonance experiments with 13C/15N-labelled Im7 in 6 M urea provided almost complete resonance assignments for the backbone nuclei, and measurement of backbone 15N relaxation parameters allowed dynamic ordering of the unfolded polypeptide chain to be investigated. Reduced spectral density mapping and fitting backbone R2 relaxation rates to a polymer dynamics model identified four clusters of interacting residues, each predicted by the average area buried upon folding for each residue. Chemical shift analyses and measurement of NOEs detected with a long mixing-time 1H-1H-15N NOESY-HSQC spectrum confirmed the formation of four clusters. Each cluster of interacting side-chains in urea-unfolded Im7 occurs in a region of the protein that forms a helix in the protein, with the largest clusters being associated with the three long helices that are formed in the on-pathway folding intermediate, whilst the smallest cluster forms a helix only in the native state. NMR studies of a Phe15Ala Im7 variant and a protein in which residues 51-56 are replaced by three glycine residues (H3G3 Im7*), indicated that the clusters do not interact with each other, possibly because they are solvated by urea, as indicated by analysis of NOEs between the protein and the solvent. Based on these data, we suggest that dilution of the chaotrope to initiate refolding will result in collapse of the clusters, leading to the formation of persistent helical structure and the generation of the three-helix folding intermediate.

Characterization of the residual structure in the unfolded state of the Δ131Δ fragment of staphylococcal nuclease

The determination of the conformational preferences in unfolded states of proteins constitutes an important challenge in structural biology. We use inter-residue distances estimated from site-directed spin-labeling NMR experimental measurements as ensemble-averaged restraints in all-atom molecular dynamics simulations to characterise the residual structure of the Delta131Delta fragment of staphylococcal nuclease under physiological conditions. Our findings indicate that Delta131Delta under these conditions shows a tendency to form transiently hydrophobic clusters similar to those present in the native state of wild-type staphylococcal nuclease. Only rarely, however, all these interactions are simultaneously realized to generate conformations with an overall native topology.

Motional properties of unfolded ubiquitin: a model for a random coil protein

The characterization of unfolded states of proteins has recently attracted considerable interest, as the residual structure present in these states may play a crucial role in determining their folding and misfolding behavior. Here, we investigated the dynamics in the denatured state of ubiquitin in 8 M urea at pH2. Under these conditions, ubiquitin does not have any detectable local residual structure, and uniform 15N relaxation rates along the sequence indicate the absence of motional restrictions caused by residual secondary structure and/or long-range interactions. A comparison of different models to predict relaxation data in unfolded proteins suggests that the subnanosecond dynamics in unfolded states depend on segmental motions only and do not show a dependence on the residue type but for proline and glycine residues.

Native state energetics of the Src SH2 domain: evidence for a partially structured state in the denatured ensemble

We have defined the free-energy profile of the Src SH2 domain using a variety of biophysical techniques. Equilibrium and kinetic experiments monitored by tryptophan fluorescence show that Src SH2 is quite stable and folds rapidly by a two-state mechanism, without populating any intermediates. Native state hydrogen-deuterium exchange confirms this two-state behavior; we detect no cooperative partially unfolded forms in equilibrium with the native conformation under any conditions. Interestingly, the apparent stability of the protein from hydrogen exchange is 2 kcal/mol greater than the stability determined by both equilibrium and kinetic studies followed by fluorescence. Native-state proteolysis demonstrates that this "super protection" does not result from a deviation from the linear extrapolation model used to fit the fluorescence data. Instead, it likely arises from a notable compaction in the unfolded state under native conditions, resulting in an ensemble of conformations with substantial solvent exposure of side chains and flexible regions sensitive to proteolysis, but backbone amides that exchange with solvent approximately 30-fold slower than would be expected for a random coil. The apparently simple behavior of Src SH2 in traditional unfolding studies masks the significant complexity present in the denatured-state ensemble.

The pH-dependence of amide chemical shift of Asp/Glu reflects its pKa in intrinsically disordered proteins with only local interactions

Detailed knowledge of the pH-dependence of ionizable residues in both folded and unfolded states of proteins is essential to understand the role of electrostatics in protein folding and stability. The reassembly of E. coli Thioredoxin (Trx) by complementation of its two disordered fragments (1-37/38-108) provides a folded heterodimer in equilibrium with its unfolded state which, based on circular dichroism and NMR spectroscopy, consists of two unfolded monomers. To gain insight into the role of electrostatics in protein folding and stability, we compared the pH-dependence of the carboxylate sidechain chemical shift of each Asp/Glu against that of its backbone amide chemical shift in the unfolded heterodimer. We monitored via C(CO)NH experiments four Asp and four Glu in fragments 38 to 108 (C37) of Trx in the pH range from 2.0 to 7.0 and compared them with results from (1)H(15)N HSQC experiments [Pujato et al., Biophys. J., 89 (2005) 3293-3302]. The (1)H(15)N HSQC analysis indicates two segments with quite distinct behavior: (A) a segment from Ala57 to Ala108 in which ionizable residues have up to three contiguous neighbors with pH-dependent backbone amide shifts, and (B) a segment of fifteen contiguous pH-dependent backbone amide shifts (Leu42 to Val56) in which two Asp and two Glu are implicated in medium range interactions. In all cases, the titration curves are simple modified sigmoidals from which a pH-midpoint (pH(m)) can be obtained by fitting. In segment A, the pH(m) of a given backbone amide of Asp/Glu mirrors within 0.15 pH-units that of its carboxylate sidechain (i.e., the pK(a)). In contrast, segment B shows significant differences with absolute values of 0.46 and 0.74 pH-units for Asp and Glu, respectively. The dispersion in the pH(m) of the backbone amide of Asp/Glu is also different in the two segments. Segment A shows a dispersion of 0.31 and 0.17 pH-units for Asp and Glu, respectively. Segment B shows a substantially larger dispersion (0.50 and 1.08 pH-units for Asp and Glu, respectively). In both segments, the dispersion in the pH(m) of its backbone amide is larger than in the pK(a) of the carboxylate sidechain (the latter is only 0.17 and 0.52 pH-units for Asp and Glu, respectively). Our results indicate that the pH(m) of the backbone amide chemical shift of Asp/Glu in a disordered polypeptide segment is a good predictor of its pK(a) whenever there are none or few neighboring backbone amides with similar pH-dependence.

Structural interpretation of hydrogen exchange protection factors in proteins: characterization of the native state fluctuations of CI2

Protection factors obtained from equilibrium hydrogen exchange experiments are an important source of structural information on both native and nonnative states of proteins. We present a method for determining ensembles of protein structures by using hydrogen exchange data as restraints in molecular dynamics simulations in conjunction with an empirical force-field. The method is applied to determine the ensemble of structures representing the native state of chymotrypsin inhibitor 2 (CI2), including the rare, large fluctuations responsible for hydrogen exchange.

Local structural preferences and dynamics restrictions in the urea-denatured state of SUMO-1: NMR characterization

We have investigated by multidimensional NMR the structural and dynamic characteristics of the urea-denatured state of activated SUMO-1, a 97-residue protein belonging to the growing family of ubiquitin-like proteins involved in post-translational modifications. Complete backbone amide and 15N resonance assignments were obtained in the denatured state by using HNN and HN(C)N experiments. These enabled other proton assignments from TOCSY-HSQC spectra. Secondary Halpha chemical shifts and 1H-1H NOE indicate that the protein chain in the denatured state has structural preferences in the broad beta-domain for many residues. Several of these are seen to populate the (phi,psi) space belonging to polyproline II structure. Although there is no evidence for any persistent structures, many contiguous stretches of three or more residues exhibit structural propensities suggesting possibilities of short-range transient structure formation. The hetero-nuclear 1H-15N NOEs are extremely weak for most residues, except for a few at the C-terminal, and the 15N relaxation rates show sequence-wise variation. Some of the regions of slow motions coincide with those of structural preferences and these are interspersed by highly flexible residues. The implications of these observations for the early folding events starting from the urea-denatured state of activated SUMO-1 have been discussed.

Human lactoferricin is partially folded in aqueous solution and is better stabilized in a membrane mimetic solvent

Lactoferricins are highly basic bioactive peptides that are released in the stomach through proteolytic cleavage of various lactoferrin proteins. Here we have determined the solution structure of human lactoferricin (LfcinH) by conventional two-dimensional nuclear magnetic resonance methods in both aqueous solution and a membrane mimetic solvent. Unlike the 25-residue bovine lactoferricin (LfcinB), which adopts a somewhat distorted antiparallel beta sheet, the longer LfcinH peptide shows a helical content from Gln14 to Lys29 in the membrane mimetic solvent but a nonexistent beta-sheet character in either the N- or C-terminal regions of the peptide. The helical characteristic of the LfcinH peptide resembles the conformation that this region adopts in the crystal structure of the intact protein. The LfcinH structure determined in aqueous solution displays a nascent helix in the form of a coiled conformation in the region from Gln14 to Lys29. Numerous hydrophobic interactions create the basis for the better-defined overall structure observed in the membrane mimetic solvent. The 49-residue LfcinH peptide isolated for these studies was found to be slightly longer than previously reported peptide preparations and was found to have an intact peptide bond between residues Ala11 and Val12. The distinct solution structures of LfcinH and LfcinB represent a novel difference in the physical properties of these two peptides, which contributes to their unique physiological activities.

pH dependence of amide chemical shifts in natively disordered polypeptides detects medium-range interactions with ionizable residues

A growing number of natively disordered proteins undergo a folding/binding process that is essential for their biological function. An interesting question is whether these proteins have incompletely solvated regions that drive the folding/binding process. Although the presence of predominantly hydrophobic buried regions can be easily ascertained by high-sensitivity differential scanning calorimetry analysis, the identification of those residues implicated in the burial requires NMR analysis. We have selected a partially solvated natively disordered fragment of Escherichia coli, thioredoxin, C37 (38-108), for full NMR spectral assignment. The secondary chemical shifts, temperature coefficients, and relaxation rates (R(1) and R(2)) of this fragment indicate the presence of a flexible backbone without a stable hydrogen bond network near neutral pH. (1)H-(15)N heteronuclear single quantum coherence analysis of the pH dependence of amide chemical shifts in fragment C37 within pH 2.0 and 7.0 suggests the presence of interactions between nonionizable residues and the carboxylate groups of four Asp and four Glu residues. The pH midpoints (pH(m)) of the amides in the ionizable residues (Asp or Glu) and, consequently, the shifts in the pH(m) (DeltapH(m)) of these residues with respect to model tetrapeptides, are sequence-dependent; and the nonionizable residues that show pH dependence cluster around the ionizable ones. The same pH dependence has been observed in two fragments: M37 (38-73) and C73 (74-108), ruling out the participation of long-range interactions. Our studies indicate the presence of a 15-residue pH-dependent segment with the highest density of ionizable sites in the disordered ensembles of fragments C37 and M37. The observed correlations between ionizable and nonionizable residues in this segment suggest the organization of the backbone and side chains through local and medium-range interactions up to nine residues apart, in contrast to only a few interactions in fragment C73. These results agree qualitatively with the predominantly hydrophobic buried surface detected only in fragments C37 and M37 by highly sensitive differential scanning calorimetry analysis. This work offers a sensitive and rapid new tool to obtain clues about local and nonlocal interactions between ionizable and nonionizable residues in the growing family of natively disordered small proteins with full NMR assignments.