Heteronuclear three-dimensional NMR spectroscopy of a partially denatured protein: the A-state of human ubiquitin

Structural Elucidation of Ubiquitin via Gas-Phase Ion/Ion Cross-Linking Reactions Using Sodium-Cationized Reagents Coupled with Infrared Multiphoton Dissociation

Accurate structural determination of proteins is critical to understanding their biological functions and the impact of structural disruption on disease progression. Gas-phase cross-linking mass spectrometry (XL-MS) via ion/ion reactions between multiply charged protein cations and singly charged cross-linker anions has previously been developed to obtain low-resolution structural information on proteins. This method significantly shortens experimental time relative to conventional solution-phase XL-MS but has several technical limitations: (1) the singly deprotonated N-hydroxysulfosuccinimide (sulfo-NHS)-based cross-linker anions are restricted to attachment at neutral amine groups of basic amino acid residues and (2) analyzing terminal cross-linked fragment ions is insufficient to unambiguously localize sites of linker attachment. Herein, we demonstrate enhanced structural information for alcohol-denatured A-state ubiquitin obtained from an alternative gas-phase XL-MS approach. Briefly, singly sodiated ethylene glycol bis(sulfosuccinimidyl succinate) (sulfo-EGS) cross-linker anions enable covalent cross-linking at both ammonium and amine groups. Additionally, covalently modified internal fragment ions, along with terminal b-/y-type counterparts, improve the determination of linker attachment sites. Molecular dynamics simulations validate experimentally obtained gas-phase conformations of denatured ubiquitin. This method has identified four cross-linking sites across 8+ ubiquitin, including two new sites in the N-terminal region of the protein that were originally inaccessible in prior gas-phase XL approaches. The two N-terminal cross-linking sites suggest that the N-terminal half of ubiquitin is more compact in gas-phase conformations. By comparison, the two C-terminal linker sites indicate the signature transformation of this region of the protein from a native to a denatured conformation. Overall, the results suggest that the solution-phase secondary structures of the A-state ubiquitin are conserved in the gas phase. This method also provides sufficient sensitivity to differentiate between two gas-phase conformers of the same charge state with subtle structural variations.

In-droplet hydrogen-deuterium exchange to examine protein/peptide solution conformer heterogeneity

RATIONALE

Many different structure analysis techniques are not capable of probing the heterogeneity of solution conformations. Here, we examine the ability of in-droplet hydrogen-deuterium exchange (HDX) to directly probe solution conformer heterogeneity of a protein with mass spectrometry (MS) detection.

METHODS

Two vibrating capillary vibrating sharp-edge spray ionization (cVSSI) devices have been arranged such that they generate microdroplet plumes of the analyte and D₂ O reagent, which coalesce to form reaction droplets where HDX takes place in the solution environment. The native HDX-MS setup has been first explored for two model peptides that have distinct structural compositions in solution. The effectiveness of the multidevice cVSSI-HDX in illustrating structural details has been further exploited to investigate coexisting solution-phase conformations of the protein ubiquitin.

RESULTS

In-droplet HDX reveals decreased backbone exchange for a model peptide that has a greater helix-forming propensity. Differences in intrinsic rates of the alanine and serine residues may account for much of the observed protection. The data allow the first estimates of backbone exchange rates for peptides undergoing in-droplet HDX. That said, the approach may hold greater potential for investigating the tertiary structure and structural transitions of proteins. For ubiquitin protein, HDX reactivity differences suggest that multiple conformers are present in native solutions. The addition of methanol to buffered aqueous solutions of ubiquitin results in increased populations of solution conformers of higher reactivity. Data analysis suggests that partially folded conformers such as the A-state of ubiquitin increase with methanol content; the native state may be preserved to a limited degree even under stronger denaturation conditions.

CONCLUSION

The deuterium uptake after in-droplet HDX has been observed to correspond to some degree with peptide backbone hydrogen protection based on differences in intrinsic rates of exchange. The presence of coexisting protein solution structures under native and denaturing solution conditions has been distinguished by the isotopic distributions of deuterated ubiquitin ions.

Protein Unfolding in Freeze Frames: Intermediate States are Revealed by Variable-Temperature Ion Mobility-Mass Spectrometry

The gas phase is an idealized laboratory for the study of protein structure, from which it is possible to examine stable and transient forms of mass-selected ions in the absence of bulk solvent. With ion mobility-mass spectrometry (IM-MS) apparatus built to operate at both cryogenic and elevated temperatures, we have examined conformational transitions that occur to the monomeric proteins: ubiquitin, lysozyme, and α-synuclein as a function of temperature and in source activation. We rationalize the experimental observations with a temperature-dependent framework model and comparison to known conformers. Data from ubiquitin show unfolding transitions that proceed through diverse and highly elongated intermediate states, which converge to more compact structures. These findings contrast with data obtained from lysozyme─a protein where (un)-folding plasticity is restricted by four disulfide linkages, although this is alleviated in its reduced form. For structured proteins, collision activation of the protein ions in-source enables subsequent "freezing" or thermal annealing of unfolding intermediates, whereas disordered proteins restructure substantially at 250 K even without activation, indicating that cold denaturation can occur without solvent. These data are presented in the context of a toy model framework that describes the relative occupancy of the available conformational space.

