Selective characterization of microsecond motions in proteins by NMR relaxation

A unique chaperoning mechanism in class A JDPs recognizes and stabilizes mutant p53

J-domain proteins (JDPs) constitute a large family of molecular chaperones that bind a broad spectrum of substrates, targeting them to Hsp70, thus determining the specificity of and activating the entire chaperone functional cycle. The malfunction of JDPs is therefore inextricably linked to myriad human disorders. Here, we uncover a unique mechanism by which chaperones recognize misfolded clients, present in human class A JDPs. Through a newly identified β-hairpin site, these chaperones detect changes in protein dynamics at the initial stages of misfolding, prior to exposure of hydrophobic regions or large structural rearrangements. The JDPs then sequester misfolding-prone proteins into large oligomeric assemblies, protecting them from aggregation. Through this mechanism, class A JDPs bind destabilized p53 mutants, preventing clearance of these oncoproteins by Hsp70-mediated degradation, thus promoting cancer progression. Removal of the β-hairpin abrogates this protective activity while minimally affecting other chaperoning functions. This suggests the class A JDP β-hairpin as a highly specific target for cancer therapeutics.

Optimization of TROSY- and anti-TROSY-based 15N CPMG relaxation dispersion experiments through phase cycling

CPMG relaxation dispersion studies of biomolecular dynamics on the μs-ms timescale can provide detailed kinetic, thermodynamic, and structural insights into function. Frequently, the 15N spin serves as the probe of choice, as uniform incorporation of the 15N isotope is facile and cost-effective, and the interpretation of the resulting data is often relatively straightforward. In conventional CPMG relaxation dispersion experiments the application of CPMG pulses with constant radiofrequency (RF) phase can lead to artifactual dispersion profiles that result from off-resonance effects, RF field inhomogeneity, and pulse miscalibration. The development of CPMG experiments with the [0013]-phase cycle has significantly reduced the impact of pulse imperfections over a greater bandwidth of frequency offsets in comparison to constant phase experiments. Application of 15N-TROSY-based CPMG schemes to studies of the dynamics of large molecules is necessary for high sensitivity, yet the correct incorporation of the [0013]-phase cycle is non-trivial. Here we present TROSY- and anti-TROSY-based 15N CPMG experiments with the [0013]-phase cycling scheme and demonstrate, through comprehensive numerical simulations and experimental validation, enhanced resistance to pulse imperfections relative to traditional schemes utilizing constant phase CPMG pulses. Notably, exchange parameters derived from the new experiments are in good agreement with those obtained using other, more established, 15N-based CPMG approaches.

Poly(ADP-ribosyl)ation enhances nucleosome dynamics and organizes DNA damage repair components within biomolecular condensates

Nucleosomes, the basic structural units of chromatin, hinder recruitment and activity of various DNA repair proteins, necessitating modifications that enhance DNA accessibility. Poly(ADP-ribosyl)ation (PARylation) of proteins near damage sites is an essential initiation step in several DNA-repair pathways; however, its effects on nucleosome structural dynamics and organization are unclear. Using NMR, cryoelectron microscopy (cryo-EM), and biochemical assays, we show that PARylation enhances motions of the histone H3 tail and DNA, leaving the configuration of the core intact while also stimulating nuclease digestion and ligation of nicked nucleosomal DNA by LIG3. PARylation disrupted interactions between nucleosomes, preventing self-association. Addition of LIG3 and XRCC1 to PARylated nucleosomes generated condensates that selectively partition DNA repair-associated proteins in a PAR- and phosphorylation-dependent manner in vitro. Our results establish that PARylation influences nucleosomes across different length scales, extending from the atom-level motions of histone tails to the mesoscale formation of condensates with selective compositions.

Correlating histone acetylation with nucleosome core particle dynamics and function

Epigenetic modifications of chromatin play a critical role in regulating the fidelity of the genetic code and in controlling the translation of genetic information into the protein components of the cell. One key posttranslational modification is acetylation of histone lysine residues. Molecular dynamics simulations, and to a smaller extent experiment, have established that lysine acetylation increases the dynamics of histone tails. However, a systematic, atomic resolution experimental investigation of how this epigenetic mark, focusing on one histone at a time, influences the structural dynamics of the nucleosome beyond the tails, and how this translates into accessibility of protein factors such as ligases and nucleases, has yet to be performed. Herein, using NMR spectroscopy of nucleosome core particles (NCPs), we evaluate the effects of acetylation of each histone on tail and core dynamics. We show that for histones H2B, H3, and H4, the histone core particle dynamics are little changed, even though the tails have increased amplitude motions. In contrast, significant increases to H2A dynamics are observed upon acetylation of this histone, with the docking domain and L1 loop particularly affected, correlating with increased susceptibility of NCPs to nuclease digestion and more robust ligation of nicked DNA. Dynamic light scattering experiments establish that acetylation decreases inter-NCP interactions in a histone-dependent manner and facilitates the development of a thermodynamic model for NCP stacking. Our data show that different acetylation patterns result in nuanced changes to NCP dynamics, modulating interactions with other protein factors, and ultimately controlling biological output.

Biophysical and Integrative Characterization of Protein Intrinsic Disorder as a Prime Target for Drug Discovery

