Protein global fold determination using site-directed spin and isotope labeling

We describe a simple experimental approach for the rapid determination of protein global folds. This strategy utilizes site-directed spin labeling (SDSL) in combination with isotope enrichment to determine long-range distance restraints between amide protons and the unpaired electron of a nitroxide spin label using the paramagnetic effect on relaxation rates. The precision and accuracy of calculating a protein global fold from only paramagnetic effects have been demonstrated on barnase, a well-characterized protein. Two monocysteine derivatives of barnase, (H102C) and (H102A/Q15C), were 15N enriched, and the paramagnetic nitroxide spin label, MTSSL, attached to the single Cys residue of each. Measurement of amide 1H longitudinal relaxation times, in both the oxidized and reduced states, allowed the determination of the paramagnetic contribution to the relaxation processes. Correlation times were obtained from the frequency dependence of these relaxation processes at 800, 600, and 500 MHz. Distances in the range of 8 to 35 A were calculated from the magnitude of the paramagnetic contribution to the relaxation processes and individual amide 1H correlation times. Distance restraints from the nitroxide spin to amide protons were used as restraints in structure calculations. Using nitroxide to amide 1H distances as long-range restraints and known secondary structure restraints, barnase global folds were calculated having backbone RMSDs <3 A from the crystal structure. This approach makes it possible to rapidly obtain the overall topology of a protein using a limited number of paramagnetic distance restraints.     

Hammond behavior versus ground state effects in protein folding: evidence for narrow free energy barriers and residual structure in unfolded states

Apparent transition state movement upon mutation or changes in solvent conditions is frequently observed in protein folding and is often interpreted in terms of Hammond behavior. This led to the conclusion that barrier regions in protein folding are broad maxima on the free energy landscape. Here, we use the concept of self-interaction and cross-interaction parameters to test experimental data of 21 well-characterized proteins for Hammond behavior. This allows us to characterize the origin of transition state movements along different reaction coordinates. Only one of the 21 proteins shows a small but coherent transition state movement in agreement with the Hammond postulate. In most proteins the structure of the transition state is insensitive to changes in protein stability. The apparent change in the position of the transition state upon mutation, which is frequently observed in phi-value analysis, is in most cases due to ground-state effects caused by structural changes in the unfolded state. This argues for significant residual structure in unfolded polypeptide chains of many proteins. Disruption of these residual interactions by mutation often leads to decreased folding rates, which implies that these interactions are still present in the transition state. The failure to detect Hammond behavior shows that the free energy barriers encountered by a folding polypeptide chain are generally rather narrow and robust maxima for all experimentally explorable reaction coordinates.     

Dynamics in the unfolded state of beta2-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.     

Mapping long-range contacts in a highly unfolded protein

Insights into the earliest events in protein folding can be obtained by analysis of the conformational propensities of unfolded or partly folded states. The structure of the acid-unfolded state of apomyoglobin has been characterized using paramagnetic spin labeling and NMR. Nitroxide side-chains, introduced by coupling to mutant cysteine residues at positions 18, 77, and 133, were used as probes of chain compaction and long-range tertiary contacts. Significant interactions are observed within and between the N and C termini, while the central region of the polypeptide chain behaves as a random polymer. Even in this highly denatured form, the protein samples transient compact states in which there are native-like contacts between the N and C-terminal regions.     

Direct characterization of the folded, unfolded and urea-denatured states of the C-terminal domain of the ribosomal protein L9

