Relaxation Dispersion Experiments Research Articles

Slow conformational exchange and overall rocking motion in ubiquitin protein crystals

Proteins perform their functions in solution but their structures are most frequently studied inside crystals. Here we probe how the crystal packing alters microsecond dynamics, using solid-state NMR measurements and multi-microsecond MD simulations of different crystal forms of ubiquitin. In particular, near-rotary-resonance relaxation dispersion (NERRD) experiments probe angular backbone motion, while Bloch–McConnell relaxation dispersion data report on fluctuations of the local electronic environment. These experiments and simulations reveal that the packing of the protein can significantly alter the thermodynamics and kinetics of local conformational exchange. Moreover, we report small-amplitude reorientational motion of protein molecules in the crystal lattice with an ~3–5° amplitude on a tens-of-microseconds time scale in one of the crystals, but not in others. An intriguing possibility arises that overall motion is to some extent coupled to local dynamics. Our study highlights the importance of considering the packing when analyzing dynamics of crystalline proteins.

Elevated μs-ms timescale backbone dynamics in the transition state analog form of arginine kinase.

Arginine kinase catalyzes reversible phosphoryl transfer between arginine and ATP. Crystal structures of arginine kinase in an open, substrate-free form and closed, transition state analog (TSA) complex indicate that the enzyme undergoes substantial domain and loop rearrangements required for substrate binding, catalysis, and product release. Nuclear magnetic resonance (NMR) has shown that substrate-free arginine kinase is rigid on the ps-ns timescale (average S2=0.84±0.08) yet quite dynamic on the µs-ms timescale (35 residues with Rex, 12%), and that movements of the N-terminal domain and the loop comprising residues I182-G209 are rate-limiting on catalysis. Here, NMR of the TSA-bound enzyme shows similar rigidity on the ps-ns timescale (average S2=0.91±0.05) and substantially increased μs-ms timescale dynamics (77 residues; 22%). Many of the residues displaying μs-ms dynamics in NMR Carr-Purcell-Meiboom-Gill (CPMG) 15N backbone relaxation dispersion experiments of the TSA complex are also dynamic in substrate-free enzyme. However, the presence of additional dynamic residues in the TSA-bound form suggests that dynamics extend through much of the C-terminal domain, which indicates that in the closed form, a larger fraction of the protein takes part in conformational transitions to the excited state(s). Conformational exchange rate constants (kex) of the TSA complex are all approximately 2500s-1, higher than any observed in the substrate-free enzyme (800-1900s-1). Elevated μs-ms timescale protein dynamics in the TSA-bound enzyme is more consistent with recently postulated catalytic networks involving multiple interconnected states at each step of the reaction, rather than a classical single stabilized transition state.

Altering the allosteric pathway in IGPS suppresses millisecond motions and catalytic activity

Imidazole glycerol phosphate synthase (IGPS) is a V-type allosteric enzyme, meaning that its catalytic rate is critically dependent on activation by its allosteric ligand, N'-[(5'-phosphoribulosyl)formimino]-5-aminoimidazole-4-carboxamide ribonucleotide (PRFAR). The allosteric mechanism of IGPS is reliant on millisecond conformational motions for efficient catalysis. We engineered four mutants of IGPS designed to disrupt millisecond motions and allosteric coupling to identify regions that are critical to IGPS function. Multiple-quantum Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments and NMR chemical shift titrations reveal diminished enzyme flexibility and a reshaping of the allosteric connectivity in each mutant construct, respectively. The functional relevance of the observed motional quenching is confirmed by significant reductions in glutaminase kinetic activity and allosteric ligand binding affinity. This work presents relevant conclusions toward the control of protein allostery and design of unique allosteric sites for potential enzyme inhibitors with regulatory or therapeutic benefit.

Protein NMR Studies of Substrate Binding to Human Blood Group A and B Glycosyltransferases.

