Two folded conformers of ubiquitin revealed by high-pressure NMR

Paramagnetic relaxation enhancement-assisted structural characterization of a partially disordered conformation of ubiquitin

Nuclear magnetic resonance (NMR) is a powerful tool to study three-dimensional structures as well as protein conformational fluctuations in solution, but it is compromised by increases in peak widths and missing signals. We previously reported that ubiquitin has two folded conformations, N₁ and N₂ and plus another folded conformation, I, in which some amide group signals of residues 33-41 almost disappeared above 3 kbar at pH 4.5 and 273 K. Thus, well-converged structural models could not be obtained for this region owing to the absence of distance restraints. Here, we reexamine the problem using the ubiquitin Q41N variant as a model for this locally disordered conformation, I. We demonstrate that the variant shows pressure-induced loss of backbone amide group signals at residues 28, 33, 36, and 39-41 like the wild-type, with a similar but smaller effect on CαH and CβH signals. In order to characterize this I structure, we measured paramagnetic relaxation enhancement (PRE) under high pressure to obtain distance restraints, and calculated the structure assisted by Bayesian inference. We conclude that the more disordered I conformation observed at pH 4.0, 278 K, and 2.5 kbar largely retained the N₂ conformation, although the amide groups at residues 33-41 have more heterogeneous conformations and more contact with water, which differ from the N₁ and N₂ states. The PRE-assisted strategy has the potential to improve structural characterization of proteins that lack NMR signals, especially for relatively more open and hydrated protein conformations.

Molecular dynamics simulation of proteins under high pressure: Structure, function and thermodynamics

BACKGROUND

Molecular dynamics (MD) simulation is well-recognized as a powerful tool to investigate protein structure, function, and thermodynamics. MD simulation is also used to investigate high pressure effects on proteins. For conducting better MD simulation under high pressure, the main issues to be addressed are: (i) protein force fields and water models were originally developed to reproduce experimental properties obtained at ambient pressure; and (ii) the timescale to observe the pressure effect is often much longer than that of conventional MD simulations.

SCOPE OF REVIEW

First, we describe recent developments in MD simulation methodologies for studying the high-pressure structure and dynamics of protein molecules. These developments include force fields for proteins and water molecules, and enhanced simulation techniques. Then, we summarize recent studies of MD simulations of proteins in water under high pressure.

MAJOR CONCLUSIONS

Recent MD simulations of proteins in solution under pressure have reproduced various phenomena identified by experiments using high pressure, such as hydration, water penetration, conformational change, helix stabilization, and molecular stiffening.

GENERAL SIGNIFICANCE

MD simulations demonstrate differences in the properties of proteins and water molecules between ambient and high-pressure conditions. Comparing the results obtained by MD calculations with those obtained experimentally could reveal the mechanism by which biological molecular machines work well in collaboration with water molecules.

Molecular simulations by generalized-ensemble algorithms in isothermal-isobaric ensemble

Generalized-ensemble algorithms are powerful techniques for investigating biomolecules such as protein, DNA, lipid membrane, and glycan. The generalized-ensemble algorithms were originally developed in the canonical ensemble. On the other hand, not only temperature but also pressure is controlled in experiments. Additionally, pressure is used as perturbation to study relationship between function and structure of biomolecules. For this reason, it is important to perform efficient conformation sampling based on the isothermal-isobaric ensemble. In this article, we review a series of the generalized-ensemble algorithms in the isothermal-isobaric ensemble: multibaric-multithermal, pressure- and temperature-simulated tempering, replica-exchange, and replica-permutation methods. These methods achieve more efficient simulation than the conventional isothermal-isobaric simulation. Furthermore, the isothermal-isobaric generalized-ensemble simulation samples conformations of biomolecules from wider range of temperature and pressure. Thus, we can estimate physical quantities more accurately at any temperature and pressure values. The applications to the biomolecular system are also presented.

Characterization of low-lying excited states of proteins by high-pressure NMR

Hydrostatic pressure alters the free energy of proteins by a few kJ mol^-1, with the amount depending on their partial molar volumes. Because the folded ground state of a protein contains cavities, it is always a state of large partial molar volume. Therefore pressure always destabilises the ground state and increases the population of partially and completely unfolded states. This is a mild and reversible conformational change, which allows the study of excited states under thermodynamic equilibrium conditions. Many of the excited states studied in this way are functionally relevant; they also seem to be very similar to kinetic folding intermediates, thus suggesting that evolution has made use of the 'natural' dynamic energy landscape of the protein fold and sculpted it to optimise function. This includes features such as ligand binding, structural change during the catalytic cycle, and dynamic allostery.

Water-Protein Interactions Coupled with Protein Conformational Transition

