Microsecond Protein Dynamics Measured by 13Cα Rotating-Frame Spin Relaxation

The A39G FF domain folds on a volcano-shaped free energy surface via separate pathways

Conformational dynamics play critical roles in protein folding, misfolding, function, misfunction, and aggregation. While detecting and studying the different conformational states populated by protein molecules on their free energy surfaces (FESs) remain a challenge, NMR spectroscopy has emerged as an invaluable experimental tool to explore the FES of a protein, as conformational dynamics can be probed at atomic resolution over a wide range of timescales. Here, we use chemical exchange saturation transfer (CEST) to detect "invisible" minor states on the energy landscape of the A39G mutant FF domain that exhibited "two-state" folding kinetics in traditional experiments. Although CEST has mostly been limited to studies of processes with rates between ∼5 to 300 s^-1 involving sparse states with populations as low as ∼1%, we show that the line broadening that is often associated with minor state dips in CEST profiles can be exploited to inform on additional conformers, with lifetimes an order of magnitude shorter and populations close to 10-fold smaller than what typically is characterized. Our analysis of CEST profiles that exploits the minor state linewidths of the 71-residue A39G FF domain establishes a folding mechanism that can be described in terms of a four-state exchange process between interconverting states spanning over two orders of magnitude in timescale from ∼100 to ∼15,000 μs. A similar folding scheme is established for the wild-type domain as well. The study shows that the folding of this small domain proceeds through a pair of sparse, partially structured intermediates via two discrete pathways on a volcano-shaped FES.

Characterizing Protein Dynamics with NMR R 1ρ Relaxation Experiments

The measurement of R_1ρ , the longitudinal relaxation rate constant in the rotating frame, is one of the few available methods to characterize the μs-ms functional dynamics of biomolecules. Here, we focus on ¹⁵N R_1ρ experiments for protein NH groups. We present protocols for both on- and off-resonance ¹⁵N R_1ρ measurements needed for relaxation dispersion studies, and describe the data analysis for extracting kinetic and thermodynamic parameters characterizing the motional processes.

Conformational heterogeneity of the calmodulin binding interface

Calmodulin (CaM) is a ubiquitous Ca(2+) sensor and a crucial signalling hub in many pathways aberrantly activated in disease. However, the mechanistic basis of its ability to bind diverse signalling molecules including G-protein-coupled receptors, ion channels and kinases remains poorly understood. Here we harness the high resolution of molecular dynamics simulations and the analytical power of Markov state models to dissect the molecular underpinnings of CaM binding diversity. Our computational model indicates that in the absence of Ca(2+), sub-states in the folded ensemble of CaM's C-terminal domain present chemically and sterically distinct topologies that may facilitate conformational selection. Furthermore, we find that local unfolding is off-pathway for the exchange process relevant for peptide binding, in contrast to prior hypotheses that unfolding might account for binding diversity. Finally, our model predicts a novel binding interface that is well-populated in the Ca(2+)-bound regime and, thus, a candidate for pharmacological intervention.

Structural Insights into the Calcium-Mediated Allosteric Transition in the C-Terminal Domain of Calmodulin from Nuclear Magnetic Resonance Measurements

Calmodulin is a two-domain signaling protein that becomes activated upon binding cooperatively two pairs of calcium ions, leading to large-scale conformational changes that expose its binding site. Despite significant advances in understanding the structural biology of calmodulin functions, the mechanistic details of the conformational transition between closed and open states have remained unclear. To investigate this transition, we used a combination of molecular dynamics simulations and nuclear magnetic resonance (NMR) experiments on the Ca(2+)-saturated E140Q C-terminal domain variant. Using chemical shift restraints in replica-averaged metadynamics simulations, we obtained a high-resolution structural ensemble consisting of two conformational states and validated such an ensemble against three independent experimental data sets, namely, interproton nuclear Overhauser enhancements, (15)N order parameters, and chemical shift differences between the exchanging states. Through a detailed analysis of this structural ensemble and of the corresponding statistical weights, we characterized a calcium-mediated conformational transition whereby the coordination of Ca(2+) by just one oxygen of the bidentate ligand E140 triggers a concerted movement of the two EF-hands that exposes the target binding site. This analysis provides atomistic insights into a possible Ca(2+)-mediated activation mechanism of calmodulin that cannot be achieved from static structures alone or from ensemble NMR measurements of the transition between conformations.