Protein intrinsic disorder is increasingly recognized for its biological and disease-driven functions. However, it represents significant challenges for biophysical studies due to its high conformational flexibility. In addressing these challenges, we highlight the complementary and distinct capabilities of a range of experimental and computational methods and further describe integrative strategies available for combining these techniques. Integrative biophysics methods provide valuable insights into the sequence–structure–function relationship of disordered proteins, setting the stage for protein intrinsic disorder to become a promising target for drug discovery. Finally, we briefly summarize recent advances in the development of new small molecule inhibitors targeting the disordered N-terminal domains of three vital transcription factors.

Surface electrostatics dictate RNA-binding protein CAPRIN1 condensate concentration and hydrodynamic properties

Biomolecular condensates concentrate proteins, nucleic acids, and small molecules and play an essential role in many biological processes. Their formation is tuned by a balance between energetically favorable and unfavorable contacts, with charge-charge interactions playing a central role in some systems. The positively charged intrinsically disordered carboxy-terminal region of the RNA-binding protein CAPRIN1 is one such example, phase separating upon addition of negatively charged ATP or high concentrations of sodium chloride (NaCl). Using solution NMR spectroscopy, we measured residue-specific near-surface electrostatic potentials (ϕENS) of CAPRIN1 along its NaCl-induced phase separation trajectory to compare with those obtained using ATP. In both cases, electrostatic shielding decreases ϕENS values, yet surface potentials of CAPRIN1 in the two condensates can be different, depending on the amount of NaCl or ATP added. Our results establish that even small differences in ϕENS can significantly affect the level of protein enrichment and the mechanical properties of the condensed phase, leading, potentially, to the regulation of biological processes.

Salt Dependence Conformational Stability of the Dimeric SAM Domain of MAPKKK Ste11 from Budding Yeast: A Native-State H/D Exchange NMR Study

The sterile α motif, also called the SAM domain, is known to form homo or heterocomplexes that modulate diverse biological functions through the regulation of specific protein-protein interactions. The MAPK pathway of budding yeast Saccharomyces cerevisiae is comprised of a three-tier kinase system akin to mammals. The MAPKKK Ste11 protein of yeast contains a homodimer SAM domain, which is critical for transmitting cues to the downstream kinases. The structural stability of the dimeric Ste11 SAM is maintained by hydrophobic and ionic interactions at the interfacial amino acids. The urea-induced equilibrium-unfolding process of the Ste11 SAM domain is cooperative without evidence of any intermediate states. The native-state H/D exchange under subdenaturing conditions is a useful method for the detection of intermediate states of proteins. In the present study, we investigated the effect of ionic strength on the conformational stability of the dimer using the H/D exchange experiments. The hydrogen exchange behavior of the Ste11 dimer under physiological salt concentrations reveals two partially unfolded metastable intermediate states, which may be generated by a sequential and cooperative unfolding of the five helices present in the domain. These intermediates appear to be significant for the reversible unfolding kinetics via hydrophobic collapse. In contrast, higher ionic concentrations eliminate this cooperative interactions that stabilize the pairs of helices.

Protein Dynamics revealed by NMR Relaxation Methods

Structural biology often focuses primarily on three-dimensional structures of biological macromolecules, deposited in the Protein Data Bank (PDB). This resource is a remarkable entity for the world-wide scientific and medical communities, as well as the general public, as it is a growing translation into three-dimensional space of the vast information in genomic databases, e.g. GENBANK. There is, however, significantly more to understanding biological function than the three-dimensional coordinate space for ground-state structures of biomolecules. The vast array of biomolecules experiences natural dynamics, interconversion between multiple conformational states, and molecular recognition and allosteric events that play out on timescales ranging from picoseconds to seconds. This wide range of timescales demands ingenious and sophisticated experimental tools to sample and interpret these motions, thus enabling clearer insight into functional annotation of the PDB. NMR spectroscopy is unique in its ability to sample this range of timescales at atomic resolution and in physiologically relevant conditions using spin relaxation methods. The field is constantly expanding to provide new creative experiments, to yield more detailed coverage of timescales, and to broaden the power of interpretation and analysis methods. This review highlights the current state of the methodology and examines the extension of analysis tools for more complex experiments and dynamic models. The future for understanding protein dynamics is bright, and these extended tools bring greater compatibility with developments in computational molecular dynamics, all of which will further our understanding of biological molecular functions. These facets place NMR as a key component in integrated structural biology.

Detection of chemical exchange in methyl groups of macromolecules

The zero- and double-quantum methyl TROSY Hahn-echo and the methyl ¹H-¹H dipole-dipole cross-correlation nuclear magnetic resonance experiments enable estimation of multiple quantum chemical exchange broadening in methyl groups in proteins. The two relaxation rate constants are established to be linearly dependent using molecular dynamics simulations and empirical analysis of experimental data. This relationship allows chemical exchange broadening to be recognized as an increase in the Hahn-echo relaxation rate constant. The approach is illustrated by analyzing relaxation data collected at three temperatures for E. coli ribonuclease HI and by analyzing relaxation data collected for different cofactor and substrate complexes of E. coli AlkB.

Reducing bias in the analysis of solution-state NMR data with dynamics detectors

