Hydrogen exchange in native and alcohol forms of ubiquitin

DMSO-Quenched H/D-Exchange 2D NMR Spectroscopy and Its Applications in Protein Science

Hydrogen/deuterium (H/D) exchange combined with two-dimensional (2D) NMR spectroscopy has been widely used for studying the structure, stability, and dynamics of proteins. When we apply the H/D-exchange method to investigate non-native states of proteins such as equilibrium and kinetic folding intermediates, H/D-exchange quenching techniques are indispensable, because the exchange reaction is usually too fast to follow by 2D NMR. In this article, we will describe the dimethylsulfoxide (DMSO)-quenched H/D-exchange method and its applications in protein science. In this method, the H/D-exchange buffer is replaced by an aprotic DMSO solution, which quenches the exchange reaction. We have improved the DMSO-quenched method by using spin desalting columns, which are used for medium exchange from the H/D-exchange buffer to the DMSO solution. This improvement has allowed us to monitor the H/D exchange of proteins at a high concentration of salts or denaturants. We describe methodological details of the improved DMSO-quenched method and present a case study using the improved method on the H/D-exchange behavior of unfolded human ubiquitin in 6 M guanidinium chloride.

Protein structure prediction using residue-resolved protection factors from hydrogen-deuterium exchange NMR

Hydrogen-deuterium exchange (HDX) measured by nuclear magnetic resonance (NMR) provides structural information for proteins relating to solvent accessibility and flexibility. While this structural information is beneficial, the data cannot be used exclusively to elucidate structures. However, the structural information provided by the HDX-NMR data can be supplemented by computational methods. In previous work, we developed an algorithm in Rosetta to predict structures using qualitative HDX-NMR data (categories of exchange rate). Here we expand on the effort, and utilize quantitative protection factors (PFs) from HDX-NMR for structure prediction. From observed correlations between PFs and solvent accessibility/flexibility measures, we present a scoring function to quantify the agreement with HDX data. Using a benchmark set of 10 proteins, an average improvement of 5.13 Å in root-mean-square deviation (RMSD) is observed for cases of inaccurate Rosetta predictions. Ultimately, seven out of 10 predictions are accurate without including HDX data, and nine out of 10 are accurate when using our PF-based HDX score.

Variable-Temperature Electrospray Ionization for Temperature-Dependent Folding/Refolding Reactions of Proteins and Ligand Binding

Stabilities and structure(s) of proteins are directly coupled to their local environment or Gibbs free energy landscape as defined by solvent, temperature, pressure, and concentration. Solution pH, ionic strength, cofactors, chemical chaperones, and osmolytes perturb the chemical potential and induce further changes in structure, stability, and function. At present, no single analytical technique can monitor these effects in a single measurement. Mass spectrometry and ion mobility-mass spectrometry play increasingly essential roles in studies of proteins, protein complexes, and even membrane protein complexes; however, with few exceptions, the effects of the solution temperature on the stability and structure(s) of analytes have not been thoroughly investigated. Here, we describe a new variable-temperature electrospray ionization (vT-ESI) source that utilizes a thermoelectric chip to cool and heat the solution contained within the static ESI emitter. This design allows for solution temperatures to be varied from ∼5 to 98 °C with short equilibration times (<2 min) between precisely controlled temperature changes. The performance of the apparatus for vT-ESI-mass spectrometry and vT-ESI-ion mobility-mass spectrometry studies of cold- and heat-folding reactions is demonstrated using ubiquitin and frataxin. Instrument performance for studies on temperature-dependent ligand binding is shown using the chaperonin GroEL.

Nonthermal excitation effects mediated by sub-terahertz radiation on hydrogen exchange in ubiquitin

Water dynamics in the hydration layers of biomolecules play crucial roles in a wide range of biological functions. A hydrated protein contains multiple components of diffusional and vibrational dynamics of water and protein, which may be coupled at ∼0.1-THz frequency (10-ps timescale) at room temperature. However, the microscopic description of biomolecular functions based on various modes of protein-water-coupled motions remains elusive. A novel approach for perturbing the hydration dynamics in the subterahertz frequency range and probing them at the atomic level is therefore warranted. In this study, we investigated the effect of klystron-based, intense 0.1-THz excitation on the slow dynamics of ubiquitin using NMR-based measurements of hydrogen-deuterium exchange. We demonstrated that the subterahertz irradiation accelerated the hydrogen-deuterium exchange of the amides located in the interior of the protein and hydrophobic surfaces while decelerating this exchange in the amides located in the surface loop and short 3₁₀ helix regions. This subterahertz-radiation-induced effect was qualitatively contradictory to the increased-temperature-induced effect. Our results suggest that the heterogeneous water dynamics occurring at the protein-water interface include components that are nonthermally excited by the subterahertz radiation. Such subterahertz-excited components may be linked to the slow function-related dynamics of the protein.

Residual Structure of Unfolded Ubiquitin as Revealed by Hydrogen/Deuterium-Exchange 2D NMR

The characterization of residual structures persistent in unfolded proteins in concentrated denaturant solution is currently an important issue in studies of protein folding because the residual structure present, if any, in the unfolded state may form a folding initiation site and guide the subsequent folding reactions. Here, we studied the hydrogen/deuterium (H/D)-exchange behavior of unfolded human ubiquitin in 6 M guanidinium chloride. We employed a dimethylsulfoxide (DMSO)-quenched H/D-exchange NMR technique with the use of spin desalting columns, which allowed us to perform a quick medium exchange from 6 M guanidinium chloride to a quenching DMSO solution. Based on the backbone resonance assignment of ubiquitin in the DMSO solution, we successfully investigated the H/D-exchange kinetics of 60 identified peptide amide groups in the ubiquitin sequence. Although a majority of these amide groups were not protected, certain amide groups involved in a middle helix (residues 23-34) and an N-terminal β-hairpin (residues 2-16) were significantly protected with a protection factor of 2.1-4.2, indicating that there were residual structures in unfolded ubiquitin and that these amide groups were more than 52% hydrogen bonded in the residual structures. We show that the hydrogen-bonded residual structures in the α-helix and the β-hairpin are formed even in 6 M guanidinium chloride, suggesting that these residual structures may function as a folding initiation site to guide the subsequent folding reactions of ubiquitin.

Structural dynamics of pentapeptide repeat proteins