Visualizing transient dark states by NMR spectroscopy

Myriad biological processes proceed through states that defy characterization by conventional atomic-resolution structural biological methods. The invisibility of these 'dark' states can arise from their transient nature, low equilibrium population, large molecular weight, and/or heterogeneity. Although they are invisible, these dark states underlie a range of processes, acting as encounter complexes between proteins and as intermediates in protein folding and aggregation. New methods have made these states accessible to high-resolution analysis by nuclear magnetic resonance (NMR) spectroscopy, as long as the dark state is in dynamic equilibrium with an NMR-visible species. These methods - paramagnetic NMR, relaxation dispersion, saturation transfer, lifetime line broadening, and hydrogen exchange - allow the exploration of otherwise invisible states in exchange with a visible species over a range of timescales, each taking advantage of some unique property of the dark state to amplify its effect on a particular NMR observable. In this review, we introduce these methods and explore two specific techniques - paramagnetic relaxation enhancement and dark state exchange saturation transfer - in greater detail.

Off-resonance rotating-frame relaxation dispersion experiment for 13C in aromatic side chains using L-optimized TROSY-selection

Protein dynamics on the microsecond-millisecond time scales often play a critical role in biological function. NMR relaxation dispersion experiments are powerful approaches for investigating biologically relevant dynamics with site-specific resolution, as shown by a growing number of publications on enzyme catalysis, protein folding, ligand binding, and allostery. To date, the majority of studies has probed the backbone amides or side-chain methyl groups, while experiments targeting other sites have been used more sparingly. Aromatic side chains are useful probes of protein dynamics, because they are over-represented in protein binding interfaces, have important catalytic roles in enzymes, and form a sizable part of the protein interior. Here we present an off-resonance R 1ρ experiment for measuring microsecond to millisecond conformational exchange of aromatic side chains in selectively (13)C labeled proteins by means of longitudinal- and transverse-relaxation optimization. Using selective excitation and inversion of the narrow component of the (13)C doublet, the experiment achieves significant sensitivity enhancement in terms of both signal intensity and the fractional contribution from exchange to transverse relaxation; additional signal enhancement is achieved by optimizing the longitudinal relaxation recovery of the covalently attached (1)H spins. We validated the L-TROSY-selected R 1ρ experiment by measuring exchange parameters for Y23 in bovine pancreatic trypsin inhibitor at a temperature of 328 K, where the ring flip is in the fast exchange regime with a mean waiting time between flips of 320 μs. The determined chemical shift difference matches perfectly with that measured from the NMR spectrum at lower temperatures, where separate peaks are observed for the two sites. We further show that potentially complicating effects of strong scalar coupling between protons (Weininger et al. in J Phys Chem B 117: 9241-9247, 2013b) can be accounted for using a simple expression, and provide recommendations for data acquisition when the studied system exhibits this behavior. The present method extends the repertoire of relaxation methods tailored for aromatic side chains by enabling studies of faster processes and improved control over artifacts due to strong coupling.

Protein conformational exchange measured by 1H R1ρ relaxation dispersion of methyl groups

Activated dynamics plays a central role in protein function, where transitions between distinct conformations often underlie the switching between active and inactive states. The characteristic time scales of these transitions typically fall in the microsecond to millisecond range, which is amenable to investigations by NMR relaxation dispersion experiments. Processes at the faster end of this range are more challenging to study, because higher RF field strengths are required to achieve refocusing of the exchanging magnetization. Here we describe a rotating-frame relaxation dispersion experiment for (1)H spins in methyl (13)CHD2 groups, which improves the characterization of fast exchange processes. The influence of (1)H-(1)H rotating-frame nuclear Overhauser effects (ROE) is shown to be negligible, based on a comparison of R 1ρ relaxation data acquired with tilt angles of 90° and 35°, in which the ROE is maximal and minimal, respectively, and on samples containing different (1)H densities surrounding the monitored methyl groups. The method was applied to ubiquitin and the apo form of calmodulin. We find that ubiquitin does not exhibit any (1)H relaxation dispersion of its methyl groups at 10 or 25 °C. By contrast, calmodulin shows significant conformational exchange of the methionine methyl groups in its C-terminal domain, as previously demonstrated by (1)H and (13)C CPMG experiments. The present R 1ρ experiment extends the relaxation dispersion profile towards higher refocusing frequencies, which improves the definition of the exchange correlation time, compared to previous results.

