Efficient and accurate determination of the overall rotational diffusion tensor of a molecule from (15)N relaxation data using computer program ROTDIF

Selective [9-15N] Guanosine for Nuclear Magnetic Resonance Studies of Large Ribonucleic Acids

RNAs regulate various cellular processes using malleable 3D structures, and understanding the factors that control RNA structure and dynamics is critical for understanding their mechanisms of action. To mitigate factors that have limited studies of large, functionally relevant RNAs by solution nuclear magnetic resonance (NMR) spectroscopy, we have extended a recently described 2H-enhanced, 1H-15N correlation approach that used uniformly 15N-labeled guanosine triphosphate (GTP) by developing a chemoenzymatic labeling technology that grafts selectively labeled [9-15N]-Guanine on to any labeled ribose to make [9-15N]-GTP. The approach exploits advantageous NMR properties of the N9 nucleus which, when combined with extensive ribose deuteration and optimized NMR pulse sequences, affords sharp signals without complications that can arise using uniform [15N]-guanine labeling. The utility of the approach for NMR signal assignment and dynamics analysis is demonstrated for three large RNAs (20-78 kDa) that play critical roles in viral replication. With this approach, NMR studies of RNAs comprising 200 nt or more should now be feasible.

Molecular dynamics as an efficient process to predict 15N chemical shift anisotropy at very high NMR magnetic fields

The emergence of very high NMR magnetic fields will certainly encourage the study of larger biological systems with their dynamics and interactions. NMR spin relaxation allows probing the dynamical properties of proteins where the 15N longitudinal (R1) and transverse (R2) relaxation rates in addition to the 1H-15N heteronuclear NOE describe the ps-ns time scale. Their analytical representation involves the chemical shift anisotropy (CSA) effect that represents the major contribution at a very high magnetic field above 18.8 T. An accurate analysis of the latter parameters in terms of model free (MF) requires considering its effect. Until now, a uniform value of -160 ppm for the CSA has been widely used to derive the backbone order parameters (S2), giving rise to a large fluctuation of its value at very high magnetic fields. Conversely, the use of a site-specific CSA improves the accurate analysis of protein dynamics but requires a cost-effective experimental multi-field approach. In the present paper, we show how the CSA mainly contributes to the relaxation parameters at 28.2 T compared to lower magnetic fields and may bias the determination of S2. We propose to replace the time-consuming measurement of spin relaxation at multiple fields by a combination of molecular dynamics (MD) and the measurement of spin relaxation at one very high magnetic field only. We applied this strategy to three well-folded proteins (ubiquitin, GB3 and ribonuclease H) to show that the determined order parameters are in good agreement with the ones obtained by means of experimental data only.

NMR Dynamic View of the Stabilization of the WW4 Domain by Neutral NaCl and Kosmotropic Na2SO4 and NaH2PO4

The Hofmeister series categorizes ions based on their effects on protein stability, yet the microscopic mechanism remains a mystery. In this series, NaCl is neutral, Na2SO4 and Na2HPO4 are kosmotropic, while GdmCl and NaSCN are chaotropic. This study employs CD and NMR to investigate the effects of NaCl, Na2SO4, and Na2HPO4 on the conformation, stability, binding, and backbone dynamics (ps-ns and µs-ms time scales) of the WW4 domain with a high stability and accessible side chains at concentrations ≤ 200 mM. The results indicated that none of the three salts altered the conformation of WW4 or showed significant binding to the four aliphatic hydrophobic side chains. NaCl had no effect on its thermal stability, while Na2SO4 and Na2HPO4 enhanced the stability by ~5 °C. Interestingly, NaCl only weakly interacted with the Arg27 amide proton, whereas Na2SO4 bound to Arg27 and Phe31 amide protons with Kd of 32.7 and 41.6 mM, respectively. Na2HPO4, however, bound in a non-saturable manner to Trp9, His24, and Asn36 amide protons. While the three salts had negligible effects on ps-ns backbone dynamics, NaCl and Na2SO4 displayed no effect while Na2HPO4 significantly increased the µs-ms backbone dynamics. These findings, combined with our recent results with GdmCl and NaSCN, suggest a microscopic mechanism for the Hofmeister series. Additionally, the data revealed a lack of simple correlation between thermodynamic stability and backbone dynamics, most likely due to enthalpy-entropy compensation. Our study rationalizes the selection of chloride and phosphate as the primary anions in extracellular and intracellular spaces, as well as polyphosphate as a primitive chaperone in certain single-cell organisms.

A Transient, Excited Species of the Autoinhibited α-State of the Bacterial Transcription Factor RfaH May Represent an Early Intermediate on the Fold-Switching Pathway

RfaH is a two-domain transcription factor in which the C-terminal domain switches fold from an α-helical hairpin to a β-roll upon binding the ops-paused RNA polymerase. To ascertain the presence of a sparsely populated excited state that may prime the autoinhibited resting state of RfaH for binding ops-paused RNA polymerase, we carried out a series of NMR-based exchange experiments to probe for conformational exchange on the millisecond time scale. Quantitative analysis of these data reveals exchange between major ground (∼95%) and sparsely populated excited (∼5%) states with an exchange lifetime of ∼3 ms involving residues at the interface between the N-terminal and C-terminal domains formed by the β3/β4 hairpin and helix α3 of the N-terminal domain and helices α4 and α5 of the C-terminal domain. The largest 15N backbone chemical shift differences are associated with the β3/β4 hairpin, leading us to suggest that the excited state may involve a rigid body lateral displacement/rotation away from the C-terminal domain to adopt a position similar to that seen in the active RNA polymerase-bound state. Such a rigid body reorientation would result in a reduction in the interface between the N- and C-terminal domains with the possible introduction of a cavity or cavities. This hypothesis is supported by the observation that the population of the excited species and the exchange rate of interconversion between ground and excited states are reduced at a high (2.5 kbar) pressure. Mechanistic implications for fold switching of the C-terminal domain in the context of RNA polymerase binding are discussed.

Insights into mechanisms of MALT1 allostery from NMR and AlphaFold dynamic analyses

Mucosa-associated lymphoid tissue lymphoma-translocation protein 1 (MALT1) is an attractive target for the development of modulatory compounds in the treatment of lymphoma and other cancers. While the three-dimensional structure of MALT1 has been previously determined through X-ray analysis, its dynamic behaviour in solution has remained unexplored. We present here dynamic analyses of the apo MALT1 form along with the E549A mutation. This investigation used NMR 15N relaxation and NOE measurements between side-chain methyl groups. Our findings confirm that MALT1 exists as a monomer in solution, and demonstrate that the domains display semi-independent movements in relation to each other. Our dynamic study, covering multiple time scales, along with the assessment of conformational populations by Molecular Dynamic simulations, Alpha Fold modelling and PCA analysis, put the side chain of residue W580 in an inward position, shedding light at potential mechanisms underlying the allosteric regulation of this enzyme.