Pentapeptide repeat proteins (PRPs) represent a large superfamily with more than 38 000 sequences in nearly 3500 species, the majority belonging to cyanobacteria but represented among all branches of life. PRPs contain at least eight consecutive pentapeptide repeats with the consensus (A/C/S/V/T/L/I)(D/N/S/K/E/I/R)(L/F)(S/T/R/E/Q/K/V/D)(G/D/E/N/R/Q/K). PRPs fold into right-handed quadrilateral β helices, also known as repeat-five-residue (Rfr)-folds, with four consecutive pentapeptide repeats comprising a single coil, the ~90° change in polypeptide direction in square-shaped coils achieved by type I, II and IV β turns, and hydrogen bonds between coils establishing β ladders on each Rfr-fold face. PRPs are broadly categorized into group 1 and 2 involved in antibiotic resistance and group 3 currently having unknown functions. Motivated by their intriguing structures, we are investigating PRP biophysical characteristics, including Rfr-fold thermal stability, β turn and β ladder hydrogen bond amide exchange rates and backbone dynamics. Here, we present analysis of 20 ns molecular dynamics (MD) simulations and all atom normal mode analysis (aaNMA) calculations for four group 1 and group 2 and four group 3 PRPs whose structures have been determined by X-ray crystallography. The MD cross-correlation matrices and aaNMA indicated strong correlated motion between adjacent coils and weak coupled motion between coils separated by one or more intervening coils. Slow anticorrelated motions were detected between adjacent coils in aaNMA modes that we hypothesize are requisite to access exchange-competent states necessary to permit solvent exchange of amide hydrogens involved in β-ladder and β-turns hydrogen bonds, which can have lifetimes on the order of months.

Ion Mobility and Gas-Phase Covalent Labeling Study of the Structure and Reactivity of Gaseous Ubiquitin Ions Electrosprayed from Aqueous and Denaturing Solutions

Gas-phase ion/ion chemistry was coupled to ion mobility/mass spectrometry analysis to correlate the structure of gaseous ubiquitin to its solution structures with selective covalent structural probes. Collision cross section (CCS) distributions were measured to ensure the ubiquitin ions were not unfolded when they were introduced to the gas phase. Aqueous solutions stabilizing the native state of ubiquitin yielded folded ubiquitin structures with CCS values consistent with previously published literature. Denaturing solutions favored several families of unfolded conformations for most of the charge states evaluated. Gas-phase covalent labeling via ion/ion reactions was followed by collision-induced dissociation of the intact, labeled protein to determine which residues were labeled. Ubiquitin 5+ and 6+ electrosprayed from aqueous conditions were covalently modified preferentially at the lysine 29 and arginine 54 positions, indicating that elements of three-dimensional structure were maintained in the gas phase. On the other hand, most ubiquitin ions produced in denaturing conditions were labeled at various other lysine residues, likely due to the availability of additional sites following methanol- and low-pH-induced unfolding. These data support the conservation of ubiquitin structural elements in the gas phase. The research presented here provides the basis for residue-specific characterization of biomolecules in the gas phase.

Insights into the Dynamic Structural Properties of a Lanthipeptide Synthetase using Hydrogen-Deuterium Exchange Mass Spectrometry

The biosynthesis of ribosomally synthesized and post-translationally modified peptides (RiPPs) proceeds via the multistep maturation of genetically encoded precursor peptides, often catalyzed by enzymes with multiple functions and iterative activities. Recent studies have suggested that, among other factors, conformational sampling of enzyme:peptide complexes likely plays a critical role in defining the kinetics and, ultimately, the set of post-translational modifications in these systems. However, detailed characterizations of these putative conformational sampling mechanisms have not yet been possible on many RiPP biosynthetic systems. In this study, we report the first comprehensive application of hydrogen-deuterium exchange mass spectrometry (HDX-MS) to study the biophysical properties of a RiPP biosynthetic enzyme. Using the well-characterized class II lanthipeptide synthetase HalM2 as a model system, we have employed HDX-MS to demonstrate that HalM2 is indeed a highly structurally dynamic enzyme. Using this HDX-MS approach, we have identified novel precursor peptide binding elements, have uncovered long-range structural communication across the enzyme that is triggered by ligand binding and ATP hydrolysis, and have detected specific interactions between the HalM2 synthetase and the leader- and core-peptide subdomains of the modular HalA2 precursor peptide substrate. The functional relevance of the dynamic HalM2 elements discovered in this study are validated with biochemical assays and kinetic analysis of a panel of HDX-MS guided variant enzymes. Overall, the data have provided a wealth of fundamentally new information on LanM systems that will inform the rational manipulation and engineering of these impressive multifunctional catalysts. Moreover, this work highlights the broad utility of the HDX-MS platform for revealing important biophysical properties and enzyme structural dynamics that likely play a widespread role in RiPP biosynthesis.

Melting proteins confined in nanodroplets with 10.6 μm light provides clues about early steps of denaturation

Ubiquitin confined within nanodroplets was irradiated with a variable-power CO2 laser. Mass spectrometry analysis shows evidence for a protein "melting"-like transition within droplets prior to solvent evaporation and ion formation. Ion mobility spectrometry reveals that structures associated with early steps of denaturation are trapped because of short droplet lifetimes.

From Compact to String-The Role of Secondary and Tertiary Structure in Charge-Induced Unzipping of Gas-Phase Proteins

In the gas phase, protein ions can adopt a broad range of structures, which have been investigated extensively in the past using ion mobility-mass spectrometry (IM-MS)-based methods. Compact ions with low number of charges undergo a Coulomb-driven transition to partially folded species when the charge increases, and finally form extended structures with presumably little or no defined structure when the charge state is high. However, with respect to the secondary structure, IM-MS methods are essentially blind. Infrared (IR) spectroscopy, on the other hand, is sensitive to such structural details and there is increasing evidence that helices as well as β-sheet-like structures can exist in the gas phase, especially for ions in low charge states. Very recently, we showed that also the fully extended form of highly charged protein ions can adopt a distinct type of secondary structure that features a characteristic C₅-type hydrogen bond pattern. Here we use a combination of IM-MS and IR spectroscopy to further investigate the influence of the initial, native conformation on the formation of these structures. Our results indicate that when intramolecular Coulomb-repulsion is large enough to overcome the stabilization energies of the genuine secondary structure, all proteins, regardless of their sequence or native conformation, form C₅-type hydrogen bond structures. Furthermore, our results suggest that in highly charged proteins the positioning of charges along the sequence is only marginally influenced by the basicity of individual residues. Graphical Abstract ᅟ.

Characterizing the Structure and Oligomerization of Major Royal Jelly Protein 1 (MRJP1) by Mass Spectrometry and Complementary Biophysical Tools

