NMR Spectroscopic Characterization of Millisecond Protein Folding by Transverse Relaxation Dispersion Measurements

Including the Ensemble of Unstructured Conformations in the Analysis of Protein's Native State by High-Pressure NMR Spectroscopy

The analysis of pressure induced changes in the chemical shift of proteins allows statements on structural fluctuations proteins exhibit at ambient pressure. The inherent issue of separating general pressure effects from structural related effects on the pressure dependence of chemical shifts has so far been addressed by considering the characteristics of random coil peptides on increasing pressure. In this work, chemically and pressure denatured states of the cold shock protein B from Bacillus subtilis (BsCspB) have been assigned in 2D 1H-15N HSQC NMR spectra and their dependence on increasing hydrostatic pressure has been evaluated. The pressure denatured polypeptide chain has been used to separate general from structural related effects on 1H and 15N chemical shifts of native BsCspB and the implications on the interpretation of pressure induced changes in the chemical shift regarding the structure of BsCspB are discussed. It has been found that the ensemble of unstructured conformations of BsCspB shows different responses to increasing pressure than random coil peptides do. Thus, the approach used for considering the general effects that arise when hydrostatic pressure increases changes the structural conclusions that are drawn from high pressure NMR spectroscopic experiments that rely on the analysis of chemical shifts.

Monitoring protein unfolding transitions by NMR-spectroscopy

NMR-spectroscopy has certain unique advantages for recording unfolding transitions of proteins compared e.g. to optical methods. It enables per-residue monitoring and separate detection of the folded and unfolded state as well as possible equilibrium intermediates. This allows a detailed view on the state and cooperativity of folding of the protein of interest and the correct interpretation of subsequent experiments. Here we summarize in detail practical and theoretical aspects of such experiments. Certain pitfalls can be avoided, and meaningful simplification can be made during the analysis. Especially a good understanding of the NMR exchange regime and relaxation properties of the system of interest is beneficial. We show by a global analysis of signals of the folded and unfolded state of GB1 how accurate values of unfolding can be extracted and what limits different NMR detection and unfolding methods. E.g. commonly used exchangeable amides can lead to a systematic under determination of the thermodynamic protein stability. We give several perspectives of how to deal with more complex proteins and how the knowledge about protein stability at residue resolution helps to understand protein properties under crowding conditions, during phase separation and under high pressure.

Conformational exchange of aromatic side chains by 1H CPMG relaxation dispersion

Aromatic side chains are attractive probes of protein dynamics on the millisecond time scale, because they are often key residues in enzyme active sites and protein binding sites. Further they allow to study specific processes, like histidine tautomerization and ring flips. Till now such processes have been studied by aromatic ¹³C CPMG relaxation dispersion experiments. Here we investigate the possibility of aromatic ¹H CPMG relaxation dispersion experiments as a complementary method. Artifact-free dispersions are possible on uniformly ¹H and ¹³C labeled samples for histidine δ2 and ε1, as well as for tryptophan δ1. The method has been validated by measuring fast folding-unfolding kinetics of the small protein CspB under native conditions. The determined rate constants and populations agree well with previous results from ¹³C CPMG relaxation dispersion experiments. The CPMG-derived chemical shift differences between the folded and unfolded states are in good agreement with those obtained directly from the spectra. In contrast, the ¹H relaxation dispersion profiles in phenylalanine, tyrosine and the six-ring moiety of tryptophan, display anomalous behavior caused by ³J ¹H-¹H couplings and, if present, strong ¹³C-¹³C couplings. Therefore they require site-selective ¹H/²H and, in case of strong couplings, ¹³C/¹²C labeling. In summary, aromatic ¹H CPMG relaxation dispersion experiments work on certain positions (His δ2, His ε1 and Trp δ1) in uniformly labeled samples, while other positions require site-selective isotope labeling.

Equilibrium and Kinetic Unfolding of GB1: Stabilization of the Native State by Pressure

NMR spectroscopy allows an all-atom view on pressure-induced protein folding, separate detection of different folding states, determination of their population, and the measurement of the folding kinetics at equilibrium. Here, we studied the folding of protein GB1 at pH 2 in a temperature and pressure dependent way. We find that the midpoints of temperature-induced unfolding increase with higher pressure. NMR relaxation dispersion experiments disclosed that the unfolding kinetics slow down at elevated pressure while the folding kinetics stay virtually the same. Therefore, pressure is stabilizing the native state of GB1. These findings extend the knowledge of the influence of pressure on protein folding kinetics, where so far typically a destabilization by increased activation volumes of folding was observed. Our findings thus point toward an exceptional section in the pressure-temperature phase diagram of protein unfolding. The stabilization of the native state could potentially be caused by a shift of p K_a values of glutamates and aspartates in favor of the negatively charged state as judged from pH sensitive chemical shifts.

