Measuring 1HN temperature coefficients in invisible protein states by relaxation dispersion NMR spectroscopy

Non-uniform sampling of similar NMR spectra and its application to studies of the interaction between alpha-synuclein and liposomes

The accelerated acquisition of multidimensional NMR spectra using sparse non-uniform sampling (NUS) has been widely adopted in recent years. The key concept in NUS is that a major part of the data is omitted during measurement, and then reconstructed using, for example, compressed sensing (CS) methods. CS requires spectra to be compressible, that is, they should contain relatively few "significant" points. The more compressible the spectrum, the fewer experimental NUS points needed in order for it to be accurately reconstructed. In this paper we show that the CS processing of similar spectra can be enhanced by reconstructing only the differences between them. Accurate reconstruction can be obtained at lower sampling levels as the difference is sparser than the spectrum itself. In many situations this method is superior to "conventional" compressed sensing. We exemplify the concept of "difference CS" with one such case-the study of alpha-synuclein binding to liposomes and its dependence on temperature. To obtain information on temperature-dependent transitions between different states, we need to acquire several dozen spectra at various temperatures, with and without the presence of liposomes. Our detailed investigation reveals that changes in the binding modes of the alpha-synuclein ensemble are not only temperature-dependent but also show non-linear behavior in their transitions. Our proposed CS processing approach dramatically reduces the number of NUS points required and thus significantly shortens the experimental time.

A fine balance of hydrophobic-electrostatic communication pathways in a pH-switching protein

Allostery is the phenomenon of coupling between distal binding sites in a protein. Such coupling is at the crux of protein function and regulation in a myriad of scenarios, yet determining the molecular mechanisms of coupling networks in proteins remains a major challenge. Here, we report mechanisms governing pH-dependent myristoyl switching in monomeric hisactophilin, whereby the myristoyl moves between a sequestered state, i.e., buried within the core of the protein, to an accessible state, in which the myristoyl has increased accessibility for membrane binding. Measurements of the pH and temperature dependence of amide chemical shifts reveal protein local structural stability and conformational heterogeneity that accompany switching. An analysis of these measurements using a thermodynamic cycle framework shows that myristoyl-proton coupling at the single-residue level exists in a fine balance and extends throughout the protein. Strikingly, small changes in the stereochemistry or size of core and surface hydrophobic residues by point mutations readily break, restore, or tune myristoyl switch energetics. Synthesizing the experimental results with those of molecular dynamics simulations illuminates atomistic details of coupling throughout the protein, featuring a large network of hydrophobic interactions that work in concert with key electrostatic interactions. The simulations were critical for discerning which of the many ionizable residues in hisactophilin are important for switching and identifying the contributions of nonnative interactions in switching. The strategy of using temperature-dependent NMR presented here offers a powerful, widely applicable way to elucidate the molecular mechanisms of allostery in proteins at high resolution.

Temperature dependence of NMR chemical shifts: Tracking and statistical analysis

Isotropic chemical shifts measured by solution nuclear magnetic resonance (NMR) spectroscopy offer extensive insights into protein structure and dynamics. Temperature dependences add a valuable dimension; notably, the temperature dependences of amide proton chemical shifts are valuable probes of hydrogen bonding, temperature‐dependent loss of structure, and exchange between distinct protein conformations. Accordingly, their uses include structural analysis of both folded and disordered proteins, and determination of the effects of mutations, binding, or solution conditions on protein energetics. Fundamentally, these temperature dependences result from changes in the local magnetic environments of nuclei, but correlations with global thermodynamic parameters measured via calorimetric methods have been observed. Although the temperature dependences of amide proton and nitrogen chemical shifts are often well approximated by a linear model, deviations from linearity are also observed and may be interpreted as evidence of fast exchange between distinct conformational states. Here, we describe computational methods, accessible via the Shift‐T web server, including an automated tracking algorithm that propagates initial (single temperature) ¹H—¹⁵N cross peak assignments to spectra collected over a range of temperatures. Amide proton and nitrogen temperature coefficients (slopes determined by fitting chemical shift vs. temperature data to a linear model) are subsequently calculated. Also included are methods for the detection of systematic, statistically significant deviation from linearity (curvature) in the temperature dependences of amide proton chemical shifts. The use and utility of these methods are illustrated by example, and the Shift‐T web server is freely available at http://meieringlab.uwaterloo.ca/shiftt.