Royal jelly (RJ) triggers the development of female honeybee larvae into queens. This effect has been attributed to the presence of major royal jelly protein 1 (MRJP1) in RJ. MRJP1 isolated from royal jelly is tightly associated with apisimin, a 54-residue α-helical peptide that promotes the noncovalent assembly of MRJP1 into multimers. No high-resolution structural data are available for these complexes, and their binding stoichiometry remains uncertain. We examined MRJP1/apisimin using a range of biophysical techniques. We also investigated the behavior of deglycosylated samples, as well as samples with reduced apisimin content. Our mass spectrometry (MS) data demonstrate that the native complexes predominantly exist in a (MRJP1₄ apisimin₄) stoichiometry. Hydrogen/deuterium exchange MS reveals that MRJP1 within these complexes is extensively disordered in the range of residues 20-265. Marginally stable secondary structure (likely antiparallel β-sheet) exists around residues 266-432. These weakly structured regions interchange with conformers that are extensively unfolded, giving rise to bimodal (EX1) isotope distributions. We propose that the native complexes have a "dimer of dimers" quaternary structure in which MRJP1 chains are bridged by apisimin. Specifically, our data suggest that apisimin acts as a linker that forms hydrophobic contacts involving the MRJP1 segment ³¹⁶VLFFGLV³²². Deglycosylation produces large soluble aggregates, highlighting the role of glycans as aggregation inhibitors. Samples with reduced apisimin content form dimeric complexes with a (MRJP1₂ apisimin₁) stoichiometry. The information uncovered in this work will help pave the way toward a better understanding of the unique physiological role played by MRJP1 during queen differentiation.

Ubiquitin utilizes an acidic surface patch to alter chromatin structure

Ubiquitylation of histone H2B, associated with gene activation, leads to chromatin decompaction through an unknown mechanism. We used a hydrogen-deuterium exchange strategy coupled with nuclear magnetic resonance spectroscopy to map the ubiquitin surface responsible for its structural effects on chromatin. Our studies revealed that a previously uncharacterized acidic patch on ubiquitin comprising residues Glu16 and Glu18 is essential for decompaction. These residues mediate promiscuous electrostatic interactions with the basic histone proteins, potentially positioning the ubiquitin moiety as a dynamic “wedge” that prevents the intimate association of neighboring nucleosomes. Using two independent cross-linking strategies and an oligomerization assay, we also showed that ubiquitin-ubiquitin contacts occur in the chromatin environment and are important for the solubilization of the chromatin polymers. Our work highlights a novel, chromatin-related aspect of the “ubiquitin code”, and sheds light on how the information rich ubiquitin modification can orchestrate different biochemical outcomes using different surface features.

Real-time HD Exchange Kinetics of Proteins from Buffered Aqueous Solution with Electrothermal Supercharging and Top-Down Tandem Mass Spectrometry

Electrothermal supercharging (ETS) with electrospray ionization produces highly charged protein ions from buffered aqueous solutions in which proteins have native folded structures. ETS increases the charge of ribonuclease A by 34%, whereas only a 6% increase in charge occurs for a reduced-alkylated form of this protein, which is unfolded and its structure is ~66% random coil in this solution. These results indicate that protein denaturation that occurs in the ESI droplets is the primary mechanism for ETS. ETS does not affect the extent of solution-phase hydrogen-deuterium exchange (HDX) that occurs for four proteins that have significantly different structures in solution, consistent with a droplet lifetime that is considerably shorter than observable rates of HDX. Rate constants for HDX of ubiquitin are obtained with a spatial resolution of ~1.3 residues with ETS and electron transfer dissociation of the 10+ charge-state using a single capillary containing a few μL of protein solution in which HDX continuously occurs. HDX protection at individual residues with ETS HDX is similar to that with reagent supercharging HDX and with solution-phase NMR, indicating that the high spray potentials required to induce ETS do not lead to HD scrambling. Graphical Abstract ᅟ.

Solution dependence of the collisional activation of ubiquitin [M + 7H](7+) ions

The solution dependence of gas-phase unfolding for ubiquitin [M + 7H](7+) ions has been studied by ion mobility spectrometry-mass spectrometry (IMS-MS). Different acidic water:methanol solutions are used to favor the native (N), more helical (A), or unfolded (U) solution states of ubiquitin. Unfolding of gas-phase ubiquitin ions is achieved by collisional heating and newly formed structures are examined by IMS. With an activation voltage of 100 V, a selected distribution of compact structures unfolds, forming three resolvable elongated states (E1-E3). The relative populations of these elongated structures depend strongly on the solution composition. Activation of compact ions from aqueous solutions known to favor N-state ubiquitin produces mostly the E1 type elongated state, whereas activation of compact ions from methanol containing solutions that populate A-state ubiquitin favors the E3 elongated state. Presumably, this difference arises because of differences in precursor ion structures emerging from solution. Thus, it appears that information about solution populations can be retained after ionization, selection, and activation to produce the elongated states. These data as well as others are discussed.

Using hydrogen-deuterium exchange to monitor protein structure in the presence of gold nanoparticles

The potential applications of protein-functionalized gold nanoparticles (AuNPs) have motivated many studies characterizing protein-AuNP interactions. However, the lack of detailed structural information has hindered our ability to understand the mechanism of protein adsorption on AuNPs. In order to determine the structural perturbations that occur during adsorption, hydrogen/deuterium exchange (HDX) of amide protons was measured for two proteins by NMR. Specifically, we measured both slow (5-300 min) and fast (10-500 ms) H/D exchange rates for GB3 and ubiquitin, two well-characterized proteins. Overall, amide exchange rates are very similar in the presence and absence of AuNPs, supporting a model where the adsorbed protein remains largely folded on the AuNP surface. Small differences in exchange rates are observed for several loop residues, suggesting that the secondary structure remains relatively rigid while loops and surface residues can experience perturbations upon binding. Strikingly, several of these residues are close to lysines, which supports a model where positive surface residues may interact favorably with AuNP-bound citrate. Because these proteins appear to remain folded on AuNP surfaces, these studies suggest that it may be possible to engineer functional AuNP-based nanoconjugates without the use of chemical linkers.

Evidence for two new solution states of ubiquitin by IMS-MS analysis

