Site-specific backbone amide (15)N chemical shift anisotropy tensors in a small protein from liquid crystal and cross-correlated relaxation measurements

Transverse relaxation optimized spectroscopy of NH2 groups in glutamine and asparagine side chains of proteins

A transverse relaxation optimized spectroscopy (TROSY) approach is described for the optimal detection of NH2 groups in asparagine and glutamine side chains of proteins. Specifically, we have developed NMR experiments for isolating the slow-relaxing 15N and 1H components of NH2 multiplets. Although even modest sensitivity gains in 2D NH2-TROSY correlation maps compared to their decoupled NH2-HSQC counterparts can be achieved only occasionally, substantial improvements in resolution of the NMR spectra are demonstrated for asparagine and glutamine NH2 sites of a buried cavity mutant, L99A, of T4 lysozyme at 5 ºC. The NH2-TROSY approach is applied to CPMG relaxation dispersion measurements at the side chain NH2 positions of the L99A T4 lysozyme mutant - a model system for studies of the role of protein dynamics in ligand binding.

Insights into Molecular Diversity within the FUS/EWS/TAF15 Protein Family: Unraveling Phase Separation of the N-Terminal Low-Complexity Domain from RNA-Binding Protein EWS

The FET protein family, comprising FUS, EWS, and TAF15, plays crucial roles in mRNA maturation, transcriptional regulation, and DNA damage response. Clinically, they are linked to Ewing family tumors and neurodegenerative diseases such as amyotrophic lateral sclerosis. The fusion protein EWS::FLI1, the causative mutation of Ewing sarcoma, arises from a genomic translocation that fuses a portion of the low-complexity domain (LCD) of EWS (EWSLCD) with the DNA binding domain of the ETS transcription factor FLI1. This fusion protein modifies transcriptional programs and disrupts native EWS functions, such as splicing. The exact role of the intrinsically disordered EWSLCD remains a topic of active investigation, but its ability to phase separate and form biomolecular condensates is believed to be central to EWS::FLI1's oncogenic properties. Here, we used paramagnetic relaxation enhancement NMR, microscopy, and all-atom molecular dynamics (MD) simulations to better understand the self-association and phase separation tendencies of the EWSLCD. Our NMR data and mutational analysis suggest that a higher density and proximity of tyrosine residues amplify the likelihood of condensate formation. MD simulations revealed that the tyrosine-rich termini exhibit compact conformations with unique contact networks and provided critical input on the relationship between contacts formed within a single molecule (intramolecular) and inside the condensed phase (intermolecular). These findings enhance our understanding of FET proteins' condensate-forming capabilities and underline differences between EWS, FUS, and TAF15.

Simulation of oriented NMR spectra: Combining molecular dynamics and chemical shift tensor calculations

Solid state NMR is widely used to study the orientation and other structural features of proteins and peptides in lipid bilayers. Using data obtained by PISEMA (Polarization Inversion Spin Exchange at Magic Angle) experiments, periodic spectral patterns arise from well-aligned α-helical molecules. Significant problems in the interpretation of PISEMA spectra may arise for systems that do not form perfectly defined secondary structures, like α-helices, or the signal pattern is disturbed by molecular motion. Here, we present a new method that combines molecular dynamics simulation with tensorial orientational constraints (MDOC) and chemical shift tensor calculations for the simulation and interpretation of PISEMA-like spectra. The calculations include the spectra arising from non α-helical molecules and molecules with non-uniform intrinsic mobility. In a first step, dipolar or quadrupolar interaction tensors drive molecular rotations and reorientations to obtain the proper mean values as observed in corresponding NMR experiments. In a second step, the coordinate snapshots of the MDOC simulations are geometry optimized with the isotropic 15 N chemical shifts as constraints using Bond Polarization Theory (BPT) to provide reliable 15 N CS tensor data. The averaged dipolar 1 H-15 N couplings and the δzz tensor components can then be combined to simulate PISEMA patterns. We apply this method to the ß-helical peptide gramicidin A (gA) and demonstrate that this method enables the assignment of most PISEMA resonances. In addition, MDOC simulations provide local order parameters for the calculated sites. These local order parameters reveal large differences in backbone mobility between L- and D-amino acids of gA.

Insights into Molecular Diversity within the FET Family: Unraveling Phase Separation of the N-Terminal Low Complexity Domain from RNA-Binding Protein EWS

The FET family proteins, which includes FUS, EWS, and TAF15, are RNA chaperones instrumental in processes such as mRNA maturation, transcriptional regulation, and the DNA damage response. These proteins have clinical significance: chromosomal rearrangements in FET proteins are implicated in Ewing family tumors and related sarcomas. Furthermore, point mutations in FUS and TAF15 are associated with neurodegenerative conditions like amyotrophic lateral sclerosis and frontotemporal lobar dementia. The fusion protein EWS::FLI1, the causative mutation of Ewing sarcoma, arises from a genomic translocation that fuses the low-complexity domain (LCD) of EWS (EWSLCD) with the DNA binding domain of the ETS transcription factor FLI1. This fusion not only alters transcriptional programs but also hinders native EWS functions like splicing. However, the precise function of the intrinsically disordered EWSLCD is still a topic of active investigation. Due to its flexible nature, EWSLCD can form transient interactions with itself and other biomolecules, leading to the formation of biomolecular condensates through phase separation - a mechanism thought to be central to the oncogenicity of EWS::FLI1. In our study, we used paramagnetic relaxation enhancement NMR, analytical ultracentrifugation, light microscopy, and all-atom molecular dynamics (MD) simulations to better understand the self-association and phase separation tendencies of EWSLCD. Our aim was to elucidate the molecular events that underpin EWSLCD-mediated biomolecular condensation. Our NMR data suggest tyrosine residues primarily drive the interactions vital for EWSLCD phase separation. Moreover, a higher density and proximity of tyrosine residues amplify the likelihood of condensate formation. Atomistic MD simulations and hydrodynamic experiments revealed that the tyrosine-rich N and C-termini tend to populate compact conformations, establishing unique contact networks, that are connected by a predominantly extended, tyrosine-depleted, linker region. MD simulations provide critical input on the relationship between contacts formed within a single molecule (intramolecular) and inside the condensed phase (intermolecular), and changes in protein conformations upon condensation. These results offer deeper insights into the condensate-forming abilities of the FET proteins and highlights unique structural and functional nuances between EWS and its counterparts, FUS and TAF15.

Intracellular environment can change protein conformational dynamics in cells through weak interactions