Multistep mutational transformation of a protein fold through structural intermediates

New protein folds may evolve from existing folds through metamorphic evolution involving a dramatic switch in structure. To mimic pathways by which amino acid sequence changes could induce a change in fold, we designed two folded hybrids of Xfaso 1 and Pfl 6, a pair of homologous Cro protein sequences with ~40% identity but different folds (all-α vs. α + β, respectively). Each hybrid, XPH1 or XPH2, is 85% identical in sequence to its parent, Xfaso 1 or Pfl 6, respectively; 55% identical to its noncognate parent; and ~70% identical to the other hybrid. XPH1 and XPH2 also feature a designed hybrid chameleon sequence corresponding to the C-terminal region, which switched from α-helical to β-sheet structure during Cro evolution. We report solution nuclear magnetic resonance (NMR) structures of XPH1 and XPH2 at 0.3 Å and 0.5 Å backbone root mean square deviation (RMSD), respectively. XPH1 retains a global fold generally similar to Xfaso 1, and XPH2 retains a fold similar to Pfl 6, as measured by TM-align scores (~0.7), DALI Z-scores (7-9), and backbone RMSD (2-3 Å RMSD for the most ordered regions). However, these scores also indicate significant deviations in structure. Most notably, XPH1 and XPH2 have different, and intermediate, secondary structure content relative to Xfaso 1 and Pfl 6. The multistep progression in sequence, from Xfaso 1 to XPH1 to XPH2 to Pfl 6, thus involves both abrupt and gradual changes in folding pattern. The plasticity of some protein folds may allow for "polymetamorphic" evolution through intermediate structures.

Using 1 HN amide temperature coefficients to define intrinsically disordered regions: An alternative NMR method

This report describes a cost-effective experimental method for determining an intrinsically disordered protein (IDP) region in a given protein sample. In this area, the most popular (and conventional) means is using the amide (¹ H^N ) NMR signal chemical shift distributed in the range of 7.5-8.5 ppm. For this study, we applied an additional step: analysis of ¹ H^N chemical shift temperature coefficients (¹ H^N -CSTCs) of the signals. We measured ¹ H-¹⁵ N two-dimensional NMR spectra of model IDP samples and ordered samples at four temperatures (288, 293, 298, and 303 K). We derived the ¹ H^N -CSTC threshold deviation, which gives the best correlation of ordered and disordered regions among the proteins examined (below -3.6 ppb/K). By combining these criteria with the newly optimized chemical shift range (7.8-8.5 ppm), the ratios of both true positive and true negative were improved by approximately 19% (62-81%) compared with the conventional "chemical shift-only" method.

Concurrent Increases and Decreases in Local Stability and Conformational Heterogeneity in Cu, Zn Superoxide Dismutase Variants Revealed by Temperature-Dependence of Amide Chemical Shifts

The chemical shifts of backbone amide protons in proteins are sensitive reporters of local structural stability and conformational heterogeneity, which can be determined from their readily measured linear and nonlinear temperature-dependences, respectively. Here we report analyses of amide proton temperature-dependences for native dimeric Cu, Zn superoxide dismutase (holo pWT SOD1) and structurally diverse mutant SOD1s associated with amyotrophic lateral sclerosis (ALS). Holo pWT SOD1 loses structure with temperature first at its periphery and, while having extremely high global stability, nevertheless exhibits extensive conformational heterogeneity, with ∼1 in 5 residues showing evidence for population of low energy alternative states. The holo G93A and E100G ALS mutants have moderately decreased global stability, whereas V148I is slightly stabilized. Comparison of the holo mutants as well as the marginally stable immature monomeric unmetalated and disulfide-reduced (apo(2SH)) pWT with holo pWT shows that changes in the local structural stability of individual amides vary greatly, with average changes corresponding to differences in global protein stability measured by differential scanning calorimetry. Mutants also exhibit altered conformational heterogeneity compared to pWT. Strikingly, substantial increases as well as decreases in local stability and conformational heterogeneity occur, in particular upon maturation and for G93A. Thus, the temperature-dependence of amide shifts for SOD1 variants is a rich source of information on the location and extent of perturbation of structure upon covalent changes and ligand binding. The implications for potential mechanisms of toxic misfolding of SOD1 in disease and for general aspects of protein energetics, including entropy-enthalpy compensation, are discussed.

