An optimized method for (15)N R(1) relaxation rate measurements in non-deuterated proteins

Dynamics of Globular Proteins when Interacting with Zwitterionic Silica Nanoparticles by Nuclear Magnetic Resonance Spin Relaxation

The many emerging applications of nanoparticles in diverse fields in chemistry and biology require the characterization of interactions between nanoparticles and surrounding biomolecules such as proteins. Nuclear magnetic resonance spin relaxation of proteins, which is highly sensitive to interactions with nanoparticles, contains rich information about protein mobility and binding kinetics. The interactions of globular proteins with silica nanoparticles differ markedly from those with liposome nanoparticles, although both are driven by electrostatic forces. For unmodified silica nanoparticles, their interactions with an internally rigid protein like ubiquitin uniformly increase the backbone amide 15N transverse R2 relaxation for most residues. In contrast, for ubiquitin-POPG liposome interactions, their characteristic transverse R2 profiles indicate that ubiquitin undergoes diffusive rotational motions on the liposome surface. Here, we show that coating silica nanoparticles with sulfobetaine siloxane zwitterionic molecules profoundly alters their interactions with proteins in a manner that closely resembles the mode of interaction observed with liposomes. 15N-R2 relaxation reveals that ubiquitin and the Ras-binding domain of B-Raf both exhibit axial reorientational motions about an axis perpendicular to the nanoparticle surface in the bound state, where the interactions involve predominantly positively charged surface regions. These findings point toward a global dynamics mechanism of proteins when interacting with organic or inorganic nanoparticles with densely charged soft surfaces. This information may help tailor the coatings of nanoparticles to adopt specific modes of interaction with proteins that can be used to control their function in vivo and in vitro.

Pitfalls in measurements of R1 relaxation rates of protein backbone 15N nuclei

The dynamics of the backbone and side-chains of protein are routinely studied by interpreting experimentally determined 15N spin relaxation rates. R1(15N), the longitudinal relaxation rate, reports on fast motions and encodes, together with the transverse relaxation R2, structural information about the shape of the molecule and the orientation of the amide bond vectors in the internal diffusion frame. Determining error-free 15N longitudinal relaxation rates remains a challenge for small, disordered, and medium-sized proteins. Here, we show that mono-exponential fitting is sufficient, with no statistical preference for bi-exponential fitting up to 800 MHz. A detailed comparison of the TROSY and HSQC techniques at medium and high fields showed no statistically significant differences. The least error-prone DD/CSA interference removal technique is the selective inversion of amide signals while avoiding water resonance. The exchange of amide with solvent deuterons appears to affect the rate R1 of solvent-exposed amides in all fields tested and in each DD/CSA interference removal technique in a statistically significant manner. In summary, the most accurate R1(15N) rates in proteins are achieved by selective amide inversion, without the addition of D2O. Importantly, at high magnetic fields stronger than 800 MHz, when non-mono-exponential decay is involved, it is advisable to consider elimination of the shortest delays (typically up to 0.32 s) or bi-exponential fitting.

The sequence-structure-function relationship of intrinsic ERα disorder

The oestrogen receptor (ER or ERα), a nuclear hormone receptor that drives most breast cancer1, is commonly activated by phosphorylation at serine 118 within its intrinsically disordered N-terminal transactivation domain2,3. Although this modification enables oestrogen-independent ER function, its mechanism has remained unclear despite ongoing clinical trials of kinase inhibitors targeting this region4-6. By integration of small-angle X-ray scattering and nuclear magnetic resonance spectroscopy with functional studies, we show that serine 118 phosphorylation triggers an unexpected expansion of the disordered domain and disrupts specific hydrophobic clustering between two aromatic-rich regions. Mutations mimicking this disruption rescue ER transcriptional activity, target-gene expression and cell growth impaired by a phosphorylation-deficient S118A mutation. These findings, driven by hydrophobic interactions, extend beyond electrostatic models and provide mechanistic insights into intrinsically disordered proteins7, with implications for other nuclear receptors8. This fundamental sequence-structure-function relationship advances our understanding of intrinsic ER disorder, crucial for developing targeted breast cancer therapeutics.

Modulation of peroxisomal import by the PEX13 SH3 domain and a proximal FxxxF binding motif