Conformational dynamics is important for protein functions, many of which are performed in cells. How the intracellular environment may affect protein conformational dynamics is largely unknown. Here, loop conformational dynamics is studied for a model protein in Escherichia coli cells by using nuclear magnetic resonance (NMR) spectroscopy. The weak interactions between the protein and surrounding macromolecules in cells hinder the protein rotational diffusion, which extends the dynamic detection timescale up to microseconds by the NMR spin relaxation method. The loop picosecond to microsecond dynamics is confirmed by nanoparticle-assisted spin relaxation and residual dipolar coupling methods. The loop interactions with the intracellular environment are perturbed through point mutation of the loop sequence. For the sequence of the protein that interacts stronger with surrounding macromolecules, the loop becomes more rigid in cells. In contrast, the mutational effect on the loop dynamics in vitro is small. This study provides direct evidence that the intracellular environment can modify protein loop conformational dynamics through weak interactions.

Fast Motions Dominate Dynamics of Intrinsically Disordered Tau Protein at High Temperatures

Reorientational dynamics of intrinsically disordered proteins (IDPs) contain multiple motions often clustered around three motional modes: ultrafast librational motions of amide groups, fast local backbone conformational fluctuations and slow chain segmental motions. This dynamic picture is mainly based on ¹⁵ N NMR relaxation studies of IDPs at relatively low temperatures where the amide-water proton exchange rates are sufficiently small. Less is known, however, about the dynamics of IDPs at more physiological temperatures. Here, we investigate protein dynamics in a 441-residue long IDP, tau protein, in the temperature range from 0-25 °C, using ¹⁵ N NMR relaxation rates and spectral density analysis. While at these temperatures relaxation rates are still better described in terms of amide group librational motions, local backbone dynamics and chain segmental motions, the temperature-dependent trend of spectral densities suggests that the timescales of fast backbone conformational fluctuations and slower chain segmental motions might become inseparable at higher temperatures. Our data demonstrate the remarkable dynamic plasticity of this prototypical IDP and highlight the need for dynamic studies of IDPs at multiple temperatures.

Bromodomain Interactions with Acetylated Histone 4 Peptides in the BRD4 Tandem Domain: Effects on Domain Dynamics and Internal Flexibility

The bromodomain and extra-terminal (BET) protein BRD4 regulates gene expression via recruitment of transcriptional regulatory complexes to acetylated chromatin. Like other BET proteins, BRD4 contains two bromodomains, BD1 and BD2, that can interact cooperatively with target proteins and designed ligands, with important implications for drug discovery. Here, we used nuclear magnetic resonance (NMR) spectroscopy to study the dynamics and interactions of the isolated bromodomains, as well as the tandem construct including both domains and the intervening linker, and investigated the effects of binding a tetra-acetylated peptide corresponding to the tail of histone 4. The peptide affinity is lower for both domains in the tandem construct than for the isolated domains. Using ¹⁵N spin relaxation, we determined the global rotational correlation times and residue-specific order parameters for BD1 and BD2. Isolated BD1 is monomeric in the apo state but apparently dimerizes upon binding the tetra-acetylated peptide. Isolated BD2 partially dimerizes in both the apo and peptide-bound states. The backbone order parameters reveal marked differences between BD1 and BD2, primarily in the acetyl-lysine binding site where the ZA loop is more flexible in BD2. Peptide binding reduces the order parameters of the ZA loop in BD1 and the ZA and BC loops in BD2. The AB loop, located distally from the binding site, shows variable dynamics that reflect the different dimerization propensities of the domains. These results provide a basis for understanding target recognition by BRD4.

Solid-State NMR Dipolar and Chemical Shift Anisotropy Recoupling Techniques for Structural and Dynamical Studies in Biological Systems

With the development of NMR methodology and technology during the past decades, solid-state NMR (ssNMR) has become a particularly important tool for investigating structure and dynamics at atomic scale in biological systems, where the recoupling techniques play pivotal roles in modern high-resolution MAS NMR. In this review, following a brief introduction on the basic theory of recoupling in ssNMR, we highlight the recent advances in dipolar and chemical shift anisotropy recoupling methods, as well as their applications in structural determination and dynamical characterization at multiple time scales (i.e., fast-, intermediate-, and slow-motion). The performances of these prevalent recoupling techniques are compared and discussed in multiple aspects, together with the representative applications in biomolecules. Given the recent emerging advances in NMR technology, new challenges for recoupling methodology development and potential opportunities for biological systems are also discussed.

Generating Five Independent Molecular Alignments for Simultaneous Protein Structure and Dynamics Determination Using Nuclear Magnetic Resonance Spectroscopy

Residual dipolar couplings (RDCs) are commonly used in NMR for protein structure and dynamics studies, but it is challenging to generate five independent RDC data sets (required for simultaneous structure and dynamics determination) for most protein molecules in the magnetic field. In this work, a reporter protein with a lanthanide tag is introduced to create five independent alignments. This reporter protein is then attached to target proteins where five independent sets of RDCs are also obtained for the target proteins. The fitting of RDCs provides important information about the structure and dynamics of the target proteins. The method is simple and effective and, in principle, can be used to generate complete sets of RDCs for different protein molecules.

Extensive tests and evaluation of the CHARMM36IDPSFF force field for intrinsically disordered proteins and folded proteins

Intrinsically disordered proteins (IDPs) have received increasing attention in recent studies due to their structural heterogeneity and critical biological functions. To fully understand the structural properties and determine accurate ensembles of IDPs, molecular dynamics (MD) simulation was widely used to sample diverse conformations and reveal the structural dynamics. However, the classical state-of-the-art force fields perform well for folded proteins while being unsatisfactory for the simulations of disordered proteins reported in many previous studies. Thus, improved force fields were developed to precisely describe both folded proteins and disordered proteins. Preliminary tests show that our newly developed CHARMM36IDPSFF (C36IDPSFF) force field can well reproduce the experimental observables of several disordered proteins, but more tests on different types of proteins are needed to further evaluate the performance of C36IDPSFF. Here, we extensively simulate short peptides, disordered proteins, and fast-folding proteins as well as folded proteins, and compare the simulated results with the experimental observables. The simulation results show that C36IDPSFF could substantially reproduce the experimental observables for most of the tested proteins but some limitations are also found in the radius of gyration of large disordered proteins and the stability of fast-folding proteins. This force field will facilitate large scale studies of protein structural dynamics and functions using MD simulations.

What Drives 15N Spin Relaxation in Disordered Proteins? Combined NMR/MD Study of the H4 Histone Tail

