Conformational Exchange Dynamics Research Articles

Nascent β Structure in the Elongated Hydrophobic Region of a Gerstmann–Sträussler–Scheinker PrP Allele

Prion diseases are neurodegenerative disorders caused by the misfolding of the cellular prion protein (PrPC). Gerstmann–Sträussler–Scheinker syndrome is an inherited prion disease with one early-onset allele (HRdup) containing an eight-amino-acid insertion; this LGGLGGYV insert is positioned after valine 129 (human PrPC sequence) in a hydrophobic tract in the natively disordered region. Here we have characterized the structure and explored the molecular motions and dynamics of HRdup PrP and a control allele. High-resolution NMR data suggest that the core of HRdup has a canonical PrPC structure, yet a nascent β-structure is observed in the flexible elongated hydrophobic region of HRdup. In addition, using mouse PrPC sequence, we observed that a methionine/valine polymorphism at codon 128 (equivalent of methionine/valine 129 in human sequence) and oligomerization caused by high protein concentration affects conformational exchange dynamics at residue G130. We hypothesize that with the β-structure at the N-terminus, the hydrophobic region of HRdup can adopt a fully extended configuration and fold back to form an extended β-sheet with the existing β-sheet. We propose that these structures are early chemical events in disease pathogenesis.

Controlling the dynamics of the Nek2 leucine zipper by engineering of "kinetic" disulphide bonds.

Nek2 is a dimeric serine/ threonine protein kinase that belongs to the family of NIMA-related kinases (Neks). Its N-terminal catalytic domain and its C-terminal regulatory region are bridged by a leucine zipper, which plays an important role in the activation of Nek2’s catalytic activity. Unusual conformational dynamics on the intermediary/slow timescale has thwarted all attempts so far to determine the structure of the Nek2 leucine zipper by means of X-ray crystallography and Nuclear Magnetic Resonance (NMR). Disulfide engineering, the strategic placement of non-native disulfide bonds into flexible regions flanking the coiled coil, was used to modulate the conformational exchange dynamics of this important dimerization domain. The resulting reduction in exchange rate leads to substantial improvements of important features in NMR spectra, such as line width, coherence transfer leakage and relaxation. These effects were comprehensively analyzed for the wild type protein, two single disulfide bond-bearing mutants and another double disulfide bonds-carrying mutant. Furthermore, exchange kinetics were measured across a wide temperature range, allowing for a detailed analysis of activation energy (ΔG‡) and maximal rate constant (k’ex). For one mutant carrying a disulfide bond at its C-terminus, a full backbone NMR assignment could be obtained for both conformers, demonstrating the benefits of the disulfide engineering. Our study demonstrates the first successful application of ‘kinetic’ disulfide bonds for the purpose of controlling the adverse effects of protein dynamics. Firstly, this provides a promising, robust platform for the full structural and functional investigation of the Nek2 leucine zipper in the future. Secondly, this work broadens the toolbox of protein engineering by disulfide bonds through the addition of a kinetic option in addition to the well-established thermodynamic uses of disulfide bonds.

Probing Conformational Exchange Dynamics in a Short-Lived Protein Folding Intermediate by Real-Time Relaxation-Dispersion NMR.

NMR spectroscopy is a powerful tool for studying molecular dynamics at atomic resolution simultaneously for a large number of nuclear sites. In this communication, we combine two powerful NMR techniques, relaxation-dispersion NMR and real-time NMR, in order to obtain unprecedented information on the conformational exchange dynamics present in short-lived excited protein states, such as those transiently accumulated during protein folding. We demonstrate the feasibility of the approach for the amyloidogenic protein β2-microglobulin that folds via an intermediate state which is believed to be responsible for the onset of the aggregation process leading to amyloid formation.

Protein Allostery and Conformational Dynamics.

The functions of many proteins are regulated through allostery, whereby effector binding at a distal site changes the functional activity (e.g., substrate binding affinity or catalytic efficiency) at the active site. Most allosteric studies have focused on thermodynamic properties, in particular, substrate binding affinity. Changes in substrate binding affinity by allosteric effectors have generally been thought to be mediated by conformational transitions of the proteins or, alternatively, by changes in the broadness of the free energy basin of the protein conformational state without shifting the basin minimum position. When effector binding changes the free energy landscape of a protein in conformational space, the change affects not only thermodynamic properties but also dynamic properties, including the amplitudes of motions on different time scales and rates of conformational transitions. Here we assess the roles of conformational dynamics in allosteric regulation. Two cases are highlighted where NMR spectroscopy and molecular dynamics simulation have been used as complementary approaches to identify residues possibly involved in allosteric communication. Perspectives on contentious issues, for example, the relationship between picosecond-nanosecond local and microsecond-millisecond conformational exchange dynamics, are presented.