Donor and acceptor substrate binding to human blood group A and B glycosyltransferases (GTA, GTB) has been studied by a variety of protein NMR experiments. Prior crystallographic studies had shown these enzymes to adopt an open conformation in the absence of substrates. Binding either of the donor substrate UDP-Gal or of UDP induces a semiclosed conformation. In the presence of both donor and acceptor substrates, the enzymes shift towards a closed conformation with ordering of an internal loop and the C-terminal residues, which then completely cover the donor-binding pocket. Chemical-shift titrations of uniformly 2 H,15 N-labeled GTA or GTB with UDP affected about 20 % of all crosspeaks in 1 H,15 N TROSY-HSQC spectra, reflecting substantial plasticity of the enzymes. On the other hand, it is this conformational flexibility that impedes NH backbone assignments. Chemical-shift-perturbation experiments with δ1-[13 C]methyl-Ile-labeled samples revealed two Ile residues-Ile123 at the bottom of the UDP binding pocket, and Ile192 as part of the internal loop-that were significantly disturbed upon stepwise addition of UDP and H-disaccharide, also revealing long-range perturbations. Finally, methyl TROSY-based relaxation dispersion experiments do not reveal micro- to millisecond timescale motions. Although this study reveals substantial conformational plasticity of GTA and GTB, the matter of how binding of substrates shifts the enzymes into catalytically competent states remains enigmatic.

The Molecular Mechanisms Underlying the Hijack of Host Proteins by the 1918 Spanish Influenza Virus.

The 1918 Spanish influenza A virus (IAV) caused one of the most serious pandemics in history. The nonstructural protein 1 (NS1) of the 1918 IAV hijacks the interaction between human CrkII and JNK1. Little is, however, known about its molecular mechanism. Here, we performed X-ray crystallography, NMR relaxation dispersion experiment, and fluorescence spectroscopy to determine the structural, kinetic, and thermodynamic mechanisms underlying the hijacking of CrkII by 1918 IAV NS1. We observed that the interaction between a proline-rich motif in NS1 and the N-terminal SH3 domain of CrkII displays strikingly rapid kinetics and exceptionally high affinity with 100-fold faster kon and 3300-fold lower Kd compared to those for the CrkII-JNK1 interaction. These results provide molecular insight into the mechanism by which 1918 IAV NS1 hijacks CrkII and disrupts its interactions with critical cellular signaling proteins.

The Michaelis Complex of Arginine Kinase Samples the Transition State at a Frequency That Matches the Catalytic Rate.

Arginine kinase (AK), which is a member of the phosphagen kinase family, serves as a model system for studying the structural and dynamic determinants of biomolecular enzyme catalysis of all major states involved of the enzymatic cycle. These states are the apo state (substrate free), the Michaelis complex analogue AK:Arg:Mg·AMPPNP (MCA), a product complex analogue AK:pAIE:Mg·ADP (PCA), and the transition state analogue AK:Arg:Mg·ADP:NO3- (TSA). The conformational dynamics of these states have been studied by NMR relaxation dispersion measurements of the methyl groups of the Ile, Leu, and Val residues at two static magnetic fields. Although all states undergo significant amounts of μs-ms time scale dynamics, only the MCA samples a dominant excited state that resembles the TSA, as evidenced by the strong correlation between the relaxation dispersion derived chemical shift differences Δω and the equilibrium chemical shift differences Δδ of these states. The average lifetime of the MCA is 36 ms and the free energy difference to the TSA-like form is 8.5 kJ/mol. It is shown that the conformational energy landscape of the Michaelis complex analogue is shaped in a way that at room temperature it channels passage to the transition state, thereby determining the rate-limiting step of the phosphorylation reaction of arginine. Conversely, relaxation dispersion experiments of the TSA reveal that it samples the structures of the Michaelis complex analogue or the apo state as its dominant excited state. This reciprocal behavior shows that the free energy of the TSA, with all ligands bound, is lower by only about 8.9 kJ/mol than that of the Michaelis or apo complex conformations with the TSA ligands present.

Probing conformational dynamics in biomolecules via chemical exchange saturation transfer: a primer.