Microfluidic Platform for Time-Resolved Characterization of Protein Higher-Order Structures and Dynamics Using Top-Down Mass Spectrometry

Characterization of protein higher-order structures and dynamics is essential for understanding the biological functions of proteins and revealing the underlying mechanisms. Top-down mass spectrometry (MS) accesses structural information at both the intact protein level and the peptide fragment level. Native top-down MS allows analysis of a protein complex's architecture and subunits' identity and modifications. Top-down hydrogen/deuterium exchange (HDX) MS offers high spatial resolution for conformational or binding interface analysis and enables conformer-specific characterization. A microfluidic chip can provide superior performance for front-end reactions useful for these MS workflows, such as flexibility in manipulating multiple reactant flows, integrating various functional modules, and automation. However, most microchip-MS devices are designed for bottom-up approaches or top-down proteomics. Here, we demonstrate a strategy for designing a microchip for top-down MS analysis of protein higher-order structures and dynamics. It is suitable for time-resolved native MS and HDX MS, with designs aiming for efficient ionization of intact protein complexes, flexible manipulation of multiple reactant flows, and precise control of reaction times over a broad range of flow rates on the submicroliter per minute scale. The performance of the prototype device is demonstrated by measurements of systems including monoclonal antibodies, antibody-antigen complexes, and coexisting protein conformers. This strategy may benefit elaborate structural analysis of biomacromolecules and inspire method development using the microchip-MS approach.

Gas-Phase Ion/Ion Chemistry for Structurally Sensitive Probes of Gaseous Protein Ion Structure: Electrostatic and Electrostatic to Covalent Cross-Linking

Intramolecular interactions within a protein are key in maintaining protein tertiary structure and understanding how proteins function. Ion mobility-mass spectrometry (IM-MS) has become a widely used approach in structural biology since it provides rapid measurements of collision cross sections (CCS), which inform on the gas-phase conformation of the biomolecule under study. Gas-phase ion/ion reactions target amino acid residues with specific chemical properties and the modified sites can be identified by MS. In this study, electrostatically reactive, gas-phase ion/ion chemistry and IM-MS are combined to characterize the structural changes between ubiquitin electrosprayed from aqueous and denaturing conditions. The electrostatic attachment of sulfo-NHS acetate to ubiquitin via ion/ion reactions and fragmentation by electron-capture dissociation (ECD) provide the identification of the most accessible protonated sites within ubiquitin as the sulfonate group forms an electrostatic complex with accessible protonated side chains. The protonated sites identified by ECD from the different solution conditions are distinct and, in some cases, reflect the disruption of interactions such as salt bridges that maintain the native protein structure. This agrees with previously published literature demonstrating that a high methanol concentration at low pH causes the structure of ubiquitin to change from a native (N) state to a more elongated A state. Results using gas-phase, electrostatic cross-linking reagents also point to similar structural changes and further confirm the role of methanol and acid in favoring a more unfolded conformation. Since cross-linking reagents have a distance constraint for the two reactive sites, the data is valuable in guiding computational structures generated by molecular dynamics. The research presented here describes a promising strategy that can detect subtle changes in the local environment of targeted amino acid residues to inform on changes in the overall protein structure.

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.

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.

Comprehensive structural and dynamical view of an unfolded protein from the combination of single-molecule FRET, NMR, and SAXS

The properties of unfolded proteins are essential both for the mechanisms of protein folding and for the function of the large group of intrinsically disordered proteins. However, the detailed structural and dynamical characterization of these highly dynamic and conformationally heterogeneous ensembles has remained challenging. Here we combine and compare three of the leading techniques for the investigation of unfolded proteins, NMR spectroscopy (NMR), small-angle X-ray scattering (SAXS), and single-molecule Förster resonance energy transfer (FRET), with the goal of quantitatively testing their consistency and complementarity and for obtaining a comprehensive view of the unfolded-state ensemble. Using unfolded ubiquitin as a test case, we find that its average dimensions derived from FRET and from structural ensembles calculated using the program X-PLOR-NIH based on NMR and SAXS restraints agree remarkably well; even the shapes of the underlying intramolecular distance distributions are in good agreement, attesting to the reliability of the approaches. The NMR-based results provide a highly sensitive way of quantifying residual structure in the unfolded state. FRET-based nanosecond fluorescence correlation spectroscopy allows long-range distances and chain dynamics to be probed in a time range inaccessible by NMR. The combined techniques thus provide a way of optimally using the complementarity of the available methods for a quantitative structural and dynamical description of unfolded proteins both at the global and the local level.

Observing a late folding intermediate of Ubiquitin at atomic resolution by NMR

