NMR spectroscopy captures the essential role of dynamics in regulating biomolecular function.

Biomolecules are in constant motion. To understand how they function, and why malfunctions can cause disease, it is necessary to describe their three-dimensional structures in terms of dynamic conformational ensembles. Here, we demonstrate how nuclear magnetic resonance (NMR) spectroscopy provides an essential, dynamic view of structural biology that captures biomolecular motions at atomic resolution. We focus on examples that emphasize the diversity of biomolecules and biochemical applications that are amenable to NMR, such as elucidating functional dynamics in large molecular machines, characterizing transient conformations implicated in the onset of disease, and obtaining atomic-level descriptions of intrinsically disordered regions that make weak interactions involved in liquid-liquid phase separation. Finally, we discuss the pivotal role that NMR has played in driving forward our understanding of the biomolecular dynamics-function paradigm.

Systematic identification of conditionally folded intrinsically disordered regions by AlphaFold2.

The AlphaFold Protein Structure Database contains predicted structures for millions of proteins. For the majority of human proteins that contain intrinsically disordered regions (IDRs), which do not adopt a stable structure, it is generally assumed that these regions have low AlphaFold2 confidence scores that reflect low-confidence structural predictions. Here, we show that AlphaFold2 assigns confident structures to nearly 15% of human IDRs. By comparison to experimental NMR data for a subset of IDRs that are known to conditionally fold (i.e., upon binding or under other specific conditions), we find that AlphaFold2 often predicts the structure of the conditionally folded state. Based on databases of IDRs that are known to conditionally fold, we estimate that AlphaFold2 can identify conditionally folding IDRs at a precision as high as 88% at a 10% false positive rate, which is remarkable considering that conditionally folded IDR structures were minimally represented in its training data. We find that human disease mutations are nearly fivefold enriched in conditionally folded IDRs over IDRs in general and that up to 80% of IDRs in prokaryotes are predicted to conditionally fold, compared to less than 20% of eukaryotic IDRs. These results indicate that a large majority of IDRs in the proteomes of human and other eukaryotes function in the absence of conditional folding, but the regions that do acquire folds are more sensitive to mutations. We emphasize that the AlphaFold2 predictions do not reveal functionally relevant structural plasticity within IDRs and cannot offer realistic ensemble representations of conditionally folded IDRs.

Activation of caspase-9 on the apoptosome as studied by methyl-TROSY NMR.

Mitochondrial apoptotic signaling cascades lead to the formation of the apoptosome, a 1.1-MDa heptameric protein scaffold that recruits and activates the caspase-9 protease. Once activated, caspase-9 cleaves and activates downstream effector caspases, triggering the onset of cell death through caspase-mediated proteolysis of cellular proteins. Failure to activate caspase-9 enables the evasion of programmed cell death, which occurs in various forms of cancer. Despite the critical apoptotic function of caspase-9, the structural mechanism by which it is activated on the apoptosome has remained elusive. Here, we used a combination of methyl-transverse relaxation-optimized NMR spectroscopy, protein engineering, and biochemical assays to study the activation of caspase-9 bound to the apoptosome. In the absence of peptide substrate, we observed that both caspase-9 and its isolated protease domain (PD) only very weakly dimerize with dissociation constants in the millimolar range. Methyl-NMR spectra of isotope-labeled caspase-9, within the 1.3-MDa native apoptosome complex or an engineered 480-kDa apoptosome mimic, reveal that the caspase-9 PD remains monomeric after recruitment to the scaffold. Binding to the apoptosome, therefore, organizes caspase-9 PDs so that they can rapidly and extensively dimerize only when substrate is present, providing an important layer in the regulation of caspase-9 activation. Our work highlights the unique role of NMR spectroscopy to structurally characterize protein domains that are flexibly tethered to large scaffolds, even in cases where the molecular targets are in excess of 1 MDa, as in the present example.