Import of proteins into peroxisomes depends on PEX5, PEX13 and PEX14. By combining biochemical methods and structural biology, we show that the C-terminal SH3 domain of PEX13 mediates intramolecular interactions with a proximal FxxxF motif. The SH3 domain also binds WxxxF peptide motifs in the import receptor PEX5, demonstrating evolutionary conservation of such interactions from yeast to human. Strikingly, intramolecular interaction of the PEX13 FxxxF motif regulates binding of PEX5 WxxxF/Y motifs to the PEX13 SH3 domain. Crystal structures reveal how FxxxF and WxxxF/Y motifs are recognized by a non-canonical surface on the SH3 domain. The PEX13 FxxxF motif also mediates binding to PEX14. Surprisingly, the potential PxxP binding surface of the SH3 domain does not recognize PEX14 PxxP motifs, distinct from its yeast ortholog. Our data show that the dynamic network of PEX13 interactions with PEX5 and PEX14, mediated by diaromatic peptide motifs, modulates peroxisomal matrix import.

Sensitivity-enhanced NMR 15N R1 and R1ρ relaxation experiments for the investigation of intrinsically disordered proteins at high magnetic fields

NMR relaxation experiments provide residue-specific insights into the structural dynamics of proteins. Here, we present an optimized set of sensitivity-enhanced 15N R1 and R1ρ relaxation experiments applicable to fully protonated proteins. The NMR pulse sequences are conceptually similar to the set of TROSY-based sequences and their HSQC counterpart (Lakomek et al., J. Biomol. NMR 2012). Instead of the TROSY read-out scheme, a sensitivity-enhanced HSQC read-out scheme is used, with improved and easier optimized water suppression. The presented pulse sequences are applied on the cytoplasmic domain of the SNARE protein Synpatobrevin-2 (Syb-2), which is intrinsically disordered in its monomeric pre-fusion state. A two-fold increase in the obtained signal-to-noise ratio is observed for this intrinsically disordered protein, therefore offering a four-fold reduction of measurement time compared to the TROSY-detected version. The inter-scan recovery delay can be shortened to two seconds. Pulse sequences were tested at 600 MHz and 1200 MHz 1H Larmor frequency, thus applicable over a wide magnetic field range. A comparison between protonated and deuterated protein samples reveals high agreement, indicating that reliable 15N R1 and R1ρ rate constants can be extracted for fully protonated and deuterated samples. The presented pulse sequences will benefit not only for IDPs but also for an entire range of low and medium-sized proteins.

Excited-state observation of active K-Ras reveals differential structural dynamics of wild-type versus oncogenic G12D and G12C mutants

Despite the prominent role of the K-Ras protein in many different types of human cancer, major gaps in atomic-level information severely limit our understanding of its functions in health and disease. Here, we report the quantitative backbone structural dynamics of K-Ras by solution nuclear magnetic resonance spectroscopy of the active state of wild-type K-Ras bound to guanosine triphosphate (GTP) nucleotide and two of its oncogenic P-loop mutants, G12D and G12C, using a new nanoparticle-assisted spin relaxation method, relaxation dispersion and chemical exchange saturation transfer experiments covering the entire range of timescales from picoseconds to milliseconds. Our combined experiments allow detection and analysis of the functionally critical Switch I and Switch II regions, which have previously remained largely unobservable by X-ray crystallography and nuclear magnetic resonance spectroscopy. Our data reveal cooperative transitions of K-Ras·GTP to a highly dynamic excited state that closely resembles the partially disordered K-Ras·GDP state. These results advance our understanding of differential GTPase activities and signaling properties of the wild type versus mutants and may thus guide new strategies for the development of therapeutics.

Elucidating the concentration-dependent effects of thiocyanate binding to carbonic anhydrase

