Transient, sparsely populated compact states of apo and calcium-loaded calmodulin probed by paramagnetic relaxation enhancement: interplay of conformational selection and induced fit

Ligand-induced protein transition state stabilization switches the binding pathway from conformational selection to induced fit

Protein-ligand complex formation is fundamental to biological function. A central question is whether proteins spontaneously adopt binding-competent conformations to which ligands bind conformational selection (CS) or whether ligands induce the binding-competent conformation induced fit (IF). Here, we resolve the CS and IF binding pathways by characterizing protein conformational dynamics over a wide range of ligand concentrations using NMR relaxation dispersion. We determined the relative flux through the two pathways using a four-state binding model that includes both CS and IF. Experiments conducted without ligand show that galectin-3 exchanges between the ground-state conformation and a high-energy conformation similar to the ligand-bound conformation, demonstrating that CS is a plausible pathway. Near-identical crystal structures of the apo and ligand-bound states suggest that the high-energy conformation in solution corresponds to the apo crystal structure. Stepwise additions of the ligand lactose induce progressive changes in the relaxation dispersions that we fit collectively to the four-state model, yielding all microscopic rate constants and binding affinities. The ligand affinity is higher for the bound-like conformation than for the ground state, as expected for CS. Nonetheless, the IF pathway contributes greater than 70% of the total flux even at low ligand concentrations. The higher flux through the IF pathway is explained by considerably higher rates of exchange between the two protein conformations in the ligand-associated state. Thus, the ligand acts to decrease the activation barrier between protein conformations in a manner reciprocal to enzymatic transition-state stabilization of reactions involving ligand transformation.

Visualizing Heterogeneous Protein Conformations with Multi-Tilt Nanoparticle-Aided Cryo-Electron Microscopy Sampling

Obtaining the heterogeneous conformation of small proteins is important for understanding their biological role, but it is still challenging. Here, we developed a multi-tilt nanoparticle-aided cryo-electron microscopy sampling (MT-NACS) technique that enables the observation of heterogeneous conformations of small proteins and applied it to calmodulin. By imaging the proteins labeled by two gold nanoparticles at multiple tilt angles and analyzing the projected positions of the nanoparticles, the distributions of 3D interparticle distances were obtained. From the measured distance distributions, the conformational changes associated with Ca²⁺ binding and salt concentration were determined. MT-NACS was also used to track the structural change accompanied by the interaction between amyloid-beta and calmodulin, which has never been observed experimentally. This work offers an alternative platform for studying the functional flexibility of small proteins.

Paramagnetic NMR restraints for the characterization of protein structural rearrangements

Mobility is a common feature of biomacromolecules, often fundamental for their function. Thus, in many cases, biomacromolecules cannot be described by a single conformation, but rather by a conformational ensemble. NMR paramagnetic data demonstrated quite informative to monitor this conformational variability, especially when used in conjunction with data from different sources. Due to their long-range nature, paramagnetic data can, for instance, i) clearly demonstrate the occurrence of conformational rearrangements, ii) reveal the presence of minor conformational states, sampled only for a short time, iii) indicate the most representative conformations within the conformational ensemble sampled by the molecule, iv) provide an upper limit to the weight of each conformation.

Single-Molecule Electrical and Spectroscopic Profiling Protein Allostery Using a Gold Plasmonic Nanopore

Direct structural and dynamic characterization of protein conformers in solution is highly desirable but currently impractical. Herein, we developed a single molecule gold plasmonic nanopore system for observation of protein allostery, enabling us to monitor translocation dynamics and conformation transition of proteins by ion current detection and SERS spectrum measurement, respectively. Allosteric transition of calmodulin (CaM) was elaborately probed by the nanopore system. Two conformers of CaM were well-resolved at a single-molecule level using both the ion current blockage signal and the SERS spectra. The collected SERS spectra provided structural evidence to confirm the interaction between CaM and the gold plasmonic nanopore, which was responsible for the different translocation behaviors of the two conformers. SERS spectra revealed the amino acid residues involved in the conformational change of CaM upon calcium binding. The results demonstrated that the excellent spectral characterization furnishes a single-molecule nanopore technique with an advanced capability of direct structure analysis.

Molecular dynamics simulations of the calmodulin-induced α-helix in the SK2 calcium-gated potassium ion channel

The family of small-conductance Ca²⁺-activated potassium ion channels (SK channels) is composed of four members (SK1, SK2, SK3, and SK4) involved in neuron-firing regulation. The gating of these channels depends on the intracellular Ca²⁺ concentration, and their sensitivity to this ion is provided by calmodulin (CaM). This protein binds to a specific region in SK channels known as the calmodulin-binding domain (CaMBD), an event which is essential for their gating. While CaMBDs are typically disordered in the absence of CaM, the SK2 channel subtype displays a small prefolded α-helical region in its CaMBD even if CaM is not present. This small helix is known to turn into a full α-helix upon CaM binding, although the molecular-level details for this conversion are not fully understood yet. In this work, we offer new insights on this physiologically relevant process by means of enhanced sampling, atomistic Hamiltonian replica exchange molecular dynamics simulations, providing a more detailed understanding of CaM binding to this target. Our results show that CaM is necessary for inducing a full α-helix along the SK2 CaMBD through hydrophobic interactions with V426 and L427. However, it is also necessary that W431 does not compete for these interactions; the role of the small prefolded α-helix in the SK2 CaMBD would be to stabilize W431 so that this is the case. In conclusion, our findings provide further insight into a key interaction between CaM and SK channels that is important for channel sensitivity to Ca²⁺.

Binding and Functional Folding (BFF): A Physiological Framework for Studying Biomolecular Interactions and Allostery

EF-hand Ca2+-binding proteins (CBPs), such as S100 proteins (S100s) and calmodulin (CaM), are signaling proteins that undergo conformational changes upon increasing intracellular Ca2+. Upon binding Ca2+, S100 proteins and CaM interact with protein targets and induce important biological responses. The Ca2+-binding affinity of CaM and most S100s in the absence of target is weak (CaKD > 1 μM). However, upon effector protein binding, the Ca2+ affinity of these proteins increases via heterotropic allostery (CaKD < 1 μM). Because of the high number and micromolar concentrations of EF-hand CBPs in a cell, at any given time, allostery is required physiologically, allowing for (i) proper Ca2+ homeostasis and (ii) strict maintenance of Ca2+-signaling within a narrow dynamic range of free Ca2+ ion concentrations, [Ca2+]free. In this review, mechanisms of allostery are coalesced into an empirical "binding and functional folding (BFF)" physiological framework. At the molecular level, folding (F), binding and folding (BF), and BFF events include all atoms in the biomolecular complex under study. The BFF framework is introduced with two straightforward BFF types for proteins (type 1, concerted; type 2, stepwise) and considers how homologous and nonhomologous amino acid residues of CBPs and their effector protein(s) evolved to provide allosteric tightening of Ca2+ and simultaneously determine how specific and relatively promiscuous CBP-target complexes form as both are needed for proper cellular function.