Conformational fluctuations of proteins are crucially important for their functions. However, changes in the location and dynamics of hydrated water in many proteins accompanied by the conformational transition have not been fully understood. Here, we used phase-modulated clean chemical exchange NMR approach to investigate pressure-induced changes in water-to-amide proton exchange occurring at sub-second time scale. With the transition of ubiquitin from its native conformation (N1) to an alternative conformation (N2) at 250 MPa, proton exchange rates of residues 32-35, 40-41, and 71, which are located at the C-terminal side of the protein, were significantly increased. These observations can be explained by the destabilization of the hydrogen bonds in the backbone and partial exposure of those amide groups to solvent in N2. We conclude that phase-modulated clean chemical exchange NMR approach coupled with pressure perturbation will be a useful tool for investigations of more open and hydrated protein structures.

Monitoring protein folding through high pressure NMR spectroscopy

High-pressure is a well-known perturbation method used to destabilize globular proteins. It is perfectly reversible, which is essential for a proper thermodynamic characterization of a protein equilibrium. In contrast to other perturbation methods such as heat or chemical denaturant that destabilize protein structures uniformly, pressure exerts local effects on regions or domains of a protein containing internal cavities. When combined with NMR spectroscopy, hydrostatic pressure offers the possibility to monitor at a residue level the structural transitions occurring upon unfolding and to determine the kinetic properties of the process. High-pressure NMR experiments can now be routinely performed, owing to the recent development of commercially available high-pressure sample cells. This review summarizes recent advances and some future directions of high-pressure NMR techniques for the characterization at atomic resolution of the energy landscape of protein folding.

Interactions Controlling the Slow Dynamic Conformational Motions of Ubiquitin

Rational mutation of proteins based on their structural and dynamic characteristics is a useful strategy for amplifying specific fluctuations in proteins. Here, we show the effects of mutation on the conformational fluctuations and thermodynamic stability of ubiquitin. In particular, we focus on the salt bridge between K11 and E34 and the hydrogen bond between I36 and Q41, which are predicted to control the fluctuation between the basic folded state, N₁, and the alternatively folded state, N₂, of the protein, using high-pressure NMR spectroscopy. The E34A mutation, which disrupts the salt bridge, did not alter picosecond–to–nanosecond, microsecond–to–millisecond dynamic motions, and stability of the protein, while the Q41N mutation, which destabilizes the hydrogen bond, specifically amplified the N₁–N₂ conformational fluctuation and decreased stability. Based on the observed thermodynamic stabilities of the various conformational states, we showed that in the Q41N mutant, the N₁ state is more significantly destabilized than the N₂ state, resulting in an increase in the relative population of N₂. Identifying the interactions controlling specific motions of a protein will facilitate molecular design to achieve functional dynamics beyond native state dynamics.

High-pressure NMR techniques for the study of protein dynamics, folding and aggregation

High-pressure is a well-known perturbation method used to destabilize globular proteins and dissociate protein complexes or aggregates. The heterogeneity of the response to pressure offers a unique opportunity to dissect the thermodynamic contributions to protein stability. In addition, pressure perturbation is generally reversible, which is essential for a proper thermodynamic characterization of a protein equilibrium. When combined with NMR spectroscopy, hydrostatic pressure offers the possibility of monitoring at an atomic resolution the structural transitions occurring upon unfolding and determining the kinetic properties of the process. The recent development of commercially available high-pressure sample cells greatly increased the potential applications for high-pressure NMR experiments that can now be routinely performed. This review summarizes the recent applications and future directions of high-pressure NMR techniques for the characterization of protein conformational fluctuations, protein folding and the stability of protein complexes and aggregates.

The unfolding effects on the protein hydration shell and partial molar volume: a computational study

In this paper we apply the computational analysis recently proposed by our group to characterize the solvation properties of a native protein in aqueous solution, and to four model aqueous solutions of globular proteins in their unfolded states thus characterizing the protein unfolded state hydration shell and quantitatively evaluating the protein unfolded state partial molar volumes.

Cavities and Excited States in Proteins

Protein cavities or voids are observed as defects in atomic packing. Cavities have long been suggested to play important roles in protein dynamics and function, but the underlying origin and mechanism remains elusive. Here, recent studies about the cavities characterized by high-pressure NMR spectroscopy have been reviewed. Analysis of the pressure-dependent chemical shifts showed both linear and nonlinear response of proteins to pressure. The linear response corresponded to compression within the native ensemble, while the nonlinear response indicated the involvement of low-lying excited states that were different from the native state. The finding of non-linear pressure shifts in various proteins suggested that the existence of the low-lying excited states was common for globular proteins. However, the absolute nonlinear coefficient values varied significantly from protein to protein, and showed a good correlation with the density of cavities. Extensive studies on hen lysozyme as a model system showed that the cavity hydration and water penetration into the interior of proteins was an origin of the conformational transition to the excited states. The importance of cavities for protein function and evolution has also been explained. In addition to these "equilibrium" cavities, there are also "transient" cavities formed in the interior of the protein structure, as manifested by the ring flip motions of aromatic rings. The significance of transient cavities, reflecting an intrinsic dynamic nature within the native state, has also been discussed.

High-Pressure NMR Spectroscopy Reveals Functional Sub-states of Ubiquitin and Ubiquitin-Like Proteins

