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.

Pan-vaccine analysis reveals innate immune endotypes predictive of antibody responses to vaccination.

Several studies have shown that the pre-vaccination immune state is associated with the antibody response to vaccination. However, the generalizability and mechanisms that underlie this association remain poorly defined. Here, we sought to identify a common pre-vaccination signature and mechanisms that could predict the immune response across 13 different vaccines. Analysis of blood transcriptional profiles across studies revealed three distinct pre-vaccination endotypes, characterized by the differential expression of genes associated with a pro-inflammatory response, cell proliferation, and metabolism alterations. Importantly, individuals whose pre-vaccination endotype was enriched in pro-inflammatory response genes known to be downstream of nuclear factor-kappa B showed significantly higher serum antibody responses 1 month after vaccination. This pro-inflammatory pre-vaccination endotype showed gene expression characteristic of the innate activation state triggered by Toll-like receptor ligands or adjuvants. These results demonstrate that wide variations in the transcriptional state of the immune system in humans can be a key determinant of responsiveness to vaccination.

Transcriptional atlas of the human immune response to 13 vaccines reveals a common predictor of vaccine-induced antibody responses.

Systems vaccinology has defined molecular signatures and mechanisms of immunity to vaccination. However, comparative analysis of immunity to different vaccines is lacking. We integrated transcriptional data of over 3,000 samples, from 820 adults across 28 studies of 13 vaccines and analyzed vaccination-induced signatures of antibody responses. Most vaccines induced signatures of innate immunity and plasmablasts at days 1 and 7, respectively, after vaccination. However, the yellow fever vaccine induced an early transient signature of T and B cell activation at day 1, followed by delayed antiviral/interferon and plasmablast signatures that peaked at days 7 and 14-21, respectively. Thus, there was no evidence for a 'universal signature' that predicted antibody response to all vaccines. However, accounting for the asynchronous nature of responses, we defined a time-adjusted signature that predicted antibody responses across vaccines. These results provide a transcriptional atlas of immunity to vaccination and define a common, time-adjusted signature of antibody responses.

Structural transitions enable interleukin-18 maturation and signaling.

Several interleukin-1 (IL-1) family members, including IL-1ß and IL-18, require processing by inflammasome-associated caspases to unleash their activities. Here, we unveil, by cryoelectron microscopy (cryo-EM), two major conformations of the complex between caspase-1 and pro-IL-18. One conformation is similar to the complex of caspase-4 and pro-IL-18, with interactions at both the active site and an exosite (closed conformation), and the other only contains interactions at the active site (open conformation). Thus, pro-IL-18 recruitment and processing by caspase-1 is less dependent on the exosite than the active site, unlike caspase-4. Structure determination by nuclear magnetic resonance uncovers a compact fold of apo pro-IL-18, which is similar to caspase-1-bound pro-IL-18 but distinct from cleaved IL-18. Binding sites for IL-18 receptor and IL-18 binding protein are only formed upon conformational changes after pro-IL-18 cleavage. These studies show how pro-IL-18 is selected as a caspase-1 substrate, and why cleavage is necessary for its inflammatory activity.

Poly(ADP-ribosyl)ation enhances nucleosome dynamics and organizes DNA damage repair components within biomolecular condensates.

Nucleosomes, the basic structural units of chromatin, hinder recruitment and activity of various DNA repair proteins, necessitating modifications that enhance DNA accessibility. Poly(ADP-ribosyl)ation (PARylation) of proteins near damage sites is an essential initiation step in several DNA-repair pathways; however, its effects on nucleosome structural dynamics and organization are unclear. Using NMR, cryoelectron microscopy (cryo-EM), and biochemical assays, we show that PARylation enhances motions of the histone H3 tail and DNA, leaving the configuration of the core intact while also stimulating nuclease digestion and ligation of nicked nucleosomal DNA by LIG3. PARylation disrupted interactions between nucleosomes, preventing self-association. Addition of LIG3 and XRCC1 to PARylated nucleosomes generated condensates that selectively partition DNA repair-associated proteins in a PAR- and phosphorylation-dependent manner in vitro. Our results establish that PARylation influences nucleosomes across different length scales, extending from the atom-level motions of histone tails to the mesoscale formation of condensates with selective compositions.