The stability of the isolated C-terminal domain of the ribosomal protein L9 (CTL9) is strongly dependent upon pH. Below pH 4.2, the folded and unfolded states are both populated significantly. Their interconversion is slow on the NMR chemical shift time-scale and separate, well-resolved resonances from each state are observed. This allows the hydrodynamic properties of both states to be studied under identical conditions by using pulse field gradient NMR experiments. Hydrodynamic radii of the folded, unfolded and urea denatured protein molecules at pD 3.8 have been derived. The acid-denatured protein has a significantly smaller hydrodynamic radius, 28.2A, compared to that of the urea-denatured protein, which is 33.6A at pD 3.8. Far-UV CD spectra show that there is more residual secondary structure retained in the acid-denatured ensemble than in the urea-denatured one. ANS binding experiments and analysis of the CD data show that this acid-denatured species is not a molten globule state. Diffusion measurements of CTL9 were conducted over the pD range from 2.1 to 7.0. The hydrodynamic radii of both the folded and the acid-unfolded protein start to increase below pD 4, with the radius of hydration of the acid-unfolded state increasing from 25.1A at pD 4.2 to 33.5A at pD 2.1. The hydrodynamic radius of the urea-denatured protein is much less sensitive to pH. The unfolded protein at pD 2.1, no urea, has almost the same hydrodynamic radius as the urea-denatured protein at pD 3.8. The CD spectra, however, show significant differences in residual secondary structure, and the acid-denatured state contains more structure.     

Residual interactions in unfolded bile acid-binding protein by 19F NMR

The folding initiation mechanism of human bile acid-binding protein (BABP) has been examined by (19) F NMR. Equilibrium unfolding studies of BABP labeled with fluorine at all eight of its phenylalanine residues showed that at least two sites experience changes in solvent exposure at high denaturant concentrations. Peak assignments were made by site-specific 4FPhe incorporation. The resonances for proteins specifically labeled at Phe17, Phe47, and Phe63 showed changes in chemical shift at denaturant concentrations at which the remaining five phenylalanine residues appear to be fully solvent-exposed. Phe17 is a helical residue that was not expected to participate in a folding initiation site. Phe47 and Phe63 form part of a hydrophobic core region that may be conserved as a site for folding initiation in the intracellular lipid-binding protein family.     

Helix formation and the unfolded state of a 52-residue helical protein

A growing class of proteins in biological processes has been found to be unfolded on isolation under normal solution conditions. We have used NMR spectroscopy to characterize the structural and dynamic properties of the unfolded and partially folded states of a 52-residue alanine-rich protein (Ala-14) at temperatures from -5 degrees C to 40 degrees C. At 40 degrees C, alanine residues in Ala-14 adopt phi and psi angles, consistent with a significant ensemble population of polyproline II conformation. Analysis of relaxation rates in the protein reveals that a series of residues, Gln 35-Ala 36-Ala 37-Lys 38-Asp 39-Asp 40-Ala 41-Ala 42, displays slow motional dynamics at both -5 degrees C and 40 degrees C. Temperature-dependent chemical shift changes indicate that this region is the site of helix initiation. The remaining N-terminal residues become increasingly dynamic as they extend from the nucleation site. The C terminus remains dynamic and changes less with temperature, indicating it is relatively unstructured. Ala-14 provides a high-resolution portrait of the unfolded state and the process of helix nucleation and propagation in the absence of tertiary contacts, information that bears on early events in protein folding.     

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 ¹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 ¹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 ¹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.     

Rapid exploration of the folding topology of helical membrane proteins using paramagnetic perturbation

An understanding of the folding states of α-helical membrane proteins in detergent systems is important for functional and structural studies of these proteins. Here, we present a rapid and simple method for identification of the folding topology and assembly of transmembrane helices using paramagnetic perturbation in nuclear magnetic resonance spectroscopy. By monitoring the perturbation of signals from glycine residues located at specific sites, the folding topology and the assembly of transmembrane helices of membrane proteins were easily identified without time-consuming backbone assignment. This method is validated with Mistic (membrane-integrating sequence for translation of integral membrane protein constructs) of known structure as a reference protein. The folding topologies of two bacterial histidine kinase membrane proteins (SCO3062 and YbdK) were investigated by this method in dodecyl phosphocholine (DPC) micelles. Combing with analytical ultracentrifugation, we identified that the transmembrane domain of YbdK is present as a parallel dimer in DPC micelle. In contrast, the interaction of transmembrane domain of SCO3062 is not maintained in DPC micelle due to disruption of native structure of the periplasmic domain by DPC micelle.     

Determination of ultrafast protein folding rates from loop formation dynamics