NMR paves the way for atomic level descriptions of sparsely populated, transiently formed biomolecular conformers

The importance of dynamics to biomolecular function is becoming increasingly clear. A description of the structure-function relationship must, therefore, include the role of motion, requiring a shift in paradigm from focus on a single static 3D picture to one where a given biomolecule is considered in terms of an ensemble of interconverting conformers, each with potentially diverse activities. In this Perspective, we describe how recent developments in solution NMR spectroscopy facilitate atomic resolution studies of sparsely populated, transiently formed biomolecular conformations that exchange with the native state. Examples of how this methodology is applied to protein folding and misfolding, ligand binding, and molecular recognition are provided as a means of illustrating both the power of the new techniques and the significant roles that conformationally excited protein states play in biology.

Isotope labeling methods for relaxation measurements

Nuclear magnetic spin relaxation has emerged as a powerful technique for probing molecular dynamics. Not only is it possible to use it for determination of time constant(s) for molecular reorientation but it can also be used to characterize internal motions on time scales from picoseconds to seconds. Traditionally, uniformly (15)N labeled samples have been used for these experiments but it is clear that this limits the applications. For instance, sensitivity for large systems is dramatically increased if dynamics is probed at methyl groups and structural characterization of low-populated states requires measurements on (13)Cα, (13)Cβ or (13)CO or (1)Hα. Unfortunately, homonuclear scalar couplings may lead to artifacts in the latter types of experiments and selective isotopic labeling schemes that only label the desired position are necessary. Both selective and uniform labeling schemes for measurements of relaxation rates for a large number of positions in proteins are discussed in this chapter.

Specific 12CβD(2)12CγD(2)S13CεHD(2) isotopomer labeling of methionine to characterize protein dynamics by 1H and 13C NMR relaxation dispersion

Protein dynamics on the micro- to millisecond time scale is increasingly found to be critical for biological function, as demonstrated by numerous NMR relaxation dispersion studies. Methyl groups are excellent probes of protein interactions and dynamics because of their favorable NMR relaxation properties, which lead to sharp signals in the ¹H and ¹³C NMR spectra. Out of the six different methyl-bearing amino acid residue types in proteins, methionine plays a special role because of its extensive side-chain flexibility and the high polarizability of the sulfur atom. Methionine is over-represented in many protein–protein recognition sites, making the methyl group of this residue type an important probe of the relationships among dynamics, interactions, and biological function. Here we present a straightforward method to label methionine residues with specific ¹³CHD₂ methyl isotopomers against a deuterated background. The resulting protein samples yield NMR spectra with improved sensitivity due to the essentially 100% population of the desired ¹³CHD₂ methyl isotopomer, which is ideal for ¹H and ¹³C spin relaxation experiments to investigate protein dynamics in general and conformational exchange in particular. We demonstrate the approach by measuring ¹H and ¹³C CPMG relaxation dispersion for the nine methionines in calcium-free calmodulin (apo-CaM). The results show that the C-terminal domain, but not the N-terminal domain, of apo-CaM undergoes fast exchange between the ground state and a high-energy state. Since target proteins are known to bind specifically to the C-terminal domain of apo-CaM, we speculate that the high-energy state might be involved in target binding through conformational selection.

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.

Conformational flexibility and the mechanisms of allosteric transitions in topologically similar proteins

Conformational flexibility plays a central role in allosteric transition of proteins. In this paper, we extend the analysis of our previous study [S. Tripathi and J. J. Portman, Proc. Natl. Acad. Sci. U.S.A. 106, 2104 (2009)] to investigate how relatively minor structural changes of the meta-stable states can significantly influence the conformational flexibility and allosteric transition mechanism. We use the allosteric transitions of the domains of calmodulin as an example system to highlight the relationship between the transition mechanism and the inter-residue contacts present in the meta-stable states. In particular, we focus on the origin of transient local unfolding (cracking), a mechanism that can lower free energy barriers of allosteric transitions, in terms of the inter-residue contacts of the meta-stable states and the pattern of local strain that develops during the transition. We find that the magnitude of the local strain in the protein is not the sole factor determining whether a region will ultimately crack during the transition. These results emphasize that the residue interactions found exclusively in one of the two meta-stable states is the key in understanding the mechanism of allosteric conformational change.