Nuclear magnetic resonance (NMR) is sensitive to dynamics on a wide range of correlation times. Recently, we have shown that analysis of relaxation rates via fitting to a correlation function with a small number of exponential terms could yield a biased characterization of molecular motion in solid-state NMR due to limited sensitivity of experimental data to certain ranges of correlation times. We introduced an alternative approach based on "detectors" in solid-state NMR, for which detector responses characterize motion for a range of correlation times and reduce potential bias resulting from the use of simple models for the motional correlation functions. Here, we show that similar bias can occur in the analysis of solution-state NMR relaxation data. We have thus adapted the detector approach to solution-state NMR, specifically separating overall tumbling motion from internal motions and accounting for contributions of chemical exchange to transverse relaxation. We demonstrate that internal protein motions can be described with detectors when the overall motion and the internal motions are statistically independent. We illustrate the detector analysis on ubiquitin with typical relaxation data sets recorded at a single high magnetic field or at multiple high magnetic fields and compare with results of model-free analysis. We also compare our methodology to LeMaster's method of dynamics analysis.

Microsecond Protein Dynamics from Combined Bloch-McConnell and Near-Rotary-Resonance R1p Relaxation-Dispersion MAS NMR

Studying protein dynamics on microsecond-to-millisecond (μs-ms) time scales can provide important insight into protein function. In magic-angle-spinning (MAS) NMR, μs dynamics can be visualized by R 1 ρ rotating-frame relaxation dispersion experiments in different regimes of radio-frequency field strengths: at low RF field strength, isotropic-chemical-shift fluctuation leads to "Bloch-McConnell-type" relaxation dispersion, while when the RF field approaches rotary resonance conditions bond angle fluctuations manifest as increased R 1 ρ rate constants ("Near-Rotary-Resonance Relaxation Dispersion", NERRD). Here we explore the joint analysis of both regimes to gain comprehensive insight into motion in terms of geometric amplitudes, chemical-shift changes, populations and exchange kinetics. We use a numerical simulation procedure to illustrate these effects and the potential of extracting exchange parameters, and apply the methodology to the study of a previously described conformational exchange process in microcrystalline ubiquitin.

Why the Energy Landscape of Barnase Is Hierarchical

We have used NMR and computational methods to characterize the dynamics of the ribonuclease barnase over a wide range of timescales in free and inhibitor-bound states. Using temperature- and denaturant-dependent measurements of chemical shift, we show that barnase undergoes frequent and highly populated hinge bending. Using relaxation dispersion, we characterize a slower and less populated motion with a rate of 750 ± 200 s⁻¹, involving residues around the lip of the active site, which occurs in both free and bound states and therefore suggests conformational selection. Normal mode calculations characterize correlated hinge bending motions on a very rapid timescale. These three measurements are combined with previous measurements and molecular dynamics calculations on barnase to characterize its dynamic landscape on timescales from picoseconds to milliseconds and length scales from 0.1 to 2.5 nm. We show that barnase has two different large-scale fluctuations: one on a timescale of 10⁻⁹−10⁻⁶ s that has no free energy barrier and is a hinge bending that is determined by the architecture of the protein; and one on a timescale of milliseconds (i.e., 750 s⁻¹) that has a significant free energy barrier and starts from a partially hinge-bent conformation. These two motions can be described as hierarchical, in that the more highly populated faster motion provides a platform for the slower (less probable) motion. The implications are discussed. The use of temperature and denaturant is suggested as a simple and general way to characterize motions on the intermediate ns-μs timescale.

Nucleosome structure and dynamics are coming of age

Since the first high-resolution structure of the nucleosome was reported in 1997, the available information on chromatin structure has increased very rapidly. Here, we review insights derived from cutting-edge biophysical and structural approaches applied to the study of nucleosome dynamics and nucleosome-binding factors, with a focus on the experimental advances driving the research. In addition, we highlight emerging challenges in nucleosome structural biology.

Exploring methods to expedite the recording of CEST datasets using selective pulse excitation

Chemical Exchange Saturation Transfer (CEST) has emerged as a powerful tool for studies of biomolecular conformational exchange involving the interconversion between a major, visible conformer and one or more minor, invisible states. Applications typically entail recording a large number of 2D datasets, each of which differs in the position of a weak radio frequency field, so as to generate a CEST profile for each nucleus from which the chemical shifts of spins in the invisible state(s) are obtained. Here we compare a number of band-selective CEST schemes for speeding up the process using either DANTE or cosine-modulated excitation approaches. We show that while both are essentially identical for applications such as ¹⁵N CEST, in cases where the probed spins are dipolar or scalar coupled to other like spins there can be advantages for the cosine-excitation scheme.

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.

X-ray crystallography is the main method for protein structure determination. Here the authors combine solid-state NMR measurements and molecular dynamics simulations and show that crystal packing alters the thermodynamics and kinetics of local conformational exchange as well as overall rocking motion of protein molecules in the crystal lattice.

Measurement of 15N longitudinal relaxation rates in 15NH4+ spin systems to characterise rotational correlation times and chemical exchange