Time-resolved solid state NMR of biomolecular processes with millisecond time resolution

We review recent efforts to develop and apply an experimental approach to the structural characterization of transient intermediate states in biomolecular processes that involve large changes in molecular conformation or assembly state. This approach depends on solid state nuclear magnetic resonance (ssNMR) measurements that are performed at very low temperatures, typically 25-30 K, with signal enhancements from dynamic nuclear polarization (DNP). This approach also involves novel technology for initiating the process of interest, either by rapid mixing of two solutions or by a rapid inverse temperature jump, and for rapid freezing to trap intermediate states. Initiation by rapid mixing or an inverse temperature jump can be accomplished in approximately-one millisecond. Freezing can be accomplished in approximately 100 microseconds. Thus, millisecond time resolution can be achieved. Recent applications to the process by which the biologically essential calcium sensor protein calmodulin forms a complex with one of its target proteins and the process by which the bee venom peptide melittin converts from an unstructured monomeric state to a helical, tetrameric state after a rapid change in pH or temperature are described briefly. Future applications of millisecond time-resolved ssNMR are also discussed briefly.

The evolution of paramagnetic NMR as a tool in structural biology

Paramagnetic NMR data contain extremely accurate long-range information on metalloprotein structures and, when used in the frame of integrative structural biology approaches, they allow for the retrieval of structural details to a resolution that is not achievable using other techniques. Paramagnetic data thus represent an extremely powerful tool to refine protein models in solution, especially when coupled to X-ray or cryoelectron microscopy data, to monitor the formation of complexes and determine the relative arrangements of their components, and to highlight the presence of conformational heterogeneity. More recently, theoretical and computational advancements in quantum chemical calculations of paramagnetic NMR observables are progressively opening new routes in structural biology, because they allow for the determination of the structure within the coordination sphere of the metal center, thus acting as a loupe on sites that are difficult to observe but very important for protein function.

Quantitative Agreement between Conformational Substates of Holo Calcium-Loaded Calmodulin Detected by Double Electron-Electron Resonance EPR and Predicted by Molecular Dynamics Simulations

Calcium-loaded calmodulin (CaM/4Ca²⁺) comprises two domains that undergo rigid body reorientation from a predominantly extended conformation to a compact one upon binding target peptides. A recent replica-exchange molecular dynamics (MD) simulation on holo CaM/4Ca²⁺ suggested the existence of distinct structural clusters (substates) along the path from extended to compact conformers in the absence of substrates. Here, we experimentally demonstrate the existence of CaM/4Ca²⁺ substates trapped in local minima by three freezing/annealing regimes (slow, 40 s; intermediate, 1.5 s; fast, 0.5 ms) using pulsed Q-band double electron-electron resonance (DEER) EPR spectroscopy to measure interdomain distances between nitroxide spin-labels positioned at A17C and A128C in the N- and C-terminal domains, respectively. The DEER echo curves were directly fit to population-optimized P(r) pairwise distance distributions calculated from the coordinates of the MD clusters and compact crystal structure. DEER data on fully deuterated CaM/4Ca²⁺ were acquired at multiple values of the second echo period (10-35 μs) and analyzed globally to eliminate instrumental and overfitting artifacts and ensure accurate populations, peak positions, and widths. The DEER data for all three freezing regimes are quantitatively accounted for within experimental error by 5-6 distinct conformers comprising a predominantly populated extended form (60-75%) and progressively more compact states whose populations decrease as the degree of compactness increases. The shortest interdomain separation is found in the compact crystal structure, which has an occupancy of 4-6%. Thus, CaM/4Ca²⁺ samples high energy local minima comprising a few discrete substates of increasing compactness in a rugged energy landscape.

Paramagnetic Chemical Probes for Studying Biological Macromolecules

Paramagnetic chemical probes have been used in electron paramagnetic resonance (EPR) and nuclear magnetic resonance (NMR) spectroscopy for more than four decades. Recent years witnessed a great increase in the variety of probes for the study of biological macromolecules (proteins, nucleic acids, and oligosaccharides). This Review aims to provide a comprehensive overview of the existing paramagnetic chemical probes, including chemical synthetic approaches, functional properties, and selected applications. Recent developments have seen, in particular, a rapid expansion of the range of lanthanoid probes with anisotropic magnetic susceptibilities for the generation of structural restraints based on residual dipolar couplings and pseudocontact shifts in solution and solid state NMR spectroscopy, mostly for protein studies. Also many new isotropic paramagnetic probes, suitable for NMR measurements of paramagnetic relaxation enhancements, as well as EPR spectroscopic studies (in particular double resonance techniques) have been developed and employed to investigate biological macromolecules. Notwithstanding the large number of reported probes, only few have found broad application and further development of probes for dedicated applications is foreseen.

DEER experiments reveal fundamental differences between calmodulin complexes with IQ and MARCKS peptides in solution

Calmodulin (CaM) is a calcium-binding protein that regulates the function of many proteins by indirectly conferring Ca²⁺ sensitivity, and it undergoes a large conformational change on partners' binding. We compared the solution binding mode of the target peptides MARCKS and IQ by double electron-electron resonance (DEER) distance measurements and paramagnetic NMR. We combined nitroxide and Gd(III) spin labels, including specific substitution of one of the Ca²⁺ ions in the CaM mutant N60D by a Gd(III) ion. The binding of MARCKS to holo-CaM resulted neither in a closed conformation nor in a unique relative orientation between the two CaM domains, in contrast with the crystal structure. Binding of IQ to holo-CaM did generate a closed conformation. Using elastic network modeling and 12 distance restraints obtained from multiple holo-CaM/IQ DEER data, we derived a model of the solution structure, which is in reasonable agreement with the crystal structure.