NMR Dynamic View of the Destabilization of WW4 Domain by Chaotropic GdmCl and NaSCN

GdmCl and NaSCN are two strong chaotropic salts commonly used in protein folding and stability studies, but their microscopic mechanisms remain enigmatic. Here, by CD and NMR, we investigated their effects on conformations, stability, binding and backbone dynamics on ps-ns and µs-ms time scales of a 39-residue but well-folded WW4 domain at salt concentrations ≤200 mM. Up to 200 mM, both denaturants did not alter the tertiary packing of WW4, but GdmCl exerted more severe destabilization than NaSCN. Intriguingly, GdmCl had only weak binding to amide protons, while NaSCN showed extensive binding to both hydrophobic side chains and amide protons. Neither denaturant significantly affected the overall ps-ns backbone dynamics, but they distinctively altered µs-ms backbone dynamics. This study unveils that GdmCl and NaSCN destabilize a protein before the global unfolding occurs with differential binding properties and µs-ms backbone dynamics, implying the absence of a simple correlation between thermodynamic stability and backbone dynamics of WW4 at both ps-ns and µs-ms time scales.

Structural basis of aggregate binding by the AAA+ disaggregase ClpG

Severe heat stress causes massive loss of essential proteins by aggregation, necessitating a cellular activity that rescues aggregated proteins. This activity is executed by ATP-dependent, ring-forming, hexameric AAA+ disaggregases. Little is known about the recognition principles of stress-induced protein aggregates. How can disaggregases specifically target aggregated proteins, while avoiding binding to soluble non-native proteins? Here, we determined by NMR spectroscopy the core structure of the aggregate-targeting N1 domain of the bacterial AAA+ disaggregase ClpG, which confers extreme heat resistance to bacteria. N1 harbors a Zn2+-coordination site that is crucial for structural integrity and disaggregase functionality. We found that conserved hydrophobic N1 residues located on a β-strand are crucial for aggregate targeting and disaggregation activity. Analysis of mixed hexamers consisting of full-length and N1-truncated subunits revealed that a minimal number of four N1 domains must be present in a AAA+ ring for high-disaggregation activity. We suggest that multiple N1 domains increase substrate affinity through avidity effects. These findings define the recognition principle of a protein aggregate by a disaggregase, involving simultaneous contacts with multiple hydrophobic substrate patches located in close vicinity on an aggregate surface. This binding mode ensures selectivity for aggregated proteins while sparing soluble, non-native protein structures from disaggregase activity.

Interaction of the bacterial division regulator MinE with lipid bicelles studied by NMR spectroscopy

The bacterial MinE and MinD division regulatory proteins form a standing wave enabling MinC, which binds MinD, to inhibit FtsZ polymerization everywhere except at the midcell, thereby assuring correct positioning of the cytokinetic septum and even distribution of contents to daughter cells. The MinE dimer undergoes major structural rearrangements between a resting six-stranded state present in the cytoplasm, a membrane-bound state, and a four-stranded active state bound to MinD on the membrane, but it is unclear which MinE motifs interact with the membrane in these different states. Using NMR, we probe the structure and global dynamics of MinE bound to disc-shaped lipid bicelles. In the bicelle-bound state, helix α1 no longer sits on top of the six-stranded β-sheet, losing any contact with the protein core, but interacts directly with the bicelle surface; the structure of the protein core remains unperturbed and also interacts with the bicelle surface via helix α2. Binding may involve a previously identified excited state of free MinE in which helix α1 is disordered, thereby allowing it to target the membrane surface. Helix α1 and the protein core undergo nanosecond rigid body motions of differing amplitudes in the plane of the bicelle surface. Global dynamics on the sub-millisecond time scale between a ground state and a sparsely populated excited state are also observed and may represent a very early intermediate on the transition path between the resting six-stranded and active four-stranded conformations. In summary, our results provide insights into MinE structural rearrangements important during bacterial cell division.

Conformational Dynamics of Phytoglobin BvPgb1.2 from Beta vulgaris ssp. vulgaris

Plant hemoglobins, often referred to as phytoglobins, play important roles in abiotic stress tolerance. Several essential small physiological metabolites can be bound to these heme proteins. In addition, phytoglobins can catalyze a range of different oxidative reactions in vivo. These proteins are often oligomeric, but the degree and relevance of subunit interactions are largely unknown. In this study, we delineate which residues are involved in dimer formation of a sugar beet phytoglobin type 1.2 (BvPgb1.2) using NMR relaxation experiments. E. coli cells harboring a phytoglobin expression vector were cultivated in isotope-labeled (²H, ¹³C and ¹⁵N) M9 medium. The triple-labeled protein was purified to homogeneity using two chromatographic steps. Two forms of BvPgb1.2 were examined, the oxy-form and the more stable cyanide-form. Using three-dimensional triple-resonance NMR experiments, sequence-specific assignments for CN-bound BvPgb1.2 were achieved for 137 backbone amide cross-peaks in the ¹H-¹⁵N TROSY spectrum, which amounts to 83% of the total number of 165 expected cross-peaks. A large proportion of the non-assigned residues are located in α-helixes G and H, which are proposed to be involved in protein dimerization. Such knowledge around dimer formation will be instrumental for developing a better understanding of phytoglobins' roles in planta.

Structure and dynamics of the von Willebrand Factor C6 domain

Von Willebrand disease (VWD) is a bleeding disorder with different levels of severity. VWD-associated mutations are located in the von Willebrand factor (VWF) gene, coding for the large multidomain plasma protein VWF with essential roles in hemostasis and thrombosis. On the one hand, a variety of mutations in the C-domains of VWF are associated with increased bleeding upon vascular injury. On the other hand, VWF gain-of-function (GOF) mutations in the C4 domain have recently been identified, which induce an increased risk of myocardial infarction. Mechanistic insights into how these mutations affect the molecular behavior of VWF are scarce and holistic approaches are challenging due to the multidomain and multimeric character of this large protein. Here, we determine the structure and dynamics of the C6 domain and the single nucleotide polymorphism (SNP) variant G2705R in C6 by combining nuclear magnetic resonance spectroscopy, molecular dynamics simulations and aggregometry. Our findings indicate that this mutation mostly destabilizes VWF by leading to a more pronounced hinging between both subdomains of C6. Hemostatic parameters of variant G2705R are close to normal under static conditions, but the missense mutation results in a gain-of-function under flow conditions, due to decreased VWF stem stability. Together with the fact that two C4 variants also exhibit GOF characteristics, our data underline the importance of the VWF stem region in VWF's hemostatic activity and the risk of mutation-associated prothrombotic properties in VWF C-domain variants due to altered stem dynamics.

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.