High-pressure nuclear magnetic resonance (NMR) spectroscopy has revealed that ubiquitin has at least two high-energy states--an alternatively folded state N2 and a locally disordered state I--between the basic folded state N1 and totally unfolded U state. The high-energy states are conserved among ubiquitin-like post-translational modifiers, ubiquitin, NEDD8, and SUMO-2, showing the E1-E2-E3 cascade reaction. It is quite intriguing that structurally similar high-energy states are evolutionally conserved in the ubiquitin-like modifiers, and the thermodynamic stabilities vary among the proteins. To investigate atomic details of the high-energy states, a Q41N mutant of ubiquitin was created as a structural model of N2, which is 71% populated even at atmospheric pressure. The convergent structure of the "pure" N2 state was obtained by nuclear Overhauser effect (NOE)-based structural analysis of the Q41N mutant at 2.5 kbar, where the N2 state is 97% populated. The N2 state of ubiquitin is closely similar to the conformation of the protein bound to the ubiquitin-activating enzyme E1. The recognition of E1 by ubiquitin is best explained by conformational selection rather than by induced-fit motion.

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.

Impact of hydrostatic pressure on an intrinsically disordered protein: a high-pressure NMR study of α-synuclein

The impact of pressure on the backbone (15) N, (1) H and (13) C chemical shifts in N-terminally acetylated α-synuclein has been evaluated over a pressure range 1-2500 bar. Even while the chemical shifts fall very close to random coil values, as expected for an intrinsically disordered protein, substantial deviations in the pressure dependence of the chemical shifts are seen relative to those in short model peptides. In particular, the nonlinear pressure response of the (1) H(N) chemical shifts, which commonly is associated with the presence of low-lying "excited states", is much larger in α-synuclein than in model peptides. The linear pressure response of (1) H(N) chemical shift, commonly linked to H-bond length change, correlates well with those in short model peptides, and is found to be anticorrelated with its temperature dependence. The pressure dependence of (13) C chemical shifts shows remarkably large variations, even when accounting for residue type, and do not point to a clear shift in population between different regions of the Ramachandran map. However, a nearly universal decrease in (3) JHN-Hα by 0.22 ± 0.05 Hz suggests a slight increase in population of the polyproline II region at 2500 bar. The first six residues of N-terminally acetylated synuclein show a transient of approximately 15% population of α-helix, which slightly diminishes at 2500 bar. The backbone dynamics of the protein is not visibly affected beyond the effect of slight increase in water viscosity at 2500 bar.

Pressure-Induced Helical Structure of a Peptide Studied by Simulated Tempering Molecular Dynamics Simulations

It is known experimentally that an AK16 peptide forms more α-helix structures with increasing pressure while proteins unfold in general. In order to understand this abnormality, molecular dynamics (MD) simulations with the simulated tempering method for the isobaric-isothermal ensemble were performed in a wide pressure range from 1.0 × 10(-4) GPa to 1.4 GPa. From the results of the simulations, it is found that the fraction of the folded state decreases once and increases after that with increasing pressure. The partial molar volume change from the folded state to unfolded state increases monotonically from a negative value to a positive value with pressure. The behavior under high pressure conditions is consistent with the experimental results. The radius of gyration of highly helical structures decreases with increasing pressure, which indicates that the helix structure shrinks with pressure. This is the reason why the fraction of the folded state increases as pressure increases.

High-pressure NMR reveals close similarity between cold and alcohol protein denaturation in ubiquitin

Proteins denature not only at high, but also at low temperature as well as high pressure. These denatured states are not easily accessible for experiment, because usually heat denaturation causes aggregation, whereas cold or pressure denaturation occurs at temperatures well below the freezing point of water or pressures above 5 kbar, respectively. Here we have obtained atomic details of the pressure-assisted, cold-denatured state of ubiquitin at 2,500 bar and 258 K by high-resolution NMR techniques. Under these conditions, a folded, native-like and a disordered state exist in slow exchange. Secondary chemical shifts show that the disordered state has structural propensities for a native-like N-terminal β-hairpin and α-helix and a nonnative C-terminal α-helix. These propensities are very similar to the previously described alcohol-denatured (A-)state. Similar to the A-state, (15)N relaxation data indicate that the secondary structure elements move as independent segments. The close similarity of pressure-assisted, cold-denatured, and alcohol-denatured states with native and nonnative secondary elements supports a hierarchical mechanism of folding and supports the notion that similar to alcohol, pressure and cold reduce the hydrophobic effect. Indeed, at nondenaturing concentrations of methanol, a complete transition from the native to the A-state can be achieved at ambient temperature by varying the pressure from 1 to 2,500 bar. The methanol-assisted pressure transition is completely reversible and can also be induced in protein G. This method should allow highly detailed studies of protein-folding transitions in a continuous and reversible manner.

Counterbalance of ligand- and self-coupled motions characterizes multispecificity of ubiquitin