Heteronuclear Adiabatic Relaxation Dispersion (HARD) for quantitative analysis of conformational dynamics in proteins

NMR relaxation methods probe biomolecular motions over a wide range of timescales. In particular, the rotating frame spin-lock R1ρ and Carr–Purcell–Meiboom–Gill (CPMG) R2 experiments are commonly used to characterize μs to ms dynamics, which play a critical role in enzyme folding and catalysis. In an effort to complement these approaches, we introduced the Heteronuclear Adiabatic Relaxation Dispersion (HARD) method, where dispersion in rotating frame relaxation rate constants (longitudinal R1ρ and transverse R2ρ) is created by modulating the shape and duration of adiabatic full passage (AFP) pulses. Previously, we showed the ability of the HARD method to detect chemical exchange dynamics in the fast exchange regime (kex∼104–105s−1). In this article, we show the sensitivity of the HARD method to slower exchange processes by measuring R1ρ and R2ρ relaxation rates for two soluble proteins (ubiquitin and 10C RNA ligase). One advantage of the HARD method is its nominal dependence on the applied radio frequency field, which can be leveraged to modulate the dispersion in the relaxation rate constants. In addition, we also include product operator simulations to define the dynamic range of adiabatic R1ρ and R2ρ that is valid under all exchange regimes. We conclude from both experimental observations and simulations that this method is complementary to CPMG-based and rotating frame spin-lock R1ρ experiments to probe conformational exchange dynamics for biomolecules. Finally, this approach is germane to several NMR-active nuclei, where relaxation rates are frequency-offset independent.

Functional Dynamics in the Enzyme Adenylate Kinase

An important question in biology is how the energy landscape of enzymes can enable efficient catalysis of chemical reactions. We have undertaken an in-depth analysis of conformational exchange dynamics in the enzyme adenylate kinase (adk). Adk is an essential enzyme in higher organisms that catalyzes the reversible interconversion of AMP and ATP into two ADP molecules. The enzyme is modular and is composed of distinct ATP and AMP binding subdomains and in addition a core subdomain (where core is responsible for global thermodynamic and thermal stability). We have previously shown that substrate binding is accompanied by rate-limiting spatial displacements of both the ATP and AMP binding motifs. Here, we present a solution state chemical shift based NMR approach to probe the native energy landscape of adenylate kinase with and without its natural substrates present. Binding of ATP induces a dynamic equilibrium in which the ATP binding subdomain populates both open and closed conformations with equal weights. A similar scenario is observed upon AMP binding in the AMP binding subdomain. These structural ensembles represent complexes that are populated transiently during the enzymatic reaction cycle. Our proposed dynamic mode of protein-ligand interaction in adk stands in contrast to the traditional view of substrate/enzyme complexes as rigid, low entropy states. Finally, by using a combination of protein engineering and hydrogen to deuterium exchange experiments we have shown that the individual subdomains in adk can fold independently of each other. Independent folding of subdomains can, in principle, be utilized to accommodate the structural change during the functional open to closed transition.

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.

Self-Threading and Dethreading Dynamics of Poly(ethylene glycol)-Substituted Cyclodextrins with Different Chain Lengths

Poly(ethylene glycol) (PEG)-substituted cyclodextrins (CDs) with different chain lengths have been synthesized. PEG-substituted CDs formed self-threading complexes in aqueous solutions. The conformation of PEG-substituted CDs changed from a self-threading form to a dethreading form in the presence of an equal molar amount of 1-adamantanecarboxylic acid (AdCA) as a competitive guest and exchanged between a self-threading form and a dethreading form in the presence of a half molar amount of AdCA. The conformational exchange dynamics of PEG-substituted CDs were studied by means of variable-temperature NMR and 2D EXSY NMR experiments, and it was revealed that the conformational exchange rates decreased exponentially with the PEG chain length of the substituent.

Self-Threading of a Poly(ethylene glycol) Chain in a Cyclodextrin-Ring: Control of the Exchange Dynamics by Chain Length

Poly(ethylene glycol) (PEG)-substituted cyclodextrins (CDs) with different chain lengths have been synthesized. PEG-substituted CDs formed self-threading complexes in aqueous solutions, and the conformational exchange dynamics between self-threading and dethreading could be regulated by its chain length.

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.

Quantitative analysis of conformational exchange contributions to 1H-15N multiple-quantum relaxation using field-dependent measurements. Time scale and structural characterization of exchange in a calmodulin C-terminal domain mutant.