Many proteins naturally carry metal centers, with a large share of them being in the active sites of several enzymes. Paramagnetic effects are a powerful source of structural information and, therefore, if the native metal is paramagnetic, or it can be functionally substituted with a paramagnetic one, paramagnetic effects can be used to study the metal sites, as well as the overall structure of the protein. One notable example is cobalt(II) substitution for zinc(II) in carbonic anhydrase. In this manuscript we investigate the effects of sodium thiocyanate on the chemical environment of the metal ion of the human carbonic anhydrase II. The electron paramagnetic resonance (EPR) titration of the cobalt(II) protein with thiocyanate shows that the EPR spectrum changes from A-type to C-type on passing from 1:1 to 1:1000-fold ligand excess. This indicates the occurrence of a change in the electronic structure, which may reflect a sizable change in the metal coordination environment in turn caused by a modification of the frozen solvent glass. However, paramagnetic nuclear magnetic resonance (NMR) data indicate that the metal coordination cage remains unperturbed even in 1:1000-fold ligand excess. This result proves that the C-type EPR spectrum observed at large ligand concentration should be ascribed to the low temperature at which EPR measurements are performed, which impacts on the structure of the protein when it is destabilized by a high concentration of a chaotropic agent.

Histone H4 Tails in Nucleosomes: a Fuzzy Interaction with DNA

The interaction of positively charged N-terminal histone tails with nucleosomal DNA plays an important role in chromatin assembly and regulation, modulating their susceptibility to post-translational modifications and recognition by chromatin-binding proteins. Here, we report residue-specific ¹⁵ N NMR relaxation rates for histone H4 tails in reconstituted nucleosomes. These data indicate that H4 tails are strongly dynamically disordered, albeit with reduced conformational flexibility compared to a free peptide with the same sequence. Remarkably, the NMR observables were successfully reproduced in a 2-μs MD trajectory of the nucleosome. This is an important step toward resolving an apparent inconsistency where prior simulations were generally at odds with experimental evidence on conformational dynamics of histone tails. Our findings indicate that histone H4 tails engage in a fuzzy interaction with nucleosomal DNA, underpinned by a variable pattern of short-lived salt bridges and hydrogen bonds, which persists at low ionic strength (0-100 mM NaCl).

Functional protein dynamics on uncharted time scales detected by nanoparticle-assisted NMR spin relaxation

Protein function depends critically on intrinsic internal dynamics, which is manifested in distinct ways, such as loop motions that regulate protein recognition and catalysis. Under physiological conditions, dynamic processes occur on a wide range of time scales from subpicoseconds to seconds. Commonly used NMR spin relaxation in solution provides valuable information on very fast and slow motions but is insensitive to the intermediate nanosecond to microsecond range that exceeds the protein tumbling correlation time. Presently, very little is known about the nature and functional role of these motions. It is demonstrated here how transverse spin relaxation becomes exquisitely sensitive to these motions at atomic resolution when studying proteins in the presence of nanoparticles. Application of this novel cross-disciplinary approach reveals large-scale dynamics of loops involved in functionally critical protein-protein interactions and protein-calcium ion recognition that were previously unobservable.

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.

TROSY pulse sequence for simultaneous measurement of the 15N R1 and {1H}-15N NOE in deuterated proteins

A TROSY-based NMR experiment is described for simultaneous measurement of the ¹⁵N longitudinal relaxation rate constant R₁ and the {¹H}-¹⁵N nuclear Overhauser enhancement. The experiment is based on the observation that the TROSY mixing pulse sequence element symmetrically exchanges ¹H and ¹⁵N magnetizations. The accuracy of the proposed technique is validated by comparison to independent measurements of both relaxation parameters for the protein ubiquitin. The simultaneous experiment is approximately 20-33% shorter than conventional sequential measurements.

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.

Rapid Determination of Fast Protein Dynamics from NMR Chemical Exchange Saturation Transfer Data

Functional motions of (15)N-labeled proteins can be monitored by solution NMR spin relaxation experiments over a broad range of timescales. These experiments however typically take of the order of several days to a week per protein. Recently, NMR chemical exchange saturation transfer (CEST) experiments have emerged to probe slow millisecond motions complementing R1ρ and CPMG-type experiments. CEST also simultaneously reports on site-specific R1 and R2 parameters. It is shown here how CEST-derived R1 and R2 relaxation parameters can be measured within a few hours at an accuracy comparable to traditional relaxation experiments. Using a "lean" version of the model-free approach S(2) order parameters can be determined that match those from the standard model-free approach applied to (15)N R1, R2 , and {(1)H}-(15)N NOE data. The new methodology, which is demonstrated for ubiquitin and arginine kinase (42 kDa), should serve as an effective screening tool of protein dynamics from picosecond-to-millisecond timescales.