The study of intermediates in the protein folding pathway provides a wealth of information about the energy landscape. The intermediates also frequently initiate pathogenic fibril formations. While observing the intermediates is difficult due to their transient nature, extreme conditions can partially unfold the proteins and provide a glimpse of the intermediate states. Here, we observe the high resolution structure of a hydrophobic core mutant of Ubiquitin at an extreme acidic pH by nuclear magnetic resonance (NMR) spectroscopy. In the structure, the native secondary and tertiary structure is conserved for a major part of the protein. However, a long loop between the beta strands β3 and β5 is partially unfolded. The altered structure is supported by fluorescence data and the difference in free energies between the native state and the intermediate is reflected in the denaturant induced melting curves. The unfolded region includes amino acids that are critical for interaction with cofactors as well as for assembly of poly-Ubiquitin chains. The structure at acidic pH resembles a late folding intermediate of Ubiquitin and indicates that upon stabilization of the protein's core, the long loop converges on the core in the final step of the folding process.

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.

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.

Microsecond molecular dynamics simulation of guanidinium chloride induced unfolding of ubiquitin

All atom molecular dynamics simulations have been used to explore the atomic detail mechanism of guanidinium induced unfolding of the protein ubiquitin.

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.

Gas-phase intramolecular protein crosslinking via ion/ion reactions: ubiquitin and a homobifunctional sulfo-NHS ester

Gas-phase intra-molecular crosslinking of protein ubiquitin cations has been demonstrated via ion/ion reactions with anions of a homobifunctional N-hydroxysulfosuccinimide (sulfo-NHS) ester reagent. The ion/ion reaction between multiply-protonated ubiquitin and crosslinker monoanions produces a stable, charge-reduced complex. Covalent crosslinking is indicated by the consecutive loss of 2 molecules of sulfo-NHS under ion trap collisional activation conditions. Covalent modification is verified by the presence of covalently crosslinked sequence ions produced by ion-trap collision-induced dissociation of the ion generated from the losses of sulfo-NHS. Analysis of the crosslinked sequence fragments allows for the localization of crosslinked primary amines, enabling proximity mapping of the gas-phase 3-D structures. The presence of two unprotonated reactive sites within the distance constraint of the crosslinker is required for successful crosslinking. The ability to covalently crosslink is, therefore, sensitive to protein charge state. As the charge state increases, fewer reactive sites are available and protein structure is more likely to become extended because of intramolecular electrostatic repulsion. At high charge states, the reagent shows little evidence for covalent crosslinking but does show evidence for 'electrostatic crosslinking' in that the binding of the sulfonate groups to the protein is sufficiently strong that backbone cleavages are favored over reagent detachment under ion trap collisional activation conditions.

Effects of Fe(II)/H2O2 oxidation on ubiquitin conformers measured by ion mobility-mass spectrometry

Oxidative modifications can have significant effects on protein structure in solution. Here, the structures and stabilities of oxidized ubiquitin ions electrosprayed from an aqueous solution (pH 2) are studied by ion mobility spectrometry-mass spectrometry (IMS-MS). IMS-MS has proven to be a valuable technique to assess gas phase and in many cases, solution structures. Herein, in vitro oxidation is performed by Fenton chemistry with Fe(II)/hydrogen peroxide. Most molecules in solution remain unmodified, whereas ∼20% of the population belongs to an M+16 Da oxidized species. Ions of low charge states (+7 and +8) show substantial variance in collision cross section distributions between unmodified and oxidized species. Novel and previously reported gaussian conformers are used to model cross section distributions for +7 and +8 oxidized ubiquitin ions, respectively, in order to correlate variances in observed gas-phase distributions to changes in populations of solution states. Based on gaussian modeling, oxidized ions of charge state +7 have an A-state conformation which is more populated for oxidized relative to unmodified ions. Oxidized ubiquitin ions of charge state +8 have a distribution of conformers arising from native-state ubiquitin and higher intensities of A- and U-state conformers relative to unmodified ions. This work provides evidence that incorporation of a single oxygen atom to ubiquitin leads to destabilization of the native state in an acidic solution (pH ∼2) and to unfolding of gas-phase compact structures.

A thermodynamic definition of protein domains

Protein domains are conspicuous structural units in globular proteins, and their identification has been a topic of intense biochemical interest dating back to the earliest crystal structures. Numerous disparate domain identification algorithms have been proposed, all involving some combination of visual intuition and/or structure-based decomposition. Instead, we present a rigorous, thermodynamically-based approach that redefines domains as cooperative chain segments. In greater detail, most small proteins fold with high cooperativity, meaning that the equilibrium population is dominated by completely folded and completely unfolded molecules, with a negligible subpopulation of partially folded intermediates. Here, we redefine structural domains in thermodynamic terms as cooperative folding units, based on m-values, which measure the cooperativity of a protein or its substructures. In our analysis, a domain is equated to a contiguous segment of the folded protein whose m-value is largely unaffected when that segment is excised from its parent structure. Defined in this way, a domain is a self-contained cooperative unit; i.e., its cooperativity depends primarily upon intrasegment interactions, not intersegment interactions. Implementing this concept computationally, the domains in a large representative set of proteins were identified; all exhibit consistency with experimental findings. Specifically, our domain divisions correspond to the experimentally determined equilibrium folding intermediates in a set of nine proteins. The approach was also proofed against a representative set of 71 additional proteins, again with confirmatory results. Our reframed interpretation of a protein domain transforms an indeterminate structural phenomenon into a quantifiable molecular property grounded in solution thermodynamics.