Quenching of the triplet state of tryptophan by contact with cysteine can be used to measure the kinetics of loop formation in unfolded proteins. Here we show that cysteine quenching dynamics also provide a novel method for measuring folding rates when the exchange between folded and unfolded states is faster than the unquenched triplet lifetime (approximately 100 micros). We use this technique to investigate folding/unfolding kinetics of the 35 residue headpiece subdomain of the protein villin, which contains a single tryptophan residue and was engineered to contain a cysteine residue at the N terminus. At intermediate concentrations of denaturant the time-course of the triplet decay consists of two relaxations, the rates and amplitudes of which reveal the fast kinetics for folding and unfolding of this protein. The folding rates extracted using a simple kinetic model are close to those reported previously from laser-induced temperature-jump experiments that employ the change in tryptophan fluorescence as a probe. However, the results differ significantly from those reported from dynamic NMR line shape analysis on a variant with methionine at the N terminus, an issue that remains to be resolved. The analysis of the triplet quenching kinetics also shows that the quenching rates in the unfolded state increase with decreasing denaturant concentration, indicating a compaction of the unfolded protein.     

Calcium structural transition of human cardiac troponin C in reconstituted muscle fibres as studied by site-directed spin labelling

The in situ structure of human cardiac troponin C (hcTnC) has been studied with site-directed, spin labelling, electron paramagnetic resonance (SDSL-EPR). Analysis of the in situ structures of hcTnC is essential for elucidating the molecular mechanism behind its Ca(2+)-sensitive regulation. We prepared two hcTnC mutants (C35S and C84S) containing one native cysteine residue (84 and 35, respectively) for spin labelling. The mutants were labelled with a methane thiosulfonate spin label (MTSSL) and the TnC was reconstituted into permeabilized muscle fibres. The mobility of Cys84-MTSSL changed markedly after addition of Ca2+, while that of the Cys35 residue did not change in the monomer state or in fibres. The rotational correlation time of Cys84-MTSSL decreased from 32ns to 13ns upon Ca(2+)-binding in the monomer state, whereas in fibres the spectrum of Cys84-MTSSL was resolved into mobile (16ns) and immobile (35ns) components and the addition of Ca2+ increased the immobile component. Moreover, the accessibility of Cys84-MTSSL to molecular oxygen increased slightly in the presence of Ca2+. These data suggest that Cys35 remains in the same location regardless of the addition of Ca2+, whereas Cys84 is located at the position that interacts with B and C helices of hcTnC and interacts with troponin I (TnI) at high concentrations of Ca2+. We determined the distances between Cys35 and Cys84 by measuring pulsed electron-electron double resonance spectra. The distances were 26.0 angstroms and 27.2 angstroms in the monomer state and in fibres, respectively, and the addition of Ca2+ decreased the distance to 23.2 angstroms in fibres but only slightly in the monomer state, showing that Ca2+ binding to the N-domain of hcTnC induced a larger structural change in muscle fibres than in the monomer state.     

Backbone structure of a small helical integral membrane protein: A unique structural characterization

The structural characterization of small integral membrane proteins pose a significant challenge for structural biology because of the multitude of molecular interactions between the protein and its heterogeneous environment. Here, the three‐dimensional backbone structure of Rv1761c from Mycobacterium tuberculosis has been characterized using solution NMR spectroscopy and dodecylphosphocholine (DPC) micelles as a membrane mimetic environment. This 127 residue single transmembrane helix protein has a significant (10 kDa) C‐terminal extramembranous domain. Five hundred and ninety distance, backbone dihedral, and orientational restraints were employed resulting in a 1.16 Å rmsd backbone structure with a transmembrane domain defined at 0.40 Å. The structure determination approach utilized residual dipolar coupling orientation data from partially aligned samples, long‐range paramagnetic relaxation enhancement derived distances, and dihedral restraints from chemical shift indices to determine the global fold. This structural model of Rv1761c displays some influences by the membrane mimetic illustrating that the structure of these membrane proteins is dictated by a combination of the amino acid sequence and the protein's environment. These results demonstrate both the efficacy of the structural approach and the necessity to consider the biophysical properties of membrane mimetics when interpreting structural data of integral membrane proteins and, in particular, small integral membrane proteins.     