Ion mobility spectrometry coupled with mass spectrometry (IMS–MS) is used to investigate the populations of different states for ubiquitin in water:methanol solutions. In these experiments, ubiquitin is electrosprayed from 20 water:methanol (100:0 to 5:95, pH = 2) solutions, ranging from native to denaturing conditions. With an increased percentage of methanol in solution, ubiquitin ions ([M + 7H]⁷⁺ to [M + 12H]¹²⁺) show substantial variations in both charge state distributions and ion mobility distributions. Analysis of these data provides evidence for the existence of five ubiquitin states in solution: the native N state, favored in solutions of 100:0 to 70:30 water:methanol for the +7 and +8 charge states; the more helical A state and a new closely related A′ state, favored in solutions of 70:30 to 5:95 water:methanol for the +9 to +12 charge states; the unfolded U state, populated in 40:60 to 5:95 water:methanol solutions for the +8 to +10 and +12 charge states; and a new low-abundance state termed the B state, observed for 100:0 to 70:30 water:methanol solutions in the +8 to +10 and +12 charge states. The relative abundances for different states in different solutions are determined. The analysis presented here provides insight into how solution structures evolve into anhydrous conformations and demonstrates the utility of IMS–MS methods as a means of characterizing populations of conformers for proteins in solution.

Time-resolved multidimensional NMR with non-uniform sampling

Time-resolved experiments demand high resolution both in spectral dimensions and in time of the studied kinetic process. The latter requirement traditionally prohibits applications of the multidimensional experiments, which, although capable of providing invaluable information about structure and dynamics and almost unlimited spectral resolution, require too lengthy data collection. Our work shows that the problem has a solution in using modern methods of NMR data collection and signal processing. A continuous fast pulsing three-dimensional experiment is acquired using non-uniform sampling during full time of the studied reaction. High sensitivity and time-resolution of a few minutes is achieved by simultaneous processing of the full data set with the multi-dimensional decomposition. The method is verified and illustrated in realistic simulations and by measuring deuterium exchange rates of amide protons in ubiquitin. We applied the method for characterizing kinetics of in vitro phosphorylation of two tyrosine residues in an intrinsically disordered cytosolic domain of the B cell receptor protein CD79b. Signals of many residues including tyrosines in both phosphorylated and unmodified forms of CD79b are found in a heavily crowded region of 2D ¹H-¹⁵N correlation spectrum and the significantly enhanced spectral resolution provided by the 3D time-resolved approach was essential for the quantitative site-specific analysis.

Conformer-specific characterization of nonnative protein states using hydrogen exchange and top-down mass spectrometry

Characterization of structure and dynamics of nonnative protein states is important for understanding molecular mechanisms of processes as diverse as folding, binding, aggregation, and enzyme catalysis to name just a few; however, selectively probing local minima within rugged energy landscapes remains a problem. Mass spectrometry (MS) coupled with hydrogen/deuterium exchange (HDX) offers a unique advantage of being able to make a distinction among multiple protein conformers that coexist in solution; however, detailed structural interrogation of such states previously remained out of reach of HDX MS. In this work, we exploited the aforementioned unique feature of HDX MS in combination with the ability of MS to isolate narrow populations of protein ions to characterize individual protein conformers coexisting in solution in equilibrium. Subsequent fragmentation of the protein ions using electron-capture dissociation allowed us to allocate the deuterium distribution along the protein backbone, yielding a backbone-amide protection map for the selected conformer unaffected by contributions from other protein states present in solution. The method was tested with the small regulatory protein ubiquitin (Ub), which is known to form nonnative intermediate states under a variety of mildly denaturing conditions. Protection maps of these intermediate states obtained at residue-level resolution provide clear evidence that they are very similar to the so-called A-state of Ub that is formed in solutions with low pH and high alcohol. Method validation was carried out by comparing the backbone-amide protection map of native Ub with those deduced from high-resolution NMR measurements.

Mapping pH-induced protein structural changes under equilibrium conditions by pulsed oxidative labeling and mass spectrometry

Mass spectrometry (MS)-based protein conformational studies are a rapidly growing field. The characterization of partially disordered conformers is of particular interest because these species are not amenable to classical high-resolution techniques. Such equilibrium intermediates can often be populated by exposure to mildly acidic pH. Hydroxyl radical (·OH) introduces oxidative modifications at solvent-accessible side chains, while buried sites are protected. ·OH can be generated by laser photolysis of H(2)O(2) (fast photochemical oxidation of proteins-FPOP). The resulting labeling pattern can be analyzed by MS. The characterization of partially disordered intermediates usually involves comparative measurements under different solvent conditions. It can be challenging to separate structurally induced labeling changes from pH-mediated "secondary" effects. The issue of secondary effects in FPOP has received little prior attention. We demonstrate that with a proper choice of conditions (e.g., in the absence of pH-dependent ·OH scavengers) such undesired phenomena can be almost completely eliminated. Using apomyoglobin as a model system, we map the structure of an intermediate that is formed at pH 4. This species retains a highly protected helix G that is surrounded by partially protected helices A, B, and H. Our results demonstrate the utility of FPOP for the structural characterization of equilibrium intermediates. The near absence of an intrinsic pH dependence represents an advantage compared to hydrogen/deuterium exchange MS.

Quantitative measurement of the solvent accessibility of histidine imidazole groups in proteins

We report a method for expressing the solvent accessibility of histidine imidazole groups in proteins. The method is based on measuring the rate of the hydrogen exchange (HX) reaction of the imidazole C(ε1)-hydrogen. The rate profile of the HX reaction as a function of pH gives a sigmoidal curve, which reaches the maximal rate constant (k(max)) on the alkaline side of the sigmoidal curve. To quantitatively describe the solvent accessibility of imidazole groups in proteins, it is necessary to compare the k(max) of the imidazole groups with their intrinsic k(max) ((i)k(max)), the maximal rate constants for the given imidazole groups when they are fully exposed to the bulk solvent. However, the mechanism of the HX reaction suggests that the (i)k(max) of an imidazole group differs depending on its pK(a), and no systematic study has been conducted to clarify how the (i)k(max) is affected by pK(a). We therefore investigated the relationship between (i)k(max) and pK(a) using four imidazole derivatives at three different temperatures. The experimentally determined pK(a)-specific (i)k(max) values allowed us to derive a general formula to estimate the (i)k(max) value of any given imidazole group exhibiting a specific pK(a) at a specific temperature. Using the formula, the protection factors (PF), the ratio of (i)k(max) to k(max), of five imidazole groups in dihydrofolate reductase were obtained and used to express the magnitude of their solvent accessibility. In this definition, the smaller the PF value, the higher the solvent accessibility, and a value of 1 indicates full exposure to the bulk solvent. The solvent accessibility expressed by the PF values agreed well with the solvent accessible surface areas obtained from the X-ray diffraction data.