Time-resolved DEER EPR and solid-state NMR afford kinetic and structural elucidation of substrate binding to Ca2+-ligated calmodulin

Complex formation between calmodulin and target proteins underlies numerous calcium signaling processes in biology, yet structural and mechanistic details, which entail major conformational changes in both calmodulin and its substrates, have been unclear. We show that a combination of time-resolved electron paramagnetic and NMR measurements can elucidate the molecular mechanism, at the quantitative kinetic and structural levels, of the binding pathway of a peptide substrate from skeletal muscle myosin light-chain kinase to calcium-loaded calmodulin. The mechanism involves coupled folding and binding and comprises a bifurcated process, with rapid, direct complex formation when the peptide interacts first with the C-terminal domain of calmodulin or a slower, two-step complex formation when the peptide interacts initially with the N-terminal domain.

Recent advances in rapid mixing and freeze quenching have opened the path for time-resolved electron paramagnetic resonance (EPR)-based double electron-electron resonance (DEER) and solid-state NMR of protein–substrate interactions. DEER, in conjunction with phase memory time filtering to quantitatively extract species populations, permits monitoring time-dependent probability distance distributions between pairs of spin labels, while solid-state NMR provides quantitative residue-specific information on the appearance of structural order and the development of intermolecular contacts between substrate and protein. Here, we demonstrate the power of these combined approaches to unravel the kinetic and structural pathways in the binding of the intrinsically disordered peptide substrate (M13) derived from myosin light-chain kinase to the universal eukaryotic calcium regulator, calmodulin. Global kinetic analysis of the data reveals coupled folding and binding of the peptide associated with large spatial rearrangements of the two domains of calmodulin. The initial binding events involve a bifurcating pathway in which the M13 peptide associates via either its N- or C-terminal regions with the C- or N-terminal domains, respectively, of calmodulin/4Ca²⁺ to yield two extended “encounter” complexes, states A and A*, without conformational ordering of M13. State A is immediately converted to the final compact complex, state C, on a timescale τ ≤ 600 μs. State A*, however, only reaches the final complex via a collapsed intermediate B (τ ∼ 1.5 to 2.5 ms), in which the peptide is only partially ordered and not all intermolecular contacts are formed. State B then undergoes a relatively slow (τ ∼ 7 to 18 ms) conformational rearrangement to state C.

Characterizing conformational ensembles of multi-domain proteins using anisotropic paramagnetic NMR restraints

It has been over two decades since paramagnetic NMR started to form part of the essential techniques for structural analysis of proteins under physiological conditions. Paramagnetic NMR has significantly expanded our understanding of the inherent flexibility of proteins, in particular, those that are formed by combinations of two or more domains. Here, we present a brief overview of techniques to characterize conformational ensembles of such multi-domain proteins using paramagnetic NMR restraints produced through anisotropic metals, with a focus on the basics of anisotropic paramagnetic effects, the general procedures of conformational ensemble reconstruction, and some representative reweighting approaches.

The mechanism of complex formation between calmodulin and voltage gated calcium channels revealed by molecular dynamics

Calmodulin, a ubiquitous eukaryotic calcium sensor responsible for the regulation of many fundamental cellular processes, is a highly flexible protein and exhibits an unusually wide range of conformations. Furthermore, CaM is known to interact with more than 300 cellular targets. Molecular dynamics (MD) simulation trajectories suggest that EF-hand loops show different magnitudes of flexibility. Therefore, the four EF-hand motifs have different affinities for Ca2+ ions, which enables CaM to function on wide range of Ca2+ ion concentrations. EF-hand loops are 2-3 times more flexible in apo CaM whereas least flexible in Ca2+/CaM-IQ motif complexes. We report a unique intermediate conformation of Ca2+/CaM while transitioning from extended to compact form. We also report the complex formation process between Ca2+/CaM and IQ CaM-binding motifs. Our results showed how IQ motif recognise its binding site on the CaM and how CaM transforms from extended to compact form upon binding to IQ motif.

Reconstruction of Coupled Intra- and Interdomain Protein Motion from Nuclear and Electron Magnetic Resonance

Proteins composed of multiple domains allow for structural heterogeneity and interdomain dynamics that may be vital for function. Intradomain structures and dynamics can influence interdomain conformations and vice versa. However, no established structure determination method is currently available that can probe the coupling of these motions. The protein Pin1 contains separate regulatory and catalytic domains that sample "extended" and "compact" states, and ligand binding changes this equilibrium. Ligand binding and interdomain distance have been shown to impact the activity of Pin1, suggesting interdomain allostery. In order to characterize the conformational equilibrium of Pin1, we describe a novel method to model the coupling between intra- and interdomain dynamics at atomic resolution using multistate ensembles. The method uses time-averaged nuclear magnetic resonance (NMR) restraints and double electron-electron resonance (DEER) data that resolve distance distributions. While the intradomain calculation is primarily driven by exact nuclear Overhauser enhancements (eNOEs), J couplings, and residual dipolar couplings (RDCs), the relative domain distribution is driven by paramagnetic relaxation enhancement (PREs), RDCs, interdomain NOEs, and DEER. Our data support a 70:30 population of the compact and extended states in apo Pin1. A multistate ensemble describes these conformations simultaneously, with distinct conformational differences located in the interdomain interface stabilizing the compact or extended states. We also describe correlated conformations between the catalytic site and interdomain interface that may explain allostery driven by interdomain contact.

Paramagnetic Relaxation Enhancement for Detecting and Characterizing Self-Associations of Intrinsically Disordered Proteins

Intrinsically disordered proteins and intrinsically disordered regions within proteins make up a large and functionally significant part of the human proteome. The highly flexible nature of these sequences allows them to form weak, long-range, and transient interactions with diverse biomolecular partners. Specific yet low-affinity interactions promote promiscuous binding and enable a single intrinsically disordered segment to interact with a multitude of target sites. Because of the transient nature of these interactions, they can be difficult to characterize by structural biology methods that rely on proteins to form a single, predominant conformation. Paramagnetic relaxation enhancement NMR is a useful tool for identifying and defining the structural underpinning of weak and transient interactions. A detailed protocol for using paramagnetic relaxation enhancement to characterize the lowly-populated encounter complexes that form between intrinsically disordered proteins and their protein, nucleic acid, or other biomolecular partners is described.

Coarse-Grained Modeling and Molecular Dynamics Simulations of Ca2+-Calmodulin