Explicit models of motions to analyze NMR relaxation data in proteins

Nuclear Magnetic Resonance (NMR) is a tool of choice to characterize molecular motions. In biological macromolecules, pico- to nano-second motions, in particular, can be probed by nuclear spin relaxation rates which depend on the time fluctuations of the orientations of spin interaction frames. For the past 40 years, relaxation rates have been successfully analyzed using the Model Free (MF) approach which makes no assumption on the nature of motions and reports on the effective amplitude and time-scale of the motions. However, obtaining a mechanistic picture of motions from this type of analysis is difficult at best, unless complemented with molecular dynamics (MD) simulations. In spite of their limited accuracy, such simulations can be used to obtain the information necessary to build explicit models of motions designed to analyze NMR relaxation data. Here, we present how to build such models, suited in particular to describe motions of methyl-bearing protein side-chains and compare them with the MF approach. We show on synthetic data that explicit models of motions are more robust in the presence of rotamer jumps which dominate the relaxation in methyl groups of protein side-chains. We expect this work to motivate the use of explicit models of motion to analyze MD and NMR data.

Structure analysis suggests Ess1 isomerizes the carboxy-terminal domain of RNA polymerase II via a bivalent anchoring mechanism

Accurate gene transcription in eukaryotes depends on isomerization of serine-proline bonds within the carboxy-terminal domain (CTD) of RNA polymerase II. Isomerization is part of the "CTD code" that regulates recruitment of proteins required for transcription and co-transcriptional RNA processing. Saccharomyces cerevisiae Ess1 and its human ortholog, Pin1, are prolyl isomerases that engage the long heptad repeat (YSPTSPS)₂₆ of the CTD by an unknown mechanism. Here, we used an integrative structural approach to decipher Ess1 interactions with the CTD. Ess1 has a rigid linker between its WW and catalytic domains that enforces a distance constraint for bivalent interaction with the ends of long CTD substrates (≥4-5 heptad repeats). Our binding results suggest that the Ess1 WW domain anchors the proximal end of the CTD substrate during isomerization, and that linker divergence may underlie evolution of substrate specificity.

Gain-of-Function Variant p.Pro2555Arg of von Willebrand Factor Increases Aggregate Size through Altering Stem Dynamics

The multimeric plasma glycoprotein (GP) von Willebrand factor (VWF) is best known for recruiting platelets to sites of injury during primary hemostasis. Generally, mutations in the VWF gene lead to loss of hemostatic activity and thus the bleeding disorder von Willebrand disease. By employing cone and platelet aggregometry and microfluidic assays, we uncovered a platelet GPIIb/IIIa-dependent prothrombotic gain of function (GOF) for variant p.Pro2555Arg, located in the C4 domain, leading to an increase in platelet aggregate size. We performed complementary biophysical and structural investigations using circular dichroism spectra, small-angle X-ray scattering, nuclear magnetic resonance spectroscopy, molecular dynamics simulations on the single C4 domain, and dimeric wild-type and p.Pro2555Arg constructs. C4-p.Pro2555Arg retained the overall structural conformation with minor populations of alternative conformations exhibiting increased hinge flexibility and slow conformational exchange. The dimeric protein becomes disordered and more flexible. Our data suggest that the GOF does not affect the binding affinity of the C4 domain for GPIIb/IIIa. Instead, the increased VWF dimer flexibility enhances temporal accessibility of platelet-binding sites. Using an interdisciplinary approach, we revealed that p.Pro2555Arg is the first VWF variant, which increases platelet aggregate size and shows a shear-dependent function of the VWF stem region, which can become hyperactive through mutations. Prothrombotic GOF variants of VWF are a novel concept of a VWF-associated pathomechanism of thromboembolic events, which is of general interest to vascular health but not yet considered in diagnostics. Thus, awareness should be raised for the risk they pose. Furthermore, our data implicate the C4 domain as a novel antithrombotic drug target.

Molecular recognition of ubiquitin and Lys63-linked diubiquitin by STAM2 UIM-SH3 dual domain: the effect of its linker length and flexibility

Multidomain proteins represent a broad spectrum of the protein landscape and are involved in various interactions. They could be considered as modular building blocks assembled in distinct fashion and connected by linkers of varying lengths and sequences. Due to their intrinsic flexibility, these linkers provide proteins a subtle way to modulate interactions and explore a wide range of conformational space. In the present study, we are seeking to understand the effect of the flexibility and dynamics of the linker involved in the STAM2 UIM-SH3 dual domain protein with respect to molecular recognition. We have engineered several constructs of UIM-SH3 with different length linkers or domain deletion. By means of SAXS and NMR experiments, we have shown that the modification of the linker modifies the flexibility and the dynamics of UIM-SH3. Indeed, the global tumbling of both the UIM and SH3 domain is different but not independent from each other while the length of the linker has an impact on the ps-ns time scale dynamics of the respective domains. Finally, the modification of the flexibility and dynamics of the linker has a drastic effect on the interaction of UIM-SH3 with Lys63-linked diubiquitin with a roughly eight-time weaker dissociation constant.

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

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

Structure and dynamics of the platelet integrin-binding C4 domain of von Willebrand factor

Von Willebrand factor (VWF) is a key player in the regulation of hemostasis by promoting recruitment of platelets to sites of vascular injury. An array of 6 C domains forms the dimeric C-terminal VWF stem. Upon shear force activation, the stem adopts an open conformation allowing the adhesion of VWF to platelets and the vessel wall. To understand the underlying molecular mechanism and associated functional perturbations in disease-related variants, knowledge of high-resolution structures and dynamics of C domains is of paramount interest. Here, we present the solution structure of the VWF C4 domain, which binds to the platelet integrin and is therefore crucial for the VWF function. In the structure, we observed 5 intra- and inter-subdomain disulfide bridges, of which 1 is unique in the C4 domain. The structure further revealed an unusually hinged 2-subdomain arrangement. The hinge is confined to a very short segment around V2547 connecting the 2 subdomains. Together with 2 nearby inter-subdomain disulfide bridges, this hinge induces slow conformational changes and positional alternations of both subdomains with respect to each other. Furthermore, the structure demonstrates that a clinical gain-of-function VWF variant (Y2561) is more likely to have an effect on the arrangement of the C4 domain with neighboring domains rather than impairing platelet integrin binding.