Backbone (¹⁵N) NMR relaxation is one of the main sources of information on dynamics of disordered proteins. Yet, we do not know very well what drives ¹⁵N relaxation in such systems, i.e., how different forms of motion contribute to the measurable relaxation rates. To address this problem, we have investigated, both experimentally and via molecular dynamics simulations, the dynamics of a 26-residue peptide imitating the N-terminal portion of the histone protein H4. One part of the peptide was found to be fully flexible, whereas the other part features some transient structure (a hairpin stabilized by hydrogen bonds). The following motional modes proved relevant for ¹⁵N relaxation. 1) Sub-picosecond librations attenuate relaxation rates according to S² ∼0.85-0.90. 2) Axial peptide-plane fluctuations along a stretch of the peptide chain contribute to relaxation-active dynamics on a fast timescale (from tens to hundreds of picoseconds). 3) φ/ψ backbone jumps contribute to relaxation-active dynamics on both fast (from tens to hundreds of picoseconds) and slow (from hundreds of picoseconds to a nanosecond) timescales. The major contribution is from polyproline II (PPII) ↔ β transitions in the Ramachandran space; in the case of glycine residues, the major contribution is from PPII ↔ (β) ↔ rPPII transitions, in which rPPII is the mirror-image (right-handed) version of the PPII geometry, whereas β geometry plays the role of an intermediate state. 4) Reorientational motion of certain (sufficiently long-lived) elements of transient structure, i.e., rotational tumbling, contributes to slow relaxation-active dynamics on ∼1-ns timescale (however, it is difficult to isolate this contribution). In conclusion, recent advances in the area of force-field development have made it possible to obtain viable Molecular Dynamics models of protein disorder. After careful validation against the experimental relaxation data, these models can provide a valuable insight into mechanistic origins of spin relaxation in disordered peptides and proteins.

Correlated motions of C'-N and Cα-Cβ pairs in protonated and per-deuterated GB3

We investigated correlated µs-ms time scale motions of neighboring ¹³C'-¹⁵N and ¹³C_α-¹³C_β nuclei in both protonated and perdeuterated samples of GB3. The techniques employed, NMR relaxation due to cross-correlated chemical shift modulations, specifically target concerted changes in the isotropic chemical shifts of the two nuclei associated with spatial fluctuations. Field-dependence of the relaxation rates permits identification of the parameters defining the chemical exchange rate constant under the assumption of a two-site exchange. The time scale of motions falls into the intermediate to fast regime (with respect to the chemical shift time scale, 100-400 s^-1 range) for the ¹³C'-¹⁵N pairs and into the slow to intermediate regime for the ¹³C_α-¹³C_β pairs (about 150 s^-1). Comparison of the results obtained for protonated and deuterated GB3 suggests that deuteration has a tendency to reduce these slow scale correlated motions, especially for the ¹³C_α-¹³C_β pairs.

NMR Methods for Characterizing the Basic Side Chains of Proteins: Electrostatic Interactions, Hydrogen Bonds, and Conformational Dynamics

NMR spectroscopy is a powerful tool for studying protein dynamics. Conventionally, NMR studies on protein dynamics have probed motions of protein backbone NH, side-chain aromatic, and CH₃ groups. Recently, there has been remarkable progress in NMR methodologies that can characterize motions of cationic groups in protein side chains. These NMR methods allow investigations of the dynamics of positively charged lysine (Lys) and arginine (Arg) side chains and their hydrogen bonds as well as their electrostatic interactions important for protein function. Here, describing various practical aspects, we provide an overview of the NMR methods for dynamics studies of Lys and Arg side chains. Some example data on protein-DNA complexes are shown. We will also explain how molecular dynamics (MD) simulations can facilitate the interpretation of the NMR data on these basic side chains. Studies combining NMR and MD have revealed the highly dynamic nature of short-range electrostatic interactions via ion pairs, especially those involving Lys side chains.

Fast evaluation of protein dynamics from deficient 15N relaxation data

Simple and convenient method of protein dynamics evaluation from the insufficient experimental ¹⁵N relaxation data is presented basing on the ratios, products, and differences of longitudinal and transverse ¹⁵N relaxation rates obtained at a single magnetic field. Firstly, the proposed approach allows evaluating overall tumbling correlation time (nanosecond time scale). Next, local parameters of the model-free approach characterizing local mobility of backbone amide N-H vectors on two different time scales, S² and R _ex , can be elucidated. The generalized order parameter, S², describes motions on the time scale faster than the overall tumbling correlation time (pico- to nanoseconds), while the chemical exchange term, R _ex , identifies processes slower than the overall tumbling correlation time (micro- to milliseconds). Advantages and disadvantages of different methods of data handling are thoroughly discussed.

Atomic resolution conformational dynamics of intrinsically disordered proteins from NMR spin relaxation

Nuclear magnetic resonance (NMR) spectroscopy is one of the most powerful experimental approaches for investigating the conformational behaviour of intrinsically disordered proteins (IDPs). IDPs represent a significant fraction of all proteomes, and, despite their importance for understanding fundamental biological processes, the molecular basis of their activity still remains largely unknown. The functional mechanisms exploited by IDPs in their interactions with other biomolecules are defined by their intrinsic dynamic modes and associated timescales, justifying the considerable interest over recent years in the development of technologies adapted to measure and describe this behaviour. NMR spin relaxation delivers information-rich, site-specific data reporting on conformational fluctuations occurring throughout the molecule. Here we review recent progress in the use of ¹⁵N relaxation to identify local backbone dynamics and long-range chain-like motions in unfolded proteins.

Intrinsic Differences in Backbone Dynamics between Wild Type and DNA-Contact Mutants of the p53 DNA Binding Domain Revealed by Nuclear Magnetic Resonance Spectroscopy

Mutations in p53's DNA binding domain (p53DBD) are associated with 50% of all cancers, making it an essential system to investigate and understand the genesis and progression of cancer. In this work, we studied the changes in the structure and dynamics of wild type p53DBD in comparison with two of its "hot-spot" DNA-contact mutants, R248Q and R273H, by analysis of backbone amide chemical shift perturbations and ¹⁵N spin relaxation measurements. The results of amide chemical shift changes indicated significantly more perturbations in the R273H mutant than in wild type and R248Q p53DBD. Analysis of ¹⁵N spin relaxation rates and the resulting nuclear magnetic resonance order parameters suggests that for most parts, the R248Q mutant exhibits limited conformational flexibility and is similar to the wild type protein. In contrast, R273H showed significant backbone dynamics extending up to its β-sandwich scaffold in addition to motions along the DNA binding interface. Furthermore, comparison of rotational correlation times between the mutants suggests that the R273H mutant, with a higher correlation time, forms an enlarged structural fold in comparison to the R248Q mutant and wild type p53DBD. Finally, we identify three regions in these proteins that show conformational flexibility to varying degrees, which suggests that the R273H mutant, in addition to being a DNA-contact mutation, exhibits properties of a conformational mutant.