Bringing dynamic molecular machines into focus by methyl-TROSY NMR.

Large macromolecular assemblies, so-called molecular machines, are critical to ensuring proper cellular function. Understanding how proper function is achieved at the atomic level is crucial to advancing multiple avenues of biomedical research. Biophysical studies often include X-ray diffraction and cryo-electron microscopy, providing detailed structural descriptions of these machines. However, their inherent flexibility has complicated an understanding of the relation between structure and function. Solution NMR spectroscopy is well suited to the study of such dynamic complexes, and continued developments have increased size boundaries; insights into function have been obtained for complexes with masses as large as 1 MDa. We highlight methyl-TROSY (transverse relaxation optimized spectroscopy) NMR, which enables the study of such large systems, and include examples of applications to several cellular machines. We show how this emerging technique contributes to an understanding of cellular function and the role of molecular plasticity in regulating an array of biochemical activities.

Enteric viruses replicate in salivary glands and infect through saliva.

Enteric viruses like norovirus, rotavirus and astrovirus have long been accepted as spreading in the population through fecal-oral transmission: viruses are shed into feces from one host and enter the oral cavity of another, bypassing salivary glands (SGs) and reaching the intestines to replicate, be shed in feces and repeat the transmission cycle1. Yet there are viruses (for example, rabies) that infect the SGs2,3, making the oral cavity one site of replication and saliva one conduit of transmission. Here we report that enteric viruses productively and persistently infect SGs, reaching titres comparable to those in the intestines. We demonstrate that enteric viruses get released into the saliva, identifying a second route of viral transmission. This is particularly significant for infected infants, whose saliva directly transmits enteric viruses to their mothers' mammary glands through backflow during suckling. This sidesteps the conventional gut-mammary axis route4 and leads to a rapid surge in maternal milk secretory IgA antibodies5,6. Lastly, we show that SG-derived spheroids7 and cell lines8 can replicate and propagate enteric viruses, generating a scalable and manageable system of production. Collectively, our research uncovers a new transmission route for enteric viruses with implications for therapeutics, diagnostics and importantly sanitation measures to prevent spread through saliva.

Design of amyloidogenic peptide traps.

Segments of proteins with high ß-strand propensity can self-associate to form amyloid fibrils implicated in many diseases. We describe a general approach to bind such segments in ß-strand and ß-hairpin conformations using de novo designed scaffolds that contain deep peptide-binding clefts. The designs bind their cognate peptides in vitro with nanomolar affinities. The crystal structure of a designed protein-peptide complex is close to the design model, and NMR characterization reveals how the peptide-binding cleft is protected in the apo state. We use the approach to design binders to the amyloid-forming proteins transthyretin, tau, serum amyloid A1 and amyloid ß1-42 (Aß42). The Aß binders block the assembly of Aß fibrils as effectively as the most potent of the clinically tested antibodies to date and protect cells from toxic Aß42 species.

The Importance of Poly(ADP-Ribose) Polymerase as a Sensor of Unligated Okazaki Fragments during DNA Replication.

Poly(ADP-ribose) is synthesized by PARP enzymes during the repair of stochastic DNA breaks. Surprisingly, however, we show that most if not all endogenous poly(ADP-ribose) is detected in normal S phase cells at sites of DNA replication. This S phase poly(ADP-ribose) does not result from damaged or misincorporated nucleotides or from DNA replication stress. Rather, perturbation of the DNA replication proteins LIG1 or FEN1 increases S phase poly(ADP-ribose) more than 10-fold, implicating unligated Okazaki fragments as the source of S phase PARP activity. Indeed, S phase PARP activity is ablated by suppressing Okazaki fragment formation with emetine, a DNA replication inhibitor that selectively inhibits lagging strand synthesis. Importantly, PARP activation during DNA replication recruits the single-strand break repair protein XRCC1, and human cells lacking PARP activity and/or XRCC1 are hypersensitive to FEN1 perturbation. Collectively, our data indicate that PARP1 is a sensor of unligated Okazaki fragments during DNA replication and facilitates their repair.