Minute Additions of DMSO Affect Protein Dynamics Measurements by NMR Relaxation Experiments through Significant Changes in Solvent Viscosity

Studies of protein-ligand binding often rely on dissolving the ligand in dimethyl sulfoxide (DMSO) to achieve sufficient solubility, and then titrating the ligand solution into the protein solution. As a result, the final protein-ligand solution contains small amounts of DMSO in the buffer. Here we report how the addition of DMSO impacts studies of protein conformational dynamics. We used 15 N NMR relaxation to compare the rotational diffusion correlation time (τC ) of proteins in aqueous buffer with and without DMSO. We found that τC scales with the viscosity of the water-DMSO mixture, which depends sensitively on the amount of DMSO and varies by a factor of 2 across the relevant concentration range. NMR relaxation studies of side chains dynamics are commonly interpreted using τC as a fixed parameter, obtained from backbone 15 N relaxation data acquired on a separate sample. Model-free calculations show that errors in τC , arising from mismatched DMSO concentration between samples, lead to significant errors in order parameters. Our results highlight the importance of determining τC for each sample or carefully matching the DMSO concentrations between samples.

Ubiquitin Modulates Liquid-Liquid Phase Separation of UBQLN2 via Disruption of Multivalent Interactions

Under stress, certain eukaryotic proteins and RNA assemble to form membraneless organelles known as stress granules. The most well-studied stress granule components are RNA-binding proteins that undergo liquid-liquid phase separation (LLPS) into protein-rich droplets mediated by intrinsically disordered low-complexity domains (LCDs). Here we show that stress granules include proteasomal shuttle factor UBQLN2, an LCD-containing protein structurally and functionally distinct from RNA-binding proteins. In vitro, UBQLN2 exhibits LLPS at physiological conditions. Deletion studies correlate oligomerization with UBQLN2's ability to phase-separate and form stress-induced cytoplasmic puncta in cells. Using nuclear magnetic resonance (NMR) spectroscopy, we mapped weak, multivalent interactions that promote UBQLN2 oligomerization and LLPS. Ubiquitin or polyubiquitin binding, obligatory for UBQLN2's biological functions, eliminates UBQLN2 LLPS, thus serving as a switch between droplet and disperse phases. We postulate that UBQLN2 LLPS enables its recruitment to stress granules, where its interactions with ubiquitinated substrates reverse LLPS to enable shuttling of clients out of stress granules.

Ab Initio Prediction of NMR Spin Relaxation Parameters from Molecular Dynamics Simulations

¹H-¹⁵N NMR spin relaxation and relaxation dispersion experiments can reveal the time scale and extent of protein motions across the ps-ms range, where the ps-ns dynamics revealed by fundamental quantities R₁, R₂, and heteronuclear NOE can be well-sampled by molecular dynamics simulations (MD). Although the principles of relaxation prediction from simulations are well-established, numerous NMR-MD comparisons have hitherto focused upon the aspect of order parameters S² due to common artifacts in the prediction of transient dynamics. We therefore summarize here all necessary components and highlight existing and proposed solutions, such as the inclusion of quantum mechanical zero-point vibrational corrections and separate MD convergence of global and local motions in coarse-grained and all-atom force fields, respectively. For the accuracy of the MD prediction to be tested, two model proteins GB3 and Ubiquitin are used to validate five atomistic force fields against published NMR data supplemented by the coarse-grained force field MARTINI+EN. In Amber and CHARMM-type force fields, quantitative agreement was achieved for structured elements with minimum adjustment of global parameters. Deviations from experiment occur in flexible loops and termini, indicating differences in both the extent and time scale of backbone motions. The lack of systematic patterns and water model dependence suggests that modeling of the local environment limits prediction accuracy. Nevertheless, qualitative accuracy in a 2 μs CHARMM36m Stam2 VHS domain simulation demonstrates the potential of MD-based interpretation in combination with NMR-measured dynamics, increasing the utility of spin relaxation in integrative structural biology.

Narrowing the gap between experimental and computational determination of methyl group dynamics in proteins

A synergistic analysis of methyl NMR relaxation data and MD simulations identifies ways to improve studies of protein dynamics.

Molecular impact of covalent modifications on nonribosomal peptide synthetase carrier protein communication

Nonribosomal peptide synthesis involves the interplay between covalent protein modifications, conformational fluctuations, catalysis, and transient protein-protein interactions. Delineating the mechanisms involved in orchestrating these various processes will deepen our understanding of domain-domain communication in nonribosomal peptide synthetases (NRPSs) and lay the groundwork for the rational reengineering of NRPSs by swapping domains handling different substrates to generate novel natural products. Although many structural and biochemical studies of NRPSs exist, few studies have focused on the energetics and dynamics governing the interactions in these systems. Here, we present detailed binding studies of an adenylation domain and its partner carrier protein in apo-, holo-, and substrate-loaded forms. Results from fluorescence anisotropy, isothermal titration calorimetry, and NMR titrations indicated that covalent modifications to a carrier protein modulate domain communication, suggesting that chemical modifications to carrier proteins during NRPS synthesis may impart directionality to sequential NRPS domain interactions. Comparison of the structure and dynamics of an apo-aryl carrier protein with those of its modified forms revealed structural fluctuations induced by post-translational modifications and mediated by modulations of protein dynamics. The results provide a comprehensive molecular description of a carrier protein throughout its life cycle and demonstrate how a network of dynamic residues can propagate the molecular impact of chemical modifications throughout a protein and influence its affinity toward partner domains.

Identification of primary and secondary UBA footprints on the surface of ubiquitin in cell-mimicking crowded solution

Despite significant advancements in our understanding of ubiquitin-mediated signaling, the influence of the intracellular environment on the formation of transient ubiquitin-partner complexes remains poorly explored. In our work, we introduce macromolecular crowding as a first level of complexity toward the imitation of a cellular environment in the study of such interactions. Using NMR spectroscopy, we find that the stereospecific complex of ubiquitin and the ubiquitin-associated domain (UBA) is minimally perturbed by the crowding agent Ficoll. However, in addition to the primary canonical recognition patch on ubiquitin, secondary patches are identified, indicating that in cell-mimicking crowded solution, UBA contacts ubiquitin at multiple sites.

A critical assessment of methods to recover information from averaged data