Using Paramagnetism to Slow Down Nuclear Relaxation in Protein NMR

Paramagnetic metal ions accelerate nuclear spin relaxation; this effect is widely used for distance measurement and called paramagnetic relaxation enhancement (PRE). Theoretical predictions established that, under special circumstances, it is also possible to achieve a reduction in nuclear relaxation rates (negative PRE). This situation would occur if the mechanism of nuclear relaxation in the diamagnetic state is counterbalanced by a paramagnetic relaxation mechanism caused by the metal ion. Here we report the first experimental evidence for such a cross-correlation effect. Using a uniformly ¹⁵N-labeled mutant of calbindin D_9k loaded with either Tm³⁺ or Tb³⁺, reduced R₁ and R₂ relaxation rates of backbone ¹⁵N spins were observed compared with the diamagnetic reference (the same protein loaded with Y³⁺). The effect arises from the compensation of the chemical shift anisotropy tensor by the anisotropic dipolar shielding generated by the unpaired electron spin.

Quantifying protein dynamics in the ps-ns time regime by NMR relaxation

Both ¹⁵N chemical shift anisotropy (CSA) and sufficiently rapid exchange linebroadening transitions exhibit relaxation contributions that are proportional to the square of the magnetic field. Deconvoluting these contributions is further complicated by residue-dependent variations in protein amide ¹⁵N CSA values which have proven difficult to accurately measure. Exploiting recently reported improvements for the implementation of T₁ and T_1ρ experiments, field strength-dependent studies have been carried out on the B3 domain of protein G (GB3) as well as on the immunophilin FKBP12 and a H87V variant of that protein in which the major conformational exchange linebroadening transition is suppressed. By applying a zero frequency spectral density rescaling analysis to the relaxation data collected at magnetic fields from 500 to 900 MHz ¹H, differential residue-specific ¹⁵N CSA values have been obtained for GB3 which correlate with those derived from solid state and liquid crystalline NMR measurements to a level similar to the correlation among those previously reported studies. Application of this analysis protocol to FKBP12 demonstrated an efficient quantitation of both weak exchange linebroadening contributions and differential residue-specific ¹⁵N CSA values. Experimental access to such differential residue-specific ¹⁵N CSA values should significantly facilitate more accurate comparisons with molecular dynamics simulations of protein motion that occurs within the timeframe of global molecular tumbling.

High-power (1)H composite pulse decoupling provides artifact free exchange-mediated saturation transfer (EST) experiments

Exchange-mediated saturation transfer (EST) provides critical information regarding dynamics of molecules. In typical applications EST is studied by either scanning a wide range of (15)N chemical shift offsets where the applied (15)N irradiation field strength is on the order of hundreds of Hertz or, scanning a narrow range of (15)N chemical shift offsets where the applied (15)N irradiation field-strength is on the order of tens of Hertz during the EST period. The (1)H decoupling during the EST delay is critical as incomplete decoupling causes broadening of the EST profile, which could possibly result in inaccuracies of the extracted kinetic parameters and transverse relaxation rates. Currently two different (1)H decoupling schemes have been employed, intermittently applied 180° pulses and composite-pulse-decoupling (CPD), for situations where a wide range, or narrow range of (15)N chemical shift offsets are scanned, respectively. We show that high-power CPD provides artifact free EST experiments, which can be universally implemented regardless of the offset range or irradiation field-strengths.

Proton-detected 3D (15)N/(1)H/(1)H isotropic/anisotropic/isotropic chemical shift correlation solid-state NMR at 70kHz MAS

Chemical shift anisotropy (CSA) tensors offer a wealth of information for structural and dynamics studies of a variety of chemical and biological systems. In particular, CSA of amide protons can provide piercing insights into hydrogen-bonding interactions that vary with the backbone conformation of a protein and dynamics. However, the narrow span of amide proton resonances makes it very difficult to measure (1)H CSAs of proteins even by using the recently proposed 2D (1)H/(1)H anisotropic/isotropic chemical shift (CSA/CS) correlation technique. Such difficulties due to overlapping proton resonances can in general be overcome by utilizing the broad span of isotropic chemical shifts of low-gamma nuclei like (15)N. In this context, we demonstrate a proton-detected 3D (15)N/(1)H/(1)H CS/CSA/CS correlation experiment at fast MAS frequency (70kHz) to measure (1)H CSA values of unresolved amide protons of N-acetyl-(15)N-l-valyl-(15)N-l-leucine (NAVL).

Identification of Dynamic Modes in an Intrinsically Disordered Protein Using Temperature-Dependent NMR Relaxation

The dynamic modes and time scales sampled by intrinsically disordered proteins (IDPs) define their function. Nuclear magnetic resonance (NMR) spin relaxation is probably the most powerful tool for investigating these motions delivering site-specific descriptions of conformational fluctuations from throughout the molecule. Despite the abundance of experimental measurement of relaxation in IDPs, the physical origin of the measured relaxation rates remains poorly understood. Here we measure an extensive range of auto- and cross-correlated spin relaxation rates at multiple magnetic field strengths on the C-terminal domain of the nucleoprotein of Sendai virus, over a large range of temperatures (268-298 K), and combine these data to describe the dynamic behavior of this archetypal IDP. An Arrhenius-type relationship is used to simultaneously analyze up to 61 relaxation rates per amino acid over the entire temperature range, allowing the measurement of local activation energies along the chain, and the assignment of physically distinct dynamic modes. Fast (τ ≤ 50 ps) components report on librational motions, a dominant mode occurs on time scales around 1 ns, apparently reporting on backbone sampling within Ramachandran substates, while a slower component (5-25 ns) reports on segmental dynamics dominated by the chain-like nature of the protein. Extending the study to three protein constructs of different lengths (59, 81, and 124 amino acids) substantiates the assignment of these contributions. The analysis is shown to be remarkably robust, accurately predicting a broad range of relaxation data measured at different magnetic field strengths and temperatures. The ability to delineate intrinsic modes and time scales from NMR spin relaxation will improve our understanding of the behavior and function of IDPs, adding a new and essential dimension to the description of this biologically important and ubiquitous class of proteins.