RNA binding and chaperone activity of the E. coli cold-shock protein CspA

Ensuring the correct folding of RNA molecules in the cell is of major importance for a large variety of biological functions. Therefore, chaperone proteins that assist RNA in adopting their functionally active states are abundant in all living organisms. An important feature of RNA chaperone proteins is that they do not require an external energy source to perform their activity, and that they interact transiently and non-specifically with their RNA targets. So far, little is known about the mechanistic details of the RNA chaperone activity of these proteins. Prominent examples of RNA chaperones are bacterial cold shock proteins (Csp) that have been reported to bind single-stranded RNA and DNA. Here, we have used advanced NMR spectroscopy techniques to investigate at atomic resolution the RNA-melting activity of CspA, the major cold shock protein of Escherichia coli, upon binding to different RNA hairpins. Real-time NMR provides detailed information on the folding kinetics and folding pathways. Finally, comparison of wild-type CspA with single-point mutants and small peptides yields insights into the complementary roles of aromatic and positively charged amino-acid side chains for the RNA chaperone activity of the protein.

Protein folding by NMR

Protein folding is a highly complex process proceeding through a number of disordered and partially folded nonnative states with various degrees of structural organization. These transiently and sparsely populated species on the protein folding energy landscape play crucial roles in driving folding toward the native conformation, yet some of these nonnative states may also serve as precursors for protein misfolding and aggregation associated with a range of devastating diseases, including neuro-degeneration, diabetes and cancer. Therefore, in vivo protein folding is often reshaped co- and post-translationally through interactions with the ribosome, molecular chaperones and/or other cellular components. Owing to developments in instrumentation and methodology, solution NMR spectroscopy has emerged as the central experimental approach for the detailed characterization of the complex protein folding processes in vitro and in vivo. NMR relaxation dispersion and saturation transfer methods provide the means for a detailed characterization of protein folding kinetics and thermodynamics under native-like conditions, as well as modeling high-resolution structures of weakly populated short-lived conformational states on the protein folding energy landscape. Continuing development of isotope labeling strategies and NMR methods to probe high molecular weight protein assemblies, along with advances of in-cell NMR, have recently allowed protein folding to be studied in the context of ribosome-nascent chain complexes and molecular chaperones, and even inside living cells. Here we review solution NMR approaches to investigate the protein folding energy landscape, and discuss selected applications of NMR methodology to studying protein folding in vitro and in vivo. Together, these examples highlight a vast potential of solution NMR in providing atomistic insights into molecular mechanisms of protein folding and homeostasis in health and disease.

Protein dynamics and function from solution state NMR spectroscopy

It is well-established that dynamics are central to protein function; their importance is implicitly acknowledged in the principles of the Monod, Wyman and Changeux model of binding cooperativity, which was originally proposed in 1965. Nowadays the concept of protein dynamics is formulated in terms of the energy landscape theory, which can be used to understand protein folding and conformational changes in proteins. Because protein dynamics are so important, a key to understanding protein function at the molecular level is to design experiments that allow their quantitative analysis. Nuclear magnetic resonance (NMR) spectroscopy is uniquely suited for this purpose because major advances in theory, hardware, and experimental methods have made it possible to characterize protein dynamics at an unprecedented level of detail. Unique features of NMR include the ability to quantify dynamics (i) under equilibrium conditions without external perturbations, (ii) using many probes simultaneously, and (iii) over large time intervals. Here we review NMR techniques for quantifying protein dynamics on fast (ps-ns), slow (μs-ms), and very slow (s-min) time scales. These techniques are discussed with reference to some major discoveries in protein science that have been made possible by NMR spectroscopy.

Solid-state NMR approaches to internal dynamics of proteins: from picoseconds to microseconds and seconds