Spatially resolved protein hydrogen exchange measured by subzero-cooled chip-based nanoelectrospray ionization tandem mass spectrometry

Mass spectrometry has become a valuable method for studying structural dynamics of proteins in solution by measuring their backbone amide hydrogen/deuterium exchange (HDX) kinetics. In a typical exchange experiment one or more proteins are incubated in deuterated buffer at physiological conditions. After a given period of deuteration, the exchange reaction is quenched by acidification (pH 2.5) and cooling (0 °C) and the deuterated protein (or a digest thereof) is analyzed by mass spectrometry. The unavoidable loss of deuterium (back-exchange) that occurs under quench conditions is undesired as it leads to loss of information. Here we describe the successful application of a chip-based nanoelectrospray ionization mass spectrometry top-down fragmentation approach based on cooling to subzero temperature (-15 °C) which reduces the back-exchange at quench conditions to very low levels. For example, only 4% and 6% deuterium loss for fully deuterated ubiquitin and β(2)-microglobulin were observed after 10 min of back-exchange. The practical value of our subzero-cooled setup for top-down fragmentation HDX analyses is demonstrated by electron-transfer dissociation of ubiquitin ions under carefully optimized mass spectrometric conditions where gas-phase hydrogen scrambling is negligible. Our results show that the known dynamic behavior of ubiquitin in solution is accurately reflected in the deuterium contents of the fragment ions.

Conformational dynamics of a membrane transport protein probed by H/D exchange and covalent labeling: the glycerol facilitator

Glycerol facilitator (GF) is a tetrameric membrane protein responsible for the selective permeation of glycerol and water. Each of the four GF subunits forms a transmembrane channel. Every subunit consists of six helices that completely span the lipid bilayer, as well as two half-helices (TM7 and TM3). X-ray crystallography has revealed that the selectivity of GF is due to its unique amphipathic channel interior. To explore the structural dynamics of GF, we employ hydrogen/deuterium exchange (HDX) and oxidative labeling with mass spectrometry (MS). HDX-MS reveals that transmembrane helices are generally more protected than extramembrane segments, consistent with data previously obtained for other membrane proteins. Interestingly, TM7 does not follow this trend. Instead, this half-helix undergoes rapid deuteration, indicative of a highly dynamic local structure. The oxidative labeling behavior of most GF residues is consistent with the static crystal structure. However, the side chains of C134 and M237 undergo labeling although they should be inaccessible according to the X-ray structure. In agreement with our HDX-MS data, this observation attests to the fact that TM7 is only marginally stable. We propose that the highly mobile nature of TM7 aids in the efficient diffusion of guest molecules through the channel ("molecular lubrication"). In the absence of such dynamics, host-guest molecular recognition would favor semipermanent binding of molecules inside the channel, thereby impeding transport. The current work highlights the complementary nature of HDX, covalent labeling, and X-ray crystallography for the characterization of membrane proteins.

Prediction of native-state hydrogen exchange from perfectly funneled energy landscapes

Simulations based on perfectly funneled energy landscapes often capture many of the kinetic features of protein folding. We examined whether simulations based on funneled energy functions can also describe fluctuations in native-state protein ensembles. We quantitatively compared the site-specific local stability determined from structure-based folding simulations, with hydrogen exchange protection factors measured experimentally for ubiquitin, chymotrypsin inhibitor 2, and staphylococcal nuclease. Different structural definitions for the open and closed states based on the number of native contacts for each residue, as well as the hydrogen-bonding state, or a combination of both criteria were evaluated. The predicted exchange patterns agree with the experiments under native conditions, indicating that protein topology indeed has a dominant effect on the exchange kinetics. Insights into the simplest mechanistic interpretation of the amide exchange process were thus obtained.

Real-time hydrogen/deuterium exchange kinetics via supercharged electrospray ionization tandem mass spectrometry

Amide hydrogen/deuterium exchange (HDX) rate constants of bovine ubiquitin in an ammonium acetate solution containing 1% of the electrospray ionization (ESI) "supercharging" reagent m-nitrobenzyl alcohol (m-NBA) were obtained using top-down, electron transfer dissociation (ETD) tandem mass spectrometry (MS). The supercharging reagent replaces the acid and temperature "quench" step in the conventional MS approach to HDX experiments by causing rapid protein denaturation to occur in the ESI droplet. The higher charge state ions that are produced with m-NBA are more unfolded, as measured by ion mobility, and result in higher fragmentation efficiency and higher sequence coverage with ETD. Single amino acid resolution was obtained for 44 of 72 exchangeable amide sites, and summed kinetic data were obtained for regions of the protein where adjacent fragment ions were not observed, resulting in an overall spatial resolution of 1.3 residues. Comparison of these results with previous values from NMR indicates that the supercharging reagent does not cause significant structural changes to the protein in the initial ESI solution and that scrambling or back-exchange is minimal. This new method for top-down HDX-MS enables real-time kinetic data measurements under physiological conditions, similar to those obtained using NMR, with comparable spatial resolution and significantly better sensitivity.

Protein vivisection reveals elusive intermediates in folding

Although most folding intermediates escape detection, their characterization is crucial to the elucidation of folding mechanisms. Here, we outline a powerful strategy to populate partially unfolded intermediates: A buried aliphatic residue is substituted with a charged residue (e.g., Leu-->Glu(-)) to destabilize and unfold a specific region of the protein. We applied this strategy to ubiquitin, reversibly trapping a folding intermediate in which the beta5-strand is unfolded. The intermediate refolds to a native-like structure upon charge neutralization under mildly acidic conditions. Characterization of the trapped intermediate using NMR and hydrogen exchange methods identifies a second folding intermediate and reveals the order and free energies of the two major folding events on the native side of the rate-limiting step. This general strategy may be combined with other methods and have broad applications in the study of protein folding and other reactions that require trapping of high-energy states.

Electrospray ionization mass spectrometry as a method for studying the high-pressure denaturation of proteins

High-pressure denaturation of proteins can provide important information concerning their folding and function. These studies require expensive and complicated equipment. In this paper, we present a new convenient method for studying high-pressure denaturation of proteins combining DHX (deuterium–hydrogen exchange) and electrospray ionization MS. Application of various values of pressure causes different degrees of protein unfolding resulting in molecules with a different number of protons available for exchange with deuterons. After decompression a protein refolds and a certain number of deuterons are trapped within the hydrophobic core of a refolded protein. Redissolving the deuterated protein in an aqueous buffer initiates the DHX of amides located on the protein surface only, which can be monitored under atmospheric pressure by MS. Depending on the degree of deuteration after high-pressure treatment, the DHX kinetics are different and indicate how many deuterons were trapped in the protein after refolding. The dependence of this number on pressure gives information on the denaturation state of a protein. The distribution of deuterium along the sequence of a high-pressure-denatured protein was studied the ECD (electron-capture-induced dissociation) on a Fourier-transform mass spectrometer, enabling the monitoring of high-pressure denaturation with single amino acid resolution.