Nucleophosmin integrates within the nucleolus via multi-modal interactions with proteins displaying R-rich linear motifs and rRNA

The nucleolus is a membrane-less organelle formed through liquid-liquid phase separation of its components from the surrounding nucleoplasm. Here, we show that nucleophosmin (NPM1) integrates within the nucleolus via a multi-modal mechanism involving multivalent interactions with proteins containing arginine-rich linear motifs (R-motifs) and ribosomal RNA (rRNA). Importantly, these R-motifs are found in canonical nucleolar localization signals. Based on a novel combination of biophysical approaches, we propose a model for the molecular organization within liquid-like droplets formed by the N-terminal domain of NPM1 and R-motif peptides, thus providing insights into the structural organization of the nucleolus. We identify multivalency of acidic tracts and folded nucleic acid binding domains, mediated by N-terminal domain oligomerization, as structural features required for phase separation of NPM1 with other nucleolar components in vitro and for localization within mammalian nucleoli. We propose that one mechanism of nucleolar localization involves phase separation of proteins within the nucleolus.

DOI:http://dx.doi.org/10.7554/eLife.13571.001

Inside cells, machines called ribosomes assemble proteins from building blocks known as amino acids. Cells can alter the numbers of ribosomes they produce to match the cell’s demand for new proteins. For instance, when cells grow they require a lot of new proteins and therefore more ribosomes are produced. However, when cells face harsh conditions that cause stress (e.g. exposure to UV radiation or a harmful chemical) they generally stop growing and therefore need fewer ribosomes.

In human and other eukaryotic cells, ribosomes are assembled in a structure called the nucleolus. However, because the nucleolus is not separated from the rest of the cell by a membrane, it was not clear how it is able to accumulate large quantities of the proteins and other molecules needed to make ribosomes. Recent work suggests that the nucleolus is formed through a process referred to as “phase separation” in which the liquid in a particular region of the cell has different physical properties to the liquid surrounding it. This is like how oil and water form separate layers when mixed.

A protein called nucleophosmin is found at high levels in the nucleolus where it interacts with many other proteins, including those involved in making ribosomes. Nucleophosmin binds to motifs within these proteins that contain multiple copies of an amino acid called arginine (referred to as R-motifs). Now, Mitrea et al. investigate how nucleophosmin binds to R-motif proteins and whether this is important for assembling the nucleolus. A search for R-motifs in a list of over a hundred proteins known to bind to nucleophosmin showed that the majority of these proteins contained multiple R-motifs. Furthermore, when high levels of nucleophosmin and the R-motif proteins were present, they underwent phase separation.

Next, Mitrea et al. examine the changes in how nucleophosmin and a ribosomal protein interact before and after phase separation. The experiments show that many molecules of nucleophosmin bind to each other and that multiple regions in nucleophosmin are able to interact with the R-motifs. Together, these interactions produce large assemblies of proteins that result in the creation of separate liquid layers. Furthermore, the experiments show that R-motif proteins and other molecules needed to make ribosomes can be brought together within the same liquid phase by nucleophosmin. Mitrea et al.’s findings provide the first insights into the role of nucleophosmin in the molecular organisation of the nucleolus. The next challenge is to understand how this organisation promotes the production of ribosomes and helps the cell to respond to stressful situations.

DOI:http://dx.doi.org/10.7554/eLife.13571.002

(15)N CSA tensors and (15)N-(1)H dipolar couplings of protein hydrophobic core residues investigated by static solid-state NMR

In this work, we assess the usefulness of static (15)N NMR techniques for the determination of the (15)N chemical shift anisotropy (CSA) tensor parameters and (15)N-(1)H dipolar splittings in powder protein samples. By using five single labeled samples of the villin headpiece subdomain protein in a hydrated lyophilized powder state, we determine the backbone (15)N CSA tensors at two temperatures, 22 and -35 °C, in order to get a snapshot of the variability across the residues and as a function of temperature. All sites probed belonged to the hydrophobic core and most of them were part of α-helical regions. The values of the anisotropy (which include the effect of the dynamics) varied between 130 and 156 ppm at 22 °C, while the values of the asymmetry were in the 0.32-0.082 range. The Leu-75 and Leu-61 backbone sites exhibited high mobility based on the values of their temperature-dependent anisotropy parameters. Under the assumption that most differences stem from dynamics, we obtained the values of the motional order parameters for the (15)N backbone sites. While a simple one-dimensional line shape experiment was used for the determination of the (15)N CSA parameters, a more advanced approach based on the "magic sandwich" SAMMY pulse sequence (Nevzorov and Opella, 2003) was employed for the determination of the (15)N-(1)H dipolar patterns, which yielded estimates of the dipolar couplings. Accordingly, the motional order parameters for the dipolar interaction were obtained. It was found that the order parameters from the CSA and dipolar measurements are highly correlated, validating that the variability between the residues is governed by the differences in dynamics. The values of the parameters obtained in this work can serve as reference values for developing more advanced magic-angle spinning recoupling techniques for multiple labeled samples.

The Quest for Simplicity: Remarks on the Free-Approach Models

Nuclear magnetic relaxation provides a powerful method giving insight into molecular motions at atomic resolution on a broad time scale. Dynamics of biological macromolecules has been widely exploited by measuring (15)N and (13)C relaxation data. Interpretation of these data relies almost exclusively on the use of the model-free approach (MFA) and its extended version (EMFA) which requires no particular physical model of motion and a small number of parameters. It is shown that EMFA is often unable to cope with three different time scales and fails to describe slow internal motions properly. In contrast to EMFA, genuine MFA with two time scales can reproduce internal motions slower than the overall tumbling. It is also shown that MFA and simplified EMFA are equivalent with respect to the values of the N-H bond length and chemical shift anisotropy. Therefore, the vast majority of (15)N relaxation data for proteins can be satisfactorily interpreted solely with MFA.

Structural dynamics and conformational equilibria of SERCA regulatory proteins in membranes by solid-state NMR restrained simulations