Many chemical and biological processes rely on the movement of monovalent cations and an understanding of such processes can therefore only be achieved by characterising the dynamics of the involved ions. It has recently been shown that 15N-ammonium can be used as a proxy for potassium to probe potassium binding in bio-molecules such as DNA quadruplexes and enzymes. Moreover, equations have been derived to describe the time-evolution of 15N-based spin density operator elements of 15NH4+ spin systems. Herein NMR pulse sequences are derived to select specific spin density matrix elements of the 15NH4+ spin system and to measure their longitudinal relaxation in order to characterise the rotational correlation time of the 15NH4+ ion as well as report on chemical exchange events of the 15NH4+ ion. Applications to 15NH4+ in acidic aqueous solutions are used to cross-validate the developed pulse sequence while measurements of spin-relaxation rates of 15NH4+ bound to a 41kDa domain of the bacterial Hsp70 homologue DnaK are presented to show the general applicability of the derived pulse sequence. The rotational correlation time obtained for 15N-ammonium bound to DnaK is similar to the correlation time that describes the rotation about the threefold axis of a methyl group. The methodology presented here provides, together with the previous theoretical framework, an important step towards characterising the motional properties of cations in macromolecular systems.

Protein dynamics and function from solution state NMR spectroscopy

It is well-established that dynamics are central to protein function; their importance is implicitly acknowledged in the principles of the Monod, Wyman and Changeux model of binding cooperativity, which was originally proposed in 1965. Nowadays the concept of protein dynamics is formulated in terms of the energy landscape theory, which can be used to understand protein folding and conformational changes in proteins. Because protein dynamics are so important, a key to understanding protein function at the molecular level is to design experiments that allow their quantitative analysis. Nuclear magnetic resonance (NMR) spectroscopy is uniquely suited for this purpose because major advances in theory, hardware, and experimental methods have made it possible to characterize protein dynamics at an unprecedented level of detail. Unique features of NMR include the ability to quantify dynamics (i) under equilibrium conditions without external perturbations, (ii) using many probes simultaneously, and (iii) over large time intervals. Here we review NMR techniques for quantifying protein dynamics on fast (ps-ns), slow (μs-ms), and very slow (s-min) time scales. These techniques are discussed with reference to some major discoveries in protein science that have been made possible by NMR spectroscopy.

Studying Dynamics by Magic-Angle Spinning Solid-State NMR Spectroscopy: Principles and Applications to Biomolecules

Magic-angle spinning solid-state NMR spectroscopy is an important technique to study molecular structure, dynamics and interactions, and is rapidly gaining importance in biomolecular sciences. Here we provide an overview of experimental approaches to study molecular dynamics by MAS solid-state NMR, with an emphasis on the underlying theoretical concepts and differences of MAS solid-state NMR compared to solution-state NMR. The theoretical foundations of nuclear spin relaxation are revisited, focusing on the particularities of spin relaxation in solid samples under magic-angle spinning. We discuss the range of validity of Redfield theory, as well as the inherent multi-exponential behavior of relaxation in solids. Experimental challenges for measuring relaxation parameters in MAS solid-state NMR and a few recently proposed relaxation approaches are discussed, which provide information about time scales and amplitudes of motions ranging from picoseconds to milliseconds. We also discuss the theoretical basis and experimental measurements of anisotropic interactions (chemical-shift anisotropies, dipolar and quadrupolar couplings), which give direct information about the amplitude of motions. The potential of combining relaxation data with such measurements of dynamically-averaged anisotropic interactions is discussed. Although the focus of this review is on the theoretical foundations of dynamics studies rather than their application, we close by discussing a small number of recent dynamics studies, where the dynamic properties of proteins in crystals are compared to those in solution.

Heteronuclear transverse and longitudinal relaxation in AX4 spin systems: application to (15)N relaxations in (15)NH4(+)

The equations that describe the time-evolution of transverse and longitudinal (15)N magnetisations in tetrahedral ammonium ions, (15)NH4(+), are derived from the Bloch-Wangsness-Redfield density operator relaxation theory. It is assumed that the relaxation of the spin-states is dominated by (1) the intra-molecular (15)N-(1)H and (1)H-(1)H dipole-dipole interactions and (2) interactions of the ammonium protons with remote spins, which also include the contribution to the relaxations that arise from the exchange of the ammonium protons with the bulk solvent. The dipole-dipole cross-correlated relaxation mechanisms between each of the (15)N-(1)H and (1)H-(1)H interactions are explicitly taken into account in the derivations. An application to (15)N-ammonium bound to a 41kDa domain of the protein DnaK is presented, where a comparison between experiments and simulations show that the ammonium ion rotates rapidly within its binding site with a local correlation time shorter than approximately 1ns. The theoretical framework provided here forms the basis for further investigations of dynamics of AX4 spin systems, with ammonium ions in solution and bound to proteins of particular interest.

NMR approaches in structure-based lead discovery: recent developments and new frontiers for targeting multi-protein complexes

Nuclear magnetic resonance (NMR) spectroscopy is a pivotal method for structure-based and fragment-based lead discovery because it is one of the most robust techniques to provide information on protein structure, dynamics and interaction at an atomic level in solution. Nowadays, in most ligand screening cascades, NMR-based methods are applied to identify and structurally validate small molecule binding. These can be high-throughput and are often used synergistically with other biophysical assays. Here, we describe current state-of-the-art in the portfolio of available NMR-based experiments that are used to aid early-stage lead discovery. We then focus on multi-protein complexes as targets and how NMR spectroscopy allows studying of interactions within the high molecular weight assemblies that make up a vast fraction of the yet untargeted proteome. Finally, we give our perspective on how currently available methods could build an improved strategy for drug discovery against such challenging targets.