Exploiting conformational dynamics to modulate the function of designed proteins.

With the recent success in calculating protein structures from amino acid sequences using artificial intelligence-based algorithms, an important next step is to decipher how dynamics is encoded by the primary protein sequence so as to better predict function. Such dynamics information is critical for protein design, where strategies could then focus not only on sequences that fold into particular structures that perform a given task, but would also include low-lying excited protein states that could influence the function of the designed protein. Herein, we illustrate the importance of dynamics in modulating the function of C34, a designed α/ß protein that captures ß-strands of target ligands and is a member of a family of proteins designed to sequester ß-strands and ß hairpins of aggregation-prone molecules that lead to a variety of pathologies. Using a strategy to "see" regions of apo C34 that are invisible to NMR spectroscopy as a result of pervasive conformational exchange, as well as a mutagenesis approach whereby C34 molecules are stabilized into a single conformer, we determine the structures of the predominant conformations that are sampled by C34 and show that these attenuate the affinity for cognate peptide. Subsequently, the observed motion is exploited to develop an allosterically regulated peptide binder whose binding affinity can be controlled through the addition of a second molecule. Our study emphasizes the unique role that NMR can play in directing the design process and in the construction of new molecules with more complex functionality.

Correlating histone acetylation with nucleosome core particle dynamics and function.

Epigenetic modifications of chromatin play a critical role in regulating the fidelity of the genetic code and in controlling the translation of genetic information into the protein components of the cell. One key posttranslational modification is acetylation of histone lysine residues. Molecular dynamics simulations, and to a smaller extent experiment, have established that lysine acetylation increases the dynamics of histone tails. However, a systematic, atomic resolution experimental investigation of how this epigenetic mark, focusing on one histone at a time, influences the structural dynamics of the nucleosome beyond the tails, and how this translates into accessibility of protein factors such as ligases and nucleases, has yet to be performed. Herein, using NMR spectroscopy of nucleosome core particles (NCPs), we evaluate the effects of acetylation of each histone on tail and core dynamics. We show that for histones H2B, H3, and H4, the histone core particle dynamics are little changed, even though the tails have increased amplitude motions. In contrast, significant increases to H2A dynamics are observed upon acetylation of this histone, with the docking domain and L1 loop particularly affected, correlating with increased susceptibility of NCPs to nuclease digestion and more robust ligation of nicked DNA. Dynamic light scattering experiments establish that acetylation decreases inter-NCP interactions in a histone-dependent manner and facilitates the development of a thermodynamic model for NCP stacking. Our data show that different acetylation patterns result in nuanced changes to NCP dynamics, modulating interactions with other protein factors, and ultimately controlling biological output.

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.

Structural basis of protein substrate processing by human mitochondrial high-temperature requirement A2 protease.