Solid-state nuclear magnetic resonance (NMR) spectroscopy has matured to the point that it is possible to determine the structure of proteins in immobilized states, such as within microcrystals or embedded in membranes. Currently, researchers continue to develop and apply NMR techniques that can deliver site-resolved dynamic information toward the goal of understanding protein function at the atomic scale. As a widely-used, natural approach, researchers have mostly measured longitudinal (T1) relaxation times, which, like in solution-state NMR, are sensitive to picosecond and nanosecond motions, and motionally averaged dipolar couplings, which provide an integral amplitude of all motions with a correlation time of up to a few microseconds. While overall Brownian tumbling in solution mostly precludes access to slower internal dynamics, dedicated solid-state NMR approaches are now emerging as powerful new options. In this Account, we give an overview of the classes of solid-state NMR experiments that have expanded the accessible range correlation times from microseconds to many milliseconds. The measurement of relaxation times in the rotating frame, T1ρ, now allows researchers to access the microsecond range. Using our recent theoretical work, researchers can now quantitatively analyze this data to distinguish relaxation due to chemical-shift anisotropy (CSA) from that due to dipole-dipole couplings. Off-resonance irradiation allows researchers to extend the frequency range of such experiments. We have built multidimensional analogues of T2-type or line shape experiments using variants of the dipolar-chemical shift correlation (DIPSHIFT) experiment that are particularly suited to extract intermediate time scale motions in the millisecond range. In addition, we have continuously improved variants of exchange experiments, mostly relying on the recoupling of anisotropic interactions to address ultraslow motions in the ms to s ranges. The NH dipolar coupling offers a useful probe of local dynamics, especially with proton-depleted samples that suppress the adverse effect of strong proton dipolar couplings. We demonstrate how these techniques have provided a concise picture of the internal dynamics in a popular model system, the SH3 domain of α-spectrin. T1-based methods have shown that large-amplitude bond orientation fluctuations in the picosecond range and slower 10 ns low-amplitude motions coexist in these structures. When we include T1ρ data, we observe that many residues undergo low amplitude motions slower than 100 ns. On the millisecond to second scale, mostly localized but potentially cooperative motions occur. Comparing different exchange experiments, we found that terminal NH2 groups in side chains can even undergo a combination of ultraslow large-angle two-site jumps accompanied by small-angle fluctuations that occur 10 times more quickly.

Structural and dynamic features of cold-shock proteins of Listeria monocytogenes, a psychrophilic bacterium

Cold-shock proteins (Csps), proteins expressed when the ambient temperature drops below the growth-supporting temperature, bind to single-stranded nucleic acids and act as RNA chaperones to regulate translation. Listeria monocytogenes is a psychrophilic food-borne pathogen that is problematic for the food industry. Structures of Csps from psychrophilic bacteria have not yet been studied. Despite dramatic differences in the thermostability of Csps of various thermophilic microorganisms, these proteins share a high degree of primary sequence homology and a high degree of three-dimensional structural similarity. Here, we investigated the structural and dynamic features as well as the thermostability of L. monocytogenes CspA (Lm-CspA). Lm-CspA has a five-stranded β-barrel structure with hydrophobic core packing and two salt bridges. When heptathymidine (dT(7)) binds, values for the heteronuclear nuclear Overhauser effect and order parameters of residues in surface loop regions near nucleic acid binding sites increase dramatically. Moreover, Carr-Purcell-Meiboom-Gill experiments showed that slow motions observed for the nucleic acid binding residues K7, W8, F15, F27, and R56 disappeared in Lm-CspA-dT(7). Lm-CspA is less thermostable than mesophilic and thermophilic Csps, with a lower melting temperature (40 °C). The structural flexibility that accompanies longer surface loops and less hydrophobic core packing and a number of salt bridges and unfavorable electrostatic repulsion are likely key factors in the low thermostability of Lm-CspA. This implies that the large conformational flexibility of psychrophilic Lm-CspA, which more easily accommodates nucleic acids at low temperature, is required for RNA chaperone function under cold-shock conditions and for the cold adaptation of L. monocytogenes.

Conformational exchange of aromatic side chains characterized by L-optimized TROSY-selected ¹³C CPMG relaxation dispersion

Protein dynamics on the millisecond time scale commonly reflect conformational transitions between distinct functional states. NMR relaxation dispersion experiments have provided important insights into biologically relevant dynamics with site-specific resolution, primarily targeting the protein backbone and methyl-bearing side chains. Aromatic side chains represent attractive probes of protein dynamics because they are over-represented in protein binding interfaces, play critical roles in enzyme catalysis, and form an important part of the core. Here we introduce a method to characterize millisecond conformational exchange of aromatic side chains in selectively (13)C labeled proteins by means of longitudinal- and transverse-relaxation optimized CPMG relaxation dispersion. By monitoring (13)C relaxation in a spin-state selective manner, significant sensitivity enhancement can be achieved in terms of both signal intensity and the relative exchange contribution to transverse relaxation. Further signal enhancement results from optimizing the longitudinal relaxation recovery of the covalently attached (1)H spins. We validated the L-TROSY-CPMG experiment by measuring fast folding-unfolding kinetics of the small protein CspB under native conditions. The determined unfolding rate matches perfectly with previous results from stopped-flow kinetics. The CPMG-derived chemical shift differences between the folded and unfolded states are in excellent agreement with those obtained by urea-dependent chemical shift analysis. The present method enables characterization of conformational exchange involving aromatic side chains and should serve as a valuable complement to methods developed for other types of protein side chains.

Probing microsecond time scale dynamics in proteins by methyl (1)H Carr-Purcell-Meiboom-Gill relaxation dispersion NMR measurements. Application to activation of the signaling protein NtrC(r)