Calmodulin (CaM) is a calcium-binding protein that transduces signals to downstream proteins through target binding upon calcium binding in a time-dependent manner. Understanding the target binding process that tunes CaM’s affinity for the calcium ions (Ca²⁺), or vice versa, may provide insight into how Ca²⁺-CaM selects its target binding proteins. However, modeling of Ca²⁺-CaM in molecular simulations is challenging because of the gross structural changes in its central linker regions while the two lobes are relatively rigid due to tight binding of the Ca²⁺ to the calcium-binding loops where the loop forms a pentagonal bipyramidal coordination geometry with Ca²⁺. This feature that underlies the reciprocal relation between Ca²⁺ binding and target binding of CaM, however, has yet to be considered in the structural modeling. Here, we presented a coarse-grained model based on the Associative memory, Water mediated, Structure, and Energy Model (AWSEM) protein force field, to investigate the salient features of CaM. Particularly, we optimized the force field of CaM and that of Ca²⁺ ions by using its coordination chemistry in the calcium-binding loops to match with experimental observations. We presented a “community model” of CaM that is capable of sampling various conformations of CaM, incorporating various calcium-binding states, and carrying the memory of binding with various targets, which sets the foundation of the reciprocal relation of target binding and Ca²⁺ binding in future studies.

An Integrative Approach to Determine 3D Protein Structures Using Sparse Paramagnetic NMR Data and Physical Modeling

Paramagnetic nuclear magnetic resonance (NMR) methods have emerged as powerful tools for structure determination of large, sparsely protonated proteins. However traditional applications face several challenges, including a need for large datasets to offset the sparsity of restraints, the difficulty in accounting for the conformational heterogeneity of the spin-label, and noisy experimental data. Here we propose an integrative approach to structure determination combining sparse paramagnetic NMR with physical modelling to infer approximate protein structural ensembles. We use calmodulin in complex with the smooth muscle myosin light chain kinase peptide as a model system. Despite acquiring data from samples labeled only at the backbone amide positions, we are able to produce an ensemble with an average RMSD of ∼2.8 Å from a reference X-ray crystal structure. Our approach requires only backbone chemical shifts and measurements of the paramagnetic relaxation enhancement and residual dipolar couplings that can be obtained from sparsely labeled samples.

Unveiling protein dynamics in solution with field-cycling NMR relaxometry

Field-cycling NMR relaxometry is a well-established technique that can give information on molecular structure and dynamics of biological systems. It provides the nuclear relaxation rates as a function of the applied magnetic field, starting from fields as low as ~ 10^-4 T up to about 1-3 T. The profiles so collected, called nuclear magnetic relaxation dispersion (NMRD) profiles, can be extended to include the relaxation rates at the largest fields achievable with high resolution NMR spectrometers. By exploiting this wide range of frequencies, the NMRD profiles can provide information on motions occurring on time scales from 10^-6 to 10^-9 s. ¹H NMRD measurements have proved very useful also for the characterization of paramagnetic proteins, because they can help characterise a number of parameters including the number, distance and residence time of water molecules coordinated to the paramagnetic center, the reorientation correlation times and the electron spin relaxation time, and the electronic structure at the metal site.

Tyrosine phosphorylation regulates hnRNPA2 granule protein partitioning and reduces neurodegeneration

mRNA transport in neurons requires formation of transport granules containing many protein components, and subsequent alterations in phosphorylation status can release transcripts for translation. Further, mutations in a structurally disordered domain of the transport granule protein hnRNPA2 increase its aggregation and cause hereditary proteinopathy of neurons, myocytes, and bone. We examine in vitro hnRNPA2 granule component phase separation, partitioning specificity, assembly/disassembly, and the link to neurodegeneration. Transport granule components hnRNPF and ch-TOG interact weakly with hnRNPA2 yet partition specifically into liquid phase droplets with the low complexity domain (LC) of hnRNPA2, but not FUS LC. In vitro hnRNPA2 tyrosine phosphorylation reduces hnRNPA2 phase separation, prevents partitioning of hnRNPF and ch-TOG into hnRNPA2 LC droplets, and decreases aggregation of hnRNPA2 disease variants. The expression of chimeric hnRNPA2 D290V in Caenorhabditis elegans results in stress-induced glutamatergic neurodegeneration; this neurodegeneration is rescued by loss of tdp-1, suggesting gain-of-function toxicity. The expression of Fyn, a tyrosine kinase that phosphorylates hnRNPA2, reduces neurodegeneration associated with chimeric hnRNPA2 D290V. These data suggest a model where phosphorylation alters LC interaction specificity, aggregation, and toxicity.

Analyzing In Silico the Relationship Between the Activation of the Edema Factor and Its Interaction With Calmodulin

Molecular dynamics (MD) simulations have been recorded on the complex between the edema factor (EF) of Bacilllus anthracis and calmodulin (CaM), starting from a structure with the orthosteric inhibitor adefovir bound in the EF catalytic site. The starting structure has been destabilized by alternately suppressing different co-factors, such as adefovir ligand or ions, revealing several long-distance correlations between the conformation of CaM, the geometry of the CaM/EF interface, the enzymatic site and the overall organization of the complex. An allosteric communication between CaM/EF interface and the EF catalytic site, highlighted by these correlations, was confirmed by several bioinformatics approaches from the literature. A network of hydrogen bonds and stacking interactions extending from the helix V of of CaM, and the residues of the switches A, B and C, and connecting to catalytic site residues, is a plausible candidate for the mediation of allosteric communication. The greatest variability in volume between the different MD conditions was also found for cavities present at the EF/CaM interface and in the EF catalytic site. The similarity between the predictions from literature and the volume variability might introduce the volume variability as new descriptor of allostery.

Millisecond Time-Resolved Solid-State NMR Reveals a Two-Stage Molecular Mechanism for Formation of Complexes between Calmodulin and a Target Peptide from Myosin Light Chain Kinase