Conformational heterogeneity is key to the function of many biomacromolecules, but only a few groups have tried to characterize it until recently. Now, thanks to the increased throughput of experimental data and the increased computational power, the problem of the characterization of protein structural variability has become more and more popular. Several groups have devoted their efforts in trying to create quantitative, reliable and accurate protocols for extracting such information from averaged data. We analyze here different approaches, discussing strengths and weaknesses of each. All approaches can roughly be clustered into two groups: those satisfying the maximum entropy principle and those recovering ensembles composed of a restricted number of molecular conformations. In the first case, the solution focuses on the features that are common to all the infinite solutions satisfying the experimental data; in the second case, the reconstructed ensemble shows the conformational regions where a large probability can be placed. The upper limits for conformational probabilities (MaxOcc) can also be calculated. We also give an overview of the mainstream experimental observables, with considerations on the assumptions underlying their usage.

Characterizing Aciniform Silk Repetitive Domain Backbone Dynamics and Hydrodynamic Modularity

Spider aciniform (wrapping) silk is a remarkable fibrillar biomaterial with outstanding mechanical properties. It is a modular protein consisting, in Argiope trifasciata, of a core repetitive domain of 200 amino acid units (W units). In solution, the W units comprise a globular folded core, with five α-helices, and disordered tails that are linked to form a ~63-residue intrinsically disordered linker in concatemers. Herein, we present nuclear magnetic resonance (NMR) spectroscopy-based ¹⁵N spin relaxation analysis, allowing characterization of backbone dynamics as a function of residue on the ps–ns timescale in the context of the single W unit (W₁) and the two unit concatemer (W₂). Unambiguous mapping of backbone dynamics throughout W₂ was made possible by segmental NMR active isotope-enrichment through split intein-mediated trans-splicing. Spectral density mapping for W₁ and W₂ reveals a striking disparity in dynamics between the folded core and the disordered linker and tail regions. These data are also consistent with rotational diffusion behaviour where each globular domain tumbles almost independently of its neighbour. At a localized level, helix 5 exhibits elevated high frequency dynamics relative to the proximal helix 4, supporting a model of fibrillogenesis where this helix unfolds as part of the transition to a mixed α-helix/β-sheet fibre.

Linkage-specific conformational ensembles of non-canonical polyubiquitin chains

Polyubiquitination is a critical protein post-translational modification involved in a variety of processes in eukaryotic cells. The molecular basis for selective recognition of the polyubiquitin signals by cellular receptors is determined by the conformations polyubiquitin chains adopt; this has been demonstrated for K48- and K63-linked chains. Recent studies of the so-called non-canonical chains (linked via K6, K11, K27, K29, or K33) suggest they play important regulatory roles in growth, development, and immune system pathways, but biophysical studies are needed to elucidate the physical/structural basis of their interactions with receptors. A first step towards this goal is characterization of the conformations these chains adopt in solution. We assembled diubiquitins (Ub2) comprised of every lysine linkage. Using solution NMR measurements, small-angle neutron scattering (SANS), and in silico ensemble generation, we determined population-weighted conformational ensembles that shed light on the structure and dynamics of the non-canonical polyubiquitin chains. We found that polyubiquitin is conformationally heterogeneous, and each chain type exhibits unique conformational ensembles. For example, K6-Ub2 and K11-Ub2 (at physiological salt concentration) are in dynamic equilibrium between at least two conformers, where one exhibits a unique Ub/Ub interface, distinct from that observed in K48-Ub2 but similar to crystal structures of these chains. Conformers for K29-Ub2 and K33-Ub2 resemble recent crystal structures in the ligand-bound state. Remarkably, a number of diubiquitins adopt conformers similar to K48-Ub2 or K63-Ub2, suggesting potential overlap of biological function among different lysine linkages. These studies highlight the potential power of determining function from elucidation of conformational states.

Solution Structure of a Nonribosomal Peptide Synthetase Carrier Protein Loaded with Its Substrate Reveals Transient, Well-Defined Contacts

Nonribosomal peptide synthetases (NRPSs) are microbial enzymes that produce a wealth of important natural products by condensing substrates in an assembly line manner. The proper sequence of substrates is obtained by tethering them to phosphopantetheinyl arms of holo carrier proteins (CPs) via a thioester bond. CPs in holo and substrate-loaded forms visit NRPS catalytic domains in a series of transient interactions. A lack of structural information on substrate-loaded carrier proteins has hindered our understanding of NRPS synthesis. Here, we present the first structure of an NRPS aryl carrier protein loaded with its substrate via a native thioester bond, together with the structure of its holo form. We also present the first quantification of NRPS CP backbone dynamics. Our results indicate that prosthetic moieties in both holo and loaded forms are in contact with the protein core, but they also sample states in which they are disordered and extend in solution. We observe that substrate loading induces a large conformational change in the phosphopantetheinyl arm, thereby modulating surfaces accessible for binding to other domains. Our results are discussed in the context of NRPS domain interactions.

Functional dynamics in replication protein A DNA binding and protein recruitment domains

Replication Protein A (RPA) is an essential scaffold for many DNA processing machines; its function relies on its modular architecture. Here, we report (15)N-nuclear magnetic resonance heteronuclear relaxation analysis to characterize the movements of single-stranded (ss) DNA binding and protein interaction modules in the RPA70 subunit. Our results provide direct evidence for coordination of the motion of the tandem RPA70AB ssDNA binding domains. Moreover, binding of ssDNA substrate is found to cause dramatic reorientation and full coupling of inter-domain motion. In contrast, the RPA70N protein interaction domain remains structurally and dynamically independent of RPA70AB regardless of binding of ssDNA. This autonomy of motion between the 70N and 70AB modules supports a model in which the two binding functions of RPA are mediated fully independently, but remain differentially coordinated depending on the length of their flexible tethers. A critical role for linkers between the globular domains in determining the functional dynamics of RPA is proposed.

Changing the topology of protein backbone: the effect of backbone cyclization on the structure and dynamics of a SH3 domain

Understanding of the effects of the backbone cyclization on the structure and dynamics of a protein is essential for using protein topology engineering to alter protein stability and function. Here we have determined, for the first time, the structure and dynamics of the linear and various circular constructs of the N-SH3 domain from protein c-Crk. These constructs differ in the length and amino acid composition of the cyclization region. The backbone cyclization was carried out using intein-mediated intramolecular chemical ligation between the juxtaposed N- and the C-termini. The structure and backbone dynamics studies were performed using solution NMR. Our data suggest that the backbone cyclization has little effect on the overall three-dimensional structure of the SH3 domain: besides the termini, only minor structural changes were found in the proximity of the cyclization region. In contrast to the structure, backbone dynamics are significantly affected by the cyclization. On the subnanosecond time scale, the backbone of all circular constructs on average appears more rigid than that of the linear SH3 domain; this effect is observed over the entire backbone and is not limited to the cyclization site. The backbone mobility of the circular constructs becomes less restricted with increasing length of the circularization loop. In addition, significant conformational exchange motions (on the sub-millisecond time scale) were found in the N-Src loop and in the adjacent β-strands in all circular constructs studied in this work. These effects of backbone cyclization on protein dynamics have potential implications for the stability of the protein fold and for ligand binding.