Conformation types of ubiquitin [M+8H]8+ Ions from water:methanol solutions: evidence for the N and A States in aqueous solution

Ion mobility and mass spectrometry measurements are used to examine the gas-phase populations of [M+8H](8+) ubiquitin ions formed upon electrospraying 20 different solutions from 100:0 to 5:95 water:methanol that are maintained at pH = 2.0. Over this range of solution conditions, mobility distributions for the +8 charge state show substantial variations. Here we develop a model that treats the combined measurements as one data set. By varying the relative abundances of a discrete set of conformation types, it is possible to represent distributions obtained from any solution. For solutions that favor the well-known A-state ubiquitin, it is possible to represent the gas-phase distributions with seven conformation types. Aqueous conditions that favor the native structure require four more structural types to represent the distribution. This analysis provides the first direct evidence for trace amounts of the A state under native conditions. The method of analysis presented here should help illuminate how solution populations evolve into new gas-phase structures as solvent is removed. Evidence for trace quantities of previously unknown states under native solution conditions may provide insight about the relationship of dynamics to protein function as well as misfolding and aggregation phenomena.

Ubiquitin fragments: their known biological activities and putative roles

Ubiquitin (Ub) is involved in many key processes of cell biology. Identification of compounds that could interfere in the ubiquitination process is of importance. It could be expected that peptides derived from the Ub-binding regions might be able to interact with Ub receptors themselves and modify an ability of the Ub receptors interactions. This review summarizes current knowledge about known Ub-derived peptides and discusses putative activity of unexplored Ub fragments. Among identified biologically active Ub-derived peptides, its decapeptide fragment of the LEDGRTLSDY sequence was found to exhibit strong immunosuppressive effects on the cellular and humoral immune responses, comparable to that of cyclosporine. Some of the Ub fragments possess strong antibacterial and antifungal potency. In the search for new peptides that could interfere in the interaction of Ub with other proteins, we investigated the pentapeptide Ub sequences present in non-ubiquitin proteins. Based on examination of the Swiss-Prot database, we postulated that sequences of some Ub fragments often exist in other protein molecules. However, some of those motives are represented more frequently than others and could be involved in regulation of cellular processes related to Ub.

Probing the folding transition state of ubiquitin mutants by temperature-jump-induced downhill unfolding

Crucial to revealing mechanistic details of protein folding is a characterization of the transition state ensemble and its structural dynamics. To probe the transition state of ubiquitin thermal unfolding, we examine unfolding dynamics and kinetics of wild-type and mutant ubiquitin using time-resolved nonlinear infrared spectroscopy after a nanosecond temperature jump. We observe spectral changes on two different time scales. A fast nonexponential microsecond phase is attributed to downhill unfolding from the transition state region, which is induced by a shift of the barrier due to the rapid temperature change. Slow millisecond changes arise from thermally activated folding and unfolding kinetics. Mutants that stabilize or destabilize beta strands III-V lead to a decreased or increased amplitude of the microsecond phase, indicating that the disruption or weakening of these strands occurs in the transition state. Unfolding features from microseconds to milliseconds can be explained by temperature-dependent changes of a two-dimensional free energy surface constructed by the native contacts between beta strands of the protein. In addition, the results support the possibility of an intermediate state in thermal unfolding.

Temperature-dependent downhill unfolding of ubiquitin. I. Nanosecond-to-millisecond resolved nonlinear infrared spectroscopy

Transient thermal unfolding of ubiquitin is investigated using nonlinear infrared spectroscopy after a nanosecond laser temperature jump (T-jump). The abrupt change in the unfolding free energy surface and the ns time resolution allow us to observe a fast response on ns to micros time-scales, which we attribute to downhill unfolding, before a cross-over to ms kinetics. The downhill unfolding by a sub-population of folded proteins is induced through a shift of the barrier toward the native state. By adjusting the T-jump width, the effect of the initial (T(i)) and final (T(f)) temperature on the unfolding dynamics can be separated. From the amplitude of the fast downhill unfolding, the fractional population prepared at the unfolding transition state is obtained. This population increases with both T(i) and with T(f). A two-state kinetic analysis of the ms refolding provides thermodynamic information about the barrier height. By a combination of the fast and slow unfolding and folding parameters, a quasi-two-state kinetic analysis is performed to calculate the time-dependent population changes of the folded state. This calculation coincides with the experimentally obtained population changes at low temperature but deviations are found in the T-jump from 67 to 78 degrees C. Using temperature-dependent barrier height changes, a temperature Phi value analysis is performed. The result shows a decreasing trend of Phi(T) with temperature, which indicates an increase of the heterogeneity of the transition state. We conclude that ubiquitin unfolds along a well-defined pathway at low temperature which expands with increasing temperature to include multiple routes.

Temperature-dependent downhill unfolding of ubiquitin. II. Modeling the free energy surface