A weakened interface in the P182L variant of HSP27 associated with severe Charcot-Marie-Tooth neuropathy causes aberrant binding to interacting proteins.

HSP27 is a human molecular chaperone that forms large, dynamic oligomers and functions in many aspects of cellular homeostasis. Mutations in HSP27 cause Charcot-Marie-Tooth (CMT) disease, the most common inherited disorder of the peripheral nervous system. A particularly severe form of CMT disease is triggered by the P182L mutation in the highly conserved IxI/V motif of the disordered C-terminal region, which interacts weakly with the structured core domain of HSP27. Here, we observed that the P182L mutation disrupts the chaperone activity and significantly increases the size of HSP27 oligomers formed in vivo, including in motor neurons differentiated from CMT patient-derived stem cells. Using NMR spectroscopy, we determined that the P182L mutation decreases the affinity of the HSP27 IxI/V motif for its own core domain, leaving this binding site more accessible for other IxI/V-containing proteins. We identified multiple IxI/V-bearing proteins that bind with higher affinity to the P182L variant due to the increased availability of the IxI/V-binding site. Our results provide a mechanistic basis for the impact of the P182L mutation on HSP27 and suggest that the IxI/V motif plays an important, regulatory role in modulating protein-protein interactions.

Study of protein folding under native conditions by rapidly switching the hydrostatic pressure inside an NMR sample cell.

In general, small proteins rapidly fold on the timescale of milliseconds or less. For proteins with a substantial volume difference between the folded and unfolded states, their thermodynamic equilibrium can be altered by varying the hydrostatic pressure. Using a pressure-sensitized mutant of ubiquitin, we demonstrate that rapidly switching the pressure within an NMR sample cell enables study of the unfolded protein under native conditions and, vice versa, study of the native protein under denaturing conditions. This approach makes it possible to record 2D and 3D NMR spectra of the unfolded protein at atmospheric pressure, providing residue-specific information on the folding process. 15N and 13C chemical shifts measured immediately after dropping the pressure from 2.5 kbar (favoring unfolding) to 1 bar (native) are close to the random-coil chemical shifts observed for a large, disordered peptide fragment of the protein. However, 15N relaxation data show evidence for rapid exchange, on a ∼100-µs timescale, between the unfolded state and unstable, structured states that can be considered as failed folding events. The NMR data also provide direct evidence for parallel folding pathways, with approximately one-half of the protein molecules efficiently folding through an on-pathway kinetic intermediate, whereas the other half fold in a single step. At protein concentrations above ∼300 µM, oligomeric off-pathway intermediates compete with folding of the native state.

Monitoring 15N Chemical Shifts During Protein Folding by Pressure-Jump NMR.

Pressure-jump hardware permits direct observation of protein NMR spectra during a cyclically repeated protein folding process. For a two-state folding protein, the change in resonance frequency will occur nearly instantaneously when the protein clears the transition state barrier, resulting in a monoexponential change of the ensemble-averaged chemical shift. However, protein folding pathways can be more complex and contain metastable intermediates. With a pseudo-3D NMR experiment that utilizes stroboscopic observation, we measure the ensemble-averaged chemical shifts, including those of exchange-broadened intermediates, during the folding process. Such measurements for a pressure-sensitized mutant of ubiquitin show an on-pathway kinetic intermediate whose 15N chemical shifts differ most from the natively folded protein for strands ß5, its preceding turn, and the two strands that pair with ß5 in the native structure.

Propensity for cis-Proline Formation in Unfolded Proteins.