The human high-temperature requirement A2 (HtrA2) protein is a trimeric protease that cleaves misfolded proteins to protect cells from stresses caused by toxic, proteinaceous aggregates, and the aberrant function of HtrA2 is closely related to the onset of neurodegenerative disorders. Our methyl-transverse relaxation optimized spectroscopy (TROSY)­based NMR studies using small-peptide ligands have previously revealed a stepwise activation mechanism involving multiple distinct conformational states. However, very little is known about how HtrA2 binds to protein substrates and if the distinct conformational states observed in previous peptide studies might be involved in the processing of protein clients. Herein, we use solution-based NMR spectroscopy to investigate the interaction between the N-terminal Src homology 3 domain from downstream of receptor kinase (drk) with an added C-terminal HtrA2-binding motif (drkN SH3-PDZbm) that exhibits marginal folding stability and serves as a mimic of a physiological protein substrate. We show that drkN SH3-PDZbm binds to HtrA2 via a two-pronged interaction, involving both its C-terminal PDZ-domain binding motif and a central hydrophobic region, with binding occurring preferentially via an unfolded ensemble of substrate molecules. Multivalent interactions between several clients and a single HtrA2 trimer significantly stimulate the catalytic activity of HtrA2, suggesting that binding avidity plays an important role in regulating substrate processing. Our results provide a thermodynamic, kinetic, and structural description of the interaction of HtrA2 with protein substrates and highlight the importance of a trimeric architecture for function as a stress-protective protease that mitigates aggregation.

Mapping the per-residue surface electrostatic potential of CAPRIN1 along its phase-separation trajectory.

Electrostatic interactions and charge balance are important for the formation of biomolecular condensates involving proteins and nucleic acids. However, a detailed, atomistic picture of the charge distribution around proteins during the phase-separation process is lacking. Here, we use solution NMR spectroscopy to measure residue-specific near-surface electrostatic potentials (ϕENS) of the positively charged carboxyl-terminal intrinsically disordered 103 residues of CAPRIN1, an RNA-binding protein localized to membraneless organelles playing an important role in messenger RNA (mRNA) storage and translation. Measured ϕENS values have been mapped along the adenosine triphosphate (ATP)-induced phase-separation trajectory. In the absence of ATP, ϕENS values for the mixed state of CAPRIN1 are positive and large and progressively decrease as ATP is added. This is coupled to increasing interchain interactions, particularly between aromatic-rich and arginine-rich regions of the protein. Upon phase separation, CAPRIN1 molecules in the condensed phase are neutral (ϕENS [Formula: see text] 0 mV), with ∼five molecules of ATP associated with each CAPRIN1 chain. Increasing the ATP concentration further inverts the CAPRIN1 electrostatic potential, so that molecules become negatively charged, especially in aromatic-rich regions, leading to re-entrance into a mixed phase. Our results collectively show that a subtle balance between electrostatic repulsion and interchain attractive interactions regulates CAPRIN1 phase separation and provides insight into how nucleotides, such as ATP, can induce formation of and subsequently dissolve protein condensates.

Cryo-EM of the Nucleosome Core Particle Bound to Ran-RCC1 Reveals a Dynamic Complex.

Ras-related nuclear protein (Ran) is a member of the Ras superfamily of small guanosine triphosphatases (GTPases) and a regulator of multiple cellular processes. In healthy cells, the GTP-bound form of Ran is concentrated at chromatin, creating a Ran•GTP gradient that provides the driving force for nucleocytoplasmic transport, mitotic spindle assembly, and nuclear envelope formation. The Ran•GTP gradient is maintained by the regulator of chromatin condensation 1 (RCC1), a guanine nucleotide exchange factor that accelerates GDP/GTP exchange in Ran. RCC1 interacts with nucleosomes, which are the fundamental repeating units of eukaryotic chromatin. Here, we present a cryo-EM analysis of a trimeric complex composed of the nucleosome core particle (NCP), RCC1, and Ran. While the contacts between RCC1 and Ran in the complex are preserved compared with a previously determined structure of RCC1-Ran, our study reveals that RCC1 and Ran interact dynamically with the NCP and undergo rocking motions on the nucleosome surface. Furthermore, the switch 1 region of Ran, which plays an important role in mediating conformational changes associated with the substitution of GDP and GTP nucleotides in Ras family members, appears to undergo disorder-order transitions and forms transient contacts with the C-terminal helix of histone H2B. Nucleotide exchange assays performed in the presence and absence of NCPs are not consistent with an active role for nucleosomes in nucleotide exchange, at least in vitro. Instead, the nucleosome stabilizes RCC1 and serves as a hub that concentrates RCC1 and Ran to promote efficient Ran•GDP to Ran•GTP conversion.