Characterization of the near native conformational states of the SAM domain of Ste11 protein by NMR spectroscopy

The sterile alpha motif or SAM domain is one of the most frequently present protein interaction modules with diverse functional attributions. SAM domain of the Ste11 protein of budding yeast plays important roles in mitogen-activated protein kinase cascades. In the current study, urea-induced, at subdenaturing concentrations, structural, and dynamical changes in the Ste11 SAM domain have been investigated by nuclear magnetic resonance spectroscopy. Our study revealed that a number of residues from Helix 1 and Helix 5 of the Ste11 SAM domain display plausible alternate conformational states and largest chemical shift perturbations at low urea concentrations. Amide proton (H/D) exchange experiments indicated that Helix 1, loop, and Helix 5 become more susceptible to solvent exchange with increased concentrations of urea. Notably, Helix 1 and Helix 5 are directly involved in binding interactions of the Ste11 SAM domain. Our data further demonstrate that the existence of alternate conformational states around the regions involved in dimeric interactions in native or near native conditions.

Proton-decoupled CPMG: a better experiment for measuring (15)N R2 relaxation in disordered proteins

(15)N R2 relaxation is one of the most informative experiments for characterization of intrinsically disordered proteins (IDPs). Small changes in nitrogen R2 rates are often used to determine how IDPs respond to various biologically relevant perturbations such as point mutations, posttranslational modifications and weak ligand interactions. However collecting high-quality (15)N relaxation data can be difficult. Of necessity, the samples of IDPs are often prepared with low protein concentration and the measurement time can be limited because of rapid sample degradation. Furthermore, due to hardware limitations standard experiments such as (15)N spin-lock and CPMG can sample the relaxation decay only to ca. 150ms. This is much shorter than (15)N T2 times in disordered proteins at or near physiological temperature. As a result, the sampling of relaxation decay profiles in these experiments is suboptimal, which further lowers the precision of the measurements. Here we report a new implementation of the proton-decoupled (PD) CPMG experiment which allows one to sample (15)N R2 relaxation decay up to ca. 0.5-1s. The new experiment has been validated through comparison with the well-established spin-lock measurement. Using dilute samples of denatured ubiquitin, we have demonstrated that PD-CPMG produces up to 3-fold improvement in the precision of the data. It is expected that for intrinsically disordered proteins the gains may be even more substantial. We have also shown that this sequence has a number of favorable properties: (i) the spectra are recorded with narrow linewidth in nitrogen dimension; (ii) (15)N offset correction is small and easy to calculate; (iii) the experiment is immune to various spurious effects arising from solvent exchange; (iv) the results are stable with respect to pulse miscalibration and rf field inhomogeneity; (v) with minimal change, the pulse sequence can also be used to measure R2 relaxation of (15)N(ε) spins in arginine side chains. We anticipate that the new experiment will be a valuable addition to the NMR toolbox for studies of IDPs.

Probing transient conformational states of proteins by solid-state R(1ρ) relaxation-dispersion NMR spectroscopy

The function of proteins depends on their ability to sample a variety of states differing in structure and free energy. Deciphering how the various thermally accessible conformations are connected, and understanding their structures and relative energies is crucial in rationalizing protein function. Many biomolecular reactions take place within microseconds to milliseconds, and this timescale is therefore of central functional importance. Here we show that R_1ρ relaxation dispersion experiments in magic-angle-spinning solid-state NMR spectroscopy make it possible to investigate the thermodynamics and kinetics of such exchange process, and gain insight into structural features of short-lived states.

NMR characterization of the near native and unfolded states of the PTB domain of Dok1: alternate conformations and residual clusters

BACKGROUND

Phosphotyrosine binding (PTB) domains are critically involved in cellular signaling and diseases. PTB domains are categorized into three distinct structural classes namely IRS-like, Shc-like and Dab-like. All PTB domains consist of a core pleckstrin homology (PH) domain with additional structural elements in Shc and Dab groups. The core PH fold of the PTB domain contains a seven stranded β-sheet and a long C-terminal helix.

PRINCIPAL FINDINGS

In this work, the PTB domain of Dok1 protein has been characterized, by use of NMR spectroscopy, in solutions containing sub-denaturing and denaturing concentrations of urea. We find that the Dok1 PTB domain displays, at sub-denaturing concentrations of urea, alternate conformational states for residues located in the C-terminal helix and in the β5 strand of the β-sheet region. The β5 strand of PTB domain has been found to be experiencing significant chemical shift perturbations in the presence of urea. Notably, many of these residues in the helix and in the β5 strand are also involved in ligand binding. Structural and dynamical analyses at 7 M urea showed that the PTB domain is unfolded with islands of motionally restricted regions in the polypeptide chain. Further, the C-terminal helix appears to be persisted in the unfolded state of the PTB domain. By contrast, residues encompassing β-sheets, loops, and the short N-terminal helix lack any preferred secondary structures. Moreover, these residues demonstrated an intimate contact with the denaturant.