The B domain of protein A retains residual structures in 6 M guanidinium chloride as revealed by hydrogen/deuterium‐exchange NMR spectroscopy

The characterization of residual structures persistent in unfolded proteins is an important issue in studies of protein folding, because the residual structures present, if any, may form a folding initiation site and guide the subsequent folding reactions. Here, we studied the residual structures of the isolated B domain (BDPA) of staphylococcal protein A in 6 M guanidinium chloride. BDPA is a small three‐helix‐bundle protein, and until recently its folding/unfolding reaction has been treated as a simple two‐state process between the native and the fully unfolded states. We employed a dimethylsulfoxide (DMSO)‐quenched hydrogen/deuterium (H/D)‐exchange 2D NMR techniques with the use of spin desalting columns, which allowed us to investigate the H/D‐exchange behavior of individually identified peptide amide (NH) protons. We obtained H/D‐exchange protection factors of the 21 NH protons that form an α‐helical hydrogen bond in the native structure, and the majority of these NH protons were significantly protected with a protection factor of 2.0–5.2 in 6 M guanidinium chloride, strongly suggesting that these weakly protected NH protons form much stronger hydrogen bonds under native folding conditions. The results can be used to deduce the structure of an early folding intermediate, when such an intermediate is shown by other methods. Among three native helical regions, the third helix in the C‐terminal side was highly protected and stabilized by side‐chain salt bridges, probably acting as the folding initiation site of BDPA. The present results are discussed in relation to previous experimental and computational findings on the folding mechanisms of BDPA.     

Dynamics of unfolded polypeptide chains as model for the earliest steps in protein folding

The rate of formation of intramolecular interactions in unfolded proteins determines how fast conformational space can be explored during folding. Characterization of the dynamics of unfolded proteins is therefore essential for the understanding of the earliest steps in protein folding. We used triplet-triplet energy transfer to measure formation of intrachain contacts in different unfolded polypeptide chains. The time constants (1/k) for contact formation over short distances are almost independent of chain length, with a maximum value of about 5 ns for flexible glycine-rich chains and of 12 ns for stiffer chains. The rates of contact formation over longer distances decrease with increasing chain length, indicating different rate-limiting steps for motions over short and long chain segments. The effect of the amino acid sequence on local chain dynamics was probed by using a series of host-guest peptides. Formation of local contacts is only sixfold slower around the stiffest amino acid (proline) compared to the most flexible amino acid (glycine). Good solvents for polypeptide chains like EtOH, GdmCl and urea were found to slow intrachain diffusion and to decrease chain stiffness. These data allow us to determine the time constants for formation of the earliest intrachain contacts during protein folding.     

Cross correlation rates between Curie spin and dipole-dipole relaxation in paramagnetic proteins: the case of cerium substituted calbindin D9k

Cross correlation rates between Curie spin relaxation and H-N dipole-dipole coupling (^CS,DD _HM,HN) have been determined for a calcium binding protein, Calbindin D_9k, in which one of the two calcium ions is substituted with cerium(III). ^CS,DD _HM,HNvalues depend on both the metal-to-proton distances and the M-H-N angles and can be used as an additional constraint in order to refine the solution structure of paramagnetic metalloproteins. For this purpose, we have implemented a new module (CCR-DYANA) in a version of the program DYANA (PARAMAGNETIC-DYANA), which can be used together with other paramagnetism-based constraints such as pseudocontact shifts, residual dipolar couplings and hyperfine based Karplus relationships. This integrated structure calculation protocol has the advantage that different paramagnetic-based constraints are treated by the same algorithm in a way that the efficiency of each class of constraints can be analyzed and compared.     

Hydrophobic clustering in nonnative states of a protein: interpretation of chemical shifts in NMR spectra of denatured states of lysozyme.