Surface electrostatics dictate RNA-binding protein CAPRIN1 condensate concentration and hydrodynamic properties.

Biomolecular condensates concentrate proteins, nucleic acids, and small molecules and play an essential role in many biological processes. Their formation is tuned by a balance between energetically favorable and unfavorable contacts, with charge-charge interactions playing a central role in some systems. The positively charged intrinsically disordered carboxy-terminal region of the RNA-binding protein CAPRIN1 is one such example, phase separating upon addition of negatively charged ATP or high concentrations of sodium chloride (NaCl). Using solution NMR spectroscopy, we measured residue-specific near-surface electrostatic potentials (ϕENS) of CAPRIN1 along its NaCl-induced phase separation trajectory to compare with those obtained using ATP. In both cases, electrostatic shielding decreases ϕENS values, yet surface potentials of CAPRIN1 in the two condensates can be different, depending on the amount of NaCl or ATP added. Our results establish that even small differences in ϕENS can significantly affect the level of protein enrichment and the mechanical properties of the condensed phase, leading, potentially, to the regulation of biological processes.

Exploring Host-Guest Interactions within a 600 kDa DegP Protease Cage Complex Using Hydrodynamics Measurements and Methyl-TROSY NMR.

The DegP protease-chaperone operates within the periplasm of Gram-negative bacteria, where it assists in the regulation of protein homeostasis, promotes virulence, and is essential to survival under stress. To carry out these tasks, DegP forms a network of preorganized apo oligomers that facilitate the capture of substrates within distributions of cage-like complexes which expand to encapsulate clients of various sizes. Although the architectures of DegP cage complexes are well understood, little is known about the structures, dynamics, and interactions of client proteins within DegP cages and the relationship between client structural dynamics and function. Here, we probe host-guest interactions within a 600 kDa DegP cage complex throughout the DegP activation cycle using a model α-helical client protein through a combination of hydrodynamics measurements, methyl-transverse relaxation optimized spectroscopy-based solution nuclear magnetic resonance studies, and proteolytic activity assays. We find that in the presence of the client, DegP cages assemble cooperatively with few intermediates. Our data further show that the N-terminal half of the bound client, which projects into the interior of the cages, is predominantly unfolded and flexible, and exchanges between multiple conformational states over a wide range of time scales. Finally, we show that a concerted structural transition of the protease domains of DegP occurs upon client engagement, leading to activation. Together, our findings support a model of DegP as a highly cooperative and dynamic molecular machine that stabilizes unfolded states of clients, primarily via interactions with their C-termini, giving rise to efficient cleavage.

Phase Separation Modulates the Thermodynamics and Kinetics of RNA Hybridization.

Biomolecular condensates can influence cellular function in a number of ways, including by changing the structural dynamics and conformational equilibria of the molecules partitioned within them. Here we use methyl transverse relaxation optimized spectroscopy (methyl-TROSY) NMR in conjunction with 2'-O-methyl labeling of RNA to characterize the thermodynamics and kinetics of RNA-RNA base pairing in condensates formed by the C-terminal intrinsically disordered region of CAPRIN1, an RNA-binding protein involved in RNA transport, translation, and stability. CAPRIN1 condensates destabilize RNA-RNA base pairing, resulting from a ∼270-fold decrease and a concomitant ∼15-fold increase in the on- and off-rates for duplex formation, respectively. The ∼30-fold slower diffusion of RNA single strands within the condensed phase partially accounts for the reduced on-rate, but the further ∼9-fold reduction likely reflects shedding of CAPRIN1 chains that are interacting with the RNA prior to hybridization. Our study emphasizes the important role of protein solvation in modulating nucleic acid recognition processes inside condensates.