Calmodulin (CaM) mediates a wide range of biological responses to changes in intracellular Ca²⁺ concentrations through its calcium-dependent binding affinities to numerous target proteins. Binding of two Ca²⁺ ions to each of the two four-helix-bundle domains of CaM results in major conformational changes that create a potential binding site for the CaM binding domain of a target protein, which also undergoes major conformational changes to form the complex with CaM. Details of the molecular mechanism of complex formation are not well established, despite numerous structural, spectroscopic, thermodynamic, and kinetic studies. Here, we report a study of the process by which the 26-residue peptide M13, which represents the CaM binding domain of skeletal muscle myosin light chain kinase, forms a complex with CaM in the presence of excess Ca²⁺ on the millisecond time scale. Our experiments use a combination of selective ¹³C labeling of CaM and M13, rapid mixing of CaM solutions with M13/Ca²⁺ solutions, rapid freeze-quenching of the mixed solutions, and low-temperature solid state nuclear magnetic resonance (ssNMR) enhanced by dynamic nuclear polarization. From measurements of the dependence of 2D ¹³C-¹³C ssNMR spectra on the time between mixing and freezing, we find that the N-terminal portion of M13 converts from a conformationally disordered state to an α-helix and develops contacts with the C-terminal domain of CaM in about 2 ms. The C-terminal portion of M13 becomes α-helical and develops contacts with the N-terminal domain of CaM more slowly, in about 8 ms. The level of structural order in the CaM/M13/Ca²⁺ complexes, indicated by ¹³C ssNMR line widths, continues to increase beyond 27 ms.

Coupled intra- and interdomain dynamics support domain cross-talk in Pin1

The functional mechanisms of multidomain proteins often exploit interdomain interactions, or "cross-talk." An example is human Pin1, an essential mitotic regulator consisting of a Trp-Trp (WW) domain flexibly tethered to a peptidyl-prolyl isomerase (PPIase) domain, resulting in interdomain interactions important for Pin1 function. Substrate binding to the WW domain alters its transient contacts with the PPIase domain via means that are only partially understood. Accordingly, we have investigated Pin1 interdomain interactions using NMR paramagnetic relaxation enhancement (PRE) and molecular dynamics (MD) simulations. The PREs show that apo-Pin1 samples interdomain contacts beyond the range suggested by previous structural studies. They further show that substrate binding to the WW domain simultaneously alters interdomain separation and the internal conformation of the WW domain. A 4.5-μs all-atom MD simulation of apo-Pin1 suggests that the fluctuations of interdomain distances are correlated with fluctuations of WW domain interresidue contacts involved in substrate binding. Thus, the interdomain/WW domain conformations sampled by apo-Pin1 may already include a range of conformations appropriate for binding Pin1's numerous substrates. The proposed coupling between intra-/interdomain conformational fluctuations is a consequence of the dynamic modular architecture of Pin1. Such modular architecture is common among cell-cycle proteins; thus, the WW-PPIase domain cross-talk mechanisms of Pin1 may be relevant for their mechanisms as well.

NMR Analysis of Apo Glutamine-Binding Protein Exposes Challenges in the Study of Interdomain Dynamics

Glutamine-binding protein (GlnBP) displays an apo, "open" and a holo, "closed" crystal form, mutually related by a rigid-body reorientation of its domains. A fundamental question about such large-scale conformational transitions, whether the closed state exists in the absence of ligand, is controversial in the case of GlnBP. NMR observations have indicated no evidence of the closed form, whereas experimentally validated computations have suggested a remarkable ca. 40 % population. Herein, a paramagnetic NMR strategy designed to detect the putative apo-closed species shows that a major population of the latter is highly improbable. Further, NMR residual dipolar couplings collected under three anisotropic conditions do not reveal differential domain alignment and establish that the average solution conformation is satisfied by the apo-open crystal structure. Our results indicate that the computational prediction of large-scale interdomain motions is not trivial and may lead to erroneous conclusions without proper experimental validation.

Integrative Approaches in Structural Biology: A More Complete Picture from the Combination of Individual Techniques

With the recent technological and computational advancements, structural biology has begun to tackle more and more difficult questions, including complex biochemical pathways and transient interactions among macromolecules. This has demonstrated that, to approach the complexity of biology, one single technique is largely insufficient and unable to yield thorough answers, whereas integrated approaches have been more and more adopted with successful results. Traditional structural techniques (X-ray crystallography and Nuclear Magnetic Resonance (NMR)) and the emerging ones (cryo-electron microscopy (cryo-EM), Small Angle X-ray Scattering (SAXS)), together with molecular modeling, have pros and cons which very nicely complement one another. In this review, three examples of synergistic approaches chosen from our previous research will be revisited. The first shows how the joint use of both solution and solid-state NMR (SSNMR), X-ray crystallography, and cryo-EM is crucial to elucidate the structure of polyethylene glycol (PEG)ylated asparaginase, which would not be obtainable through any of the techniques taken alone. The second deals with the integrated use of NMR, X-ray crystallography, and SAXS in order to elucidate the catalytic mechanism of an enzyme that is based on the flexibility of the enzyme itself. The third one shows how it is possible to put together experimental data from X-ray crystallography and NMR restraints in order to refine a protein model in order to obtain a structure which simultaneously satisfies both experimental datasets and is therefore closer to the ‘real structure’.

Consequences of Hydrophobic Nanotube Binding on the Functional Dynamics of Signaling Protein Calmodulin

The wide applications of nanomaterials in industry and our daily life have raised growing concerns on their toxicity to human body. Increasing evidence links the cytotoxicity of nanoparticles to the disruption of cellular signaling pathways. Here, we report a computational study on the mechanisms of the cytotoxicity of carbon nanotubes (CNTs) by investigating the direct impacts of CNTs on the functional motions of calmodulin (CaM), which is one of the most important signaling proteins in a cell, and its signaling function relies on the Ca2+ binding-coupled conformational switching. Computational simulations with a coarse-grained model showed that binding of CNTs modifies the conformational equilibrium of CaM and induces the closed-to-open conformational transition, leading to the loss of its Ca2+-sensing ability. In addition, the binding of CNTs drastically increases the calcium affinity of CaM, which may disrupt the Ca2+ homeostasis in a cell. These results suggest that the binding of hydrophobic nanotubes not only inhibits the signaling function of CaM as a calcium sensor but also renders CaM to toxic species through sequestering Ca2+ from other competing calcium-binding proteins, suggesting a new physical mechanism of the cytotoxicity of nanoparticles.

High Accuracy Protein Structures from Minimal Sparse Paramagnetic Solid-State NMR Restraints