To study microsecond processes by relaxation dispersion NMR spectroscopy, low power deposition and short pulses are crucial and encourage the development of experiments that employ (1)H Carr-Purcell-Meiboom-Gill (CPMG) pulse trains. Herein, a method is described for the comprehensive study of microsecond to millisecond time scale dynamics of methyl groups in proteins, exploiting their high abundance and favorable relaxation properties. In our approach, protein samples are produced using [(1)H, (13)C]-d-glucose in ∼100% D(2)O, which yields CHD(2) methyl groups for alanine, valine, threonine, isoleucine, leucine, and methionine residues with high abundance, in an otherwise largely deuterated background. Methyl groups in such samples can be sequence-specifically assigned to near completion, using (13)C TOCSY NMR spectroscopy, as was recently demonstrated (Otten, R.; et al. J. Am. Chem. Soc. 2010, 132, 2952-2960). In this Article, NMR pulse schemes are presented to measure (1)H CPMG relaxation dispersion profiles for CHD(2) methyl groups, in a vein similar to that of backbone relaxation experiments. Because of the high deuteration level of methyl-bearing side chains, artifacts arising from proton scalar coupling during the CPMG pulse train are negligible, with the exception of Ile-δ1 and Thr-γ2 methyl groups, and a pulse scheme is described to remove the artifacts for those residues. Strong (13)C scalar coupling effects, observed for several leucine residues, are removed by alternative biochemical and NMR approaches. The methodology is applied to the transcriptional activator NtrC(r), for which an inactive/active state transition was previously measured and the motions in the microsecond time range were estimated through a combination of backbone (15)N CPMG dispersion NMR spectroscopy and a collection of experiments to determine the exchange-free component to the transverse relaxation rate. Exchange contributions to the (1)H line width were detected for 21 methyl groups, and these probes were found to collectively report on a local structural rearrangement around the phosphorylation site, with a rate constant of (15.5 ± 0.5) × 10(3) per second (i.e., τ(ex) = 64.7 ± 1.9 μs). The affected methyl groups indicate that, already before phosphorylation, a substantial, transient rearrangement takes place between helices 3 and 4 and strands 4 and 5. This conformational equilibrium allows the protein to gain access to the active, signaling state in the absence of covalent modification through a shift in a pre-existing dynamic equilibrium. Moreover, the conformational switching maps exactly to the regions that differ between the solution NMR structures of the fully inactive and active states. These results demonstrate that a cost-effective and quantitative study of protein methyl group dynamics by (1)H CPMG relaxation dispersion NMR spectroscopy is possible and can be applied to study functional motions on the microsecond time scale that cannot be accessed by backbone (15)N relaxation dispersion NMR. The use of methyl groups as dynamics probes extends such applications also to larger proteins.

Mutational investigation of protein folding transition states by Phi-value analysis and beyond: lessons from SH3 domain folding

Understanding how proteins adopt their unique native structures requires a complete structural characterization of the rate-limiting transition state(s) along the folding pathway. By definition, transition states are not significantly populated and are only accessible via folding kinetics studies. In this respect, interpreting the kinetic effects of amino acid substitutions (especially to Ala) via Phi-value analysis is the most common method to probe the structure of these transient, yet important states. A critical review of the key assumptions required for rigorous interpretation of Phi values reveals that a multiple substitution strategy in which a position of interest is mutated to a variety of amino acids, and not exclusively to Ala, provides the best means to characterize folding transition states. This approach has proven useful in revealing non-native interactions and (or) conformations in folding transition states. Moreover, by simultaneously examining the folding kinetics of multiple substitutions made at a single surface-exposed position using the Brønsted analysis the backbone conformation in a folding transition state can be investigated. For folding equilibria with exchange rates on the order of milliseconds, the kinetic parameters for Phi-value analysis can be obtained from NMR relaxation dispersion experiments, under fully native conditions, along with a wealth of high-resolution structural information about the states in exchange (native, denatured, and intermediate states that populate the pathway). This additional structural information, which is not readily obtained through stopped-flow based methods, can significantly facilitate the interpretation of Phi values because it often reports on the validity of the assumptions required for a rigorous interpretation of Phi values.

Kicked by Mos and tuned by MPF-the initiation of the MAPK cascade in Xenopus oocytes

The mitogen-activated protein kinase (MAPK) cascade is a paradigmatic signaling cascade, which plays a crucial role in many aspects of cellular events. The main initiator of the cascade in Xenopus oocytes is the oncoprotein Mos. After activation of the cascade, Mos activity is stabilized by MAPK via a feedback loop. Mos concentration levels are, however, not controlled by MAPK alone. In this paper we show, by imposing either a sustained or a peaked activity of M-phase promoting factor (MPF) (Cdc2-cyclin B), how the latter regulates the dynamics of Mos. Our experiments are supported by a detailed kinetic model for the Mos-MPF-MAPK network, which takes into account the three different phosphorylation states of Mos and, as a consequence, allows us to determine the time evolution of Mos under control of MPF. Our work opens a path toward a more complete and biologically realistic quantitative understanding of the dynamic interdependence of Mos and MPF in Xenopus oocytes.