Date hub proteins are a type of proteins that show multispecificity in a time-dependent manner. To understand dynamic aspects of such multispecificity we studied Ubiquitin as a typical example of a date hub protein. Here we analyzed 9 biologically relevant Ubiquitin-protein (ligand) heterodimer structures by using normal mode analysis based on an elastic network model. Our result showed that the self-coupled motion of Ubiquitin in the complex, rather than its ligand-coupled motion, is similar to the motion of Ubiquitin in the unbound condition. The ligand-coupled motions are correlated to the conformational change between the unbound and bound conditions of Ubiquitin. Moreover, ligand-coupled motions favor the formation of the bound states, due to its in-phase movements of the contacting atoms at the interface. The self-coupled motions at the interface indicated loss of conformational entropy due to binding. Therefore, such motions disfavor the formation of the bound state. We observed that the ligand-coupled motions are embedded in the motions of unbound Ubiquitin. In conclusion, multispecificity of Ubiquitin can be characterized by an intricate balance of the ligand- and self-coupled motions, both of which are embedded in the motions of the unbound form.

Remodeling of the folding free energy landscape of staphylococcal nuclease by cavity-creating mutations

The folding of staphylococcal nuclease (SNase) is known to proceed via a major intermediate in which the central OB subdomain is folded and the C-terminal helical subdomain is disordered. To identify the structural and energetic determinants of this folding free energy landscape, we have examined in detail, using high-pressure NMR, the consequences of cavity creating mutations in each of the two subdomains of an ultrastable SNase, Δ+PHS. The stabilizing mutations of Δ+PHS enhanced the population of the major folding intermediate. Cavity creation in two different regions of the Δ+PHS reference protein, despite equivalent effects on global stability, had very distinct consequences on the complexity of the folding free energy landscape. The L125A substitution in the C-terminal helix of Δ+PHS slightly suppressed the major intermediate and promoted an additional excited state involving disorder in the N-terminus, but otherwise decreased landscape heterogeneity with respect to the Δ+PHS background protein. The I92A substitution, located in the hydrophobic OB-fold core, had a much more profound effect, resulting in a significant increase in the number of intermediate states and implicating the entire protein structure. Denaturant (GuHCl) had very subtle and specific effects on the landscape, suppressing some states and favoring others, depending upon the mutational context. These results demonstrate that disrupting interactions in a region of the protein with highly cooperative, unfrustrated folding has very profound effects on the roughness of the folding landscape, whereas the effects are less pronounced for an energetically equivalent substitution in an already frustrated region.

Extensively hydrated but folded: a novel state of globular proteins stabilized at high pressure and low temperature

We studied conformational fluctuations of the transcription factor c-Myb R2 subdomain (52 residues with three Trp) at high pressure and low temperature (5°C) using two different spectroscopic methods, Trp fluorescence and (1)H NMR, on its chemically stable mutant C130I (pseudo-wild-type (WT(S))), which has a large internal cavity. As pressure was increased from 3 to 300 MPa, the Trp fluorescence λ(max) of WT(S) shifted from 342 to ∼355 nm, clearly showing that the three Trp rings become fully exposed to the polar environment, which usually is taken to indicate that the protein underwent unfolding. In contrast, as pressure was increased from 3 to 300 MPa, the high-field-shifted (1)H NMR signals characteristic of the folded state showed a still higher-field shift, but no significant changes in their intensity. The last result unequivocally shows that the protein remains largely folded at 300 MPa. The apparent discrepancy between the two predictions would only be solved if one were to postulate the existence of an extensively hydrated but folded state in WT(S). Intriguingly, such a state was not found in a cavity-filling mutant of WT(S), C130I/V103L, suggesting that this state is mediated by cavity hydration. The generality and significance of this state in proteins are discussed.

High-pressure protein crystallography and NMR to explore protein conformations

High-pressure methods for solving protein structures by X-ray crystallography and NMR are maturing. These techniques are beginning to impact our understanding of thermodynamic and structural features that define not only the protein's native conformation, but also the higher free energy conformations. The ability of high-pressure methods to visualize these mostly unexplored conformations provides new insight into protein function and dynamics. In this review, we begin with a historical discussion of high-pressure structural studies, with an eye toward early results that paved the way to mapping the multiple conformations of proteins. This is followed by an examination of several recent studies that emphasize different strengths and uses of high-pressure structural studies, ranging from basic thermodynamics to the suggestion of high-pressure structural methods as a tool for protein engineering.

Computing the stability diagram of the Trp-cage miniprotein

We report molecular dynamics simulations of the equilibrium folding/unfolding thermodynamics of an all-atom model of the Trp-cage miniprotein in explicit solvent. Simulations are used to sample the folding/unfolding free energy difference and its derivatives along 2 isochores. We model the DeltaG(u)(P,T) landscape using the simulation data and propose a stability diagram model for Trp-cage. We find the proposed diagram to exhibit features similar to globular proteins with increasing hydrostatic pressure destabilizing the native fold. The observed energy differences DeltaE(u) are roughly linearly temperature-dependent and approach DeltaE(u) = 0 with decreasing temperature, suggesting that the system approached the region of cold denaturation. In the low-temperature denatured state, the native helical secondary structure elements are largely preserved, whereas the protein conformation changes to an "open-clamp" configuration. A tighter packing of water around nonpolar sites, accompanied by an increasing solvent-accessible surface area of the unfolded ensemble, seems to stabilize the unfolded state at elevated pressures.

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.