Beyond slow two-state protein conformational exchange using CEST: applications to three-state protein interconversion on the millisecond timescale.

Although NMR spectroscopy is routinely used to study the conformational dynamics of biomolecules, robust analyses of the data are challenged in cases where exchange is more complex than two-state, such as when a 'visible' major conformer exchanges with two 'invisible' minor states on the millisecond timescale. It is becoming increasingly clear that chemical exchange saturation transfer (CEST) NMR experiments that were initially developed to study systems undergoing slow interconversion are also sensitive to intermediate-fast timescale biomolecular conformational exchange. Here we investigate the utility of the amide 15N CEST experiment to characterise protein three-state exchange occurring on the millisecond timescale by studying the interconversion between the folded (F) state of the FF domain from human HYPA/FBP11 (WT FF) and two of its folding intermediates I1 and I2. Although 15N CPMG experiments are consistent with the F state interconverting with a single minor state on the millisecond timescale, 15N CEST data clearly establish an exchange process between F and a pair of minor states. A unique three-state exchange model cannot be obtained by analysis of 15N CEST data recorded at a single temperature. However, including the relative sign of the difference in the chemical shifts of the two minor states based on a simple two-state analysis of CEST data recorded at multiple temperatures, results in a robust three-state model in which the F, I1 and I2 states interconvert with each other on the millisecond timescale ( k e x , F I 1 ~ 550 s-1, k e x , F I 2 ~ 1200 s-1, k e x , I 1 I 2 ~ 5000 s-1), with I1 and I2 sparsely populated at ~ 0.15% and ~ 0.35%, respectively, at 15 °C. A computationally demanding grid-search of exchange parameter space is not required to extract the best-fit exchange parameters from the CEST data. The utility of the CEST experiment, thus, extends well beyond studies of conformers in slow exchange on the NMR chemical shift timescale, to include systems with interconversion rates on the order of thousands/second.

Predictors at 6 and 12 Months for Social Participation Outcome at 24 Months in the Adult Burn Injury Population: A Burn Model System National Database Study.

OBJECTIVE: To identify clinical factors (physical and psychological symptoms and post-traumatic growth) that predict social participation outcome at 24-month after burn injury. DESIGN: A prospective cohort study based on Burn Model System National Database. SETTING: Burn Model System centers. PARTICIPANTS: 181 adult participants less than 2 years after burn injury (N=181). INTERVENTIONS: Not applicable. MAIN OUTCOME MEASURES: Demographic and injury variables were collected at discharge. Predictor variables were assessed at 6 and 12 months: Post-Traumatic Growth Inventory Short Form (PTGI-SF), Post-Traumatic Stress Disorder Checklist Civilian Version (PCL-C), Patient-Reported Outcomes Measurement Information System (PROMIS-29) Depression, Anxiety, Sleep Disturbance, Fatigue, and Pain Interference short forms, and self-reported Heat Intolerance. Social participation was measured at 24 months using the Life Impact Burn Recovery Evaluation (LIBRE) Social Interactions and Social Activities short forms. RESULTS: Linear and multivariable regression models were used to examine predictor variables for social participation outcomes, controlling for demographic and injury variables. For LIBRE Social Interactions, significant predictors included the PCL-C total score at 6 months (ß=-0.27, P<.001) and 12 months (ß=-0.39, P<.001), and PROMIS-29 Pain Interference at 6 months (ß=-0.20, P<.01). For LIBRE Social Activities, significant predictors consisted of the PROMIS-29 Depression at 6 months (ß=-0.37, P<.001) and 12 months (ß=-0.37, P<.001), PROMIS-29 Pain Interference at 6 months (ß=-0.40, P<.001) and 12 months (ß=-0.37, P<.001), and Heat Intolerance at 12 months (ß=-4.55, P<.01). CONCLUSIONS: Post-traumatic stress and pain predicted social interactions outcomes, while depression, pain and heat intolerance predicted social activities outcomes in people with burn injury.