Isotope labeling methods for studies of excited protein states by relaxation dispersion NMR spectroscopy

The utility of nuclear magnetic resonance (NMR) spectroscopy as a tool for the study of biomolecular structure and dynamics has benefited from the development of facile labeling methods that incorporate NMR active probes at key positions in the molecule. Here we describe a protocol for the labeling of proteins that facilitates their study using a technique that is sensitive to millisecond conformational exchange processes. The samples necessary for an analysis of exchange dynamics are discussed, using the Abp1p SH3 domain from Saccharomyces cerevisiae as an example. For this system, the time frame for production of each sample, including in vitro refolding, is about 80 h. The samples so produced facilitate the measurement of accurate chemical shifts of low populated, invisible conformers that are part of the exchange pathway. The accuracy of the methodology has been established experimentally and the chemical shifts that are obtained provide important restraints in structure calculations of the excited state.

Measuring 13Cbeta chemical shifts of invisible excited states in proteins by relaxation dispersion NMR spectroscopy

A labeling scheme is introduced that facilitates the measurement of accurate (13)C(beta) chemical shifts of invisible, excited states of proteins by relaxation dispersion NMR spectroscopy. The approach makes use of protein over-expression in a strain of E. coli in which the TCA cycle enzyme succinate dehydrogenase is knocked out, leading to the production of samples with high levels of (13)C enrichment (30-40%) at C(beta) side-chain carbon positions for 15 of the amino acids with little (13)C label at positions one bond removed (approximately 5%). A pair of samples are produced using [1-(13)C]-glucose/NaH(12)CO(3) or [2-(13)C]-glucose as carbon sources with isolated and enriched (>30%) (13)C(beta) positions for 11 and 4 residues, respectively. The efficacy of the labeling procedure is established by NMR spectroscopy. The utility of such samples for measurement of (13)C(beta) chemical shifts of invisible, excited states in exchange with visible, ground conformations is confirmed by relaxation dispersion studies of a protein-ligand binding exchange reaction in which the extracted chemical shift differences from dispersion profiles compare favorably with those obtained directly from measurements on ligand free and fully bound protein samples.

Accurate measurement of alpha proton chemical shifts of excited protein states by relaxation dispersion NMR spectroscopy

Carr-Purcell-Meiboom-Gill relaxation dispersion NMR spectroscopy can provide detailed information about low populated, invisible states of protein molecules, including backbone chemical shifts of the invisible conformer and bond vector orientations that can be used as structural constraints. Notably, the measurement of 1Halpha chemical shifts in excited protein states has not been possible to date because, in the absence of suitable labeling, the homonuclear proton scalar coupling network in side chains of proteins leads to a significant degradation in the performance of proton-based relaxation dispersion experiments. Here we have overcome this problem through a labeling scheme in which proteins are prepared with U-2H glucose and 50% D2O/50% H2O that results in deuteration levels of between 50-88% at the Cbeta carbon. Effects from residual 1Halpha-1Hbeta scalar couplings can be suppressed through a new NMR experiment that is presented here. The utility of the methodology is demonstrated on a ligand binding exchanging system and it is shown that 1Halpha chemical shifts extracted from dispersion profiles are, on average, accurate to 0.03 ppm, an order of magnitude better than they can be predicted from structure using a database approach. The ability to measure 1Halpha chemical shifts of invisible conformers is particularly important because such shifts are sensitive to both secondary and tertiary structure. Thus, the methodology presented is a valuable addition to a growing list of experiments for characterizing excited protein states that are difficult to study using the traditional techniques of structural biology.

Measurement of carbonyl chemical shifts of excited protein states by relaxation dispersion NMR spectroscopy: comparison between uniformly and selectively (13)C labeled samples

Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion nuclear magnetic resonance (NMR) spectroscopy has emerged as a powerful method for quantifying chemical shifts of excited protein states. For many applications of the technique that involve the measurement of relaxation rates of carbon magnetization it is necessary to prepare samples with isolated (13)C spins so that experiments do not suffer from magnetization transfer between coupled carbon spins that would otherwise occur during the CPMG pulse train. In the case of (13)CO experiments however the large separation between (13)CO and (13)C(alpha) chemical shifts offers hope that robust (13)CO dispersion profiles can be recorded on uniformly (13)C labeled samples, leading to the extraction of accurate (13)CO chemical shifts of the invisible, excited state. Here we compare such chemical shifts recorded on samples that are selectively labeled, prepared using [1-(13)C]-pyruvate and NaH(13)CO(3,) or uniformly labeled, generated from (13)C-glucose. Very similar (13)CO chemical shifts are obtained from analysis of CPMG experiments recorded on both samples, and comparison with chemical shifts measured using a second approach establishes that the shifts measured from relaxation dispersion are very accurate.