Pressure effects on the ensemble dynamics of ubiquitin inspected with molecular dynamics simulations and isotropic reorientational eigenmode dynamics

According to NMR chemical shift data, the ensemble of ubiquitin is a mixture of "open" and "closed" conformations at rapid equilibrium. Pressure perturbations provide the means to study the transition between the two conformers by imposing an additional constraint on the system's partial molar volume. Here we use nanosecond-timescale molecular dynamics simulations to characterize the network of correlated motions accessible to the conformers at low- and high-pressure conditions. Using the isotropic reorientational eigenmode dynamics formalism to analyze our simulation trajectories, we reproduce NMR relaxation data without fitting any parameters of our model. Comparative analysis of our results suggests that the two conformations behave very differently. The dynamics of the "closed" conformation are almost unaffected by pressure and are dominated by large-amplitude correlated motions of residues 23-34 in the extended alpha-helix. The "open" conformation under conditions of normal pressure displays increased mobility, focused on the loop residues 17-20, 46-55, and 58-59 at the bottom of the core of the structure, as well as the C-terminal residues 69-76, that directly participate in key protein-protein interactions. For the same conformation, a pressure increase induces a loss of separability between molecular tumbling and internal dynamics, while motions between different backbone sites become uncorrelated.

Implications from a network-based topological analysis of ubiquitin unfolding simulations

Background

The architectural organization of protein structures has been the focus of intense research since it can hopefully lead to an understanding of how proteins fold. In earlier works we had attempted to identify the inherent structural organization in proteins through a study of protein topology. We obtained a modular partitioning of protein structures with the modules correlating well with experimental evidence of early folding units or “foldons”. Residues that connect different modules were shown to be those that were protected during the transition phase of folding.

Methodology/Principal Findings

In this work, we follow the topological path of ubiquitin through molecular dynamics unfolding simulations. We observed that the use of recurrence quantification analysis (RQA) could lead to the identification of the transition state during unfolding. Additionally, our earlier contention that the modules uncovered through our graph partitioning approach correlated well with early folding units was vindicated through our simulations. Moreover, residues identified from native structure as connector hubs and which had been shown to be those that were protected during the transition phase of folding were indeed more stable (less flexible) well beyond the transition state. Further analysis of the topological pathway suggests that the all pairs shortest path in a protein is minimized during folding.

Conclusions

We observed that treating a protein native structure as a network by having amino acid residues as nodes and the non-covalent interactions among them as links allows for the rationalization of many aspects of the folding process. The possibility to derive this information directly from 3D structure opens the way to the prediction of important residues in proteins, while the confirmation of the minimization of APSP for folding allows for the establishment of a potentially useful proxy for kinetic optimality in the validation of sequence-structure predictions.

Nonuniform isotope patterns produced by collision-induced dissociation of homogeneously labeled ubiquitin: implications for spatially resolved hydrogen/deuterium exchange ESI-MS studies

There is an ongoing debate whether collision-induced dissociation (CID) of electrosprayed proteins after solution-phase hydrogen/deuterium exchange (HDX) is a viable approach for determining spatially resolved deuteration patterns. This work explores the use of two methods, source-CID and hexapole tandem mass spectrometry (MS/MS) on a quadrupole time-of-flight (Q-TOF) mass spectrometer, for measuring the fragment deuteration levels of regioselectively labeled ubiquitin. Both methods reveal that b-ions exhibit HDX levels significantly below that of the intact protein, whereas several y'' fragments are labeled to a much greater extent. These results are consistent with earlier source-CID data (Akashi, S.; Naito, Y.; Takio, K. Anal. Chem. 1999, 71, 4974-4980). However, the measured b-ion deuteration levels are in disagreement with the known solution-phase behavior of ubiquitin. Partial agreement is observed for y''-ions. Control experiments on homogeneously labeled ubiquitin (having the same average deuteration level at every exchangeable site) result in highly nonuniform fragment HDX levels. In particular, b-ions exhibit deuteration levels significantly below that of intact ubiquitin, thereby mimicking the behavior seen for the regioselectively labeled protein. This effect is likely caused by isotope fractionation during collisional activation, facilitated by the high mobility of charge carriers (scrambling) in the gas phase. The observation that the b-ion labeling behavior is largely independent of the spatial isotope distribution within solution-phase ubiquitin invalidates these ions as reporters of the protein deuteration pattern. This work questions the common practice of interpreting any nonuniformities in fragment deuteration as being indicative of regioselective solution-phase labeling. Artifactual deuterium enrichment or depletion during collisional activation may have contributed to the current lack of consensus as to whether HDX/CID represents a potentially viable tool for measuring solution-phase deuteration patterns.

Characterization of the molten globule conformation of V26A ubiquitin by far-UV circular dichroic spectroscopy and amide hydrogen/deuterium exchange

The molten globular conformation of V26A ubiquitin (valine to alanine mutation at residue 26) was studied by nuclear magnetic resonance spectroscopy in conjunction with amide hydrogen/deuterium exchange. Most of the amide protons that are involved in the native secondary structures were observed to be protected in the molten globule state with the protection factors from 1.2 to 6.7. These protection factors are about 2 to 6 orders of magnitude smaller than those of the native state. These observations indicate that V26A molten globule has native-like backbone structure with marginal stability. The comparison of amide protection factors of V26A ubiquitin molten globule state with those of initial collapsed state of the wild type ubiquitin suggests that V26A ubiquitin molten globule state is located close to unfolded state in the folding reaction coordinate. It is considered that V26A ubiquitin molten globule is useful model to study early events in protein folding reaction.

Molecular dynamics simulations of the native and partially folded states of ubiquitin: influence of methanol cosolvent, pH, and temperature on the protein structure and dynamics