There is a pressing need for new computational tools to integrate data from diverse experimental approaches in structural biology. We present a strategy that combines sparse paramagnetic solid-state NMR restraints with physics-based atomistic simulations. Our approach explicitly accounts for uncertainty in the interpretation of experimental data through the use of a semi-quantitative mapping between the data and the restraint energy that is calibrated by extensive simulations. We apply our approach to solid-state NMR data for the model protein GB1 labeled with Cu²⁺ -EDTA at six different sites. We are able to determine the structure to 0.9 Å accuracy within a single day of computation on a GPU cluster. We further show that in some cases, the data from only a single paramagnetic tag are sufficient for accurate folding.

Assessing protein conformational landscapes: integration of DEER data in Maximum Occurrence analysis

The properties of the conformational landscape of a biomolecule are of capital importance to understand its function. It is widely accepted that a statistical ensemble is far more representative than a single structure, especially for proteins with disordered regions. While experimental data provide the most important handle on the conformational variability that the system is experiencing, they usually report on either time or ensemble averages. Since the available conformations largely outnumber the (independent) available experimental data, the latter can be equally well reproduced by a variety of ensembles. We have proposed the Maximum Occurrence (MaxOcc) approach to provide an upper bound of the statistical weight of each conformation. This method is expected to converge towards the true statistical weights by increasing the number of independent experimental datasets. In this paper we explore the ability of DEER (Double Electron Electron Resonance) data, which report on the distance distribution between two spin labels attached to a biomolecule, to restrain the MaxOcc values and its complementarity to previously introduced experimental techniques such as NMR and Small-Angle X-ray Scattering. We here present the case of Ca2+ bound calmodulin (CaM) as a test case and show that DEER data impose a sizeable reduction of the conformational space described by high MaxOcc conformations.

Arrhythmia mutations in calmodulin cause conformational changes that affect interactions with the cardiac voltage-gated calcium channel

Calmodulin (CaM) represents one of the most conserved proteins among eukaryotes and is known to bind and modulate more than a 100 targets. Recently, several disease-associated mutations have been identified in the CALM genes that are causative of severe cardiac arrhythmia syndromes. Although several mutations have been shown to affect the function of various cardiac ion channels, direct structural insights into any CaM disease mutation have been lacking. Here we report a crystallographic and NMR investigation of several disease mutant CaMs, linked to long-QT syndrome, in complex with the IQ domain of the cardiac voltage-gated calcium channel (Ca_V1.2). Surprisingly, two mutants (D95V, N97I) cause a major distortion of the C-terminal lobe, resulting in a pathological conformation not reported before. These structural changes result in altered interactions with the Ca_V1.2 IQ domain. Another mutation (N97S) reduces the affinity for Ca²⁺ by introducing strain in EF hand 3. A fourth mutant (F141L) shows structural changes in the Ca²⁺-free state that increase the affinity for the IQ domain. These results thus show that different mechanisms underlie the ability of CaM disease mutations to affect Ca²⁺-dependent inactivation of the voltage-gated calcium channel.

Chemical shift-based methods in NMR structure determination

Chemical shifts are highly sensitive probes harnessed by NMR spectroscopists and structural biologists as conformational parameters to characterize a range of biological molecules. Traditionally, assignment of chemical shifts has been a labor-intensive process requiring numerous samples and a suite of multidimensional experiments. Over the past two decades, the development of complementary computational approaches has bolstered the analysis, interpretation and utilization of chemical shifts for elucidation of high resolution protein and nucleic acid structures. Here, we review the development and application of chemical shift-based methods for structure determination with a focus on ab initio fragment assembly, comparative modeling, oligomeric systems, and automated assignment methods. Throughout our discussion, we point out practical uses, as well as advantages and caveats, of using chemical shifts in structure modeling. We additionally highlight (i) hybrid methods that employ chemical shifts with other types of NMR restraints (residual dipolar couplings, paramagnetic relaxation enhancements and pseudocontact shifts) that allow for improved accuracy and resolution of generated 3D structures, (ii) the utilization of chemical shifts to model the structures of sparsely populated excited states, and (iii) modeling of sidechain conformations. Finally, we briefly discuss the advantages of contemporary methods that employ sparse NMR data recorded using site-specific isotope labeling schemes for chemical shift-driven structure determination of larger molecules. With this review, we aim to emphasize the accessibility and versatility of chemical shifts for structure determination of challenging biological systems, and to point out emerging areas of development that lead us towards the next generation of tools.

Gold nanocrystal labels provide a sequence-to-3D structure map in SAXS reconstructions

Small-angle x-ray scattering (SAXS) is a powerful technique to probe the structure of biological macromolecules and their complexes under virtually arbitrary solution conditions, without the need for crystallization. While it is possible to reconstruct molecular shapes from SAXS data ab initio, the resulting electron density maps have a resolution of ~1 nm and are often insufficient to reliably assign secondary structure elements or domains. We show that SAXS data of gold-labeled samples significantly enhance the information content of SAXS measurements, allowing the unambiguous assignment of macromolecular sequence motifs to specific locations within a SAXS structure. We first demonstrate our approach for site-specifically internally and end-labeled DNA and an RNA motif. In addition, we present a protocol for highly uniform and site-specific labeling of proteins with small (~1.4 nm diameter) gold particles and apply our method to the signaling protein calmodulin. In all cases, the position of the small gold probes can be reliably identified in low-resolution electron density maps. Enhancing low-resolution measurements by site-selective gold labeling provides an attractive approach to aid modeling of a large range of macromolecular systems.

Theory and practice of using solvent paramagnetic relaxation enhancement to characterize protein conformational dynamics

Paramagnetic relaxation enhancement (PRE) has been established as a powerful tool in NMR for investigating protein structure and dynamics. The PRE is usually measured with a paramagnetic probe covalently attached at a specific site of an otherwise diamagnetic protein. The present work provides the numerical formulation for probing protein structure and conformational dynamics based on the solvent PRE (sPRE) measurement, using two alternative approaches. An inert paramagnetic cosolute randomly collides with the protein, and the resulting sPRE manifests the relative solvent exposure of protein nuclei. To make the back-calculated sPRE values most consistent with the observed values, the protein structure is either refined against the sPRE, or an ensemble of conformers is selected from a pre-generated library using a Monte Carlo algorithm. The ensemble structure comprises either N conformers of equal occupancy, or two conformers with different relative populations. We demonstrate the sPRE method using GB1, a structurally rigid protein, and calmodulin, a protein comprising two domains and existing in open and closed states. The sPRE can be computed with a stand-alone program for rapid evaluation, or with the invocation of a module in the latest release of the structure calculation software Xplor-NIH. As a label-free method, the sPRE measurement can be readily integrated with other biophysical techniques. The current limitations of the sPRE method are also discussed, regarding accurate measurement and theoretical calculation, model selection and suitable timescale.