SIGNIFICANCE

This study implicates existence of alternate conformational states around the ligand binding pocket of the PTB domain either in the native or in the near native conditions. Further, the current study demonstrates that the C-terminal helical region of PTB domain may be considered as a potential site for the initiation of folding.

Conformational dynamics control ubiquitin-deubiquitinase interactions and influence in vivo signaling

Ubiquitin is a highly conserved eukaryotic protein that interacts with a diverse set of partners to act as a cellular signaling hub. Ubiquitin's conformational flexibility has been postulated to underlie its multifaceted recognition. Here we use computational and library-based means to interrogate core mutations that modulate the conformational dynamics of human ubiquitin. These ubiquitin variants exhibit increased affinity for the USP14 deubiquitinase, with concomitantly reduced affinity for other deubiquitinases. Strikingly, the kinetics of conformational motion are dramatically slowed in these variants without a detectable change in either the ground state fold or excited state population. These variants can be ligated into substrate-linked chains in vitro and in vivo but cannot solely support growth in eukaryotic cells. Proteomic analyses reveal nearly identical interaction profiles between WT ubiquitin and the variants but identify a small subset of altered interactions. Taken together, these results show that conformational dynamics are critical for ubiquitin-deubiquitinase interactions and imply that the fine tuning of motion has played a key role in the evolution of ubiquitin as a signaling hub.

Conformational stabilization of ubiquitin yields potent and selective inhibitors of USP7

Protein conformation and function are often inextricably linked, such that the states a protein adopts define its enzymatic activity or its affinity for various partners. Here we combine computational design with macromolecular display to isolate functional conformations of ubiquitin that tightly bind the catalytic core of the oncogenic ubiquitin-specific protease 7 (USP7) deubiquitinase. Structural and biochemical characterization of these ubiquitin variants suggest that remodeled backbone conformations and core packing poise these molecules for stronger interactions, leading to potent and specific inhibition of enzymatic activity. A ubiquitin variant expressed in human tumor cell lines binds and inhibits endogenous USP7, thereby enhancing Mdm2 proteasomal turnover and stabilizing p53. In sum, we have developed an approach to rationally target macromolecular libraries toward the remodeling of protein conformation, shown that engineering of ubiquitin conformation can greatly increase its interaction with deubiquitinases and developed powerful tools to probe the cellular role of USP7.

Conformational and dynamic changes at the interface contribute to ligand binding by ubiquitin

Ubiquitin interacts with numerous domains and motifs in its lifetime that vary in structure but bind the same hydrophobic patch. To identify the structural features of ubiquitin that make it an exceptional protein-protein interaction partner, we have studied the interaction of ubiquitin with the signal transducing adaptor molecule-1 ubiquitin interacting motif (UIM) using nuclear magnetic resonance. Our studies bring to light the role of the inherent backbone flexibility of ubiquitin in its interactions with a large array of binding partners, revealed from the changes in C(α) chemical shifts, backbone dynamics, and hydrogen bond lengths upon UIM binding. The crystal structures of ubiquitin complexes lend further support to our findings, underscoring the importance of the unique and flexible hydrogen bond network within ubiquitin and simultaneously providing insights into the nature of the slow motions. Taken together, our studies provide an in-depth view of the molecular changes associated with ligand recognition by ubiquitin.

Site-resolved measurement of microsecond-to-millisecond conformational-exchange processes in proteins by solid-state NMR spectroscopy

We demonstrate that conformational exchange processes in proteins on microsecond-to-millisecond time scales can be detected and quantified by solid-state NMR spectroscopy. We show two independent approaches that measure the effect of conformational exchange on transverse relaxation parameters, namely Carr-Purcell-Meiboom-Gill relaxation-dispersion experiments and measurement of differential multiple-quantum coherence decay. Long coherence lifetimes, as required for these experiments, are achieved by the use of highly deuterated samples and fast magic-angle spinning. The usefulness of the approaches is demonstrated by application to microcrystalline ubiquitin. We detect a conformational exchange process in a region of the protein for which dynamics have also been observed in solution. Interestingly, quantitative analysis of the data reveals that the exchange process is more than 1 order of magnitude slower than in solution, and this points to the impact of the crystalline environment on free energy barriers.

Local folding and misfolding in the PBX homeodomain from a three-state analysis of CPMG relaxation dispersion NMR data

NMR Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments represent a powerful approach for characterizing protein internal motions and for gaining insight into fundamental biological processes such as protein folding, catalysis, and allostery. In most cases, CPMG data are analyzed assuming that the protein exchanges between two different conformational states. Systems exchanging among more than two states are far more challenging to characterize by CPMG NMR. For example, in the case of three-state exchange in the fast time scale regime, it is difficult to uniquely connect the parameters extracted from CPMG analyses with the physical parameters of most interest, intercoversion rates, populations, and chemical shift differences for exchanging states. We have developed a grid search selection procedure that allows these physical parameters to be uniquely determined from CPMG data, based on additional information, which in this study comprises ligand-induced chemical shift perturbations. We applied this approach to the PBX homeodomain (PBX-HD), a three-helix protein with a C-terminal extension that folds into a fourth helix upon binding to DNA. We recently showed that the C-terminal extension transiently folds, even in the absence DNA, in a process that is likely tied to the cooperative binding of PBX-HD to DNA and other homeodomains. Using the grid search selection procedure, we found that PBX-HD undergoes exchange between three different conformational states, a major form in which the C-terminal extension is unfolded, the previously identified state in which the C-terminal extension forms a fourth helix, and an additional state in which the C-terminal extension is misfolded.