To provide evidence for the interpretation of temperature-dependent unfolding kinetics and the downhill unfolding scenario presented in the accompanying experimental article (Part I), the free energy surface of ubiquitin unfolding is calculated using statistical mechanical models of the Muñoz-Eaton (ME) form. The models allow only two states for each amino acid residue, folded or unfolded, and permutations of these states generate an ensemble of microstates. One-dimensional free energy curves are calculated using the number of folded residues as a reaction coordinate. The proposed sequential unfolding of ubiquitin's beta-sheet is tested by mapping the free energy onto two reaction coordinates inspired by the experiment as follows: the number of folded residues in ubiquitin's stable beta-strands I and II and those of the less stable strands III-V. Although the original ME model successfully captures folding features of zipper-like one-dimensional folders, it misses important tertiary interactions between residues that are far from each other in primary sequence. To take tertiary contacts into account, partially folded microstates based on a spherical growth model are included in the calculation and compared with the original model. By calculating the folding probability of each residue for a given point on the free energy surface, the unfolding pathway of ubiquitin is visualized. At low temperature, thermal unfolding occurs along a sequential unfolding pathway as follows: disruption of the beta-strands III-V followed by unfolding of the strands I and II. At high temperature, multiple unfolding routes are formed. The heterogeneity of the transition state explains the global nonexponential unfolding observed in the T-jump experiment at high temperature. The calculation also reports a high stability for the alpha-helix of ubiquitin.

Transient 2D IR spectroscopy of ubiquitin unfolding dynamics

Transient two-dimensional infrared (2D IR) spectroscopy is used as a probe of protein unfolding dynamics in a direct comparison of fast unfolding experiments with molecular dynamics simulations. In the experiments, the unfolding of ubiquitin is initiated by a laser temperature jump, and protein structural evolution from nanoseconds to milliseconds is probed using amide I 2D IR spectroscopy. The temperature jump prepares a subensemble near the unfolding transition state, leading to quasi-barrierless unfolding (the "burst phase") before the millisecond activated unfolding kinetics. The burst phase unfolding of ubiquitin is characterized by a loss of the coupling between vibrations of the beta-sheet, a process that manifests itself in the 2D IR spectrum as a frequency blue-shift and intensity decrease of the diagonal and cross-peaks of the sheet's two IR active modes. As the sheet unfolds, increased fluctuations and solvent exposure of the beta-sheet amide groups are also characterized by increases in homogeneous linewidth. Experimental spectra are compared with 2D IR spectra calculated from the time-evolving structures in a molecular dynamics simulation of ubiquitin unfolding. Unfolding is described as a sequential unfolding of strands in ubiquitin's beta-sheet, using two collective coordinates of the sheet: (i) the native interstrand contacts between adjacent beta-strands I and II and (ii) the remaining beta-strand contacts within the sheet. The methods used illustrate the general principles by which 2D IR spectroscopy can be used for detailed dynamical comparisons of experiment and simulation.

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.

Characterisation of disulfide-bond dynamics in non-native states of lysozyme and its disulfide deletion mutants by NMR

This report describes NMR-spectroscopic investigations of the conformational dynamics of disulfide bonds in hen-egg-white lysozyme substitution mutants. The following four systems have been investigated: 2SS(alpha), a lysozyme variant that contains C64A, C76A, C80A and C94A substitutions, was studied in water at pH 2 and 3.8 and in urea (8 M, pH 2); 2SS(beta) lysozyme, which has C6S, C30A, C115A and C127A substitutions, was studied in water (pH 2) and urea (8 M, pH 2). The NMR analysis of heteronuclear 15N-relaxation rates shows that the barrier to disulfide-bond isomerisation can vary substantially in different lysozyme mutants and depends on the residual structure present in these states. The investigations reveal cooperativity in the modulation of micro- to millisecond dynamics that is due to the presence of multiple disulfide bridges in lysozyme. Mutation of cysteines in one of the two structural domains substantially diminishes the barrier to rotational isomerisation in the other domain. However, the interactions between hydrophobic clusters within and across the domains remains intact.

Unfolding dynamics of the protein ubiquitin: insight from simulation

The temperature-induced unfolding pathway of ubiquitin has been investigated by molecular dynamics simulation at four different temperatures. It has been observed that the sequences of the unfolding events are same at all the temperatures. However, the time scale of the dynamics at different temperatures are different. The transition states at various temperatures also possess similar secondary structural elements. The intermediate conformations visited by the protein at different temperatures can help determination of the transition states. The well known " state" of ubiquitin, hitherto found to be stable only in methanol water mixture, have been observed to be a common transient intermediate conformation in the unfolding path of the protein in water. Our observation about the similarities of the unfolding process at different temperatures strongly recommend for a defined pathway for the unfolding process.

Influence of Fluoro, Chloro and Alkyl Alcohols on the Folding Pathway of Human Serum Albumin