Theoretical study of the partial molar volume change associated with the pressure-induced structural transition of ubiquitin

The partial molar volume (PMV) change associated with the pressure-induced structural transition of ubiquitin is analyzed by the three-dimensional reference interaction site model (3D-RISM) theory of molecular solvation. The theory predicts that the PMV decreases upon the structural transition, which is consistent with the experimental observation. The volume decomposition analysis demonstrates that the PMV reduction is primarily caused by the decrease in the volume of structural voids in the protein, which is partially canceled by the volume expansion due to the hydration effects. It is found from further analysis that the PMV reduction is ascribed substantially to the penetration of water molecules into a specific part of the protein. Based on the thermodynamic relation, this result implies that the water penetration causes the pressure-induced structural transition. It supports the water penetration model of pressure denaturation of proteins proposed earlier.

A general algorithm for peak-tracking in multi-dimensional NMR experiments

We present an algorithmic method allowing automatic tracking of NMR peaks in a series of spectra. It consists in a two phase analysis. The first phase is a local modeling of the peak displacement between two consecutive experiments using distance matrices. Then, from the coefficients of these matrices, a value graph containing the a priori set of possible paths used by these peaks is generated. On this set, the minimization under constraint of the target function by a heuristic approach provides a solution to the peak-tracking problem. This approach has been named GAPT, standing for General Algorithm for NMR Peak Tracking. It has been validated in numerous simulations resembling those encountered in NMR spectroscopy. We show the robustness and limits of the method for situations with many peak-picking errors, and presenting a high local density of peaks. It is then applied to the case of a temperature study of the NMR spectrum of the Lipid Transfer Protein (LTP).

The transition state of the ras binding domain of Raf is structurally polarized based on Phi-values but is energetically diffuse

The ras binding domain (RBD) of the Ser/Thr kinase c-Raf/Raf-1 spans 78 residues and adopts a structure characteristic of the beta-grasp ubiquitin-like topology. Recently, the primary sequence of Raf RBD has been nearly exhaustively mutated experimentally by insertion of stretches of degenerate codons, which revealed sequence conservation and hydrophobic core organization similar to that found in an alignment of beta-grasp ubiquitin-like proteins. These results now allow us to examine the relationship between sequence conservation and the folding process, particularly viewed through the analysis of transition state (TS) structure. Specifically, we present herein a protein engineering study combining classic truncation (Ala/Gly) and atypical mutants to predict folding TS ensemble properties. Based on classical Phi-value analysis, Raf RBD TS structure is particularly polarized around the N-terminal beta-hairpin. However, all residues constituting the inner layer of the hydrophobic core are involved in TS stabilization, although they are clearly found in a less native-like environment. The TS structure can also be probed by a direct measure of its destabilization upon mutation, DeltaDeltaG(U-++). Viewed through this analysis, Raf RBD TS is a more diffuse structure, in which all residues of the hydrophobic core including beta-strands 1, 2, 3 and 5 and the major alpha-helix play similar roles in TS stabilization. In addition, Phi-values and DeltaDeltaG(U-++) reveal striking similarities in the TS of Raf RBD and ubiquitin, a structural analogue displaying insignificant sequence identity (<12%). However, ubiquitin TS appears more denatured-like and polarized around the N-terminal beta-hairpin. We suggest that analysis of Phi-values should also consider the direct impact of mutations on differences in free energy between the unfolded and TS (DeltaDeltaG(U-++)) to ensure that the description of TS properties is accurate. Finally, the impact of these findings on the modeling of protein folding is discussed.

Massive sequence perturbation of the Raf ras binding domain reveals relationships between sequence conservation, secondary structure propensity, hydrophobic core organization and stability

The contributions of specific residues to the delicate balance between function, stability and folding rates could be determined, in part by [corrected] comparing the sequences of structures having identical folds, but insignificant sequence homology. Recently, we have devised an experimental strategy to thoroughly explore residue substitutions consistent with a specific class of structure. Using this approach, the amino acids tolerated at virtually all residues of the c-Raf/Raf1 ras binding domain (Raf RBD), an exemplar of the common beta-grasp ubiquitin-like topology, were obtained and used to define the sequence determinants of this fold. Herein, we present analyses suggesting that more subtle sequence selection pressure, including propensity for secondary structure, the hydrophobic core organization and charge distribution are imposed on the Raf RBD sequence. Secondly, using the Gibbs free energies (DeltaG(F-U)) obtained for 51 mutants of Raf RBD, we demonstrate a strong correlation between amino acid conservation and the destabilization induced by truncating mutants. In addition, four mutants are shown to significantly stabilize Raf RBD native structure. Two of these mutations, including the well-studied R89L, are known to severely compromise binding affinity for ras. Another stabilized mutant consisted of a deletion of amino acid residues E104-K106. This deletion naturally occurs in the homologues a-Raf and b-Raf and could indicate functional divergence. Finally, the combination of mutations affecting five of 78 residues of Raf RBD results in stabilization of the structure by approximately 12 kJ mol(-1) (DeltaG(F-U) is -22 and -34 kJ mol(-1) for wt and mutant, respectively). The sequence perturbation approach combined with sequence/structure analysis of the ubiquitin-like fold provide a basis for the identification of sequence-specific requirements for function, stability and folding rate of the Raf RBD and structural analogues, highlighting the utility of conservation profiles as predictive tools of structural organization.