Using relaxation dispersion NMR spectroscopy to determine structures of excited, invisible protein states

Currently the main focus of structural biology is the determination of static three-dimensional representations of biomolecules that for the most part correspond to low energy (ground state) conformations. However, it is becoming increasingly well recognized that higher energy structures often play important roles in function as well. Because these conformers are populated to only low levels and are often only transiently formed their study is not amenable to many of the tools of structural biology. In this perspective we discuss the role of CPMG-based relaxation dispersion NMR spectroscopy in characterizing these low populated, invisible states. It is shown that robust methods for measuring both backbone chemical shifts and residual anisotropic interactions in the excited state are in place and that these data provide valuable restraints for structural studies of invisible conformers.

Probing chemical shifts of invisible states of proteins with relaxation dispersion NMR spectroscopy: how well can we do?

Carr-Purcell-Meiboom-Gill relaxation dispersion NMR spectroscopy has evolved into a powerful approach for the study of low populated, invisible conformations of biological molecules. One of the powerful features of the experiment is that chemical shift differences between the exchanging conformers can be obtained, providing structural information about invisible excited states. Through the development of new labeling approaches and NMR experiments it is now possible to measure backbone 13C(alpha) and 13CO relaxation dispersion profiles in proteins without complications from 13C-13C couplings. Such measurements are presented here, along with those that probe exchange using 15N and 1HN nuclei. A key experimental design has been the choice of an exchanging system where excited-state chemical shifts were known from independent measurement. Thus it is possible to evaluate quantitatively the accuracy of chemical shift differences obtained in dispersion experiments and to establish that in general very accurate values can be obtained. The experimental work is supplemented by computations that suggest that similarly accurate shifts can be measured in many cases for systems with exchange rates and populations that fall within the range of those that can be quantified by relaxation dispersion. The accuracy of the extracted chemical shifts opens up the possibility of obtaining quantitative structural information of invisible states of the sort that is now available from chemical shifts recorded on ground states of proteins.

Protein folding and unfolding by all-atom molecular dynamics simulations

Computational protein folding can be classified into pathway and sampling approaches. Here, we use the AMBER simulation package as an example to illustrate the protocols for all-atom molecular simulations of protein folding, including system setup, simulation, and analysis. We introduced two traditional pathway approaches: ab inito folding and high-temperature unfolding. The popular replica exchange method was chosen to represent sampling approaches. Our emphasis is placed on the analysis of the simulation trajectories, and some in-depth discussions are provided for commonly encountered problems.

Separating degenerate (1)H transitions in methyl group probes for single-quantum (1)H-CPMG relaxation dispersion NMR spectroscopy

A relaxation dispersion pulse scheme is presented for quantifying chemical exchange processes in proteins that exploits 1H chemical shifts as probes of changes in conformation. The experiment selects 1H single-quantum magnetization from the I = 1/2 manifolds of the methyl group, which behave like AX spin systems, while suppressing coherences that derive from the 3/2 manifold that are extremely sensitive to pulse imperfections and that would otherwise severely compromise the accuracy of the experiment. The utility of the sequence is first demonstrated with an application to a protein system that is known not to undergo chemical exchange and flat dispersion profiles are obtained. Subsequently, the methodology is applied to study the folding of a G48M mutant of the Fyn SH3 domain that has been shown previously to undergo exchange between folded and unfolded states on the millisecond time scale.

The folding pathway of an FF domain: characterization of an on-pathway intermediate state under folding conditions by (15)N, (13)C(alpha) and (13)C-methyl relaxation dispersion and (1)H/(2)H-exchange NMR spectroscopy

The FF domain from the human protein HYPA/FBP11 folds via a low-energy on-pathway intermediate (I). Elucidation of the structure of such folding intermediates and denatured states under conditions that favour folding are difficult tasks. Here, we investigated the millisecond time-scale equilibrium folding transition of the 71-residue four-helix bundle wild-type protein by (15)N, (13)C(alpha) and methyl(13)C Carr-Purcell-Meiboom-Gill (CPMG) NMR relaxation dispersion experiments and by (1)H/(2)H-exchange measurements. The relaxation data for the wild-type protein fitted a simple two-site exchange process between the folded state (F) and I. Destabilization of F in mutants A17G and Q19G allowed the detection of the unfolded state U by (15)N CPMG relaxation dispersion. The dispersion data for these mutants fitted a three-site exchange scheme, U<-->I<-->F, with I populated higher than U. The kinetics and thermodynamics of the folding reaction were obtained via temperature and urea-dependent relaxation dispersion experiments, along with structural information on I from backbone (15)N, (13)C(alpha) and side-chain methyl (13)C chemical shifts, with further information from protection factors for the backbone amide groups from (1)H/(2)H-exchange. Notably, helices H1-H3 are at least partially formed in I, while helix H4 is largely disordered. Chemical shift differences for the methyl (13)C nuclei suggest a paucity of stable, native-like hydrophobic interactions in I. These data are consistent with Phi-analysis of the rate-limiting transition state between I and F. The combination of relaxation dispersion and Phi data can elucidate whole experimental folding pathways.