In unfolded proteins, peptide bonds involving Pro residues exist in equilibrium between the minor cis and major trans conformations. Folded proteins predominantly contain trans-Pro bonds, and slow cis-trans Pro isomerization in the unfolded state is often found to be a rate-limiting step in protein folding. Moreover, kinases and phosphatases that act upon Ser/Thr-Pro motifs exhibit preferential recognition of either the cis- or trans-Pro conformer. Here, NMR spectra obtained at both atmospheric and high pressures indicate that the population of cis-Pro falls well below previous estimates, an effect attributed to the use of short peptides with charged termini in most prior model studies. For the intrinsically disordered protein α-synuclein, cis-Pro populations at all of its five X-Pro bonds are less than 5 %, with only modest ionic strength dependence and no detectable effect of the previously demonstrated interaction between the N- and C-terminal halves of the protein. Comparison to small peptides with the same amino-acid sequence indicates that peptides, particularly those with unblocked, oppositely charged amino and carboxyl end groups, strongly overestimate the amount of cis-Pro.

Monitoring Hydrogen Exchange During Protein Folding by Fast Pressure Jump NMR Spectroscopy.

A method is introduced that permits direct observation of the rates at which backbone amide hydrogens become protected from solvent exchange after rapidly dropping the hydrostatic pressure inside the NMR sample cell from denaturing (2.5 kbar) to native (1 bar) conditions. The method is demonstrated for a pressure-sensitized ubiquitin variant that contains two Val to Ala mutations. Increased protection against hydrogen exchange with solvent is monitored as a function of time during the folding process. Results for 53 backbone amides show narrow clustering with protection occurring with a time constant of ca. 85 ms, but slower protection is observed around a reverse turn near the C-terminus of the protein. Remarkably, the native NMR spectrum returns with this slower time constant of ca. 150 ms, indicating that the almost fully folded protein retains molten globule characteristics with severe NMR line broadening until the final hydrogen bonds are formed. Prior to crossing the transition state barrier, hydrogen exchange protection factors are close to unity, but with slightly elevated values in the ß1-ß2 hairpin, previously shown to be already lowly populated in the urea-denatured state.

Automatic Assignment of Methyl-NMR Spectra of Supramolecular Machines Using Graph Theory.

Methyl groups are powerful probes for the analysis of structure, dynamics and function of supramolecular assemblies, using both solution- and solid-state NMR. Widespread application of the methodology has been limited due to the challenges associated with assigning spectral resonances to specific locations within a biomolecule. Here, we present Methyl Assignment by Graph Matching (MAGMA), for the automatic assignment of methyl resonances. A graph matching protocol examines all possibilities for each resonance in order to determine an exact assignment that includes a complete description of any ambiguity. MAGMA gives 100% accuracy in confident assignments when tested against both synthetic data, and 9 cross-validated examples using both solution- and solid-state NMR data. We show that this remarkable accuracy enables a user to distinguish between alternative protein structures. In a drug discovery application on HSP90, we show the method can rapidly and efficiently distinguish between possible ligand binding modes. By providing an exact and robust solution to methyl resonance assignment, MAGMA can facilitate significantly accelerated studies of supramolecular machines using methyl-based NMR spectroscopy.

Tangled web of interactions among proteins involved in iron-sulfur cluster assembly as unraveled by NMR, SAXS, chemical crosslinking, and functional studies.

Proteins containing iron-sulfur (Fe-S) clusters arose early in evolution and are essential to life. Organisms have evolved machinery consisting of specialized proteins that operate together to assemble Fe-S clusters efficiently so as to minimize cellular exposure to their toxic constituents: iron and sulfide ions. To date, the best studied system is the iron-sulfur cluster (isc) operon of Escherichia coli, and the eight ISC proteins it encodes. Our investigations over the past five years have identified two functional conformational states for the scaffold protein (IscU) and have shown that the other ISC proteins that interact with IscU prefer to bind one conformational state or the other. From analyses of the NMR spectroscopy-derived network of interactions of ISC proteins, small-angle X-ray scattering (SAXS) data, chemical crosslinking experiments, and functional assays, we have constructed working models for Fe-S cluster assembly and delivery. Future work is needed to validate and refine what has been learned about the E. coli system and to extend these findings to the homologous Fe-S cluster biosynthetic machinery of yeast and human mitochondria. This article is part of a Special Issue entitled: Fe/S proteins: Analysis, structure, function, biogenesis and diseases.

