Interpretation of biomolecular NMR spin relaxation parameters

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.

Dynamic 15N{1H} NOE measurements: a tool for studying protein dynamics

Intramolecular motions in proteins are one of the important factors that determine their biological activity and interactions with molecules of biological importance. Magnetic relaxation of 15N amide nuclei allows one to monitor motions of protein backbone over a wide range of time scales. 15N{1H} nuclear Overhauser effect is essential for the identification of fast backbone motions in proteins. Therefore, exact measurements of NOE values and their accuracies are critical for determining the picosecond time scale of protein backbone. Measurement of dynamic NOE allows for the determination of NOE values and their probable errors defined by any sound criterion of nonlinear regression methods. The dynamic NOE measurements can be readily applied for non-deuterated or deuterated proteins in both HSQC and TROSY-type experiments. Comparison of the dynamic NOE method with commonly implied steady-state NOE is presented in measurements performed at three magnetic field strengths. It is also shown that improperly set NOE measurement cannot be restored with correction factors reported in the literature.

The dormancy-specific regulator, SutA, is intrinsically disordered and modulates transcription initiation in Pseudomonas aeruginosa

Though most bacteria in nature are nutritionally limited and grow slowly, our understanding of core processes like transcription comes largely from studies in model organisms doubling rapidly. We previously identified a small protein of unknown function, SutA, in a screen of proteins synthesized in Pseudomonas aeruginosa during dormancy. SutA binds RNA polymerase (RNAP), causing widespread changes in gene expression, including upregulation of the ribosomal RNA genes. Here, using biochemical and structural methods, we examine how SutA interacts with RNAP and the functional consequences of these interactions. We show that SutA comprises a central α-helix with unstructured N- and C-terminal tails, and binds to the β1 domain of RNAP. It activates transcription from the rrn promoter by both the housekeeping sigma factor holoenzyme (Eσ⁷⁰ ) and the stress sigma factor holoenzyme (Eσ^S ) in vitro, but has a greater impact on Eσ^S . In both cases, SutA appears to affect intermediates in the open complex formation and its N-terminal tail is required for activation. The small magnitudes of in vitro effects are consistent with a role in maintaining activity required for homeostasis during dormancy. Our results add SutA to a growing list of transcription regulators that use their intrinsically disordered regions to remodel transcription complexes.

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.

Preserved Transmembrane Segment Topology, Structure, and Dynamics in Disparate Micellar Environments

Detergent micelles are frequently employed as membrane mimetics for solution-state membrane protein nuclear magnetic resonance spectroscopy. Here we compare topology, structure, ps-ns time-scale dynamics, and hydrodynamics of a model protein with one transmembrane (TM) segment (residues 1-55 of the apelin receptor, APJ, a G-protein-coupled receptor) in three distinct, commonly used micellar environments. In each environment, two solvent-protected helical segments connected by a solvent-exposed kink were observed. The break in helical character at the kink was maintained in a helix-stabilizing fluorinated alcohol environment, implying that this structural feature is inherent. Molecular dynamics simulations also substantiate favorable self-assembly of compact protein-micelle complexes with a more dynamic, solvent-exposed kink. Despite the observed similarity in TM segment behavior, micelle-dependent differences were clear in the structure, dynamics, and compactness of the 30-residue, extramembrane N-terminal tail of the protein. This would affect intermolecular interactions and, correspondingly, the functional state of the membrane protein.

General order parameter based correlation analysis of protein backbone motions between experimental NMR relaxation measurements and molecular dynamics simulations

Internal backbone dynamic motions are essential for different protein functions and occur on a wide range of time scales, from femtoseconds to seconds. Molecular dynamic (MD) simulations and nuclear magnetic resonance (NMR) spin relaxation measurements are valuable tools to gain access to fast (nanosecond) internal motions. However, there exist few reports on correlation analysis between MD and NMR relaxation data. Here, backbone relaxation measurements of (15)N-labeled SH3 (Src homology 3) domain proteins in aqueous buffer were used to generate general order parameters (S(2)) using a model-free approach. Simultaneously, 80 ns MD simulations of SH3 domain proteins in a defined hydrated box at neutral pH were conducted and the general order parameters (S(2)) were derived from the MD trajectory. Correlation analysis using the Gromos force field indicated that S(2) values from NMR relaxation measurements and MD simulations were significantly different. MD simulations were performed on models with different charge states for three histidine residues, and with different water models, which were SPC (simple point charge) water model and SPC/E (extended simple point charge) water model. S(2) parameters from MD simulations with charges for all three histidines and with the SPC/E water model correlated well with S(2) calculated from the experimental NMR relaxation measurements, in a site-specific manner.