Urea-induced equilibrium unfolding of human serum albumin (HSA) when studied by mean residue ellipticity at 222 nm (MRE(222)) or intrinsic fluorescence measurements showed a two-step, three-state transition with a stable intermediate around 4.6-5.2 M urea. The presence of 2,2,2-trifluoroethanol (TFE) resulted in a single-step, two-state transition with a significant shift towards higher urea concentration, suggesting the stabilizing effect of TFE. The free energy of stabilization (DeltaDeltaG(D)(H(2)O)) in the presence of 3.0 M TFE was determined to be 2.68 and 2.72 kcal/mol by MRE(222) and fluorescence measurements, respectively. The stabilizing potential of other alcohols on the refolding behavior of HSA at 5.0 M urea (where the intermediate exists) as studied by MRE(222) and intrinsic fluorescence measurements showed the following order: 1,1,1,3,3,3-hexafluoroisopropanol (HFIP) > TFE > 2-chloroethanol > tert-butanol > iso-propanol > ethanol > methanol. Further, the extent of refolding at the highest concentration of alcohol was similar in all cases. The stabilizing effect of TFE on guanidine hydrochloride (GdnHCl)-induced unfolding of HSA was nearly equal to that found for urea denaturation, as reflected in the DeltaDeltaG(D)(H(2)O) value (2.38 kcal/mol). Taken together, these results suggest that the stabilizing effect of TFE and other alcohols on urea/GdnHCl-induced unfolding of HSA is higher for alcohols that contain bulky groups or fluorine atoms.

Conformational changes during the nanosecond-to-millisecond unfolding of ubiquitin

Steady-state and transient conformational changes upon the thermal unfolding of ubiquitin were investigated with nonlinear IR spectroscopy of the amide I vibrations. Equilibrium temperature-dependent 2D IR spectroscopy reveals the unfolding of the beta-sheet of ubiquitin through the loss of cross peaks formed between transitions arising from delocalized vibrations of the beta-sheet. Transient unfolding after a nanosecond temperature jump is monitored with dispersed vibrational echo spectroscopy, a projection of the 2D IR spectrum. Whereas the equilibrium study follows a simple two-state unfolding, the transient experiments observe complex relaxation behavior that differs for various spectral components and spans 6 decades in time. The transient behavior can be separated into fast and slow time scales. From 100 ns to 0.5 ms, the spectral features associated with beta-sheet unfolding relax in a sequential, nonexponential manner, with time constants of 3 micros and 80 micros. By modeling the amide I vibrations of ubiquitin, this observation is explained as unfolding of the less stable strands III-V of the beta-sheet before unfolding of the hairpin that forms part of the hydrophobic core. This downhill unfolding is followed by exponential barrier-crossing kinetics on a 3-ms time scale.

Efficient RMSD measures for the comparison of two molecular ensembles. Root-mean-square deviation

Quantitative measures are presented for comparing the conformations of two molecular ensembles. The measures are based on Kabsch's formula for the root-mean-square deviation (RMSD) and the covariance matrix of atomic positions of isotropically distributed ensembles (IDE). By using a Taylor series expansion, it is shown that the RMSD can be expressed solely in terms of the IDE matrices. A fast approximate method is introduced for the pairwise RMSD determination whose computational cost scales linearly with the number of structures. A similarity measure for two structural ensembles that is based on the trace metric of the differences of powers of the IDE matrices is presented. The measures are illustrated for conformational ensembles generated by a molecular dynamics computer simulation of a partially folded A-state analog of ubiquitin.

Dynamic and structural analysis of isotropically distributed molecular ensembles

An efficient new method is presented for the characterization of motional correlations derived from a set of protein structures without requiring the separation of overall and internal motion. In this method, termed isotropically distributed ensemble (IDE) analysis, each structure is represented by an ensemble of isotropically distributed replicas corresponding to the situation found in an isotropic protein solution. This leads to a covariance matrix of the cartesian atomic positions with elements proportional to the ensemble average of scalar products of the position vectors with respect to the center of mass. Diagonalization of the covariance matrix yields eigenmodes and amplitudes that describe concerted motions of atoms, including overall rotational and intramolecular dynamics. It is demonstrated that this covariance matrix naturally distinguishes between "rigid" and "mobile" parts without necessitating a priori selection of a reference structure and an atom set for the orientational alignment process. The method was applied to the analysis of a 5-ns molecular dynamics trajectory of native ubiquitin and a 40-ns trajectory of a partially folded state of ubiquitin. The results were compared with essential dynamics analysis. By taking advantage of the spherical symmetry of the IDE covariance matrix, more than a 10-fold speed up is achieved for the computation of eigenmodes and mode amplitudes. IDE analysis is particularly suitable for studying the correlated dynamics of flexible and large molecules.

Structure and properties of a dimeric N-terminal fragment of human ubiquitin