Chemical shifts of resonances of specific protons in the 1H NMR spectrum of thermally denatured hen lysozyme have been determined by exchange correlation with assigned native state resonances in 2D NOESY spectra obtained under conditions where the two states are interconverting. There are subtle but widespread deviations of the measured shifts from the values which would be anticipated for a random coil; in the case of side chain protons these are virtually all net upfield shifts and it is shown that this may be the averaged effect of interactions with aromatic rings in a partially collapsed denatured state. In a very few cases, notably that of two sequential tryptophan residues, it is possible to interpret these effects in terms of specific, local interresidue interactions. Generally, however, there is no correlation with either native state shift perturbations or with sequence proximity to aromatic groups. Diminution of most of the residual shift perturbations on reduction of the disulfide cross-links confirms that they are not simply effects of residues adjacent in the sequence. Similar effects of chemical denaturants, with the disulfides intact, demonstrate that the shift perturbations reflect an enhanced tendency to side chain clustering in the thermally denatured state. The temperature dependences of the shift perturbations suggest that this clustering is noncooperative and is driven by small, favorable enthalpy changes. While the extent of conformational averaging is clearly much greater than that observed for a homologous protein, alpha-lactalbumin, in its partially folded "molten globule" state, the results clearly show that thermally denatured lysozyme differs substantially from a random coil, principally in that it is partially hydrophobically collapsed.     

Kinetics of intramolecular contact formation in a denatured protein

Quenching of the triplet state of tryptophan by cysteine has provided a new tool for measuring the rate of forming a specific intramolecular contact in disordered polypeptides. Here, we use this technique to investigate contact formation in the denatured state of CspTm, a small cold-shock protein from Thermotoga maritima, engineered to contain a single tryptophan residue (W29) and a single cysteine residue at the C terminus (C67). At all concentrations of denaturant, the decay rate of the W29 triplet of the unfolded protein is more than tenfold faster than the rate observed for the native protein ( approximately 10(4)s(-1)). Experiments on the unfolded protein without the added C-terminal cysteine residue show that this faster rate results entirely from contact quenching by C67. The quenching rate in the unfolded state by C67 increases at concentrations of denaturant that favor folding, indicating a compaction of the unfolded protein as observed previously in single-molecule F?rster resonance energy transfer (FRET) experiments.     

Relaxation data in NMR structure determination: model calculations for the lysozyme-Gd3+ complex

The effect of including paramagnetic relaxation data as additional restraints in the determination of protein tertiary structures from NMR data has been explored by a systematic series of model calculations. The system used for testing the method was the 2.0 A resolution tetragonal crystal structure of hen egg white lysozyme (129 amino acid residues) and structures were generated using a version of the hybrid "distance geometry-dynamic simulated annealing" procedure. A limited set of 769 NOEs was used as restraints in all the calculations; the strengths of these were categorized into three classes on the basis of distances observed in the crystal structure. The values of 50 phi angles were also restrained on the basis of amide-alpha coupling constants calculated from the X-ray structure. Five sets of 12 structures were determined using differing sets of paramagnetic relaxation data as restraints additional to those involving the NOE and coupling constant data. The paramagnetic relaxation data were modeled on the basis of the distances of defined protons from the crystallographic binding site of Gd3+ in lysozyme. Analysis of the results showed that the relaxation data significantly improved the correspondence between the set of generated structures and the crystal structure, and that the more well defined the relaxation data, the more significant the improvement in the quality of the structures. The results suggest that the inclusion of paramagnetic relaxation restraints could be of significant value for the experimental determination of protein structures from NMR data.     

Polyproline II helical structure in protein unfolded states: lysine peptides revisited

The left-handed polyproline II (PPII) helix gives rise to a circular dichroism spectrum that is remarkably similar to that of unfolded proteins. This similarity has been used as the basis for the hypothesis that unfolded proteins possess considerable PPII helical content. It has long been known that homopolymers of lysine adopt the PPII helical conformation at neutral pH, presumably a result of electrostatic repulsion between side chains. It is shown here that a seven-residue lysine peptide also adopts the PPII conformation. In contrast with homopolymers of lysine, this short peptide is shown to retain PPII helical character under conditions in which side-chain charges are heavily screened or even neutralized. The most plausible explanation for these observations is that the peptide backbone favors the PPII conformation to maximize favorable interactions with solvent. These data are evidence that unfolded proteins do indeed possess PPII content, indicating that the ensemble of unfolded states is significantly smaller than is commonly assumed.