Mechanism of Protein Denaturation: Partial Unfolding of the P22 Coat Protein I-Domain by Urea Binding

The I-domain is an insertion domain of the bacteriophage P22 coat protein that drives rapid folding and accounts for over half of the stability of the full-length protein. We sought to determine the role of hydrogen bonds (H-bonds) in the unfolding of the I-domain by examining (3)JNC' couplings transmitted through H-bonds, the temperature and urea-concentration dependence of (1)HN and (15)N chemical shifts, and native-state hydrogen exchange at urea concentrations where the domain is predominantly folded. The native-state hydrogen-exchange data suggest that the six-stranded β-barrel core of the I-domain is more stable against unfolding than a smaller subdomain comprised of a short α-helix and three-stranded β-sheet. H-bonds, separately determined from solvent protection and (3)JNC' H-bond couplings, are identified with an accuracy of 90% by (1)HN temperature coefficients. The accuracy is improved to 95% when (15)N temperature coefficients are also included. In contrast, the urea dependence of (1)HN and (15)N chemical shifts is unrelated to H-bonding. The protein segments with the largest chemical-shift changes in the presence of urea show curved or sigmoidal titration curves suggestive of direct urea binding. Nuclear Overhauser effects to urea for these segments are also consistent with specific urea-binding sites in the I-domain. Taken together, the results support a mechanism of urea unfolding in which denaturant binds to distinct sites in the I-domain. Disordered segments bind urea more readily than regions in stable secondary structure. The locations of the putative urea-binding sites correlate with the lower stability of the structure against solvent exchange, suggesting that partial unfolding of the structure is related to urea accessibility.

A role for indels in the evolution of Cro protein folds

Insertions and deletions in protein sequences, or indels, can disrupt structure and may result in changes in protein folds during evolution or in association with alternative splicing. Pfl 6 and Xfaso 1 are two proteins in the Cro family that share a common ancestor but have different folds. Sequence alignments of the two proteins show two gaps, one at the N terminus, where the sequence of Xfaso 1 is two residues shorter, and one near the center of the sequence, where the sequence of Pfl 6 is five residues shorter. To test the potential importance of indels in Cro protein evolution, we generated hybrid variants of Pfl 6 and Xfaso 1 with indels in one or both regions, chosen according to several plausible sequence alignments. All but one deletion variant completely unfolded both proteins, showing that a longer N-terminal sequence was critical for Pfl 6 folding and a longer central region sequence was critical for Xfaso 1 folding. By contrast, Xfaso 1 tolerated a longer N-terminal sequence with little destabilization, and Pfl 6 tolerated central region insertions, albeit with substantial effects on thermal stability and some perturbation of the surrounding structure. None of the mutations appeared to convert one stable fold into the other. On the basis of this two-protein comparison, short insertion and deletion mutations probably played a role in evolutionary fold change in the Cro family, but were also not the only factors.

Increasing the exchange time-scale that can be probed by CPMG relaxation dispersion NMR

Carr-Purcell-Meiboom-Gill relaxation dispersion NMR spectroscopy has emerged as a valuable tool to characterize conformational exchange between major and minor states in a large variety of biomolecules. The window of exchange that is amenable for study, corresponding to rates on the order of 2000 s(-1) or less, is limiting, however. Here we show that a combined analysis of both amide (15)N and (1)H(N) CPMG profiles and major state exchange induced (15)N chemical shift changes leads to significant increases in the exchange time scale for which accurate exchange parameters and chemical shift differences between the interconverting states can be obtained. The utility of the approach is illustrated with examples involving a pair of protein systems that are in the moderately fast exchange regime. In these cases the analysis of dispersion profiles alone is not sufficient to obtain robust measures of exchange parameters and chemical shift differences. Inclusion of major state exchange induced (15)N chemical shift changes measured in ((15)N-(1)H(N)) HMQC and HSQC data sets in addition to the (15)N and (1)H(N) dispersion profiles in the analysis "breaks" the correlation in parameters, allowing accurate values to be obtained. The approach is straightforward to implement and makes use of HMQC/HSQC data sets that are recorded as a matter of routine to obtain chemical shifts of the excited state. It promises to increase the range of exchanging systems involving low populated, transiently formed excited states that can be studied by relaxation dispersion NMR.