Previous peptide dissection and kinetic experiments have indicated that in vitro folding of ubiquitin may proceed via transient species in which native-like structure has been acquired in the first 45 residues. A peptide fragment, UQ(1-51), encompassing residues 1 to 51 of ubiquitin was produced in order to test whether this portion has propensity for independent self-assembly. Surprisingly, the construct formed a folded symmetrical dimer that was stabilised by 0.8 M sodium sulphate at 298 K (the S state). The solution structure of the UQ(1-51) dimer was determined by multinuclear NMR spectroscopy. Each subunit of UQ(1-51) consists of an N-terminal beta-hairpin followed by an alpha-helix and a final beta-strand, with orientations similar to intact ubiquitin. The dimer is formed by the third beta-strand of one subunit interleaving between the hairpin and third strand of the other to give a six-stranded beta-sheet, with the two alpha-helices sitting on top. The helix-helix and strand portions of the dimer interface also mimic related features in the structure of ubiquitin. The structural specificity of the UQ(1-51) peptide is tuneable: as the concentration of sodium sulphate is decreased, near-native alternative conformations are populated in slow chemical exchange. Magnetization transfer experiments were performed to characterize the various species present in 0.35 M sodium sulphate, namely the S state and two minor forms. Chemical shift differences suggest that one minor form is very similar to the S state, while the other experiences a significant conformational change in the third strand. A segmental rearrangement of the third strand in one subunit of the S state would render the dimer asymmetric, accounting for most of our results. Similar small-scale transitions in proteins are often invoked to explain solvent exchange at backbone amide proton sites that have an intermediate level of protection.

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.

Identification and ab initio simulations of early folding units in proteins

The location of protein subunits that form early during folding, constituted of consecutive secondary structure elements with some intrinsic stability and favorable tertiary interactions, is predicted using a combination of threading algorithms and local structure prediction methods. Two folding units are selected among the candidates identified in a database of known protein structures: the fragment 15-55 of 434 cro, an all-alpha protein, and the fragment 1-35 of ubiquitin, an alpha/beta protein. These units are further analyzed by means of Monte Carlo simulated annealing using several database-derived potentials describing different types of interactions. Our results suggest that the local interactions along the chain dominate in the first folding steps of both fragments, and that the formation of some of the secondary structures necessarily occurs before structure compaction. These findings led us to define a prediction protocol, which is efficient to improve the accuracy of the predicted structures. It involves a first simulation with a local interaction potential only, whose final conformation is used as a starting structure of a second simulation that uses a combination of local interaction and distance potentials. The root mean square deviations between the coordinates of predicted and native structures are as low as 2-4 A in most trials. The possibility of extending this protocol to the prediction of full proteins is discussed. Proteins 2001;42:164-176.

Structural characterization of a mutant peptide derived from ubiquitin: implications for protein folding

The formation of the N-terminal beta-hairpin of ubiquitin is thought to be an early event in the folding of this small protein. Previously, we have shown that a peptide corresponding to residues 1-17 of ubiquitin folds autonomously and is likely to have a native-like hairpin register. To investigate the causes of the stability of this fold, we have made mutations in the amino acids at the apex of the turn. We find that in a peptide where Thr9 is replaced by Asp, U(1-17)T9D, the native conformation is stabilized with respect to the wild-type sequence, so much so that we are able to characterize the structure of the mutant peptide fully by NMR spectroscopy. The data indicate that U(1-17)T9D peptide does indeed form a hairpin with a native-like register and a type I turn with a G1 beta-bulge, as in the full-length protein. The reason for the greater stability of the U(1-17)T9D mutant remains uncertain, but there are nuclear Overhauser effects between the side chains of Asp9 and Lys 11, which may indicate that a charge-charge interaction between these residues is responsible.

Solution structure and dynamics of a designed hydrophobic core variant of ubiquitin

BACKGROUND

The recent merger of computation and protein design has resulted in a burst of success in the generation of novel proteins with native-like properties. A critical component of this coupling between theory and experiment is a detailed analysis of the structures and stabilities of designed proteins to assess and improve the accuracy of design algorithms.

RESULTS

Here we report the solution structure of a hydrophobic core variant of ubiquitin, referred to as 1D7, which was designed with the core-repacking algorithm ROC. As a measure of conformational specificity, we also present amide exchange protection factors and backbone and sidechain dynamics. The results indicate that 1D7 is similar to wild-type (WT) ubiquitin in backbone structure and degree of conformational specificity. We also observe a good correlation between experimentally determined sidechain structures and those predicted by ROC. However, evaluation of the core sidechain conformations indicates that, in general, 1D7 has more sidechains in less statistically favorable conformations than WT.

CONCLUSIONS

Our results provide an explanation for the lower stability of 1D7 compared to WT, and suggest modifications to design algorithms that may improve the accuracy with which structure and stability are predicted. The results also demonstrate that core packing can affect conformational flexibility in subtle ways that are likely to be important for the design of function and protein-ligand interactions.

The hydrogen exchange core and protein folding