Solid-state NMR spectroscopy is emerging as a powerful approach to determine structure, topology, and conformational dynamics of membrane proteins at the atomic level. Conformational dynamics are often inferred and quantified from the motional averaging of the NMR parameters. However, the nature of these motions is difficult to envision based only on spectroscopic data. Here, we utilized restrained molecular dynamics simulations to probe the structural dynamics, topology and conformational transitions of regulatory membrane proteins of the calcium ATPase SERCA, namely sarcolipin and phospholamban, in explicit lipid bilayers. Specifically, we employed oriented solid-state NMR data, such as dipolar couplings and chemical shift anisotropy measured in lipid bicelles, to refine the conformational ensemble of these proteins in lipid membranes. The samplings accurately reproduced the orientations of transmembrane helices and showed a significant degree of convergence with all of the NMR parameters. Unlike the unrestrained simulations, the resulting sarcolipin structures are in agreement with distances and angles for hydrogen bonds in ideal helices. In the case of phospholamban, the restrained ensemble sampled the conformational interconversion between T (helical) and R (unfolded) states for the cytoplasmic region that could not be observed using standard structural refinements with the same experimental data set. This study underscores the importance of implementing NMR data in molecular dynamics protocols to better describe the conformational landscapes of membrane proteins embedded in realistic lipid membranes.

NMR studies of dynamic biomolecular conformational ensembles

Multidimensional heteronuclear NMR approaches can provide nearly complete sequential signal assignments of isotopically enriched biomolecules. The availability of assignments together with measurements of spin relaxation rates, residual spin interactions, J-couplings and chemical shifts provides information at atomic resolution about internal dynamics on timescales ranging from ps to ms, both in solution and in the solid state. However, due to the complexity of biomolecules, it is not possible to extract a unique atomic-resolution description of biomolecular motions even from extensive NMR data when many conformations are sampled on multiple timescales. For this reason, powerful computational approaches are increasingly applied to large NMR data sets to elucidate conformational ensembles sampled by biomolecules. In the past decade, considerable attention has been directed at an important class of biomolecules that function by binding to a wide variety of target molecules. Questions of current interest are: "Does the free biomolecule sample a conformational ensemble that encompasses the conformations found when it binds to various targets; and if so, on what time scale is the ensemble sampled?" This article reviews recent efforts to answer these questions, with a focus on comparing ensembles obtained for the same biomolecules by different investigators. A detailed comparison of results obtained is provided for three biomolecules: ubiquitin, calmodulin and the HIV-1 trans-activation response RNA.

Solid-state nuclear magnetic resonance spectroscopy for membrane protein structure determination

Solid-state NMR (ssNMR) is a versatile technique that can provide high-resolution (sub-angstrom) structural data for integral membrane proteins embedded in native and model membrane environments. The methodologies for a priori structure determination have for the most part been developed using samples with crystalline and fibrous morphologies. However, the techniques are now being applied to large, polytopic membrane proteins including receptors, ion channels, and porins. ssNMR data may be used to annotate and refine existing structures in regions of the protein not fully resolved by crystallography (including ligand-binding sites and mobile solvent accessible loop regions). This review describes the spectroscopic experiments and data analysis methods (including assignment) used to generate high-resolution structural data for membrane proteins. We also consider the range of sample morphologies that are appropriate for study by this method.

Time-shared experiments for efficient assignment of triple-selectively labeled proteins

Combinatorial triple-selective labeling facilitates the NMR assignment process for proteins that are subject to signal overlap and insufficient signal-to-noise in standard triple-resonance experiments. Aiming at maximum amino-acid type and sequence-specific information, the method represents a trade-off between the number of selectively labeled samples that have to be prepared and the number of spectra to be recorded per sample. In order to address the demand of long measurement times, we here propose pulse sequences in which individual phase-shifted transients are stored separately and recombined later to produce several 2D HN(CX) type spectra that are usually acquired sequentially. Sign encoding by the phases of (13)C 90° pulses allows to either select or discriminate against (13)C' or (13)C(α) spins coupled to (15)N. As a result, (1)H-(15)N correlation maps of the various isotopomeric species present in triple-selectively labeled proteins are deconvoluted which in turn reduces problems due to spectral overlap. The new methods are demonstrated with four different membrane proteins with rotational correlation times ranging from 18 to 52 ns.

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.

Spectral density mapping protocols for analysis of molecular motions in disordered proteins

Spectral density mapping represents the method of choice for investigations of molecular motions of intrinsically disordered proteins (IDPs). However, the current methodology has been developed for well-folded proteins. In order to find conditions for a reliable analysis of relaxation of IDPs, accuracy of the current reduced spectral density mapping protocols applied to IDPs was examined and new spectral density mapping methods employing cross-correlated relaxation rates have been designed. Various sources of possible systematic errors were analyzed theoretically and the presented approaches were tested on a partially disordered protein, delta subunit of bacterial RNA polymerase. Results showed that the proposed protocols provide unbiased description of molecular motions of IDPs and allow to separate slow exchange from fast dynamics.

Quantum Mechanical Study of Vicinal J Spin-Spin Coupling Constants for the Protein Backbone