NMR-based dynamics of free glycosaminoglycans in solution

Dynamical behaviors of glycosaminoglycans, as here illustrated with a hyaluronan oligosaccharide, are key regulators of biological functions.

Structure and dynamic properties of membrane proteins using NMR

Integral membrane proteins are one of the most challenging groups of macromolecules despite their apparent conformational simplicity. They manage and drive transport, circulate information, and participate in cellular movements via interactions with other proteins and through intricate conformational changes. Their structural and functional decoding is challenging and has imposed demanding experimental development. Solution nuclear magnetic resonance (NMR) spectroscopy is one of the techniques providing the capacity to make a significant difference in the deciphering of the membrane protein structure-function paradigm. The method has evolved dramatically during the last decade resulting in a plethora of new experiments leading to a significant increase in the scientific repertoire for studying membrane proteins. Besides solving the three-dimensional structures using state-of-the-art approaches, a large variety of developments of well-established techniques are available providing insight into membrane protein flexibility, dynamics, and interactions. Inspired by the speed of development in the application of new strategies, by invention of methods to measure solvent accessibility and describe low-populated states, this review seeks to introduce the vast possibilities solution NMR can offer to the study of membrane protein structure-function analyses with special focus on applicability.

The hepatitis B virus preS1 domain hijacks host trafficking proteins by motif mimicry

Hepatitis B virus (HBV) is an infectious, potentially lethal human pathogen. However, there are no effective therapies for chronic HBV infections. Antiviral development is hampered by the lack of high-resolution structures for essential HBV protein-protein interactions. The interaction between preS1, an HBV surface-protein domain, and its human binding partner, γ2-adaptin, subverts the membrane-trafficking apparatus to mediate virion export. This interaction is a putative drug target. We report here atomic-resolution descriptions of the binding thermodynamics and structural biology of the interaction between preS1 and the EAR domain of γ2-adaptin. NMR, protein engineering, X-ray crystallography and MS showed that preS1 contains multiple γ2-EAR-binding motifs that mimic the membrane-trafficking motifs (and binding modes) of host proteins. These motifs localize together to a relatively rigid, functionally important region of preS1, an intrinsically disordered protein. The preS1-γ2-EAR interaction was relatively weak and efficiently outcompeted by a synthetic peptide. Our data provide the structural road map for developing peptidomimetic antivirals targeting the γ2-EAR-preS1 interaction.

A monomeric TIM-barrel structure from Pyrococcus furiosus is optimized for extreme temperatures

The structure of phosphoribosyl anthranilate isomerase (TrpF) from the hyperthermophilic archaeon Pyrococcus furiosus (PfTrpF) has been determined at 1.75 Å resolution. The PfTrpF structure has a monomeric TIM-barrel fold which differs from the dimeric structures of two other known thermophilic TrpF proteins. A comparison of the PfTrpF structure with the two known bacterial thermophilic TrpF structures and the structure of a related mesophilic protein from Escherichia coli (EcTrpF) is presented. The thermophilic TrpF structures contain a higher proportion of ion pairs and charged residues compared with the mesophilic EcTrpF. These residues contribute to the closure of the central barrel and the stabilization of the barrel and the surrounding α-helices. In the monomeric PfTrpF conserved structural water molecules are mostly absent; instead, the structural waters are replaced by direct side-chain-main-chain interactions. As a consequence of these combined mechanisms, the P. furiosus enzyme is a thermodynamically stable and entropically optimized monomeric TIM-barrel enzyme which defines a good framework for further protein engineering for industrial applications.

Guanine nucleotides differentially modulate backbone dynamics of the STAS domain of the SulP/SLC26 transport protein Rv1739c of Mycobacterium tuberculosis