pH dependence of conformational fluctuations of the protein backbone

Proton binding equilibria (pK(a) values) of ionizable groups in proteins are exquisitely sensitive to their microenvironments. Apparent pK(a) values measured for individual ionizable residues with NMR spectroscopy are actually population-weighted averages of the pK(a) in different conformational microstates. NMR spectroscopy experiments with staphylococcal nuclease were used to test the hypothesis that pK(a) values of surface Glu and Asp residues are affected by pH-sensitive fluctuations of the backbone between folded and locally unfolded conformations. (15)N spin relaxation studies showed that as the pH decreases from the neutral into the acidic range the amplitudes of backbone fluctuations in the ps-ns timescale increase near carboxylic residues. Hydrogen exchange experiments suggested that backbone conformational fluctuations promoted by decreasing pH also reflect slower local or sub-global unfolding near carboxylic groups. This study has implications for structure-based pKa calculations: (1) The timescale of the backbone's response to ionization events in proteins can range from ps to ms, and even longer; (2) pH-sensitive fluctuations of the backbone can be localized to both the segment the ionizable residue is attached to or the one that occludes the ionizable group; (3) Structural perturbations are not necessarily propagated through Coulomb interactions; instead, local fluctuations appear to be coupled through the co-operativity inherent to elements of secondary structure and to networks of hydrogen bonds. These results are consistent with the idea that local conformational fluctuations and stabilities are important determinants of apparent pK(a) values of ionizable residues in proteins.

The C113D mutation in human Pin1 causes allosteric structural changes in the phosphate binding pocket of the PPIase domain through the tug of war in the dual-histidine motif

Pin1 peptidyl-prolyl isomerase (PPIase) catalyzes specifically the pSer/pThr-Pro motif. The cis-trans isomerization mechanism has been studied by various approaches, including X-ray crystallography, site-directed mutagenesis, and the kinetic isotope effect on isomerization. However, a complete picture of the reaction mechanism remains elusive. On the basis of the X-ray structure of Pin1, residue C113 was proposed to play a nucleophile attacker to catalyze the isomerization. The controversial result that the C113D Pin1 mutant retains the activity, albeit at a reduced level, challenges the importance of C113 as a catalyst. To facilitate our understanding of the Pin1 isomerization process, we compared the structures and dynamics of the wild type with those of the C113D mutant Pin1 PPIase domains (residues 51-163). We found the C113D mutation disturbed the hydrogen bonds between the conserved histidine residues, H59 and H157 ("dual-histidine motif"); H59 imidazole forms a stable hydrogen bond to H157 in the wild type, whereas it has a strong hydrogen bond to D113 with weakened bonding to H157 in the C113D mutant. The C113D mutation unbalanced the hydrogen bonding tug of war for H59 between C113/D113 and H157 and destabilized the catalytic site structure, which eventually resulted in an altered conformation of the basic triad (K63, R68, and R69) that binds to the phosphate group in a substrate. The change in the basic triad structure could explain the severely weakened substrate binding ability of the C113D mutant. Overall, this work demonstrated that C113 plays a role in keeping the catalytic site in an active fold, which has never before been described.

Unusual structural features revealed by the solution NMR structure of the NLRC5 caspase recruitment domain

The cytosolic nucleotide-binding domain and leucine-rich repeat-containing receptors (NLRs) are key sensors for bacterial and viral invaders and endogenous stress signals. NLRs contain a varying N-terminal effector domain that regulates the downstream signaling events upon its activation and determines the subclass to which a NLR member belongs. NLRC5 contains an unclassified N-terminal effector domain that has been reported to interact downstream with the tandem caspase recruitment domain (CARD) of retinoic acid-inducible gene I (RIG-I). Here we report the solution structure of the N-terminal effector domain of NLRC5 and in vitro interaction experiments with the tandem CARD of RIG-I. The N-terminal effector domain of NLRC5 adopts a six α-helix bundle with a general death fold, though it displays specific structural features that are strikingly different from the CARD. Notably, α-helix 3 is replaced by an ordered loop, and α-helix 1 is devoid of the characteristic interruption. Detailed structural alignments between the N-terminal effector domains of NLRC5 with a representative of each death-fold subfamily showed that NLRC5 fits best to the CARD subfamily and can be called an atypical CARD. Due to the specific structural features, the atypical CARD also displays a different electrostatic surface. Because the shape and charge of the surface is crucial for the establishment of a homotypic CARD-CARD interaction, these specific structural features seem to have a significant effect on the interaction between the atypical CARD of NLRC5 and the tandem RIG-I CARD.

Unique structural, dynamical, and functional properties of k11-linked polyubiquitin chains

K11-linked polyubiquitin chains play important signaling and regulatory roles in both degradative and nonproteolytic pathways in eukaryotes. To understand the structural basis of how these chains are recognized and distinguished from other polyubiquitins, we determined solution structures of K11-linked diubiquitin (K11-Ub2) in the absence and presence of salt. These structures reveal that K11-Ub2 adopts conformations distinct from those of K48-linked or K63-linked chains. Importantly, our solution NMR and SANS data are inconsistent with published crystal structures of K11-Ub2. We found that increasing salt concentration compacts K11-Ub2 and strengthens interactions between the two Ub units. Binding studies indicate that K11-Ub2 interacts with ubiquitin-receptor proteins from both proteasomal and nonproteasomal pathways but with intermediate affinity and different binding modes than either K48-linked or K63-linked diubiquitin. Our data support the hypothesis that polyubiquitin chains of different linkages possess unique conformational and dynamical properties, allowing them to be recognized differently by downstream receptor proteins.

Alanine scan of core positions in ubiquitin reveals links between dynamics, stability, and function