FRET-FCS detection of intralobe dynamics in calmodulin

Fluorescence correlation spectroscopy (FCS) can be coupled with Förster resonance energy transfer (FRET) to detect intramolecular dynamics of proteins on the microsecond time scale. Here we describe application of FRET-FCS to detect fluctuations within the N-terminal and C-terminal domains of the Ca(2+)-signaling protein calmodulin. Intramolecular fluctuations were resolved by global fitting of the two fluorescence autocorrelation functions (green-green and red-red) together with the two cross-correlation functions (green-red and red-green). To match the Förster radius for FRET to the dimensions of the N-terminal and C-terminal domains, a near-infrared acceptor fluorophore (Atto 740) was coupled with a green-emitting donor (Alexa Fluor 488). Fluctuations were detected in both domains on the time scale of 30 to 40 μs. In the N-terminal domain, the amplitude of the fluctuations was dependent on occupancy of Ca(2+) binding sites. A high amplitude of dynamics in apo-calmodulin (in the absence of Ca(2+)) was nearly abolished at a high Ca(2+) concentration. For the C-terminal domain, the dynamic amplitude changed little with Ca(2+) concentration. The Ca(2+) dependence of dynamics for the N-terminal domain suggests that the fluctuations detected by FCS in the N-terminal domain are coupled to the opening and closing of the EF-hand Ca(2+)-binding loops.

Toward flexibility-activity relationships by NMR spectroscopy: dynamics of Pin1 ligands

Drug design involves iterative ligand modifications. For flexible ligands, these modifications often entail restricting conformational flexibility. However, defining optimal restriction strategies can be challenging if the relationship between ligand flexibility and biological activity is unclear. Here, we describe an approach for ligand flexibility-activity studies using Nuclear Magnetic Resonance (NMR) spin relaxation. Specifically, we use (13)C relaxation dispersion measurements to compare site-specific changes in ligand flexibility for a series of related ligands that bind a common macromolecular receptor. The flexibility changes reflect conformational reorganization resulting from formation of the receptor-ligand complex. We demonstrate this approach on three structurally similar but flexibly differentiated ligands of human Pin1, a peptidyl-prolyl isomerase. The approach is able to map the ligand dynamics relevant for activity and expose changes in those dynamics caused by conformational locking. Thus, NMR flexibility-activity studies can provide information to guide strategic ligand rigidification. As such, they help establish an experimental basis for developing flexibility-activity relationships (FAR) to complement traditional structure-activity relationships (SAR) in molecular design.

Monitoring conformational dynamics with solid-state R 1rho experiments

A new application of solid-state rotating frame (R(1rho)) relaxation experiments to observe conformational dynamics is presented. Studies on a model compound, dimethyl sulfone (DMS), show that R(1rho) relaxation due to reorientation of a chemical shift anisotropy (CSA) tensor undergoing chemical exchange can be used to monitor slow-to-intermediate timescale conformational exchange processes. Control experiments used d ( 6 ) -DMS and alanine to confirm that the technique is monitoring reorientation of the CSA tensor rather than dipolar interactions or methyl group rotation. The application of this method to proteins could represent a new site-specific probe of conformational dynamics.

Differential responses of the backbone and side-chain conformational dynamics in FKBP12 upon binding the transition-state analog FK506: implications for transition-state stabilization and target protein recognition