Enzymatic catalysis and protein signaling are dynamic processes that involve local and/or global conformational changes occurring across a broad range of time scales. (1) H-(15) N relaxation NMR provides a comprehensive understanding of protein backbone dynamics both in the apo (unliganded) and ligand-bound conformations, enabling both fast and slow internal motions of individual amino acid residues to be observed. We recently reported the structure and nucleotide binding properties of the sulfate transporter and anti-sigma factor antagonist (STAS) domain of Rv1739c, a SulP anion transporter protein of Mycobacterium tuberculosis. In the present study, we report (1) H-(15) N NMR backbone dynamics measurements [longitudinal (T(1) ), transverse (T(2) ) and steady-state ({(1) H}-(15) N) heteronuclear NOE] of the Rv1739c STAS domain, in the absence and presence of saturating concentrations of GTP and GDP. Analysis of measured relaxation data and estimated dynamic parameters indicated distinct features differentiating the binding of GTP and GDP to Rv1739c STAS. The 9.55 ns overall rotational correlation time of Rv1739c STAS increased to 10.48 ns in the presence of GTP, and to 13.25 ns in the presence of GDP, indicating significant nucleotide-induced conformational changes. These conformational changes were accompanied by slow time scale (μs to ms) motions in discrete regions of the protein, as reflected by guanine nucleotide-induced changes in relaxation parameters. The observed nucleotide-specific alterations in the relaxation properties of individual STAS residues reflect an increased molecular anisotropy and/or the emergence of conformational equilibria governing functional properties of the STAS domain.

Structure and dynamics of Mycobacterium tuberculosis truncated hemoglobin N: insights from NMR spectroscopy and molecular dynamics simulations

The potent nitric oxide dioxygenase (NOD) activity (trHbN-Fe²⁺-O₂ + (•)NO → trHbN-Fe³⁺-OH₂ + NO₃⁻) of Mycobacterium tuberculosis truncated hemoglobin N (trHbN) protects aerobic respiration from inhibition by (•)NO. The high activity of trHbN has been attributed in part to the presence of numerous short-lived hydrophobic cavities that allow partition and diffusion of the gaseous substrates (•)NO and O₂ to the active site. We investigated the relation between these cavities and the dynamics of the protein using solution NMR spectroscopy and molecular dynamics (MD). Results from both approaches indicate that the protein is mainly rigid with very limited motions of the backbone N-H bond vectors on the picoseconds-nanoseconds time scale, indicating that substrate diffusion and partition within trHbN may be controlled by side-chains movements. Model-free analysis also revealed the presence of slow motions (microseconds-milliseconds), not observed in MD simulations, for many residues located in helices B and G including the distal heme pocket Tyr33(B10). All currently known crystal structures and molecular dynamics data of truncated hemoglobins with the so-called pre-A N-terminal extension suggest a stable α-helical conformation that extends in solution. Moreover, a recent study attributed a crucial role to the pre-A helix for NOD activity. However, solution NMR data clearly show that in near-physiological conditions these residues do not adopt an α-helical conformation and are significantly disordered and that the helical conformation seen in crystal structures is likely induced by crystal contacts. Although this lack of order for the pre-A does not disagree with an important functional role for these residues, our data show that one should not assume an helical conformation for these residues in any functional interpretation. Moreover, future molecular dynamics simulations should not use an initial α-helical conformation for these residues in order to avoid a bias based on an erroneous initial structure for the N-termini residues. This work constitutes the first study of a truncated hemoglobin dynamics performed by solution heteronuclear relaxation NMR spectroscopy.

Increased hydrophobicity and decreased backbone flexibility explain the lower solubility of a cataract-linked mutant of γD-crystallin

A number of point mutations in γD-crystallin are associated with human cataract. The Pro23-to-Thr (P23T) mutation is perhaps the most common, is geographically widespread, and presents itself in a variety of phenotypes. It is therefore important to understand the molecular basis of lens opacity due to this mutation. In our earlier studies, we noted that P23T shows retrograde and sharply lowered solubility, most likely due to the emergence of hydrophobic patches involved in protein aggregation. Binding of 4,4'-dianilino-1,1'-binaphthyl-5,5'-disulfonate (Bis-ANS) dye (a probe commonly used for detecting surface hydrophobicity) competed with aggregation, suggesting that the residues involved in Bis-ANS binding are also involved in protein aggregation. Here, using NMR spectroscopy in conjunction with Bis-ANS binding, we identify three residues (Y16, D21, and Y50) in P23T that are involved in binding the dye. Furthermore, using (15)N NMR relaxation experiments, we show that, in the mutant protein, backbone fluctuations are restricted to the picosecond-to-nanosecond and microsecond timescales relative to the wild type. Our present studies specify the residues involved in these two pivotal characteristics of the mutant protein, namely increased surface hydrophobicity and restricted mobility of the protein backbone, which can explain the nucleation and further propagation of protein aggregates. Thus, we have now identified the residues in the P23T mutant that give rise to novel hydrophobic surfaces, as well as those regions of the protein backbone where fluctuations in different timescales are restricted, providing a comprehensive understanding of how lens opacity could result from this mutation.

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.