Intermediates: ubiquitous species on folding energy landscapes?

Although intermediates have long been recognised as fascinating species that form during the folding of large proteins, the role that intermediates play in the folding of small, single-domain proteins has been widely debated. Recent discoveries using new, sensitive methods of detection and studies combining simulation and experiment have now converged on a common vision for folding, involving intermediates as ubiquitous stepping stones en route to the native state. The results suggest that the folding energy landscapes of even the smallest proteins possess significant ruggedness in which intermediates stabilized by both native and non-native interactions are common features.

Modification and optimization of the united-residue (UNRES) potential energy function for canonical simulations. I. Temperature dependence of the effective energy function and tests of the optimization method with single training proteins

We report the modification and parametrization of the united-residue (UNRES) force field for energy-based protein structure prediction and protein folding simulations. We tested the approach on three training proteins separately: 1E0L (beta), 1GAB (alpha), and 1E0G (alpha + beta). Heretofore, the UNRES force field had been designed and parametrized to locate native-like structures of proteins as global minima of their effective potential energy surfaces, which largely neglected the conformational entropy because decoys composed of only lowest-energy conformations were used to optimize the force field. Recently, we developed a mesoscopic dynamics procedure for UNRES and applied it with success to simulate protein folding pathways. However, the force field turned out to be largely biased toward -helical structures in canonical simulations because the conformational entropy had been neglected in the parametrization. We applied the hierarchical optimization method, developed in our earlier work, to optimize the force field; in this method, the conformational space of a training protein is divided into levels, each corresponding to a certain degree of native-likeness. The levels are ordered according to increasing native-likeness; level 0 corresponds to structures with no native-like elements, and the highest level corresponds to the fully native-like structures. The aim of optimization is to achieve the order of the free energies of levels, decreasing as their native-likeness increases. The procedure is iterative, and decoys of the training protein(s) generated with the energy function parameters of the preceding iteration are used to optimize the force field in a current iteration. We applied the multiplexing replica-exchange molecular dynamics (MREMD) method, recently implemented in UNRES, to generate decoys; with this modification, conformational entropy is taken into account. Moreover, we optimized the free-energy gaps between levels at temperatures corresponding to a predominance of folded or unfolded structures, as well as to structures at the putative folding-transition temperature, changing the sign of the gaps at the transition temperature. This enabled us to obtain force fields characterized by a single peak in the heat capacity at the transition temperature. Furthermore, we introduced temperature dependence to the UNRES force field; this is consistent with the fact that it is a free-energy and not a potential energy function. beta

A Thiol Labelling Competition Experiment as a Probe for Sidechain Packing in the Kinetic Folding Intermediate of N-PGK

Protein folding is directed by the sequence of sidechains along the polypeptide backbone, but despite this the developement of sidechain interactions during folding is not well understood. Here, the thiol-active reagent, dithio-nitrobenzoic acid (DTNB), is used to probe the exposure of the cysteine sidechain thiols in the kinetic folding intermediates of the N-terminal domain of phosphoglycerate kinase (N-PGK) and a number of conservative (I-, L-, or V-to-C) single cysteine variants. Rapid dilution of chemically denatured protein into folding conditions in the presence of DTNB allowed the degree of sidechain protection in any rapidly formed intermediate to be determined through the analysis of the kinetics of labelling. The protection factors derived for the intermediate(s) were generally small (<25), indicating only partial burial of the sidechains. The distribution of protection parallels the previously reported backbone amide protection for the folding intermediate of N-PGK. These observations are consistent with the hypothesis that such intermediates resemble molten globule states; i.e. with native-like backbone hydrogen bonding and overall tertiary structure, but with the sidechains that make up the hydrophobic protein core dynamic and intermittently solvent exposed. The success of the competition technique in characterizing this kinetic intermediate invites application to other model systems.

Abp1p and Fyn SH3 domains fold through similar low-populated intermediate states