Although Chemical Exchange Saturation Transfer (CEST) type NMR experiments have been used to study chemical exchange processes in molecules since the early 1960s, there has been renewed interest in the past several years in using this approach to study biomolecular conformational dynamics. The methodology is particularly powerful for the study of sparsely populated, transiently formed conformers that are recalcitrant to investigation using traditional biophysical tools, and it is complementary to relaxation dispersion and magnetization transfer experiments that have traditionally been used to study chemical exchange processes. Here we discuss the concepts behind the CEST experiment, focusing on practical aspects as well, we review some of the pulse sequences that have been developed to characterize protein and RNA conformational dynamics, and we discuss a number of examples where the CEST methodology has provided important insights into the role of dynamics in biomolecular function.

Application of geometric approximation to the CPMG experiment: Two- and three-site exchange

The Carr–Purcell–Meiboom–Gill (CPMG) experiment is one of the most classical and well-known relaxation dispersion experiments in NMR spectroscopy, and it has been successfully applied to characterize biologically relevant conformational dynamics in many cases. Although the data analysis of the CPMG experiment for the 2-site exchange model can be facilitated by analytical solutions, the data analysis in a more complex exchange model generally requires computationally-intensive numerical analysis. Recently, a powerful computational strategy, geometric approximation, has been proposed to provide approximate numerical solutions for the adiabatic relaxation dispersion experiments where analytical solutions are neither available nor feasible. Here, we demonstrate the general potential of geometric approximation by providing a data analysis solution of the CPMG experiment for both the traditional 2-site model and a linear 3-site exchange model. The approximate numerical solution deviates less than 0.5% from the numerical solution on average, and the new approach is computationally 60,000-fold more efficient than the numerical approach. Moreover, we find that accurate dynamic parameters can be determined in most cases, and, for a range of experimental conditions, the relaxation can be assumed to follow mono-exponential decay. The method is general and applicable to any CPMG RD experiment (e.g. N, C′, Cα, Hα, etc.) The approach forms a foundation of building solution surfaces to analyze the CPMG experiment for different models of 3-site exchange. Thus, the geometric approximation is a general strategy to analyze relaxation dispersion data in any system (biological or chemical) if the appropriate library can be built in a physically meaningful domain.

Global Disordering in Stereo-Specific Protein Association

Protein-protein recognition is of fundamental importance for a myriad of biological processes and is ultimately a prerequisite for life as we know it. There exist several established mechanisms that promote formation of stereo-specific protein complexes. Many of these mechanism involve conformational changes of one or both proteins in dimeric assemblies as observed in “conformational selection” and “coupled folding and binding” scenarios. In “coupled folding and binding” events, at least one of the proteins undergoes a global ordering event. By using an integrated computational and experimental approach we have discovered that also global disordering can be a productive route for formation of a stereo-specific protein complex. This mechanism was observed for the chaperone binding domain of the Yersinia effector protein YopH upon binding to its specific chaperone SycH. These two proteins are crucial for type III secretion system mediated infectivity by Yersinia and several other gram negative pathogens. NMR relaxation dispersion experiments demonstrated that the otherwise well folded YopH protein dynamically samples an expanded high-energy state that corresponds to the SycH binding competent conformation. A structure of the protein complex determined from a hybrid SAXS and computational approach revealed that YopH wraps around SycH in a horse shoe like conformation. The binding model was validated by site specific YopH mutations that promoted the disordering event and at the same time displayed improved binding affinity towards SycH. Taken together the data illustrates a tight coupling between a proteins unfolding and functional free energy landscapes and add valuable mechanistic insight into protein-protein recognition.

Entropy Transfer between Residue Pairs and Allostery in Proteins: Quantifying Allosteric Communication in Ubiquitin

It has recently been proposed by Gunasakaran et al. that allostery may be an intrinsic property of all proteins. Here, we develop a computational method that can determine and quantify allosteric activity in any given protein. Based on Schreiber's transfer entropy formulation, our approach leads to an information transfer landscape for the protein that shows the presence of entropy sinks and sources and explains how pairs of residues communicate with each other using entropy transfer. The model can identify the residues that drive the fluctuations of others. We apply the model to Ubiquitin, whose allosteric activity has not been emphasized until recently, and show that there are indeed systematic pathways of entropy and information transfer between residues that correlate well with the activities of the protein. We use 600 nanosecond molecular dynamics trajectories for Ubiquitin and its complex with human polymerase iota and evaluate entropy transfer between all pairs of residues of Ubiquitin and quantify the binding susceptibility changes upon complex formation. We explain the complex formation propensities of Ubiquitin in terms of entropy transfer. Important residues taking part in allosteric communication in Ubiquitin predicted by our approach are in agreement with results of NMR relaxation dispersion experiments. Finally, we show that time delayed correlation of fluctuations of two interacting residues possesses an intrinsic causality that tells which residue controls the interaction and which one is controlled. Our work shows that time delayed correlations, entropy transfer and causality are the required new concepts for explaining allosteric communication in proteins.