Mutations at solvent-inaccessible core positions in proteins can impact function through many biophysical mechanisms including alterations to thermodynamic stability and protein dynamics. As these properties of proteins are difficult to investigate, the impacts of core mutations on protein function are poorly understood for most systems. Here, we determined the effects of alanine mutations at all 15 core positions in ubiquitin on function in yeast. The majority (13 of 15) of alanine substitutions supported yeast growth as the sole ubiquitin. Both the two null mutants (I30A and L43A) were less stable to temperature-induced unfolding in vitro than wild type (WT) but were well folded at physiological temperatures. Heteronuclear NMR studies indicated that the L43A mutation reduces temperature stability while retaining a ground-state structure similar to WT. This structure enables L43A to bind to common ubiquitin receptors in vitro. Many of the core alanine ubiquitin mutants, including one of the null variants (I30A), exhibited an increased accumulation of high-molecular-weight species, suggesting that these mutants caused a defect in the processing of ubiquitin-substrate conjugates. In contrast, L43A exhibited a unique accumulation pattern with reduced levels of high-molecular-weight species and undetectable levels of free ubiquitin. When conjugation to other proteins was blocked, L43A ubiquitin accumulated as free ubiquitin in yeast. Based on these findings, we speculate that ubiquitin's stability to unfolding may be required for efficient recycling during proteasome-mediated substrate degradation.

Deriving quantitative dynamics information for proteins and RNAs using ROTDIF with a graphical user interface

To facilitate rigorous analysis of molecular motions in proteins, DNA, and RNA, we present a new version of ROTDIF, a program for determining the overall rotational diffusion tensor from single- or multiple-field nuclear magnetic resonance relaxation data. We introduce four major features that expand the program's versatility and usability. The first feature is the ability to analyze, separately or together, (13)C and/or (15)N relaxation data collected at a single or multiple fields. A significant improvement in the accuracy compared to direct analysis of R2/R1 ratios, especially critical for analysis of (13)C relaxation data, is achieved by subtracting high-frequency contributions to relaxation rates. The second new feature is an improved method for computing the rotational diffusion tensor in the presence of biased errors, such as large conformational exchange contributions, that significantly enhances the accuracy of the computation. The third new feature is the integration of the domain alignment and docking module for relaxation-based structure determination of multi-domain systems. Finally, to improve accessibility to all the program features, we introduced a graphical user interface that simplifies and speeds up the analysis of the data. Written in Java, the new ROTDIF can run on virtually any computer platform. In addition, the new ROTDIF achieves an order of magnitude speedup over the previous version by implementing a more efficient deterministic minimization algorithm. We not only demonstrate the improvement in accuracy and speed of the new algorithm for synthetic and experimental (13)C and (15)N relaxation data for several proteins and nucleic acids, but also show that careful analysis required especially for characterizing RNA dynamics allowed us to uncover subtle conformational changes in RNA as a function of temperature that were opaque to previous analysis.

Nanosecond time scale motions in proteins revealed by high-resolution NMR relaxometry

Understanding the molecular determinants underlying protein function requires the characterization of both structure and dynamics at atomic resolution. Nuclear relaxation rates allow a precise characterization of protein dynamics at the Larmor frequencies of spins. This usually limits the sampling of motions to a narrow range of frequencies corresponding to high magnetic fields. At lower fields one cannot achieve sufficient sensitivity and resolution in NMR. Here, we use a fast shuttle device where the polarization builds up and the signals are detected at high field, while longitudinal relaxation takes place at low fields 0.5 < B₀ < 14.1 T. The sample is propelled over a distance up to 50 cm by a blowgun-like system in about 50 ms. The analysis of nitrogen-15 relaxation in the protein ubiquitin over such a wide range of magnetic fields offers unprecedented insights into molecular dynamics. Some key regions of the protein feature structural fluctuations on nanosecond time scales, which have so far been overlooked in high-field relaxation studies. Nanosecond motions in proteins may have been underestimated by traditional high-field approaches, and slower supra-τ_c motions that have no effect on relaxation may have been overestimated. High-resolution relaxometry thus opens the way to a quantitative characterization of nanosecond motions in proteins.

Equivalence between Euler angle conventions for the description of tensorial interactions in liquid NMR: application to different software programs

Long-range orientational restraints derived from alignment or rotational diffusion tensors have greatly contributed to the expansion of applications in biomolecular NMR. The orientation of the principal axis system of these tensors is usually described by the so-called Euler angles. However, no clear consensus has emerged concerning the convention of the associated orthogonal rotations. As a result, the different programs that derive or predict them have adopted different conventions, which make comparison between their results difficult. Moreover, the rotation schemes are seldom completely described. Here, we summarize the different conventions, determine which ones are adopted by commonly used software packages, and establish the formal equivalencies between the different calculated Euler angles.

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.

Analysis of 15N-1H NMR relaxation in proteins by a combined experimental and molecular dynamics simulation approach: picosecond-nanosecond dynamics of the Rho GTPase binding domain of plexin-B1 in the dimeric state indicates allosteric pathways

We investigate picosecond–nanosecond dynamics of the Rho-GTPase Binding Domain (RBD) of plexin-B1, which plays a key role in plexin-mediated cell signaling. Backbone 15N relaxation data of the dimeric RBD are analyzed with the model-free (MF) method, and with the slowly relaxing local structure/molecular dynamics (SRLS-MD) approach. Independent analysis of the MD trajectories, based on the MF paradigm, is also carried out. MF is a widely popular and simple method, SRLS is a general approach, and SRLS-MD is an integrated approach we developed recently. Corresponding parameters from the RBD dimer, a previously studied RBD monomer mutant, and the previously studied complex of the latter with the GTPase Rac1, are compared. The L2, L3, and L4 loops of the plexin-B1 RBD are involved in interactions with other plexin domains, GTPase binding, and RBD dimerization, respectively. Peptide groups in the loops of both the monomeric and dimeric RBD are found to experience weak and moderately asymmetric local ordering centered approximately at the C(i–1)(α)–C(i)(α) axes, and nanosecond backbone motion. Peptide groups in the α-helices and the β-strands of the dimer (the β-strands of the monomer) experience strong and highly asymmetric local ordering centered approximately at the C(i–1)(α)–C(i)(α) axes (N–H bonds). N–H fluctuations occur on the picosecond time scale. An allosteric pathway for GTPase binding, providing new insights into plexin function, is delineated.

A small molecule agonist of EphA2 receptor tyrosine kinase inhibits tumor cell migration in vitro and prostate cancer metastasis in vivo