The specialized Hsp70 (HscA) interdomain linker binds to its nucleotide-binding domain and stimulates ATP hydrolysis in both cis and trans configurations.

Proteins from the isc operon of Escherichia coli constitute the machinery used to synthesize iron-sulfur (Fe-S) clusters for delivery to recipient apoproteins. Efficient and rapid [2Fe-2S] cluster transfer from the holo-scaffold protein IscU depends on ATP hydrolysis in the nucleotide-binding domain (NBD) of HscA, a specialized Hsp70-type molecular chaperone with low intrinsic ATPase activity (0.02 min(-1) at 25 °C, henceforth reported in units of min(-1)). HscB, an Hsp40-type cochaperone, binds to HscA and stimulates ATP hydrolysis to promote cluster transfer, yet while the interactions between HscA and HscB have been investigated, the role of HscA's interdomain linker in modulating ATPase activity has not been explored. To address this issue, we created three variants of the 40 kDa NBD of HscA: NBD alone (HscA386), NBD with a partial linker (HscA389), and NBD with the full linker (HscA395). We found that the rate of ATP hydrolysis of HscA395 (0.45 min(-1)) is nearly 15-fold higher than that of HscA386 (0.035 min(-1)), although their apparent affinities for ATP are equivalent. HscA395, which contains the full covalently linked linker peptide, exhibited intrinsic tryptophan fluorescence emission and basal thermostability that were higher than those of HscA386. Furthermore, HscA395 displayed narrower (1)H(N) line widths in its two-dimensional (1)H-(15)N TROSY-HSQC spectrum in comparison to HscA386, indicating that the peptide in the cis configuration binds to and stabilizes the structure of the NBD. The addition to HscA386 of a synthetic peptide with a sequence identical to that of the interdomain linker (L(387)LLDVIPLS(395)) stimulated its ATPase activity and induced widespread NMR chemical shift perturbations indicative of a binding interaction in the trans configuration.

Nucleotide-dependent interactions within a specialized Hsp70/Hsp40 complex involved in Fe-S cluster biogenesis.

The structural mechanism by which Hsp70-type chaperones interact with Hsp40-type co-chaperones has been of great interest, yet still remains a matter of debate. Here, we used solution NMR spectroscopy to investigate the ATP-/ADP-dependent interactions between Escherichia coli HscA and HscB, the specialized Hsp70/Hsp40 molecular chaperones that mediate iron-sulfur cluster transfer. We observed that NMR signals assigned to amino acid residues in the J-domain and its "HPD" motif of HscB broadened severely upon the addition of ATP-bound HscA, but these signals were not similarly broadened by ADP-bound HscA or the isolated nucleotide binding domain of HscA complexed with either ATP or ADP. An HscB variant with an altered HPD motif, HscB(H32A,P33A,D34A), failed to manifest WT-like NMR signal perturbations and also abolished WT-like stimulation of ATP hydrolysis by HscA. In addition, residues 153-171 in the C-terminal region of HscB exhibited NMR signal perturbations upon interaction with HscA, alone or complexed with ADP or ATP. These results demonstrate that the HPD motif in the J-domain of HscB directly interacts with ATP-bound HscA and suggest that a second, less nucleotide-dependent binding site for HscA resides in the C-terminal region of HscB.

PARP1 condensates differentially partition DNA repair proteins and enhance DNA ligation.