Kinetic Insights into the Binding between the nSH3 Domain of CrkII and Proline-Rich Motifs in cAbl

The interaction between CrkII and cAbl is implicated in diverse cellular processes. This interaction starts with the binding of the N-terminal Src homology 3 (nSH3) domain of CrkII to the proline-rich motifs of cAbl (PRMscAbl). Despite its critical importance, the detailed binding mechanism between the nSH3 domain and PRMs remains elusive. In this study, we used nuclear magnetic resonance Carr-Purcell-Meiboom-Gill relaxation dispersion experiment to study the binding kinetics between the nSH3 domain of CrkII and PRMscAbl. Our results highlight that the nSH3 domain binds to three PRMscAbl with very high on- and off-rate constants, indicating the transient nature of the binding. To further characterize the binding transition state, we conducted the Eyring and linear free energy relationship analyses using temperature-dependent kinetic data. These data indicate that the binding transition state of the nSH3 domain and PRM is accompanied by small activation enthalpy, owing to partial desolvation of the transition state. These results also highlight the similarity between the transition and free states, in terms of structure and energetics. Although the binding of the nSH3 domain and PRM displays the features consistent with a diffusion-limited process within our experimental conditions, further tests are necessary to determine if the binding is a true diffusion-limited process.

Coexistence of multiple minor states of fatty acid binding protein and their functional relevance.

Proteins are dynamic over a wide range of timescales, but determining the number of distinct dynamic processes and identifying functionally relevant dynamics are still challenging. Here we present the study on human intestinal fatty acid binding protein (hIFABP) using a novel analysis of 15N relaxation dispersion (RD) and chemical shift saturation transfer (CEST) experiments. Through combined analysis of the two types of experiments, we found that hIFABP exists in a four-state equilibrium in which three minor states interconvert directly with the major state. According to conversion rates from the major “closed” state to minor states, these minor states are irrelevant to the function of fatty acid transport. Based on chemical shifts of the minor states which could not be determined from RD data alone but were extracted from a combined analysis of RD and CEST data, we found that all the minor states are native-like. This conclusion is further supported by hydrogen-deuterium exchange experiments. Direct conversions between the native state and native-like intermediate states may suggest parallel multitrack unfolding/folding pathways of hIFABP. Moreover, hydrogen-deuterium exchange data indicate the existence of another locally unfolded minor state that is relevant to the fatty acid entry process.

13C-NMR studies on disulfide bond isomerization in bovine pancreatic trypsin inhibitor (BPTI)

Conformational isomerization of disulfide bonds is associated with the dynamics and thus the functional aspects of proteins. However, our understanding of the isomerization is limited by experimental difficulties in probing it. We explored the disulfide conformational isomerization of the Cys14-Cys38 disulfide bond in bovine pancreatic trypsin inhibitor (BPTI), by performing an NMR line-shape analysis of its Cys carbon peaks. In this approach, 1D (13)C spectra were recorded at small temperature intervals for BPTI samples selectively labeled with site-specifically (13)C-enriched Cys, and the recorded peaks were displayed in the order of the temperature after the spectral scales were normalized to a carbon peak. Over the profile of the line-shape, exchange broadening that altered with temperature was manifested for the carbon peaks of Cys14 and Cys38. The Cys14-Cys38 disulfide bond reportedly exists in equilibrium between a high-populated (M) and two low-populated states (m c14 and m c38). Consistent with the three-site exchange model, biphasic exchange broadening arising from the two processes was observed for the peak of the Cys14 α-carbon. As the exchange broadening is maximized when the exchange rate equals the chemical shift difference in Hz between equilibrating sites, semi-quantitative information that was useful for establishing conditions for (13)C relaxation dispersion experiments was obtained through the carbon line-shape profile. With respect to the m c38 isomerization, the (1)H-(13)C signals at the β-position of the minor state were resolved from the major peaks and detected by exchange experiments at a low temperature.