During tumor progression, EphA2 receptor can gain ligand-independent pro-oncogenic functions due to Akt activation and reduced ephrin-A ligand engagement. The effects can be reversed by ligand stimulation, which triggers the intrinsic tumor suppressive signaling pathways of EphA2 including inhibition of PI3/Akt and Ras/ERK pathways. These observations argue for development of small molecule agonists for EphA2 as potential tumor intervention agents. Through virtual screening and cell-based assays, we report here the identification and characterization of doxazosin as a novel small molecule agonist for EphA2 and EphA4, but not for other Eph receptors tested. NMR studies revealed extensive contacts of doxazosin with EphA2/A4, recapitulating both hydrophobic and electrostatic interactions recently found in the EphA2/ephrin-A1 complex. Clinically used as an α1-adrenoreceptor antagonist (Cardura®) for treating hypertension and benign prostate hyperplasia, doxazosin activated EphA2 independent of α1-adrenoreceptor. Similar to ephrin-A1, doxazosin inhibited Akt and ERK kinase activities in an EphA2-dependent manner. Treatment with doxazosin triggered EphA2 receptor internalization, and suppressed haptotactic and chemotactic migration of prostate cancer, breast cancer, and glioma cells. Moreover, in an orthotopic xenograft model, doxazosin reduced distal metastasis of human prostate cancer cells and prolonged survival in recipient mice. To our knowledge, doxazosin is the first small molecule agonist of a receptor tyrosine kinase that is capable of inhibiting malignant behaviors in vitro and in vivo.

Measurement of ¹⁵N relaxation rates in perdeuterated proteins by TROSY-based methods

While extracting dynamics parameters from backbone (15)N relaxation measurements in proteins has become routine over the past two decades, it is increasingly recognized that accurate quantitative analysis can remain limited by the potential presence of systematic errors associated with the measurement of (15)N R(1) and R(2) or R(1ρ) relaxation rates as well as heteronuclear (15)N-{(1)H} NOE values. We show that systematic errors in such measurements can be far larger than the statistical error derived from either the observed signal-to-noise ratio, or from the reproducibility of the measurement. Unless special precautions are taken, the problem of systematic errors is shown to be particularly acute in perdeuterated systems, and even more so when TROSY instead of HSQC elements are used to read out the (15)N magnetization through the NMR-sensitive (1)H nucleus. A discussion of the most common sources of systematic errors is presented, as well as TROSY-based pulse schemes that appear free of systematic errors to the level of <1 %. Application to the small perdeuterated protein GB3, which yields exceptionally high S/N and therefore is an ideal test molecule for detection of systematic errors, yields relaxation rates that show considerably less residue by residue variation than previous measurements. Measured R(2)'/R(1)' ratios fit an axially symmetric diffusion tensor with a Pearson's correlation coefficient of 0.97, comparable to fits obtained for backbone amide RDCs to the Saupe matrix.

Combining NMR and molecular dynamics studies for insights into the allostery of small GTPase-protein interactions

Combinations of experimentally derived data from nuclear magnetic resonance spectroscopy and analyses of molecular dynamics trajectories increasingly allow us to obtain a detailed description of the molecular mechanisms by which proteins function in signal transduction. This chapter provides an introduction into these two methodologies, illustrated by example of a small GTPase-effector interaction. It is increasingly becoming clear that new insights are provided by the combination of experimental and computational methods. Understanding the structural and protein dynamical contributions to allostery will be useful for the engineering of new binding interfaces and protein functions, as well as for the design/in silico screening of chemical agents that can manipulate the function of small GTPase-protein interactions in diseases such as cancer.

Determining protein dynamics from ¹⁵N relaxation data by using DYNAMICS

Motions are essential for protein function, and knowledge of protein dynamics is a key to our understanding the mechanisms underlying protein folding and stability, ligand recognition, allostery, and catalysis. In the last two decades, NMR relaxation measurements have become a powerful tool for characterizing backbone and side chain dynamics in complex biological macromolecules such as proteins and nucleic acids. Accurate analysis of the experimental data in terms of motional parameters is an essential prerequisite for developing physical models of motions to paint an adequate picture of protein dynamics. Here, I describe in detail how to use the software package DYNAMICS that was developed for accurate characterization of the overall tumbling and local dynamics in a protein from nuclear spin-relaxation rates measured by NMR. Step-by-step instructions are provided and illustrated through an analysis of (15)N relaxation data for protein G.

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.

Fast approximations of the rotational diffusion tensor and their application to structural assembly of molecular complexes

We present and evaluate a rigid-body, deterministic, molecular docking method, called ELMDOCK, that relies solely on the three-dimensional structure of the individual components and the overall rotational diffusion tensor of the complex, obtained from nuclear spin-relaxation measurements. We also introduce a docking method, called ELMPATIDOCK, derived from ELMDOCK and based on the new concept of combining the shape-related restraints from rotational diffusion with those from residual dipolar couplings, along with ambiguous contact/interface-related restraints obtained from chemical shift perturbations. ELMDOCK and ELMPATIDOCK use two novel approximations of the molecular rotational diffusion tensor that allow computationally efficient docking. We show that these approximations are accurate enough to properly dock the two components of a complex without the need to recompute the diffusion tensor at each iteration step. We analyze the accuracy, robustness, and efficiency of these methods using synthetic relaxation data for a large variety of protein-protein complexes. We also test our method on three protein systems for which the structure of the complex and experimental relaxation data are available, and analyze the effect of flexible unstructured tails on the outcome of docking. Additionally, we describe a method for integrating the new approximation methods into the existing docking approaches that use the rotational diffusion tensor as a restraint. The results show that the proposed docking method is robust against experimental errors in the relaxation data or structural rearrangements upon complex formation and is computationally more efficient than current methods. The developed approximations are accurate enough to be used in structure refinement protocols.

NMR reveals a different mode of binding of the Stam2 VHS domain to ubiquitin and diubiquitin

The VHS domain of the Stam2 protein is a ubiquitin binding domain involved in the recognition of ubiquitinated proteins committed to lysosomal degradation. Among all VHS domains, the VHS domain of Stam proteins is the strongest binder to monoubiqiuitin and exhibits preferences for K63-linked chains. In the present paper, we report the solution NMR structure of the Stam2-VHS domain in complex with monoubiquitin by means of chemical shift perturbations, spin relaxation, and paramagnetic relaxation enhancements. We also characterize the interaction of Stam2-VHS with K48- and K63-linked diubiquitin chains and report the first evidence that VHS binds differently to these two chains. Our data reveal that VHS enters the hydrophobic pocket of K48-linked diubiquitin and binds the two ubiquitin subunits with different affinities. In contrast, VHS interacts with K63-linked diubiquitin in a mode similar to its interaction with monoubiquitin. We also suggest possible structural models for both K48- and K63-linked diubiquitin in interaction with VHS. Our results, which demonstrate a different mode of binding of VHS for K48- and K63-linked diubiquitin, may explain the preference of VHS for K63- over K48-linked diubiquitin chains and monoubiquitin.