Multiple-quantum spin relaxation is a sensitive probe for correlated conformational exchange dynamics on microsecond to millisecond time scales in biomolecules. We measured differential 1H-15N multiple-quantum relaxation rates for the backbone amide groups of the E140Q mutant of the C-terminal domain of calmodulin at three static magnetic field strengths. The differential multiple-quantum relaxation rates range between -88.7 and 92.7 s(-1), and the mean and standard deviation are 7.0 +/- 24 s(-1), at a static magnetic field strength of 14.1 T. Together with values of the 1H and 15N chemical shift anisotropies (CSA) determined separately, the field-dependent data enable separation of the different contributions from dipolar-dipolar, CSA-CSA, and conformational exchange cross-correlated relaxation mechanisms to the differential multiple-quantum relaxation rates. The procedure yields precise quantitative information on the dominant conformational exchange contributions observed in this protein. The field-dependent differences between double- and zero-quantum relaxation rates directly benchmark the rates of conformational exchange, showing that these are fast on the chemical shift time scale for the large majority of residues in the protein. Further analysis of the differential 1H-15N multiple-quantum relaxation rates using previously determined exchange rate constants and populations, obtained from 15N off-resonance rotating-frame relaxation data, enables extraction of the product of the chemical shift differences between the resonance frequencies of the 1H and 15N spins in the exchanging conformations, deltasigma(H)deltasigma(N). Thus, information on the 1H chemical shift differences is obtained, while circumventing complications associated with direct measurements of conformational exchange effects on 1H single-quantum coherences in nondeuterated proteins. The method significantly increases the information content available for structural interpretation of the conformational exchange process, partly because deltasigma(H)deltasigma(N) is a signed quantity, and partly because two chemical shifts are probed simultaneously. The present results support the hypothesis that the exchange in the calcium-loaded state of the E140Q mutant involves conformations similar to those of the wild-type apo (closed) and calcium-loaded (open) states.

Dynamics of the Transition between Open and Closed Conformations in a Calmodulin C-Terminal Domain Mutant

Background: Calmodulin is a ubiquitous Ca 2+-activated regulator of cellular processes in eukaryotes. The structures of the Ca 2+-free (apo) and Ca 2+-loaded states of calmodulin have revealed that Ca 2+ binding is associated with a transition in each of the two domains from a closed to an open conformation that is central to target recognition. However, little is known about the dynamics of this conformational switch. Results: The dynamics of the transition between closed and open conformations in the Ca 2+-loaded state of the E140Q mutant of the calmodulin C-terminal domain were characterized under equilibrium conditions. The exchange time constants (τ ex) measured for 42 residues range from 13 to 46 μs, with a mean of 21 ± 3 μs. The results suggest that τ ex varies significantly between different groups of residues and that residues with similar values exhibit spatial proximity in the structures of apo and/or Ca 2+-saturated wild-type calmodulin. Using data for one of these groups, we obtained an open population of p o = 0.50 ± 0.17 and a closed → open rate constant of k o = 2.7 ± 1.0 × 10 4 s −1 . Conclusions: The conformational exchange dynamics appear to involve locally collective processes that depend on the structural topology. Comparisons with previous results indicate that similar processes occur in the wild-type protein. The measured rates match the estimated Ca 2+ off rate, suggesting that Ca 2+ release may be gated by the conformational dynamics. Structural interpretation of estimated chemical shifts suggests a mechanism for ion release.

The influence of DNA binding on the backbone dynamics of the yeast cell-cycle protein Mbp1.

Mbp1 is a transcription factor involved in the regulation of the cell cycle in yeast. The N-terminus of this protein contains a DNA binding domain that includes a winged helix-turn-helix motif. The C-terminal 24 residues of this domain (the 'tail') are disordered in the crystal state, but are important for DNA binding. We have measured 15N NMR relaxation rates at 11.75 and 14.1 T to determine the dynamics of the free protein and in its complex with a specific DNA duplex. The dynamics data were quantitatively analysed using both spectral density mapping and the Lipari-Szabo formalism including the effects of chemical exchange and rotational anisotropy. A detailed analysis has been made of the effect of anisotropy, exchange and experimental precision on the recovered motional parameters. The backbone NH relaxation is affected by motions on a variety of time scales from millisecond to tens of picoseconds. The relaxation data show a structured core of 100 residues corresponding to that observed in the crystal state. Within the core of the protein, two regions on either side of the putative recognition helix (helix B) show slow (ca. 0.2 ms) conformational exchange dynamics that are quenched upon DNA binding. The C-terminal 24 residues are generally more dynamic than in the core. However, in the free protein, a stretch of approximately 8 residues in the middle of the tail show relaxation behaviour similar to that in the core, indicating a structured region. NOEs between Ala 114 in this structured part of the tail and residues in the N-terminal beta strand of the core of the protein demonstrate that the tail folds back onto the core of the protein. In the complex with DNA, the structured part of the tail extends by ca. 3 residues. These data provide a framework for understanding the biochemical data on the mechanism and specificity of DNA binding.