A database of hydrogen-deuterium exchange results has been compiled for proteins for which there are published rates of out-exchange in the native state, protection against exchange during folding, and out-exchange in partially folded forms. The question of whether the slow exchange core is the folding core (Woodward C, 1993, Trends Biochem Sci 18:359-360) is reexamined in a detailed comparison of the specific amide protons (NHs) and the elements of secondary structure on which they are located. For each pulsed exchange or competition experiment, probe NHs are shown explicitly; the large number and broad distribution of probe NHs support the validity of comparing out-exchange with pulsed-exchange/competition experiments. There is a strong tendency for the same elements of secondary structure to carry NHs most protected in the native state, NHs first protected during folding, and NHs most protected in partially folded species. There is not a one-to-one correspondence of individual NHs. Proteins for which there are published data for native state out-exchange and theta values are also reviewed. The elements of secondary structure containing the slowest exchanging NHs in native proteins tend to contain side chains with high theta values or be connected to a turn/loop with high theta values. A definition for a protein core is proposed, and the implications for protein folding are discussed. Apparently, during folding and in the native state, nonlocal interactions between core sequences are favored more than other possible nonlocal interactions. Other studies of partially folded bovine pancreatic trypsin inhibitor (Barbar E, Barany G, Woodward C, 1995, Biochemistry 34:11423-11434; Barber E, Hare M, Daragan V, Barany G, Woodward C, 1998, Biochemistry 37:7822-7833), suggest that developing cores have site-specific energy barriers between microstates, one disordered, and the other(s) more ordered.

Autonomous folding of a peptide corresponding to the N-terminal beta-hairpin from ubiquitin

The N-terminal 17 residues of ubiquitin have been shown by 1H NMR to fold autonomously into a beta-hairpin structure in aqueous solution. This structure has a specific, native-like register, though side-chain contacts differ in detail from those observed in the intact protein. An autonomously folding hairpin has previously been identified in the case of streptococcal protein G, which is structurally homologous with ubiquitin, but remarkably, the two are not in topologically equivalent positions in the fold. This suggests that the organization of folding may be quite different for proteins sharing similar tertiary structures. Two smaller peptides have also been studied, corresponding to the isolated arms of the N-terminal hairpin of ubiquitin, and significant differences from simple random coil predictions observed in the spectra of these subfragments, suggestive of significant limitation of the backbone conformational space sampled, presumably as a consequence of the strongly beta-structure favoring composition of the sequences. This illustrates the ability of local sequence elements to express a propensity for beta-structure even in the absence of actual sheet formation. Attempts were made to estimate the population of the folded state of the hairpin, in terms of a simple two-state folding model. Using published "random coil" values to model the unfolded state, and values derived from native ubiquitin for the putative unique, folded state, it was found that the apparent population varied widely for different residues and with different NMR parameters. Use of the spectra of the subfragment peptides to provide a more realistic model of the unfolded state led to better agreement in the estimates that could be obtained from chemical shift and coupling constant measurements, while making it clear that some other approaches to population estimation could not give meaningful results, because of the tendency to populate the beta-region of conformational space even in the absence of the hairpin structure.

Unfolding of proteins monitored by electrospray ionization mass spectrometry: a comparison of positive and negative ion modes

Electrospray ionization (ESI) mass spectrometry (MS) in both the positive and negative ion mode has been used to study protein unfolding transitions of lysozyme, cytochrome c (cyt c), and ubiquitin in solution. As expected, ESI of unfolded lysozyme leads to the formation of substantially higher charge states than the tightly folded protein in both modes of operation. Surprisingly, the acid-induced unfolding of cyt c as well as the acid and the base-induced unfolding of ubiquitin show different behavior: In these three cases protein unfolding only leads to marginal changes in the negative ion charge state distributions, whereas in the positive ion mode pronounced shifts to higher charge states are observed. This shows that ESI MS in the negative ion mode as a method for probing conformational changes of proteins in solution should be treated with caution. The data presented in this work provide further evidence that the conformation of a protein in solution not its charge state is the predominant factor for determining the ESI charge state distribution in the positive ion mode. Furthermore, these data support the hypothesis of a recent study (Konermann and Douglas, Biochemistry 1997, 36, 12,296-12,302) which suggested that ESI in the positive ion mode is not sensitive to changes in the secondary structure of proteins but only to changes in the tertiary structure.

A strategy for detecting the conservation of folding-nucleus residues in protein superfamilies

BACKGROUND

Nucleation-growth theory predicts that fast-folding peptide sequences fold to their native structure via structures in a transition-state ensemble that share a small number of native contacts (the folding nucleus). Experimental and theoretical studies of proteins suggest that residues participating in folding nuclei are conserved among homologs. We attempted to determine if this is true in proteins with highly diverged sequences but identical folds (superfamilies).

RESULTS

We describe a strategy based on comparisons of residue conservation in natural superfamily sequences with simulated sequences (generated with a Monte-Carlo sequence design strategy) for the same proteins. The basic assumptions of the strategy were that natural sequences will conserve residues needed for folding and stability plus function, the simulated sequences contain no functional conservation, and nucleus residues make native contacts with each other. Based on these assumptions, we identified seven potential nucleus residues in ubiquitin superfamily members. Non-nucleus conserved residues were also identified; these are proposed to be involved in stabilizing native interactions. We found that all superfamily members conserved the same potential nucleus residue positions, except those for which the structural topology is significantly different.

CONCLUSIONS

Our results suggest that the conservation of the nucleus of a specific fold can be predicted by comparing designed simulated sequences with natural highly diverged sequences that fold to the same structure. We suggest that such a strategy could be used to help plan protein folding and design experiments, to identify new superfamily members, and to subdivide superfamilies further into classes having a similar folding mechanism.