Evolutionally conserved intermediates between ubiquitin and NEDD8

The investigation of common structural motifs provides additional information on why proteins conserve similar topologies yet may have non-conserved amino acid sequences. Proteins containing the ubiquitin superfold have similar topologies, although the sequence conservation is rather poor. Here, we present novel similarities and differences between the proteins ubiquitin and NEDD8. They have 57% identical sequence, almost identical backbone topology and similar functional strategy, although their physiological functions are mutually different. Using variable pressure NMR spectroscopy, we found that the two proteins have similar conformational fluctuation in the evolutionary conserved enzyme-binding region and contain a structurally similar locally disordered conformer (I) in equilibrium with the basic folded conformer (N). A notable difference between the two proteins is that the equilibrium population of I is far greater for NEDD8 (DeltaG(0)(NI)<5 kJ/mol) than for ubiquitin (DeltaG(0)(NI)=15.2(+/-1.0) kJ/mol), and that the tendency for overall unfolding (U) is also far higher for NEDD8 (DeltaG(0)(NU)=11.0(+/-1.5) kJ/mol) than for ubiquitin (DeltaG(0)(NU)=31.3(+/-4.7) kJ/mol). These results suggest that the marked differences in thermodynamic stabilities of the locally disordered conformer (I) and the overall unfolding species (U) are a key to determine the functional differences of the two structurally similar proteins in physiology.

Cold denaturation of ubiquitin at high pressure