We have performed densisty functional theory (DFT) calculations of vicinal J coupling constants involving the backbone torsional angle for the protein GB3 using our recently developed automatic fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach (Xiao He et al. J. Phys. Chem. B 2009, 113, 10380-10388). Interestingly, the calculated values based on an NMR structure are more accurate than those based on a high-resolution X-ray strucure because the NMR structure was refined using a large number of residual dipolar couplings (RDCs) whereas the hydrogen atoms were added into the X-ray structure in idealized positions, confirming that the postioning of the hydrogen atoms relative to the backbone atoms is important to the accuracy of J coupling constant prediction. By comparing three Karplus equations, our results have demonstrated that hydrogen bonding, substituent and electrostatic effects could have significant impacts on vicinal J couplings even though they depend mostly on the intervening dihedral angles. The root-mean-square deviations (RMSDs) of the calculated (3)J(H(N),H(α)), (3)J(H(N),C(β)), (3)J(H(N),C') values based on the NMR structure are 0.52, 0.25, and 0.35 Hz, respectively, after taking the dynamic effect into consideration. The excellent accuracy demonstrates that our AF-QM/MM approach is a useful tool to study the relationship between J coupling constants and the structure and dynamics of proteins.

A 13C labeling strategy reveals a range of aromatic side chain motion in calmodulin

NMR relaxation experiments often require site-specific isotopic enrichment schemes in order to allow for quantitative interpretation. Here we describe a new labeling scheme for site-specific (13)C-(1)H enrichment of a single ortho position of aromatic amino acid side chains in an otherwise perdeuterated background by employing a combination of [4-(13)C]erythrose and deuterated pyruvate during growth on deuterium oxide. This labeling scheme largely eliminates undesired contributions to (13)C relaxation and greatly simplifies the fitting of relaxation data using the Lipari-Szabo model-free formalism. This approach is illustrated with calcium-saturated vertebrate calmodulin and oxidized flavodoxin from Cyanobacterium anabaena . Analysis of (13)C relaxation in the aromatic groups of calcium-saturated calmodulin indicates a wide range of motion in the subnanosecond time regime.

Structural plasticity in human heterochromatin protein 1β

As essential components of the molecular machine assembling heterochromatin in eukaryotes, HP1 (Heterochromatin Protein 1) proteins are key regulators of genome function. While several high-resolution structures of the two globular regions of HP1, chromo and chromoshadow domains, in their free form or in complex with recognition-motif peptides are available, less is known about the conformational behavior of the full-length protein. Here, we used NMR spectroscopy in combination with small angle X-ray scattering and dynamic light scattering to characterize the dynamic and structural properties of full-length human HP1β (hHP1β) in solution. We show that the hinge region is highly flexible and enables a largely unrestricted spatial search by the two globular domains for their binding partners. In addition, the binding pockets within the chromo and chromoshadow domains experience internal dynamics that can be useful for the versatile recognition of different binding partners. In particular, we provide evidence for the presence of a distinct structural propensity in free hHP1β that prepares a binding-competent interface for the formation of the intermolecular β-sheet with methylated histone H3. The structural plasticity of hHP1β supports its ability to bind and connect a wide variety of binding partners in epigenetic processes.

Quantum chemical calculations of amide-15N chemical shift anisotropy tensors for a membrane-bound cytochrome-b5

There is considerable interest in determining amide-(15)N chemical shift anisotropy (CSA) tensors from biomolecules and understanding their variation for structural and dynamics studies using solution and solid-state NMR spectroscopy and also by quantum chemical calculations. Due to the difficulties associated with the measurement of CSA tensors from membrane proteins, NMR-based structural studies heavily relied on the CSA tensors determined from model systems, typically single crystals of model peptides. In the present study, the principal components of backbone amide-(15)N CSA tensors have been determined using density functional theory for a 16.7 kDa membrane-bound paramagnetic heme containing protein, cytochrome-b(5) (cytb(5)). All the calculations were performed by taking residues within 5 Å distance from the backbone amide-(15)N nucleus of interest. The calculated amide-(15)N CSA spans agree less well with our solution NMR data determined for an effective internuclear distance r(N-H) = 1.023 Å and a constant angle β = 18° that the least shielded component (δ(11)) makes with the N-H bond. The variation of amide-(15)N CSA span obtained using quantum chemical calculations is found to be smaller than that obtained from solution NMR measurements, whereas the trends of the variations are found to be in close agreement. We believe that the results reported in this study will be useful in studying the structure and dynamics of membrane proteins and heme-containing proteins, and also membrane-bound protein-protein complexes such as cytochromes-b5-P450.

Starting-structure dependence of nanosecond timescale intersubstate transitions and reproducibility of MD-derived order parameters

Factors affecting the accuracy of molecular dynamics (MD) simulations are investigated by comparing generalized order parameters for backbone NH vectors of the B3 immunoglobulin-binding domain of streptococcal protein G (GB3) derived from simulations with values obtained from NMR spin relaxation (Yao L, Grishaev A, Cornilescu G, Bax A, J Am Chem Soc 2010;132:4295-4309.). Choices for many parameters of the simulations, such as buffer volume, water model, or salt concentration, have only minor influences on the resulting order parameters. In contrast, seemingly minor conformational differences in starting structures, such as orientations of sidechain hydroxyl groups, resulting from applying different protonation algorithms to the same structure, have major effects on backbone dynamics. Some, but not all, of these effects are mitigated by increased sampling in simulations. Most discrepancies between simulated and experimental results occur for residues located at the ends of secondary structures and involve large amplitude nanosecond timescale transitions between distinct conformational substates. These transitions result in autocorrelation functions for bond vector reorientation that do not converge when calculated over individual simulation blocks, typically of length similar to the overall rotational diffusion time. A test for convergence before averaging the order parameters from different blocks results in better agreement between order parameters calculated from different sets of simulations and with NMR-derived order parameters. Thus, MD-derived order parameters are more strongly affected by transitions between conformational substates than by fluctuations within individual substates themselves, while conformational differences in the starting structures affect the frequency and scale of such transitions.

Determination of accurate 1H positions of an alanine tripeptide with anti-parallel and parallel β-sheet structures by high resolution 1H solid state NMR and GIPAW chemical shift calculation

The accurate (1)H positions of alanine tripeptide, A(3), with anti-parallel and parallel β-sheet structures could be determined by highly resolved (1)H DQMAS solid-state NMR spectra and (1)H chemical shift calculation with gauge-including projector augmented wave calculations.

The impact of influenza hemagglutinin fusion peptide length and viral subtype on its structure and dynamics

A peptide comprising no fewer than the first 20 residues of the influenza hemagglutinin HA2 subunit suffices to induce lipid mixing between the membranes of different unilamellar vesicles. This 20-residue peptide was previously reported to adopt an open "boomerang" structure that differs significantly from the closed helical-hairpin structure of a fusion peptide consisting of the first 23 residues of the HA2 sequence. This study investigates the structural and dynamic features of fusion peptides of different length and subtype. Lacking key interactions that stabilize the closed, helical-hairpin structure, the 20-residue peptide is in a dynamic equilibrium between closed and open states, adopting a ca. 11% population of the former when solubilized by DPC micelles. Peptides shorter than 20 residues would have even fewer interactions to stabilize a helical hairpin fold, resulting in a vanishing hairpin population. Considering the conserved nature of hairpin-stabilizing interactions across all serotypes, and the minimum of 20 residues needed for fusion, we postulate that the closed state plays an essential role in the fusion process. However, opening of this hairpin structure may be essential to the formation of a membrane pore at the final stage of the fusion process.

Very large residual dipolar couplings from deuterated ubiquitin

Main-chain (1)H(N)-(15)N residual dipolar couplings (RDCs) ranging from approximately -200 to 200 Hz have been measured for ubiquitin under strong alignment conditions in Pf1 phage. This represents a ten-fold increase in the degree of alignment over the typical weakly aligned samples. The measurements are made possible by extensive proton-dilution of the sample, achieved by deuteration of the protein with partial back-substitution of labile protons from 25 % H(2)O / 75 % D(2)O buffer. The spectral quality is further improved by application of deuterium decoupling. Since standard experiments using fixed-delay INEPT elements cannot accommodate a broad range of couplings, the measurements were conducted using J-resolved and J-modulated versions of the HSQC and TROSY sequences. Due to unusually large variations in dipolar couplings, the trosy (sharp) and anti-trosy (broad) signals are often found to be interchanged in the TROSY spectra. To distinguish between the two, we have relied on their respective (15)N linewidths. This strategy ultimately allowed us to determine the signs of RDCs. The fitting of the measured RDC values to the crystallographic coordinates of ubiquitin yields the quality factor Q = 0.16, which confirms the perturbation-free character of the Pf1 alignment. Our results demonstrate that RDC data can be successfully acquired not only in dilute liquid crystals, but also in more concentrated ones. As a general rule, the increase in liquid crystal concentration improves the stability of alignment media and makes them more tolerant to variations in sample conditions. The technical ability to measure RDCs under moderately strong alignment conditions may open the door for development of alternative alignment media, including new types of media that mimic biologically relevant systems.

Determination of 15N chemical shift anisotropy from a membrane-bound protein by NMR spectroscopy

Chemical shift anisotropy (CSA) tensors are essential in the structural and dynamic studies of proteins using NMR spectroscopy. Results from relaxation studies in biomolecular solution and solid-state NMR experiments on aligned samples are routinely interpreted using well-characterized CSA tensors determined from model compounds. Since CSA tensors, particularly the (15)N CSA, highly depend on a number of parameters including secondary structure, electrostatic interaction, and the amino acid sequence, there is a need for accurately determined CSA tensors from proteins. In this study, we report the backbone amide-(15)N CSA tensors for a 16.7-kDa membrane-bound and paramagnetic-heme containing protein, rabbit Cytochrome b(5) (cytb(5)), determined using the (15)N CSA/(15)N-(1)H dipolar transverse cross-correlation rates. The mean values of (15)N CSA determined for residues in helical, sheet, and turn regions are -187.9, -166.0, and -161.1 ppm, respectively, with an overall average value of -171.7 ppm. While the average CSA value determined from this study is in good agreement with previous solution NMR experiments on small globular proteins, the CSA value determined for residues in helical conformation is slightly larger, which may be attributed to the paramagnetic effect from Fe(III) of the heme unit in cytb(5). However, like in previous solution NMR studies, the CSA values reported in this study are larger than the values measured from solid-state NMR experiments. We believe that the CSA parameters reported in this study will be useful in determining the structure, dynamics, and orientation of proteins, including membrane proteins, using NMR spectroscopy.

Coupled motion in proteins revealed by pressure perturbation

The cooperative nature of protein substructure and internal motion is a critical aspect of their functional competence about which little is known experimentally. NMR relaxation is used here to monitor the effects of high pressure on fast internal motion in the protein ubiquitin. In contrast to the main chain, the motions of the methyl-bearing side chains have a large and variable pressure dependence. Within the core, this pressure sensitivity correlates with the magnitude of motion at ambient pressure. Spatial clustering of the dynamic response to applied hydrostatic pressure is also seen, indicating localized cooperativity of motion on the sub-nanosecond time scale and suggesting regions of variable compressibility. These and other features indicate that the native ensemble contains a significant fraction of members with characteristics ascribed to the recently postulated "dry molten globule". The accompanying variable side-chain conformational entropy helps complete our view of the thermodynamic architecture underlying protein stability, folding, and function.

Calculation of chemical shift anisotropy in proteins

Individual peptide groups in proteins must exhibit some variation in the chemical shift anisotropy (CSA) of their constituent atoms, but not much is known about the extent or origins of this dispersion. Direct spectroscopic measurement of CSA remains technically challenging, and theoretical methods can help to overcome these limitations by estimating shielding tensors for arbitrary structures. Here we use an automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach to compute (15)N, (13)C' and (1)H chemical shift tensors for human ubiquitin and the GB1 and GB3 fragments of staphylococcal protein G. The average and range of variation of the anisotropies is in good agreement with experimental estimates from solid-state NMR, and the variation among residues is somewhat smaller than that estimated from solution-state measurements. Hydrogen-bond effects account for much of the variation, both between helix and sheet regions, and within elements of secondary structure, but other effects (including variations in torsion angles) may play a role as well.

Ultrahigh resolution protein structures using NMR chemical shift tensors

NMR chemical shift tensors (CSTs) in proteins, as well as their orientations, represent an important new restraint class for protein structure refinement and determination. Here, we present the first determination of both CST magnitudes and orientations for (13)Cα and (15)N (peptide backbone) groups in a protein, the β1 IgG binding domain of protein G from Streptococcus spp., GB1. Site-specific (13)Cα and (15)N CSTs were measured using synchronously evolved recoupling experiments in which (13)C and (15)N tensors were projected onto the (1)H-(13)C and (1)H-(15)N vectors, respectively, and onto the (15)N-(13)C vector in the case of (13)Cα. The orientations of the (13)Cα CSTs to the (1)H-(13)C and (13)C-(15)N vectors agreed well with the results of ab initio calculations, with an rmsd of approximately 8°. In addition, the measured (15)N tensors exhibited larger reduced anisotropies in α-helical versus β-sheet regions, with very limited variation (18 ± 4°) in the orientation of the z-axis of the (15)N CST with respect to the (1)H-(15)N vector. Incorporation of the (13)Cα CST restraints into structure calculations, in combination with isotropic chemical shifts, transferred echo double resonance (13)C-(15)N distances and vector angle restraints, improved the backbone rmsd to 0.16 Å (PDB ID code 2LGI) and is consistent with existing X-ray structures (0.51 Å agreement with PDB ID code 2QMT). These results demonstrate that chemical shift tensors have considerable utility in protein structure refinement, with the best structures comparable to 1.0-Å crystal structures, based upon empirical metrics such as Ramachandran geometries and χ(1)/χ(2) distributions, providing solid-state NMR with a powerful tool for de novo structure determination.

Dynamics of biomolecules from picoseconds to seconds at atomic resolution

Although biomolecular dynamics has been investigated using NMR for at least 40 years, only in the past 20 years have internal motions been characterized at atomic resolution throughout proteins and nucleic acids. This development was made possible by multidimensional heteronuclear NMR approaches that provide near complete sequential signal assignments of uniformly labeled biomolecules. Recent methodological advances have enabled characterization of internal dynamics on timescales ranging from picoseconds to seconds, both in solution and in the solid state. The size, complexity and functional significance of biomolecules investigated by NMR continue to grow, as do the insights that have been obtained about function. In this article I review a number of recent advances that have made such studies possible, and provide a few examples of where NMR either by itself or in combination with other approaches has paved the way to a better understanding of the complex relationship between dynamics and biomolecular function. Finally, I discuss prospects for further advances in this field.