Excited States of Nucleic Acids Probed by Proton Relaxation Dispersion NMR Spectroscopy.

In this work an improved stable isotope labeling protocol for nucleic acids is introduced. The novel building blocks eliminate/minimize homonuclear 13C and 1H scalar couplings thus allowing proton relaxation dispersion (RD) experiments to report accurately on the chemical exchange of nucleic acids. Using site‐specific 2H and 13C labeling, spin topologies are introduced into DNA and RNA that make 1H relaxation dispersion experiments applicable in a straightforward manner. The novel RNA/DNA building blocks were successfully incorporated into two nucleic acids. The A‐site RNA was previously shown to undergo a two site exchange process in the micro‐ to millisecond time regime. Using proton relaxation dispersion experiments the exchange parameters determined earlier could be recapitulated, thus validating the proposed approach. We further investigated the dynamics of the cTAR DNA, a DNA transcript that is involved in the viral replication cycle of HIV‐1. Again, an exchange process could be characterized and quantified. This shows the general applicablility of the novel labeling scheme for 1H RD experiments of nucleic acids.

Excited States of Nucleic Acids Probed by Proton Relaxation Dispersion NMR Spectroscopy

AbstractIn this work an improved stable isotope labeling protocol for nucleic acids is introduced. The novel building blocks eliminate/minimize homonuclear 13C and 1H scalar couplings thus allowing proton relaxation dispersion (RD) experiments to report accurately on the chemical exchange of nucleic acids. Using site‐specific 2H and 13C labeling, spin topologies are introduced into DNA and RNA that make 1H relaxation dispersion experiments applicable in a straightforward manner. The novel RNA/DNA building blocks were successfully incorporated into two nucleic acids. The A‐site RNA was previously shown to undergo a two site exchange process in the micro‐ to millisecond time regime. Using proton relaxation dispersion experiments the exchange parameters determined earlier could be recapitulated, thus validating the proposed approach. We further investigated the dynamics of the cTAR DNA, a DNA transcript that is involved in the viral replication cycle of HIV‐1. Again, an exchange process could be characterized and quantified. This shows the general applicablility of the novel labeling scheme for 1H RD experiments of nucleic acids.

Enhancing the Sensitivity of CPMG Relaxation Dispersion to Conformational Exchange Processes by Multiple-Quantum Spectroscopy.

A triple-quantum (1) H Carr-Purcell-Meiboom-Gill NMR relaxation dispersion experiment is presented that uses methyl group probes as reporters of conformational exchange in highly deuterated, methyl-protonated proteins. Significantly larger dispersion profiles, by as much as a factor of nine, can be obtained relative to single-quantum approaches, thus offering very significant advantages in applications involving interconverting conformers with only small changes in structure or in studies of rare states that are at very low populations. Applications to a number of protein systems are presented where the utility of the method, including its improved sensitivity to chemical exchange processes, is established.

Ligand Binding Enhances Millisecond Conformational Exchange in Xylanase B2 from Streptomyces lividans.