Multicolor protein FRET with tryptophan, selective coumarin-cysteine labeling, and genetic acridonylalanine encoding

Site-specific fluorescence probes can be used to measure distances within proteins when used as part of a Förster resonance energy transfer (FRET) pair. Here we report the synthesis of a coumarin maleimide (Mcm-Mal) that is fluorogenic upon reaction with cysteine. We demonstrate that cysteine, acridonylalanine (Acd) double mutant proteins can be produced by unnatural amino acid mutagenesis and reacted with Mcm-Mal to generate Mcm/Acd labeled proteins for FRET studies. The Mcm/Acd FRET pair is minimally-perturbing, easy to install, and well-suited to studying protein distances in the 15-40 Å range. Furthermore, Mcm/Acd labeling can be combined with tryptophan fluorescence in three color FRET to monitor multiple interactions in one experiment.

Conserved conformational selection mechanism of Hsp70 chaperone-substrate interactions

Molecular recognition is integral to biological function and frequently involves preferred binding of a molecule to one of several exchanging ligand conformations in solution. In such a process the bound structure can be selected from the ensemble of interconverting ligands a priori (conformational selection, CS) or may form once the ligand is bound (induced fit, IF). Here we focus on the ubiquitous and conserved Hsp70 chaperone which oversees the integrity of the cellular proteome through its ATP-dependent interaction with client proteins. We directly quantify the flux along CS and IF pathways using solution NMR spectroscopy that exploits a methyl TROSY effect and selective isotope-labeling methodologies. Our measurements establish that both bacterial and human Hsp70 chaperones interact with clients by selecting the unfolded state from a pre-existing array of interconverting structures, suggesting a conserved mode of client recognition among Hsp70s and highlighting the importance of molecular dynamics in this recognition event.

Proteins are the workhorses of a cell and are involved in almost all biological processes. Newly made proteins need to ‘fold’ into precise three-dimensional shapes in order to carry out their roles. However, proteins sometimes fold incorrectly or unfold. These protein forms are not able to work effectively and in some cases may even cause diseases.

Chaperone proteins help other proteins to fold correctly and are found in living organisms ranging in complexity from bacteria to humans. There are many different types of chaperones that play different roles inside cells. One, called Hsp70, binds to proteins that are incorrectly folded to help them to mature into their correct structures. However, it was not clear whether Hsp70 can also associate with the mature, correctly folded form of the proteins.

A technique called Nuclear Magnetic Resonance (NMR) spectroscopy can distinguish between mature, unfolded and chaperone-bound forms of the same protein. Sekhar et al. therefore used NMR to investigate which forms of a protein Hsp70 binds to. This revealed that both the bacterial and human versions of the Hsp70 chaperone interact only with unfolded proteins.

The results presented by Sekhar et al. also explain why Hsp70 does not disrupt the routine workings of the cell: because it does not bind to mature forms of proteins. These observations extend our understanding of how chaperones assist in folding proteins, and fit into a broader research theme exploring how proteins recognize one another. It will now be interesting to see whether the same mechanism holds for more complex forms of proteins, such as aggregates, or larger protein structures with regions of both folded and unfolded elements.

Maintaining and Enhancing Diversity of Sampled Protein Conformations in Robotics-Inspired Methods

The ability to efficiently sample structurally diverse protein conformations allows one to gain a high-level view of a protein's energy landscape. Algorithms from robot motion planning have been used for conformational sampling, and several of these algorithms promote diversity by keeping track of "coverage" in conformational space based on the local sampling density. However, large proteins present special challenges. In particular, larger systems require running many concurrent instances of these algorithms, but these algorithms can quickly become memory intensive because they typically keep previously sampled conformations in memory to maintain coverage estimates. In addition, robotics-inspired algorithms depend on defining useful perturbation strategies for exploring the conformational space, which is a difficult task for large proteins because such systems are typically more constrained and exhibit complex motions. In this article, we introduce two methodologies for maintaining and enhancing diversity in robotics-inspired conformational sampling. The first method addresses algorithms based on coverage estimates and leverages the use of a low-dimensional projection to define a global coverage grid that maintains coverage across concurrent runs of sampling. The second method is an automatic definition of a perturbation strategy through readily available flexibility information derived from B-factors, secondary structure, and rigidity analysis. Our results show a significant increase in the diversity of the conformations sampled for proteins consisting of up to 500 residues when applied to a specific robotics-inspired algorithm for conformational sampling. The methodologies presented in this article may be vital components for the scalability of robotics-inspired approaches.

Utilizing tagged paramagnetic shift reagents to monitor protein dynamics by NMR

Calmodulin is a ubiquitous calcium sensor protein, known to serve as a critical interaction hub with a wide range of signaling partners. While the holo form of calmodulin (CaM-4Ca²⁺) has a well-defined ground state structure, it has been shown to undergo exchange, on a millisecond timescale, to a conformation resembling that of the peptide bound state. Tagged paramagnetic relaxation agents have been previously used to identify long-range dipolar interactions through relaxation effects on nuclear spins of interest. In the case of calmodulin, this lead to the determination of the relative orientation of the N- and C-terminal domains and the presence of a weakly populated peptide bound like state. Here, we make use of pseudocontact shifts from a tagged paramagnetic shift reagent which allows us to define minor states both in ¹³C and ¹⁵N NMR spectra and through ¹³C- and ¹⁵N-edited ¹H-CPMG relaxation dispersion measurements. This is validated by pulsed EPR (DEER) spectroscopy which reveals an ensemble consisting of a compact peptide-bound like conformer, an intermediate peptide-bound like conformer, and a (dumbbell-like) extended ground state conformer of CaM-4Ca²⁺, where addition of the MLCK peptide increases the population of the peptide-bound conformers. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

Opposing Intermolecular Tuning of Ca2+ Affinity for Calmodulin by Neurogranin and CaMKII Peptides