Poly(ADP-ribose) polymerase 1 (PARP1) is one of the first responders to DNA damage and plays crucial roles in recruiting DNA repair proteins through its activity - poly(ADP-ribosyl)ation (PARylation). The enrichment of DNA repair proteins at sites of DNA damage has been described as the formation of a biomolecular condensate. However, it is not understood how PARP1 and PARylation contribute to the formation and organization of DNA repair condensates. Using recombinant human PARP1 in vitro, we find that PARP1 readily forms viscous biomolecular condensates in a DNA-dependent manner and that this depends on its three zinc finger (ZnF) domains. PARylation enhances PARP1 condensation in a PAR chain-length dependent manner and increases the internal dynamics of PARP1 condensates. DNA and single-strand break repair proteins XRCC1, LigIII, Polß, and FUS partition in PARP1 condensates, although in different patterns. While Polß and FUS are both homogeneously mixed within PARP1 condensates, FUS enrichment is greatly enhanced upon PARylation whereas Polß partitioning is not. XRCC1 and LigIII display an inhomogeneous organization within PARP1 condensates; their enrichment in these multiphase condensates is enhanced by PARylation. Functionally, PARP1 condensates concentrate short DNA fragments and facilitate compaction of long DNA and bridge DNA ends. Furthermore, the presence of PARP1 condensates significantly promotes DNA ligation upon PARylation. These findings provide insight into how PARP1 condensation and PARylation regulate the assembly and biochemical activities in DNA repair foci, which may inform on how PARPs function in other PAR-driven condensates.

FOXO transcription factors differ in their dynamics and intra/intermolecular interactions.

Transcription factors play key roles in orchestrating a plethora of cellular mechanisms and controlling cellular homeostasis. Transcription factors share distinct DNA binding domains, which allows to group them into protein families. Among them, the Forkhead box O (FOXO) family contains transcription factors crucial for cellular homeostasis, longevity and response to stress. The dysregulation of FOXO signaling is linked to drug resistance in cancer therapy or cellular senescence, however, selective drugs targeting FOXOs are limited, thus knowledge about structure and dynamics of FOXO proteins is essential. Here, we provide an extensive study of structure and dynamics of all FOXO family members. We identify residues accounting for different dynamic and structural features. Furthermore, we show that the auto-inhibition of FOXO proteins by their C-terminal trans-activation domain is conserved throughout the family and that these interactions are not only possible intra-, but also inter-molecularly. This indicates a model in which FOXO transcription factors would modulate their activities by interacting mutually.

Automated assignment of methyl NMR spectra from large proteins.

As structural biology trends towards larger and more complex biomolecular targets, a detailed understanding of their interactions and underlying structures and dynamics is required. The development of methyl-TROSY has enabled NMR spectroscopy to provide atomic-resolution insight into the mechanisms of large molecular assemblies in solution. However, the applicability of methyl-TROSY has been hindered by the laborious and time-consuming resonance assignment process, typically performed with domain fragmentation, site-directed mutagenesis, and analysis of NOE data in the context of a crystal structure. In response, several structure-based automatic methyl assignment strategies have been developed over the past decade. Here, we present a comprehensive analysis of all available methods and compare their input data requirements, algorithmic strategies, and reported performance. In general, the methods fall into two categories: those that primarily rely on inter-methyl NOEs, and those that utilize methyl PRE- and PCS-based restraints. We discuss their advantages and limitations, and highlight the potential benefits from standardizing and combining different methods.

Unveiling invisible protein states with NMR spectroscopy.

Proteins interconvert between multiple conformations, including sparsely populated and transiently formed states that are difficult to characterize in structural detail using standard biophysical methods. In some cases, changes to the dynamical equilibria between conformations can lead to pathological protein aggregation and to the disruption of cellular homeostasis. The detection and characterization of lowly populated conformers is therefore crucial for understanding the basis of protein misfolding. NMR spectroscopy is exquisitely sensitive to the conformational dynamics of biomolecules and can be used to study sparsely populated states at the atomic level. Here, we review recent progress toward understanding the roles of sparsely populated, otherwise 'invisible' states present in protein folding and misfolding, where NMR has provided unique insight into folding intermediates, transiently misfolded states, and soluble oligomers that precede amyloid fibril formation.