Structure and dynamics of a conformationally constrained nitroxide side chain and applications in EPR spectroscopy

A disulfide-linked nitroxide side chain (R1) is the most widely used spin label for determining protein topology, mapping structural changes, and characterizing nanosecond backbone motions by site-directed spin labeling. Although the internal motion of R1 and the number of preferred rotamers are limited, translating interspin distance measurements and spatial orientation information into structural constraints is challenging. Here, we introduce a highly constrained nitroxide side chain designated RX as an alternative to R1 for these applications. RX is formed by a facile cross-linking reaction of a bifunctional methanethiosulfonate reagent with pairs of cysteine residues at i and i + 3 or i and i + 4 in an α-helix, at i and i + 2 in a β-strand, or with cysteine residues in adjacent strands in a β-sheet. Analysis of EPR spectra, a crystal structure of RX in T4 lysozyme, and pulsed electron-electron double resonance (ELDOR) spectroscopy on an immobilized protein containing RX all reveal a highly constrained internal motion of the side chain. Consistent with the constrained geometry, interspin distance distributions between pairs of RX side chains are narrower than those from analogous R1 pairs. As an important consequence of the constrained internal motion of RX, spectral diffusion detected with ELDOR reveals microsecond internal motions of the protein. Collectively, the data suggest that the RX side chain will be useful for distance mapping by EPR spectroscopy, determining spatial orientation of helical segments in oriented specimens, and measuring structural fluctuations on the microsecond time scale.

Active site dynamics in NADH oxidase from Thermus thermophilus studied by NMR spin relaxation

We have characterized the backbone dynamics of NADH oxidase from Thermus thermophilus (NOX) using a recently-developed suite of NMR experiments designed to isolate exchange broadening, together with (15)N R (1), R (1ρ ), and {(1)H}-(15)N steady-state NOE relaxation measurements performed at 11.7 and 18.8 T. NOX is a 54 kDa homodimeric enzyme that belongs to a family of structurally homologous flavin reductases and nitroreductases with many potential biotechnology applications. Prior studies have suggested that flexibility is involved in the catalytic mechanism of the enzyme. The active site residue W47 was previously identified as being particularly important, as its level of solvent exposure correlates with enzyme activity, and it was observed to undergo "gating" motions in computer simulations. The NMR data are consistent with these findings. Signals from W47 are dynamically broadened beyond detection and several other residues in the active site have significant R ( ex ) contributions to transverse relaxation rates. In addition, the backbone of S193, whose side chain hydroxyl proton hydrogen bonds directly with the FMN cofactor, exhibits extensive mobility on the ns-ps timescale. We hypothesize that these motions may facilitate structural rearrangements of the active site that allow NOX to accept both FMN and FAD as cofactors.

A hydrogen bond regulates slow motions in ubiquitin by modulating a β-turn flip

Proteins exist as conformational ensembles composed of multiple interchanging substates separated by kinetic barriers. Interconverting conformations are often difficult to probe, owing to their sparse population and transient nature. Here, we report the identification and characterization of a subset of conformations in ubiquitin that participate in microsecond-to-millisecond motions in the amides of Ile23, Asn25, and Thr55. A novel side chain to the backbone hydrogen bond that regulates these motions has also been identified. Combining our NMR studies with the available X-ray data, we have unearthed the physical process underlying slow motions-the interconversion of a type I into a type II β-turn flip at residues Glu51 through Arg54. Interestingly, the dominant conformer of wild-type ubiquitin observed in solution near neutral pH is only represented by about 22% of the crystal structures. The conformers generated as a result of the dynamics of the hydrogen bond appear to be correlated to ligand recognition by ubiquitin.

Recent developments in (15)N NMR relaxation studies that probe protein backbone dynamics

Nuclear Magnetic Resonance (NMR) relaxation is a powerful technique that provides information about internal dynamics associated with configurational energetics in proteins, as well as site-specific information involved in conformational equilibria. In particular, (15)N relaxation is a useful probe to characterize overall and internal backbone dynamics of proteins because the relaxation mainly reflects reorientational motion of the N-H bond vector. Over the past 20 years, experiments and protocols for analysis of (15)N R (1), R 2, and the heteronuclear (15)N-{(1)H} NOE data have been well established. The development of these methods has kept pace with the increase in the available static-magnetic field strength, providing dynamic parameters optimized from data fitting at multiple field strengths. Using these methodological advances, correlation times for global tumbling and order parameters and correlation times for internal motions of many proteins have been determined. More recently, transverse relaxation dispersion experiments have extended the range of NMR relaxation studies to the milli- to microsecond time scale, and have provided quantitative information about functional conformational exchange in proteins. Here, we present an overview of recent advances in (15)N relaxation experiments to characterize protein backbone dynamics.

Mechanism for recognition of polyubiquitin chains: balancing affinity through interplay between multivalent binding and dynamics