FKBP12 serves a dual role as a peptidyl-prolyl cis-trans isomerase and as a modulator of several cell signaling pathways. The macrolide FK506 is a transition-state analog of the catalyzed reaction and displaces FKBP12 from its natural target proteins. We compared the conformational exchange dynamics of the backbone and methyl-bearing side chains of FKBP12 in the free and FK506-bound states using NMR relaxation-dispersion experiments. Our results show that the free enzyme exchanges between the ground state and an excited state that resembles the ligand-bound state or Michaelis complex. In FK506-bound FKBP12, the backbone is confined to a single conformation, while conformational exchange prevails for many methyl groups. The residual side-chain dynamics in the transition-state analog-bound state suggests that the transition-state ensemble involves multiple conformations, a finding that challenges the long-standing concept of conformational restriction in the transition-state complex. Furthermore, exchange between alternative conformations is observed in the bound state for an extended network of methyl groups that includes locations remote from the active site. Several of these locations are known to be important for interactions with cellular target proteins, including calcineurin and the ryanodine receptor, suggesting that the conformational heterogeneity might play a role in the promiscuous binding of FKBP12 to different targets.

Dynamics of ligand binding from 13C NMR relaxation dispersion at natural abundance

We show that Carr-Purcell-Meiboom-Gill (CPMG) 13Calpha NMR relaxation dispersion measurements are a viable means for profiling mus-ms ligand dynamics involved in receptor binding. Critically, the dispersion is at natural 13C abundance; this matches typical pharmaceutical research settings in which ligand isotope-labeling is often impractical. The dispersion reveals ligand 13Calpha nuclei that experience mus-ms modulation of their chemical shifts due to binding. 13Calpha shifts are dominated by local torsion angles , psi, chi1; hence, these experiments identify flexible torsion angles that may assist complex formation. Since the experiments detect the ligand, they are viable even in the absence of a receptor structure. The mus-ms dynamic information gained helps establish flexibility-activity relationships. We apply these experiments to study the binding of a phospho-peptide substrate ligand to the peptidyl-prolyl isomerase Pin1.

CEESY: characterizing the conformation of unobservable protein states

Protein conformations that are only marginally populated often play important roles as intermediate states in many processes such as ligand binding, enzyme catalysis, allostery, and protein folding. An NMR method is presented that can give valuable information about the structure of these "excited states" by measuring the relative position of exchanging excited- and ground-state resonances using a single 2D spectrum. This new approach can be applied to any nucleus, which will facilitate a complete structural characterization of these states.

Biosynthetic 13C labeling of aromatic side chains in proteins for NMR relaxation measurements

Site-specific 13C labeling offers a desirable means of eliminating unwanted relaxation pathways and coherent magnetization transfer in NMR relaxation experiments. Here we use [1-13C]-glucose as the sole carbon source in the growth media for protein overexpression in Escherichia coli. The approach results in specific incorporation of 13C at isolated positions in the side chains of aromatic amino acids, which greatly simplifies the measurements and interpretation of 13C relaxation rates in these spin systems. The method is well suited for characterization of chemical exchange by CPMG or spin-lock relaxation methods. We validated the method by acquiring 13C rotating-frame relaxation dispersion data on the E140Q mutant of the C-terminal domain of calmodulin, which reveal conformational exchange dynamics with a time constant of 71 mus for Y138.

Correlated dynamics of consecutive residues reveal transient and cooperative unfolding of secondary structure in proteins

Nuclear spin relaxation is a powerful method for studying molecular dynamics at atomic resolution. Recent methods development in biomolecular NMR spectroscopy has enabled detailed investigations of molecular dynamics that are critical for biological function, with prominent examples addressing allostery, enzyme catalysis, and protein folding. Dynamic processes with similar correlation times are often detected in multiple locations of the molecule, raising the question of whether the underlying motions are correlated (corresponding to concerted fluctuations involving many atoms distributed across extended regions of the molecule) or uncorrelated (corresponding to independent fluctuations involving few atoms in localized regions). Here, we have used (13)C(alpha)(i - 1)/(13)C(alpha)(i) differential multiple-quantum spin relaxation to provide direct evidence for correlated dynamics of consecutive amino acid residues in the protein sequence. By monitoring overlapping pairs of residues (i - 1 and i, i and i + 1, etc.), we identified correlated motions that extend through continuous segments of the sequence. We detected significant correlated conformational transitions in the native state of the E140Q mutant of the calmodulin C-terminal domain. Previous work has shown that this domain exchanges between two major conformational states that resemble the functionally relevant open and closed states of the WT protein, with a mean correlation time of approximately 20 micros. The present results reveal that an entire alpha-helix undergoes partial unraveling in a transient and cooperative manner.