Src homology 3 (SH3) domains are small modules that are thought to fold via a two-state mechanism, without the accumulation of significant populations of intermediate states. Relaxation dispersion NMR studies of the folding of G48V and G48M mutants of the Fyn SH3 domain have established that, at least for these modules, folding proceeds through the formation of a transient on-pathway intermediate with an equilibrium population of 1-2% that can be readily detected [Korzhnev, D. M., et al. (2004) Nature 430, 586-590]. To investigate the generality of this result, we present an (15)N relaxation dispersion NMR study of a pair of additional SH3 domains, including a G48V mutant of a stabilized Abp1p SH3 domain that shares 36% sequence identity with the Fyn SH3 module, and a A39V/N53P/V55L mutant Fyn SH3 domain. A transient folding intermediate is detected for both of the proteins studied here, and the dispersion data are well fit to a folding model of the form F <--> I <--> U, where F, I, and U correspond to folded, intermediate, and unfolded states, respectively. The temperature dependencies of the folding/unfolding rate constants were obtained so that the thermodynamic properties of each of F, I, and U could be established. The detection of I states in folding pathways of all SH3 domains examined to date via relaxation dispersion NMR spectroscopy indicates that such intermediates may well be a conserved feature in the folding of such domains in general but that their transient nature along with their low population makes detection difficult using more well-established approaches to the study of folding.

Recognition of T-rich single-stranded DNA by the cold shock protein Bs-CspB in solution

Cold shock proteins (CSP) belong to the family of single-stranded nucleic acid binding proteins with OB-fold. CSP are believed to function as 'RNA chaperones' and during anti-termination. We determined the solution structure of Bs-CspB bound to the single-stranded DNA (ssDNA) fragment heptathymidine (dT7) by NMR spectroscopy. Bs-CspB reveals an almost invariant conformation when bound to dT7 with only minor reorientations in loop beta1-beta2 and beta3-beta4 and of few aromatic side chains involved in base stacking. Binding studies of protein variants and mutated ssDNA demonstrated that Bs-CspB associates with ssDNA at almost diffusion controlled rates and low sequence specificity consistent with its biological function. A variation of the ssDNA affinity is accomplished solely by changes of the dissociation rate. 15N NMR relaxation and H/D exchange experiments revealed that binding of dT7 increases the stability of Bs-CspB and reduces the sub-nanosecond dynamics of the entire protein and especially of loop beta3-beta4.

Similarity and difference in the unfolding of thermophilic and mesophilic cold shock proteins studied by molecular dynamics simulations

Molecular dynamics simulations were performed to unfold a homologous pair of thermophilic and mesophilic cold shock proteins at high temperatures. The two proteins differ in just 11 of 66 residues and have very similar structures with a closed five-stranded antiparallel beta-barrel. A long flexible loop connects the N-terminal side of the barrel, formed by three strands (beta1-beta3), with the C-terminal side, formed by two strands (beta4-beta5). The two proteins were found to follow the same unfolding pathway, but with the thermophilic protein showing much slower unfolding. Unfolding started with the melting of C-terminal strands, leading to exposure of the hydrophobic core. Subsequent melting of beta3 and the beta-hairpin formed by the first two strands then resulted in unfolding of the whole protein. The slower unfolding of the thermophilic protein could be attributed to ion pair formation of Arg-3 with Glu-46, Glu-21, and the C-terminal. These ion pairs were also found to be important for the difference in folding stability between the pair of proteins. Thus electrostatic interactions appear to play similar roles in the difference in folding stability and kinetics between the pair of proteins.

The inverted chevron plot measured by NMR relaxation reveals a native-like unfolding intermediate in acyl-CoA binding protein

The folding kinetics of bovine acyl-CoA binding protein was studied by 15N relaxation dispersion measurements under equilibrium conditions. Relaxation dispersion profiles were measured at several concentrations of guanidine hydrochloride (GuHCl). The unfolding rate constant (k(u)) was determined under conditions favoring folding, for which the folding rate constant (k(f)) dominates the relaxation in stopped-flow kinetic measurements. Conversely, k(f) was determined under conditions favoring unfolding, for which k(u) dominates stopped-flow data. The rates determined by NMR therefore complement those from stopped-flow kinetics and define an "inverted chevron" plot. The combination of NMR relaxation and stopped-flow kinetic measurements allowed determination of k(f) and k(u) in the range from 0.48 M GuHCl to 1.28 M GuHCl. Individually, the stopped-flow and NMR data fit two-state models for folding. However, although the values of k(f) determined by the two methods agree, the values of k(u) do not. As a result, a combined analysis of all data does not comply with a two-state model but indicates that an unfolding intermediate exists on the native side of the dominant energy barrier. The denaturant and temperature dependencies of the chemical shifts and k(u) indicate that the intermediate state is structurally similar to the native state. Equilibrium unfolding monitored by optical spectroscopy corroborate these conclusions. The temperature dependence of the chemical shifts identifies regions of the protein that are selectively destabilized in the intermediate. These results illustrate the power of combining stopped-flow kinetics and NMR spectroscopy to analyze protein folding.