Cold-induced conformational transition of ubiquitin was studied at pH 4.5 under a constant pressure of 2 kbar using variable pressure one-dimensional 1H and two-dimensional 15N/1H NMR spectroscopy as well as IR spectroscopy. Although a tendency for preferential stabilization of a peculiar locally disordered and partially hydrated conformer I, identical with that previously found with variable-pressure NMR at 0 degrees C, is recognized, the transition of the folded conformer N to the unfolded conformer U occurs largely cooperatively with decreasing temperature, reaching near completion at - 21 degrees C. NMR spectral features as well as the analysis of NMR relaxation parameters indicate that the polypeptide chain is almost fully unfolded, fairly well-hydrated and floppy at - 21 degrees C, whereas the IR spectrum shows a substantial decrease of the beta-sheet. The Gibbs energy change from the folded state (a mixture of N and I) to the unfolded state at 2 kbar obtained from the 1H NMR data is fitted well with a single DeltaCp value of 2.43 +/- 0.13 (kJ/K mol) for the entire temperature range between - 21 and 90 degrees C, covering both the cold denaturation and heat denaturation, showing that the two denatured states actually belong to a single thermodynamic phase of the protein. The DeltaCp value determined at 2 kbar is substantially smaller than the DeltaCp determined at 1 bar (3.8-5.8 (kJ/Kmol), which is consistent with the fact that the denaturation takes place from a mixture of N and I at 2 kbar rather than from pure N at 1 bar.

Conformational fluctuations of proteins revealed by variable pressure NMR

With the high-resolution variable-pressure NMR spectroscopy, one can study conformational fluctuations of proteins in a much wider conformational space than hitherto explored by NMR and other spectroscopic techniques. This is because a protein in solution generally exists as a dynamic mixture of conformers mutually differing in partial molar volume, and pressure can select the population of a conformer according to its relative volume. In this review, we describe how variable-pressure NMR can be used to probe conformational fluctuations of proteins in a wide conformational space from the folded to the fully unfolded structures, with actual examples. Furthermore, the newly emerging technique "NMR snapshots" expresses amply fluctuating protein structures as changes in atomic coordinates. Finally, the concept of conformational fluctuation is extended to include intermolecular association leading to amyloidosis.

Combined NMR-observation of cold denaturation in supercooled water and heat denaturation enables accurate measurement of deltaC(p) of protein unfolding

Cold and heat denaturation of the double mutant Arg 3-->Glu/Leu 66-->Glu of cold shock protein Csp of Bacillus caldolyticus was monitored using 1D (1)H NMR spectroscopy in the temperature range from -12 degrees C in supercooled water up to +70 degrees C. The fraction of unfolded protein, f (u), was determined as a function of the temperature. The data characterizing the unfolding transitions could be consistently interpreted in the framework of two-state models: cold and heat denaturation temperatures were determined to be -11 degrees C and 39 degrees C, respectively. A joint fit to both cold and heat transition data enabled the accurate spectroscopic determination of the heat capacity difference between native and denatured state, DeltaC(p) of unfolding. The approach described in this letter, or a variant thereof, is generally applicable and promises to be of value for routine studies of protein folding.

High Resolution1H NMR Structural Studies of Sucrose Octaacetate in Supercritical Carbon Dioxide

High pressure (HP), high resolution (HR), proton nuclear magnetic resonance (1H NMR) spectroscopy has been utilized for the first time to investigate the solution structure of a carbohydrate based system, sucrose octaacetate (SOA), in supercritical CO2. The studies indicate that the average solution state conformation of the alpha-D-Glucopyranosyl ring of SOA in scCO2 medium is consistent with the 4C1 chair form, while the beta-D-fructofuranosyl ring adopts an envelope conformation. The investigations also suggest that scCO2 is a promising medium to study the solution structure and conformation of acetylated sugar systems. Spectral manifestations of a specific interaction between the acetate methyl protons and CO2 molecules are also presented.

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.

Molecular dynamics simulations of HPr under hydrostatic pressure

The histidine-containing protein (HPr) plays an important role in the phosphotransferase system (PTS). The deformations induced on the protein structure at high hydrostatic pressure values (4, 50, 100, 150, and 200 MPa) were previously (H. Kalbitzer, A. Görler, H. Li, P. Dubovskii, A. Hengstenberg, C. Kowolik, H. Yamada, and K. Akasaka, Protein Science 2000, Vol. 9, pp. 693-703) analyzed by NMR experiments: the nonlinear variations of the amide chemical shifts at high pressure values were supposed to arise from induced shifts in the protein conformational equilibrium. Molecular dynamics (MD) simulations are here performed, to analyze the protein internal mobility at 0.1 MPa, and to relate the nonlinear variations of chemical shifts observed at high pressure, to variations in conformational equilibrium. The global features of the protein structure are only slightly modified along the pressure. Nevertheless, the values of the Voronoi residues volumes show that the residues of alpha-helices are more compressed that those belonging to the beta-sheet. The alpha-helices are also displaying the largest internal mobility and deformation in the simulations. The nonlinearity of the 1H chemical shifts, computed from the MD simulation snapshots, is in qualitative agreement with the nonlinearity of the experimentally observed chemical shifts.

How large are the volume changes accompanying protein transitions and binding?

We present a simple model to describe volume changes accompanying protein folding and binding events. The model enables one to resolve the changes in volume accompanying conformational transitions of proteins as well as association of proteins with other molecules in terms of the intrinsic, thermal and interaction (hydration) contributions. The thermal contribution to protein volume results from thermally activated mutual vibrational motions of contacting solute and solvent molecules. Our calculations suggest that near zero volume changes accompanying protein folding and binding events reflect compensation between significant changes in the intrinsic, thermal and interaction terms. We have quantitatively estimated these terms as a function of the protein's molecular weight and degree of its unfolding. Results described in this work lay foundation for more reliable and physically justified interpretations of volumetric data on protein folding and binding events. We also discuss potential ways of extending applications of our model to analyzing other macromolecular systems and events, including drug-DNA and protein-DNA interactions and helix-to-helix and helix-to-coil transitions of nucleic acids.

Close identity of a pressure-stabilized intermediate with a kinetic intermediate in protein folding

Atomic detailed structural study of a transiently existing folding intermediate is severely limited because of its short life. In ubiquitin, we found that a pressure-stabilized equilibrium conformer shares a common structural feature with the proline-trapped kinetic intermediate found in a pulse-labeling (1)H(2)H exchange NMR study [Briggs, M. S. & Roder, H. (1992) Proc. Natl. Acad. Sci. USA 89, 2017-2021]. The conformer is locally unfolded in the entire segment from residues 33 to 42 and in C-terminal residues 70-76. The close structural identity of an equilibrium intermediate stabilized under pressure with a transiently observed folding intermediate is likely to be general in terms of a folding funnel common to both experiments.

Ras-effector interactions: after one decade

Ras effectors have convergently developed a common subdomain in their otherwise unrelated protein body for their interaction with Ras. Structural analysis revealed that the mode of interaction is highly similar for all Ras effectors, but is completely different from that of effectors of other subfamilies of small GTPases. Whereas the molecular mechanism of effector activation is still elusive, detailed knowledge about the thermodynamics and dynamics of the interaction with Ras has accumulated.

High-pressure 1H NMR study of pressure-induced structural changes in the heme environments of metcyanomyoglobins

The effect of pressure on the heme environment structure of sperm whale and horse heart metcyanomyoglobins was investigated up to 300 MPa by high-pressure (1)H NMR spectroscopy. Pressure-induced changes in the distances between the observed protons and the heme iron atom were estimated from changes in the dipolar shift due to the paramagnetic effect on the protons. The changes showed that the heme peripheral structure as a whole was compressed by pressure; the movements of the protons in the heme peripheral residues were in the range of +0.16 to -0.54 A/300 MPa. One-dimensional compressibilities for the protons, excluding the protons of the distal His residue, were in the range of 1.0 x 10(-4) to 6.1 x 10(-4)/MPa. The movements of the protons induced by pressure correlated well with the distance between the protons and cavities in the protein. The distal His residue (His 64) moved toward the outside of the heme pocket, but remained in the pocket even at 300 MPa. This movement was driven dominantly by a change in the dihedral angle around the C(alpha)-C(beta) rotational bond of the residue. Comparative work on horse heart metcyanomyoglobin implied that the conformational change of the His 64 imidazole ring was larger in the horse heart metcyanomyoglobin than in the sperm whale metcyanomyoglobin.

Investigations into the chromatographic behavior of a doxorubicin-peptide conjugate

HPLC impurity profile method development for a doxorubicin-heptapeptide conjugate included significant changes of the separation profile with diluent, eluent and pH. These separation variables were also temperature-dependent with a shift in retention from 35 to 45 degrees C. There was also a direct relationship of temperature with LC retention, and a pH minimum at 5.9. Atypical dependence of the impurity profile on diluent at a k' of 18 led to further investigation. A large change in retention by several minutes was a function of both the organic eluent composition and temperature between 15 and 30 degrees C. Several Van't Hoff temperature studies from 5 to 65 degrees C on several column types resulted in non-linear plots. Analysis of the molecular subunits suggested that the peptide portion of the analyte influenced the non-linear retention behavior. The stationary phase type was not a significant factor causing non-linearity. Circular dichroism-temperature studies indicated a notable transition in ellipticity for the amine regions (198-202 nm) that occurred between 39 and 44 degrees C. This transition temperature range coincided with the results of the Van't Hoff analysis, between 35 and 44 degrees C, to indicate that these effects were not primarily stationary phase induced.

Equilibrium and pressure-jump relaxation studies of the conformational transitions of P13MTCP1

The conformational transitions of a small oncogene product, p13(MTCP1), have been studied by high-pressure fluorescence of the intrinsic tryptophan emission and high-pressure 1D and 2D 1H-15N NMR. While the unfolding transition monitored by fluorescence is cooperative, two kinds of NMR spectral changes were observed, depending on the pressure range. Below approximately 200 MPa, pressure caused continuous, non-linear shifts of many of the 15N and 1H signals, suggesting the presence of an alternate folded conformer(s) in rapid equilibrium (tau<<ms) with the basic native structure. Above approximately 200 MPa, pressure caused a sharp decrease in the intensity of the folded proteins signals, while the peaks corresponding to disordered structures increased, yielding a free energy of unfolding change of 6.0 kcal/mol and associated volume change of -100 ml/mol, in agreement with the fluorescence result. Differential scanning calorimetry also reveals two transitions between 21 and 65 degrees C, confirming the existence of an additional species under mildly denaturing conditions. We report here a real-time observation of pressure-jump unfolding kinetics by 2D NMR spectroscopy on P13MTCP1 made possible due to its very long relaxation times at high pressure revealed by fluorescence studies. Within the dead-time after the pressure-jump, the NMR spectra of the native conformer changed to those of the transient conformational species, identified in the equilibrium studies, demonstrating the equivalence between a transient species and an equilibrium excited state. After these rapid spectral changes, the intensities of all of the individual 15N-1H cross-peaks decreased gradually, and those of the disordered structure increased, consistent with the slow relaxation to the unfolded form at this pressure. Rate constants of unfolding monitored at individual amide sites within the beta-barrel were similar to those obtained from fluorescence and from side-chain protons in the hydrophobic core region, consistent with nearly cooperative unfolding. However, some heterogeneity in the apparent unfolding rate constants is apparent across the sequence and can be understood as non-uniform effects of pressure on the unfolding rate constant due to non-uniform hydration.

High pressure NMR reveals that apomyoglobin is an equilibrium mixture from the native to the unfolded

Pressure-induced reversible conformational changes of sperm whale apomyoglobin have been studied between 30 bar and 3000 bar on individual residue basis by utilizing 1H/15N hetero nuclear single-quantum coherence two-dimensional NMR spectroscopy at pH 6.0 and 35 degrees C. Apomyoglobin showed a series of pressure-dependent NMR spectra as a function of pressure, assignable to the native (N), intermediates (I), molten globule (MG) and unfolded (U) conformers. At 30 bar, the native fold (N) shows disorder only in the F helix. Between 500 bar and 1200 bar, a series of locally disordered conformers I are produced, in which local disorder occurs in the C helix, the CD loop, the G helix and part of the H helix. At 2000 bar, most cross-peaks exhibit severe line-broadening, suggesting the formation of a molten globule, but at 3000 bar all the cross-peaks reappear, showing that the molten globule turns into a well-hydrated, mobile unfolded conformation U. Since all the spectral changes were reversible with pressure, apomyoglobin is considered to exist as an equilibrium mixture of the N, I, MG and U conformers at all pressures. MG is situated at 2.4+/-(0.1) kcal/mol above N at 1 bar and the unfolding transition from the combined N-I state to MG is accompanied by a loss of partial molar volume by 75+/-(3) ml/mol. On the basis of these observations, we postulate a theorem that the partial molar volume of a protein decreases in parallel with the loss of its conformational order.