We investigated the impact of bound calmodulin (CaM)-target compound structure on the affinity of calcium (Ca²⁺) by integrating coarse-grained models and all-atomistic simulations with nonequilibrium physics. We focused on binding between CaM and two specific targets, Ca²⁺/CaM-dependent protein kinase II (CaMKII) and neurogranin (Ng), as they both regulate CaM-dependent Ca²⁺ signaling pathways in neurons. It was shown experimentally that Ca²⁺/CaM (holoCaM) binds to the CaMKII peptide with overwhelmingly higher affinity than Ca²⁺-free CaM (apoCaM); the binding of CaMKII peptide to CaM in return increases the Ca²⁺ affinity for CaM. However, this reciprocal relation was not observed in the Ng peptide (Ng_13–49), which binds to apoCaM or holoCaM with binding affinities of the same order of magnitude. Unlike the holoCaM-CaMKII peptide, whose structure can be determined by crystallography, the structural description of the apoCaM-Ng_13–49 is unknown due to low binding affinity, therefore we computationally generated an ensemble of apoCaM-Ng_13–49 structures by matching the changes in the chemical shifts of CaM upon Ng_13–49 binding from nuclear magnetic resonance experiments. Next, we computed the changes in Ca²⁺ affinity for CaM with and without binding targets in atomistic models using Jarzynski’s equality. We discovered the molecular underpinnings of lowered affinity of Ca²⁺ for CaM in the presence of Ng_13–49 by showing that the N-terminal acidic region of Ng peptide pries open the β-sheet structure between the Ca²⁺ binding loops particularly at C-domain of CaM, enabling Ca²⁺ release. In contrast, CaMKII peptide increases Ca²⁺ affinity for the C-domain of CaM by stabilizing the two Ca²⁺ binding loops. We speculate that the distinctive structural difference in the bound complexes of apoCaM-Ng_13–49 and holoCaM-CaMKII delineates the importance of CaM’s progressive mechanism of target binding on its Ca²⁺ binding affinities.

Cobalt-based paramagnetic probe to study RNA-protein interactions by NMR

Paramagnetic resonance enhancement (PRE) is an NMR technique that allows studying three-dimensional structures of RNA-protein complexes in solution. RNA strands are typically spin labeled using nitroxide reagents, which provide minimal perturbation to the native structure. The current work describes an alternative approach, which is based on a Co2+-based probe that can be covalently attached to RNA in the vicinity of the protein's binding site using 'click' chemistry. Similar to nitroxide spin labels, the transition metal based probe is capable of attenuating NMR signal intensities from protein residues localized <40Å away. The extent of attenuation is related to the probe's distance, thus allowing for construction of the protein's contact surface map. This new paradigm has been applied to study binding of HIV-1 nucleocapsid protein 7, NCp7, to a model RNA pentanucleotide.

Molecular mechanism of multispecific recognition of Calmodulin through conformational changes

Calmodulin (CaM) is found to have the capability to bind multiple targets. Investigations on the association mechanism of CaM to its targets are crucial for understanding protein-protein binding and recognition. Here, we developed a structure-based model to explore the binding process between CaM and skMLCK binding peptide. We found the cooperation between nonnative electrostatic interaction and nonnative hydrophobic interaction plays an important role in nonspecific recognition between CaM and its target. We also found that the conserved hydrophobic anchors of skMLCK and binding patches of CaM are crucial for the transition from high affinity to high specificity. Furthermore, this association process involves simultaneously both local conformational change of CaM and global conformational changes of the skMLCK binding peptide. We found a landscape with a mixture of the atypical "induced fit," the atypical "conformational selection," and "simultaneously binding-folding," depending on the synchronization of folding and binding. Finally, we extend our discussions on multispecific binding between CaM and its targets. These association characteristics proposed for CaM and skMLCK can provide insights into multispecific binding of CaM.

Cu(II)-Based Paramagnetic Probe to Study RNA-Protein Interactions by NMR

Paramagnetic NMR techniques allow for studying three-dimensional structures of RNA-protein complexes. In particular, paramagnetic relaxation enhancement (PRE) data can provide valuable information about long-range distances between different structural components. For PRE NMR experiments, oligonucleotides are typically spin-labeled using nitroxide reagents. The current work describes an alternative approach involving a Cu(II) cyclen-based probe that can be covalently attached to an RNA strand in the vicinity of the protein's binding site using "click" chemistry. The approach has been applied to study binding of HIV-1 nucleocapsid protein 7 (NCp7) to a model RNA pentanucleotide, 5'-ACGCU-3'. Coordination of the paramagnetic metal to glutamic acid residue of NCp7 reduced flexibility of the probe, thus simplifying interpretation of the PRE data. NMR experiments showed attenuation of signal intensities from protein residues localized in proximity to the paramagnetic probe as the result of RNA-protein interactions. The extent of the attenuation was related to the probe's proximity allowing us to construct the protein's contact surface map.

Modeling Protein Excited-state Structures from "Over-length" Chemical Cross-links

Chemical cross-linking coupled with mass spectroscopy (CXMS) provides proximity information for the cross-linked residues and is used increasingly for modeling protein structures. However, experimentally identified cross-links are sometimes incompatible with the known structure of a protein, as the distance calculated between the cross-linked residues far exceeds the maximum length of the cross-linker. The discrepancies may persist even after eliminating potentially false cross-links and excluding intermolecular ones. Thus the "over-length" cross-links may arise from alternative excited-state conformation of the protein. Here we present a method and associated software DynaXL for visualizing the ensemble structures of multidomain proteins based on intramolecular cross-links identified by mass spectrometry with high confidence. Representing the cross-linkers and cross-linking reactions explicitly, we show that the protein excited-state structure can be modeled with as few as two over-length cross-links. We demonstrate the generality of our method with three systems: calmodulin, enzyme I, and glutamine-binding protein, and we show that these proteins alternate between different conformations for interacting with other proteins and ligands. Taken together, the over-length chemical cross-links contain valuable information about protein dynamics, and our findings here illustrate the relationship between dynamic domain movement and protein function.

Separating Dipolar and Chemical Exchange Magnetization Transfer Processes in 1 H-CEST

An amide ¹ H-Chemical Exchange Saturation Transfer (CEST) experiment is presented for studies of conformational exchange in proteins. The approach, exploiting spin-state-selective magnetization transfer, completely suppresses undesired NOE-based dips in CEST profiles so that chemical exchange processes can be studied. The methodology is demonstrated with applications involving proteins that interconvert on the millisecond timescale between major and invisible minor states, and accurate amide ¹ H chemical shifts of the minor conformer are obtained in each case. The spin-state-selective magnetization transfer approach offers unique possibilities for quantitative studies of protein exchange through ¹ H-CEST.