Protein structural changes characterized by high-pressure, pulsed field gradient diffusion NMR spectroscopy.

Pulsed-field gradient NMR spectroscopy is widely used to measure the translational diffusion and hydrodynamic radius (Rh) of biomolecules in solution. For unfolded proteins, the Rh provides a sensitive reporter on the ensemble-averaged conformation and the extent of polypeptide chain expansion as a function of added denaturant. Hydrostatic pressure is a convenient and reversible alternative to chemical denaturants for the study of protein folding, and enables NMR measurements to be performed on a single sample. While the impact of pressure on the viscosity of water is well known, and our water diffusivity measurements agree closely with theoretical expectations, we find that elevated pressures increase the Rh of dioxane and other small molecules by amounts that correlate with their hydrophobicity, with parallel increases in rotational friction indicated by 13C longitudinal relaxation times. These data point to a tighter coupling with water for hydrophobic surfaces at elevated pressures. Translational diffusion measurement of the unfolded state of a pressure-sensitized ubiquitin mutant (VA2-ubiquitin) as a function of hydrostatic pressure or urea concentration shows that Rh values of both the folded and the unfolded states remain nearly invariant. At ca 23 Å, the Rh of the fully pressure-denatured state is essentially indistinguishable from the urea-denatured state, and close to the value expected for an idealized random coil of 76 residues. The intrinsically disordered protein (IDP) α-synuclein shows slight compaction at pressures above 2 kbar. Diffusion of unfolded ubiquitin and α-synuclein is significantly impacted by sample concentration, indicating that quantitative measurements need to be carried out under dilute conditions.

Conditional Disorder in Small Heat-shock Proteins.

Small heat-shock proteins (sHSPs) are molecular chaperones that respond to cellular stresses to combat protein aggregation. HSP27 is a critical human sHSP that forms large, dynamic oligomers whose quaternary structures and chaperone activities depend on environmental factors. Upon exposure to cellular stresses, such as heat shock or acidosis, HSP27 oligomers can dissociate into dimers and monomers, which leads to significantly enhanced chaperone activity. The structured core of the protein, the α-crystallin domain (ACD), forms dimers and can prevent the aggregation of substrate proteins to a similar degree as the full-length protein. When the ACD dimer dissociates into monomers, it partially unfolds and exhibits enhanced activity. Here, we used solution-state NMR spectroscopy to characterize the structure and dynamics of the HSP27 ACD monomer. Web show that the monomer is stabilized at low pH and that its backbone chemical shifts, 15N relaxation rates, and 1H-15N residual dipolar couplings suggest structural changes and rapid motions in the region responsible for dimerization. By analyzing the solvent accessible and buried surface areas of sHSP structures in the context of a database of dimers that are known to dissociate into disordered monomers, we predict that ACD dimers from sHSPs across all kingdoms of life may partially unfold upon dissociation. We propose a general model in which conditional disorder-the partial unfolding of ACDs upon monomerization-is a common mechanism for sHSP activity.

Automatic structure-based NMR methyl resonance assignment in large proteins.

Isotopically labeled methyl groups provide NMR probes in large, otherwise deuterated proteins. However, the resonance assignment constitutes a bottleneck for broader applicability of methyl-based NMR. Here, we present the automated MethylFLYA method for the assignment of methyl groups that is based on methyl-methyl nuclear Overhauser effect spectroscopy (NOESY) peak lists. MethylFLYA is applied to five proteins (28-358 kDa) comprising a total of 708 isotope-labeled methyl groups, of which 612 contribute NOESY cross peaks. MethylFLYA confidently assigns 488 methyl groups, i.e. 80% of those with NOESY data. Of these, 459 agree with the reference, 6 were different, and 23 were without reference assignment. MethylFLYA assigns significantly more methyl groups than alternative algorithms, has an average error rate of 1%, modest runtimes of 0.4-1.2 h, and can handle arbitrary isotope labeling patterns and data from other types of NMR spectra.