Xylanases catalyze the hydrolysis of xylan, an abundant carbon and energy source with important commercial ramifications. Despite tremendous efforts devoted to the catalytic improvement of xylanases, success remains limited because of our relatively poor understanding of their molecular properties. Previous reports suggested the potential role of atomic-scale residue dynamics in modulating the catalytic activity of GH11 xylanases; however, dynamics in these studies was probed on time scales orders of magnitude faster than the catalytic time frame. Here, we used nuclear magnetic resonance titration and relaxation dispersion experiments ((15)N-CPMG) in combination with X-ray crystallography and computational simulations to probe conformational motions occurring on the catalytically relevant millisecond time frame in xylanase B2 (XlnB2) and its catalytically impaired mutant E87A from Streptomyces lividans 66. Our results show distinct dynamical properties for the apo and ligand-bound states of the enzymes. The apo form of XlnB2 experiences conformational exchange for residues in the fingers and palm regions of the catalytic cleft, while the catalytically impaired E87A variant displays millisecond dynamics only in the fingers, demonstrating the long-range effect of the mutation on flexibility. Ligand binding induces enhanced conformational exchange of residues interacting with the ligand in the fingers and thumb loop regions, emphasizing the potential role of residue motions in the fingers and thumb loop regions for recognition, positioning, processivity, and/or stabilization of ligands in XlnB2. To the best of our knowledge, this work represents the first experimental characterization of millisecond dynamics in a GH11 xylanase family member. These results offer new insights into the potential role of conformational exchange in GH11 enzymes, providing essential dynamic information to help improve protein engineering and design applications.

The S. pombe mRNA decapping complex recruits cofactors and an Edc1-like activator through a single dynamic surface.

The removal of the 5′ 7-methylguanosine mRNA cap structure (decapping) is a central step in the 5′–3′ mRNA degradation pathway and is performed by the Dcp1:Dcp2 decapping complex. The activity of this complex is tightly regulated to prevent premature degradation of the transcript. Here, we establish that the aromatic groove of the EVH1 domain of Schizosaccharomyces pombe Dcp1 can interact with proline-rich sequences in the exonuclease Xrn1, the scaffolding protein Pat1, the helicase Dhh1, and the C-terminal disordered region of Dcp2. We show that this region of Dcp1 can also recruit a previously unidentified enhancer of decapping protein (Edc1) and solved the crystal structure of the complex. NMR relaxation dispersion experiments reveal that the Dcp1 binding site can adopt multiple conformations, thus providing the plasticity that is required to accommodate different ligands. We show that the activator Edc1 makes additional contacts with the regulatory domain of Dcp2 and that an activation motif in Edc1 increases the RNA affinity of Dcp1:Dcp2. Our data support a model where Edc1 stabilizes the RNA in the active site, which results in enhanced decapping rates. In summary, we show that multiple decapping factors, including the Dcp2 C-terminal region, compete with Edc1 for Dcp1 binding. Our data thus reveal a network of interactions that can fine-tune the catalytic activity of the decapping complex.

Kinetics of the Antibody Recognition Site in the Third IgG‐Binding Domain of Protein G

AbstractProtein dynamics occurring on a wide range of timescales play a crucial role in governing protein function. Particularly, motions between the globular rotational correlation time ( ) and 40 μs (supra‐ window), strongly influence molecular recognition. This supra‐ window was previously hidden, owing to a lack of experimental methods. Recently, we have developed a high‐power relaxation dispersion (RD) experiment for measuring kinetics as fast as 4 μs. For the first time, this method, performed under super‐cooled conditions, enabled us to detect a global motion in the first β‐turn of the third IgG‐binding domain of protein G (GB3), which was extrapolated to 371±115 ns at 310 K. Furthermore, the same residues show the plasticity in the model‐free residual dipolar coupling (RDC) order parameters and in an ensemble encoding the supra‐ dynamics. This β‐turn is involved in antibody binding, exhibiting the potential link of the observed supra‐ motion with molecular recognition.

Kinetics of the Antibody Recognition Site in the Third IgG-Binding Domain of Protein G.

Protein dynamics occurring on a wide range of timescales play a crucial role in governing protein function. Particularly, motions between the globular rotational correlation time (τc ) and 40 μs (supra-τc window), strongly influence molecular recognition. This supra-τc window was previously hidden, owing to a lack of experimental methods. Recently, we have developed a high-power relaxation dispersion (RD) experiment for measuring kinetics as fast as 4 μs. For the first time, this method, performed under super-cooled conditions, enabled us to detect a global motion in the first β-turn of the third IgG-binding domain of protein G (GB3), which was extrapolated to 371±115 ns at 310 K. Furthermore, the same residues show the plasticity in the model-free residual dipolar coupling (RDC) order parameters and in an ensemble encoding the supra-τc dynamics. This β-turn is involved in antibody binding, exhibiting the potential link of the observed supra-τc motion with molecular recognition.