RAP80 plays a key role in signal transduction in the DNA damage response by recruiting proteins to DNA damage foci by binding K63-polyubiquitin chains with two tandem ubiquitin-interacting motifs (tUIM). It is generally recognized that the typically weak interaction between ubiquitin (Ub) and various recognition motifs is intensified by themes such as tandem recognition motifs and Ub polymerization to achieve biological relevance. However, it remains an intricate problem to develop a detailed molecular mechanism to describe the process that leads to amplification of the Ub signal. A battery of solution-state NMR methods and molecular dynamics simulations were used to demonstrate that RAP80-tUIM employs mono- and multivalent interactions with polyUb chains to achieve enhanced affinity in comparison to monoUb interactions for signal amplification. The enhanced affinity is balanced by unfavorable entropic effects that include partial quenching of rapid reorientation between individual UIM domains and individual Ub domains in the bound state. For the RAP80-tUIM-polyUb interaction, increases in affinity with increasing chain length are a result of increased numbers of mono- and multivalent binding sites in the longer polyUb chains. The mono- and multivalent interactions are characterized by intrinsically weak binding and fast off-rates; these weak interactions with fast kinetics may be an important factor underlying the transient nature of protein-protein interactions that comprise DNA damage foci.

An introduction to NMR-based approaches for measuring protein dynamics

Proteins are inherently flexible at ambient temperature. At equilibrium, they are characterized by a set of conformations that undergo continuous exchange within a hierarchy of spatial and temporal scales ranging from nanometers to micrometers and femtoseconds to hours. Dynamic properties of proteins are essential for describing the structural bases of their biological functions including catalysis, binding, regulation and cellular structure. Nuclear magnetic resonance (NMR) spectroscopy represents a powerful technique for measuring these essential features of proteins. Here we provide an introduction to NMR-based approaches for studying protein dynamics, highlighting eight distinct methods with recent examples, contextualized within a common experimental and analytical framework. The selected methods are (1) Real-time NMR, (2) Exchange spectroscopy, (3) Lineshape analysis, (4) CPMG relaxation dispersion, (5) Rotating frame relaxation dispersion, (6) Nuclear spin relaxation, (7) Residual dipolar coupling, (8) Paramagnetic relaxation enhancement. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Protein flexibility and conformational entropy in ligand design targeting the carbohydrate recognition domain of galectin-3

Rational drug design is predicated on knowledge of the three-dimensional structure of the protein-ligand complex and the thermodynamics of ligand binding. Despite the fundamental importance of both enthalpy and entropy in driving ligand binding, the role of conformational entropy is rarely addressed in drug design. In this work, we have probed the conformational entropy and its relative contribution to the free energy of ligand binding to the carbohydrate recognition domain of galectin-3. Using a combination of NMR spectroscopy, isothermal titration calorimetry, and X-ray crystallography, we characterized the binding of three ligands with dissociation constants ranging over 2 orders of magnitude. (15)N and (2)H spin relaxation measurements showed that the protein backbone and side chains respond to ligand binding by increased conformational fluctuations, on average, that differ among the three ligand-bound states. Variability in the response to ligand binding is prominent in the hydrophobic core, where a distal cluster of methyl groups becomes more rigid, whereas methyl groups closer to the binding site become more flexible. The results reveal an intricate interplay between structure and conformational fluctuations in the different complexes that fine-tunes the affinity. The estimated change in conformational entropy is comparable in magnitude to the binding enthalpy, demonstrating that it contributes favorably and significantly to ligand binding. We speculate that the relatively weak inherent protein-carbohydrate interactions and limited hydrophobic effect associated with oligosaccharide binding might have exerted evolutionary pressure on carbohydrate-binding proteins to increase the affinity by means of conformational entropy.

Dynamics of protein folding: probing the kinetic network of folding-unfolding transitions with experiment and theory

The problem of spontaneous folding of amino acid chains into highly organized, biologically functional three-dimensional protein structures continues to challenge the modern science. Understanding how proteins fold requires characterization of the underlying energy landscapes as well as the dynamics of the polypeptide chains in all stages of the folding process. In recent years, important advances toward these goals have been achieved owing to the rapidly growing interdisciplinary interest and significant progress in both experimental techniques and theoretical methods. Improvements in the experimental time resolution led to determination of the timescales of the important elementary events in folding, such as formation of secondary structure and tertiary contacts. Sensitive single molecule methods made possible probing the distributions of the unfolded and folded states and following the folding reaction of individual protein molecules. Discovery of proteins that fold in microseconds opened the possibility of atomic-level theoretical simulations of folding and their direct comparisons with experimental data, as well as of direct experimental observation of the barrier-less folding transition. The ultra-fast folding also brought new questions, concerning the intrinsic limits of the folding rates and experimental signatures of barrier-less "downhill" folding. These problems will require novel approaches for even more detailed experimental investigations of the folding dynamics as well as for the analysis of the folding kinetic data. For theoretical simulations of folding, a main challenge is how to extract the relevant information from overwhelmingly detailed atomistic trajectories. New theoretical methods have been devised to allow a systematic approach towards a quantitative analysis of the kinetic network of folding-unfolding transitions between various configuration states of a protein, revealing the transition states and the associated folding pathways at multiple levels, from atomistic to coarse-grained representations. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.