A series of explicit-solvent molecular dynamics simulations of the protein ubiquitin are reported, which investigate the effect of environmental factors (presence of methanol cosolvent in the aqueous solution, neutral or low pH value, room or elevated temperature) on the structure, stability, and dynamics of the protein. The simulations are initiated either from the native structure of the protein or from a model of a partially folded state (A-state) that is known to exist at low pH in methanol-water mixtures. The main results of the simulations are: (1) The ubiquitin native structure is remarkably stable at neutral pH in water; (2) the addition of the methanol cosolvent enhances the stability of the secondary structure but weakens tertiary interactions within the protein; (3) this influence of methanol on the protein structure is enhanced at low pH, while the effect of lowering the pH in pure water is limited; and (4) the A-state of ubiquitin can be described as a set of relatively rigid secondary structure elements (a native-like beta-sheet and native-like alpha-helix plus two nonnative alpha-helices) connected by flexible linkers.

Residue-specific NH exchange rates studied by NMR diffusion experiments

We present a novel approach to the investigation of rapid (>2s(-1)) NH exchange rates in proteins, based on residue-specific diffusion measurements. (1)H, (15)N-DOSY-HSQC spectra are recorded in order to observe resolved amide proton signals for most residues of the protein. Human ubiquitin was used to demonstrate the proposed method. Exchange rates are derived directly from the decay data of the diffusion experiment by applying a model deduced from the assumption of a two-site exchange with water and the "pure" diffusion coefficients of water and protein. The "pure" diffusion coefficient of the protein is determined in an experiment with selective excitation of the amide protons in order to suppress the influence of magnetization transfer from water to amide protons on the decay data. For rapidly exchanging residues a comparison of our results with the exchange rates obtained in a MEXICO experiment showed good agreement. Molecular dynamics (MD) and quantum mechanical calculations were performed to find molecular parameters correlating with the exchangeability of the NH protons. The RMS fluctuations of the amide protons, obtained from the MD simulations, together with the NH coupling constants provide a bilinear model which shows a good correlation with the experimental NH exchange rates.

Protein folding and unfolding studied at atomic resolution by fast two-dimensional NMR spectroscopy

Atom-resolved real-time studies of kinetic processes in proteins have been hampered in the past by the lack of experimental techniques that yield sufficient temporal and atomic resolution. Here we present band-selective optimized flip-angle short transient (SOFAST) real-time 2D NMR spectroscopy, a method that allows simultaneous observation of reaction kinetics for a large number of nuclear sites along the polypeptide chain of a protein with an unprecedented time resolution of a few seconds. SOFAST real-time 2D NMR spectroscopy combines fast NMR data acquisition techniques with rapid sample mixing inside the NMR magnet to initiate the kinetic event. We demonstrate the use of SOFAST real-time 2D NMR to monitor the conformational transition of alpha-lactalbumin from a molten globular to the native state for a large number of amide sites along the polypeptide chain. The kinetic behavior observed for the disappearance of the molten globule and the appearance of the native state is monoexponential and uniform along the polypeptide chain. This observation confirms previous findings that a single transition state ensemble controls folding of alpha-lactalbumin from the molten globule to the native state. In a second application, the spontaneous unfolding of native ubiquitin under nondenaturing conditions is characterized by amide hydrogen exchange rate constants measured at high pH by using SOFAST real-time 2D NMR. Our data reveal that ubiquitin unfolds in a gradual manner with distinct unfolding regimes.

Real-time monitoring of chemical transformations by ultrafast 2D NMR spectroscopy

An approach enabling the acquisition of 2D nuclear magnetic resonance (NMR) spectra within a single scan has been recently proposed. A promising application opened up by this "ultrafast" data acquisition format concerns the monitoring of chemical transformations as they happen, in real time. The present paper illustrates some of this potential with two examples: (i) following an H/D exchange process that occurs upon dissolving a protonated protein in D2O, and (ii) real-time in situ tracking of a transient Meisenheimer complex that forms upon rapidly mixing two organic reactants inside the NMR observation tube. The first of these measurements involved acquiring a train of 2D 1H-15N HSQC NMR spectra separated by ca. 4 s; following an initial dead time, this allowed us to monitor the kinetics of hydrogen exchange in ubiquitin at a site-resolved level. The second approach enabled us to observe, within ca. 2 s after the triggering of the reaction, a competition between thermodynamic and kinetic controls via changes in a series of 2D TOCSY patterns. The real-time dynamic experiments hereby introduced thus add to an increasing family of fast characterization techniques based on 2D NMR; their potential and limitations are briefly discussed.

Structural and dynamic characteristics of a partially folded state of ubiquitin revealed by hydrogen exchange mass spectrometry

Structural and dynamic properties of a partially folded conformation (A-state) of ubiquitin are studied using amide hydrogen exchange in solution (HDX) and mass spectrometric detection. A clear distinction between the native state of the protein and the A-state can be made when HDX is carried out in a semicorrelated regime. Convoluted exchange patterns are interpreted with the aid of HDX simulations in a three-state system (highly structured, partially unstructured, and fully unstructured states). The data clearly indicate a highly dynamic character of the non-native state. Furthermore, combination of HDX and protein ion fragmentation in the gas phase [by means of collision-induced dissociation (CAD)] is used to evaluate the conformational stability of various protein segments specifically in the molten globular state. Chain flexibility appears to be distributed very unevenly in this non-native conformation. The highest degree of structural disorder is displayed by the C-terminal segment (Gly(53)-Gly(76)), which was previously suggested to form a transient alpha-helix. The least dynamic segment of ubiquitin in the A-state is Thr(9)-Glu(18) (which was previously suggested to form a stable nativelike beta-strand), with the adjacent segments exhibiting somewhat diminished conformational stability. The study also demonstrates the power of mass spectrometry as a tool in providing conformer-specific information about the structure and dynamics of both native and non-native protein states coexisting in solution under equilibrium.

NMR snapshots of a fluctuating protein structure: ubiquitin at 30 bar-3 kbar

Conformational fluctuation plays a key role in protein function, but we know little about the associated structural changes. Here we present a general method for elucidating, at the atomic level, a large-scale shape change of a protein molecule in solution undergoing conformational fluctuation. The method utilizes the intimate relationship between conformation and partial molar volume and determines three-dimensional structures of a protein at different pressures using variable pressure NMR technique, whereby NOE distance and torsion angle constraints are used to create average coordinates. Ubiquitin (pH 4.6 at 20 degrees C) was chosen as the first target, for which structures were determined at 30 bar and at 3 kbar, giving "NMR snapshots" of a fluctuating protein structure at atomic resolution. The result reveals that the helix swings in and out by >3 angstroms with a simultaneous reorientation of the C-terminal segment, providing an "open" conformer suitable for enzyme recognition. Spin relaxation analysis indicates that this fluctuation occurs in the ten microsecond time range with activation volumes -4.2(+/-3.2) and 18.5(+/-3.0) ml/mol for the "closed-to-open" and the "open-to-closed" transitions, respectively.