Conformational Dynamics and Molecular Recognition: Backbone Dynamics of the Estrogen Receptor DNA-binding Domain

We examined the internal mobility of the estrogen receptor DNA-binding domain (ER DBD) using NMR15N relaxation measurements and compared it to that of the glucocorticoid receptor DNA-binding domain (GR DBD). The studied protein fragments consist of residues Arg183-His267 of the human ER and residues Lys438-Gln520 of the rat GR. The15N longitudinal (R1) and transverse (R2) relaxation rates and steady state {1H}-15N nuclear Overhauser enhancements (NOEs) were measured at 30°C at1H NMR frequencies of 500 and 600 MHz. The NOE versus sequence profile and calculated order parameters for ER DBD backbone motions indicate enhanced internal dynamics on pico- to nanosecond time-scales in two regions of the core DBD. These are the extended strand which links the DNA recognition helix to the second zinc domain and the larger loop region of the second zinc domain. The mobility of the corresponding regions of the GR DBD, in particular that of the second zinc domain, is more limited. In addition, we find large differences between the ER and GR DBDs in the extent of conformational exchange mobility on micro- to millisecond time-scales. Based on measurements of R2as a function of the15N refocusing (CPMG) delay and quantitative (Lipari-Szabo-type) analysis, we conclude that conformational exchange occurs in the loop of the first zinc domain and throughout most of the second zinc domain of the ER DBD. The conformational exchange dynamics in GR DBD is less extensive and localized to two sites in the second zinc domain. The different dynamical features seen in the two proteins is consistent with previous studies of the free state structures in which the second zinc domain in the ER DBD was concluded to be disordered whereas the corresponding region of the GR DBD adopts a stable fold. Moreover, the regions of the ER DBD that undergo conformational dynamics on the micro- to millisecond time-scales in the free state are involved in intermolecular protein-DNA and protein-protein interactions in the dimeric bound state. Based on the present data and the previously published dynamical and DNA binding properties of a GR DBD triple mutant which recognize an ER binding site on DNA, we argue that the free state dynamical properties of the nuclear receptor DBDs is an important element in molecular recognition upon DNA binding.

Dynamics study on the anti-human immunodeficiency virus chemokine viral macrophage-inflammatory protein-II (VMIP-II) reveals a fully monomeric protein.

Encoded by Kaposi's sarcoma-associated herpesvirus, viral macrophage-inflammatory protein-II (VMIP-II) is unique among CC chemokines in that it has been shown to bind to the CXC chemokine receptor CXCR4 as well as to a variety of CC chemokine receptors. This unique binding ability allows vMIP-II to block infection by a wide range of human immunodeficiency virus type I (HIV-1) strains, but the structural and dynamic basis for this broad range of binding is not known. 15N T1, T2 and 15N[-HN] nuclear Overhauser effect (NOE) values of vMIP-II, determined through a series of heteronuclear multidimensional nuclear magnetic resonance (NMR) experiments, were used to obtain information about the backbone dynamics of the protein. Whereas almost all chemokine structures reveal a dimer or multimer, vMIP-II has a rotational correlation time (tauc) of 4.7 +/- 0.3 ns, which is consistent with a monomeric chemokine. The rotational diffusion anisotropy, D parallel/D perpendicular, is approximately 1.5 +/- 0.1. The conformation of vMIP-II is quite similar to other known chemokines, containing an unstructured N-terminus followed by an ordered turn, three beta-strands arranged in an antiparallel fashion, and one C-terminal alpha-helix that lies across the beta-strands. Most of the protein is well-ordered on a picosecond time scale, with an average order parameter S2 (excluding the N-terminal 13 amino acids) of 0.83 +/- 0. 09, and with even greater order in regions of secondary structure. The NMR data reveal that the N-terminus, which in other chemokines has been implicated in receptor binding, extends like a flexible tail in solution and possesses no secondary structure. The region of the ordered turn, including residues 25-28, experiences conformational exchange dynamics. The implications of these NMR data to the broad receptor binding capability of vMIP-II are discussed.