Discerning the Structure and Energy of Multiple Transition States in Protein Folding using ψ-Analysis

We quantify the degree to which folding occurs along a complex landscape with structurally distinct pathways using psi-analysis in combination with a protein engineering method that identifies native, non-covalent polypeptide interactions and their relative populations at the rate-limiting step. By probing the proximity of two specific partners, this method is extremely well-suited for comparison to theoretical simulations. Using ubiquitin as a model system, we detect individual pathways with site-resolved resolution, demonstrating that the protein folds through a native-like transition state ensemble with a common nucleus that contains heterogeneous features on its periphery. The consensus transition state topology has part of the major helix docked against four properly aligned beta-strands. However, structural heterogeneity exists in the transition state ensemble, wherein peripheral regions are differentially populated according to their relative stability. Pathway diversity reflects the variable order of formation of these peripheral elements, which radiate outward from the common nucleus. These results, which show only moderate agreement with traditional mutational phi-analysis, provide an extraordinarily detailed and quantitative description of protein folding.

Identifying protein folding cores from the evolution of flexible regions during unfolding

The unfolding of a protein can be described as a transition from a predominantly rigid, folded structure to an ensemble of denatured states. During unfolding, the hydrogen bonds and salt bridges break, destabilizing the secondary and tertiary structure. Our previous work shows that the network of covalent bonds, salt bridges, hydrogen bonds, and hydrophobic interactions forms constraints that define which regions of the native protein are flexible or rigid (structurally stable). Here, we test the hypothesis that information about the folding pathway is encoded in the energetic hierarchy of non-covalent interactions in the native-state structure. The incremental thermal denaturation of protein structures is simulated by diluting the network of salt bridges and hydrogen bonds, breaking them one by one, from weakest to strongest. The structurally stable and flexible regions are identified at each step, providing information about the evolution of flexible regions during denaturation. The folding core, or center of structure formation during folding, is predicted as the region formed by two or more secondary structures having the greatest stability against denaturation. For 10 proteins with different architectures, we show that the predicted folding cores from this flexibility/stability analysis are in good agreement with those identified by native-state hydrogen-deuterium exchange experiments.

NMR of hydrogen bonding in cold-shock protein A and an analysis of the influence of crystallographic resolution on comparisons of hydrogen bond lengths

Hydrogen bonding in cold-shock protein A of Escherichia coli has been investigated using long-range HNCO spectroscopy. Nearly half of the amide protons involved in hydrogen bonds in solution show no measurable protection from exchange in water, cautioning against a direct correspondence between hydrogen bonding and hydrogen exchange protection. The N to O atom distance across a hydrogen bond, R(NO), is related to the size of the (3h)J(NC') trans hydrogen bond coupling constant and the amide proton chemical shift. Both NMR parameters show poorer agreement with the 2.0-A resolution X-ray structure of the cold-shock protein studied by NMR than with a 1.2-A resolution X-ray structure of a homologous cold-shock protein from the thermophile B. caldolyticus. The influence of crystallographic resolution on comparisons of hydrogen bond lengths was further investigated using a database of 33 X-ray structures of ribonuclease A. For highly similar structures, both hydrogen bond R(NO) distance and Calpha coordinate root mean square deviations (RMSD) show systematic increases as the resolution of the X-ray structure used for comparison decreases. As structures diverge, the effects of coordinate errors on R(NO) distance and Calpha coordinate root mean square deviations become progressively smaller. The results of this study are discussed with regard to the influence of data precision on establishing structure similarity relationships between proteins.

The methanol-induced conformational transitions of beta-lactoglobulin, cytochrome c, and ubiquitin at low pH: a study by electrospray ionization mass spectrometry

The methanol-induced conformational transitions under acidic conditions for beta-lactoglobulin, cytochrome c, and ubiquitin, representing three different classes of proteins with beta-sheets, alpha-helices, and both alpha-helices and beta-sheets, respectively, are studied under equilibrium conditions by electrospray ionization mass spectrometry (ESI-MS). The folding states of proteins in solution are monitored by the charge state distributions that they produce during ESI and by hydrogen/deuterium (H/D) exchange followed by ESI-MS. The changes in charge state distributions are correlated with earlier studies by optical and other methods which have shown that, in methanol, these proteins form partially unfolded intermediates with induced alpha-helix structure. Intermediate states formed at about 35% methanol concentration are found to give bimodal charge state distributions. The same rate of H/D exchange is shown by the two contributions to the bimodal distributions. This suggests the intermediates are highly flexible and may consist of a mixture of two or more rapidly interconverting conformers. H/D exchange of proteins followed by ESI-MS shows that helical denatured states, populated at around 50% methanol concentration, transform into more protected structures with further increases in methanol concentration, consistent with previous circular dicroism studies. These more protected structures still produce high charge states in ESI, similar to those of the fully denatured proteins.

Characterization of NMR relaxation-active motions of a partially folded A-state analogue of ubiquitin

The dominant dynamics of a partially folded A-state analogue of ubiquitin that give rise to NMR 15N spin relaxation have been investigated using molecular dynamics (MD) computer simulations and reorientational quasiharmonic analysis. Starting from the X-ray structure of native ubiquitin with a protonation state corresponding to a low pH, the A-state analogue was generated by a MD simulation of a total length of 33 ns in a 60%/40% methanol/water mixture using a variable temperature scheme to control and speed up the structural transformation. The N-terminal half of the A-state analogue consists of loosely coupled native-like secondary structural elements, while the C-terminal half is mostly irregular in structure. Analysis of dipolar N-H backbone correlation functions reveals reorientational amplitudes and time-scale distributions that are comparable to those observed experimentally. Thus, the trajectory provides a realistic picture of a partially folded protein that can be used for gaining a better understanding of the various types of reorientational motions that are manifested in spin-relaxation parameters of partially folded systems. For this purpose, a reorientational quasiharmonic reorientational analysis was performed on the final 5 ns of the trajectory of the A-state analogue, and for comparison on a 5 ns trajectory of native ubiquitin. The largest amplitude reorientational modes show a markedly distinct behavior for the two states. While for native ubiquitin, such motions have a more local character involving loops and the C-terminal end of the polypeptide chain, the A-state analogue shows highly collective motions in the nanosecond time-scale range corresponding to larger-scale movements between different segments. Changes in reorientational backbone entropy between the A-state analogue and the native state of ubiquitin, which were computed from the reorientational quasiharmonic analyses, are found to depend significantly on motional correlation effects.