Visualizing protein breathing motions associated with aromatic ring flipping.

Aromatic residues cluster in the core of folded proteins, where they stabilize the structure through multiple interactions. Nuclear magnetic resonance (NMR) studies in the 1970s showed that aromatic side chains can undergo ring flips-that is, 180° rotations-despite their role in maintaining the protein fold1-3. It was suggested that large-scale 'breathing' motions of the surrounding protein environment would be necessary to accommodate these ring flipping events1. However, the structural details of these motions have remained unclear. Here we uncover the structural rearrangements that accompany ring flipping of a buried tyrosine residue in an SH3 domain. Using NMR, we show that the tyrosine side chain flips to a low-populated, minor state and, through a proteome-wide sequence analysis, we design mutants that stabilize this state, which allows us to capture its high-resolution structure by X-ray crystallography. A void volume is generated around the tyrosine ring during the structural transition between the major and minor state, and this allows fast flipping to take place. Our results provide structural insights into the protein breathing motions that are associated with ring flipping. More generally, our study has implications for protein design and structure prediction by showing how the local protein environment influences amino acid side chain conformations and vice versa.

NMR Provides Unique Insight into the Functional Dynamics and Interactions of Intrinsically Disordered Proteins.

Intrinsically disordered proteins are ubiquitous throughout all known proteomes, playing essential roles in all aspects of cellular and extracellular biochemistry. To understand their function, it is necessary to determine their structural and dynamic behavior and to describe the physical chemistry of their interaction trajectories. Nuclear magnetic resonance is perfectly adapted to this task, providing ensemble averaged structural and dynamic parameters that report on each assigned resonance in the molecule, unveiling otherwise inaccessible insight into the reaction kinetics and thermodynamics that are essential for function. In this review, we describe recent applications of NMR-based approaches to understanding the conformational energy landscape, the nature and time scales of local and long-range dynamics and how they depend on the environment, even in the cell. Finally, we illustrate the ability of NMR to uncover the mechanistic basis of functional disordered molecular assemblies that are important for human health.

Liquid-Liquid Phase Separation Modifies the Dynamic Properties of Intrinsically Disordered Proteins.

Liquid-liquid phase separation of flexible biomolecules has been identified as a ubiquitous phenomenon underlying the formation of membraneless organelles that harbor a multitude of essential cellular processes. We use nuclear magnetic resonance (NMR) spectroscopy to compare the dynamic properties of an intrinsically disordered protein (measles virus NTAIL) in the dilute and dense phases at atomic resolution. By measuring 15N NMR relaxation at different magnetic field strengths, we are able to characterize the dynamics of the protein in dilute and crowded conditions and to compare the amplitude and timescale of the different motional modes to those present in the membraneless organelle. Although the local backbone conformational sampling appears to be largely retained, dynamics occurring on all detectable timescales, including librational, backbone dihedral angle dynamics and segmental, chainlike motions, are considerably slowed down. Their relative amplitudes are also drastically modified, with slower, chain-like motions dominating the dynamic profile. In order to provide additional mechanistic insight, we performed extensive molecular dynamics simulations of the protein under self-crowding conditions at concentrations comparable to those found in the dense liquid phase. Simulation broadly reproduces the impact of formation of the condensed phase on both the free energy landscape and the kinetic interconversion between states. In particular, the experimentally observed reduction in the amplitude of the fastest component of backbone dynamics correlates with higher levels of intermolecular contacts or entanglement observed in simulations, reducing the conformational space available to this mode under strongly self-crowding conditions.

Convergent views on disordered protein dynamics from NMR and computational approaches.

Intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs) is a class of biologically important proteins exhibiting specific biophysical characteristics. They lack a hydrophobic core, and their conformational behavior is strongly influenced by electrostatic interactions. IDPs and IDRs are highly dynamic, and a characterization of the motions of IDPs and IDRs is essential for their physically correct description. NMR together with molecular dynamics simulations are the methods best suited to such a task because they provide information about dynamics of proteins with atomistic resolution. Here, we present a study of motions of a disordered C-terminal domain of the delta subunit of RNA polymerase from Bacillus subtilis. Positively and negatively charged residues in the studied domain form transient electrostatic contacts critical for the biological function. Our study is focused on investigation of ps-ns dynamics of backbone of the delta subunit based on analysis of amide 15N NMR relaxation data and molecular dynamics simulations. In order to extend an informational content of NMR data to lower frequencies, which are more sensitive to slower motions, we combined standard (high-field) NMR relaxation experiments with high-resolution relaxometry. Altogether, we collected data reporting the relaxation at 12 different magnetic fields, resulting in an unprecedented data set. Our results document that the analysis of such data provides a consistent description of dynamics and confirms the validity of so far used protocols of the analysis of dynamics of IDPs also for a partially folded protein. In addition, the potential to access detailed description of motions at the timescale of tens of ns with the help of relaxometry data is discussed. Interestingly, in our case, it appears to be mostly relevant for a region involved in the formation of temporary contacts within the disordered region, which was previously proven to be biologically important.

Intrinsically Disordered Tardigrade Proteins Self-Assemble into Fibrous Gels in Response to Environmental Stress.

Tardigrades are remarkable for their ability to survive harsh stress conditions as diverse as extreme temperature and desiccation. The molecular mechanisms that confer this unusual resistance to physical stress remain unknown. Recently, tardigrade-unique intrinsically disordered proteins have been shown to play an essential role in tardigrade anhydrobiosis. Here, we characterize the conformational and physical behaviour of CAHS-8 from Hypsibius exemplaris. NMR spectroscopy reveals that the protein comprises an extended central helical domain flanked by disordered termini. Upon concentration, the protein is shown to successively form oligomers, long fibres, and finally gels constituted of fibres in a strongly temperature-dependent manner. The helical domain forms the core of the fibrillar structure, with the disordered termini remaining highly dynamic within the gel. Soluble proteins can be encapsulated within cavities in the gel, maintaining their functional form. The ability to reversibly form fibrous gels may be associated with the enhanced protective properties of these proteins.

A Unified Description of Intrinsically Disordered Protein Dynamics under Physiological Conditions Using NMR Spectroscopy.

Intrinsically disordered proteins (IDPs) are flexible biomolecules whose essential functions are defined by their dynamic nature. Nuclear magnetic resonance (NMR) spectroscopy is ideally suited to the investigation of this behavior at atomic resolution. NMR relaxation is increasingly used to detect conformational dynamics in free and bound forms of IDPs under conditions approaching physiological, although a general framework providing a quantitative interpretation of these exquisitely sensitive probes as a function of experimental conditions is still lacking. Here, measuring an extensive set of relaxation rates sampling multiple-time-scale dynamics over a broad range of crowding conditions, we develop and test an integrated analytical description that accurately portrays the motion of IDPs as a function of the intrinsic properties of the crowded molecular environment. In particular we observe a strong dependence of both short-range and long-range motional time scales of the protein on the friction of the solvent. This tight coupling between the dynamic behavior of the IDP and its environment allows us to develop analytical expressions for protein motions and NMR relaxation properties that can be accurately applied over a vast range of experimental conditions. This unified dynamic description provides new insight into the physical behavior of IDPs, extending our ability to quantitatively investigate their conformational dynamics under complex environmental conditions, and accurately predicting relaxation rates reporting on motions on time scales up to tens of nanoseconds, both in vitro and in cellulo.

Deciphering the Dynamic Interaction Profile of an Intrinsically Disordered Protein by NMR Exchange Spectroscopy.

Intrinsically disordered proteins (IDPs) display a large number of interaction modes including folding-upon-binding, binding without major structural transitions, or binding through highly dynamic, so-called fuzzy, complexes. The vast majority of experimental information about IDP binding modes have been inferred from crystal structures of proteins in complex with short peptides of IDPs. However, crystal structures provide a mainly static view of the complexes and do not give information about the conformational dynamics experienced by the IDP in the bound state. Knowledge of the dynamics of IDP complexes is of fundamental importance to understand how IDPs engage in highly specific interactions without concomitantly high binding affinity. Here, we combine rotating-frame R1ρ, Carr-Purcell-Meiboom Gill relaxation dispersion as well as chemical exchange saturation transfer to decipher the dynamic interaction profile of an IDP in complex with its partner. We apply the approach to the dynamic signaling complex formed between the mitogen-activated protein kinase (MAPK) p38α and the intrinsically disordered regulatory domain of the MAPK kinase MKK4. Our study demonstrates that MKK4 employs a subtle combination of interaction modes in order to bind to p38α, leading to a complex displaying significantly different dynamics across the bound regions.

The membrane anchor of the transcriptional activator SREBP is characterized by intrinsic conformational flexibility.

Regulated intramembrane proteolysis (RIP) is a conserved mechanism crucial for numerous cellular processes, including signaling, transcriptional regulation, axon guidance, cell adhesion, cellular stress responses, and transmembrane protein fragment degradation. Importantly, it is relevant in various diseases including Alzheimer's disease, cardiovascular diseases, and cancers. Even though a number of structures of different intramembrane proteases have been solved recently, fundamental questions concerning mechanistic underpinnings of RIP and therapeutic interventions remain. In particular, this includes substrate recognition, what properties render a given substrate amenable for RIP, and how the lipid environment affects the substrate cleavage. Members of the sterol regulatory element-binding protein (SREBP) family of transcription factors are critical regulators of genes involved in cholesterol/lipid homeostasis. After site-1 protease cleavage of the inactive SREBP transmembrane precursor protein, RIP of the anchor intermediate by site-2 protease generates the mature transcription factor. In this work, we have investigated the labile anchor intermediate of SREBP-1 using NMR spectroscopy. Surprisingly, NMR chemical shifts, site-resolved solvent exposure, and relaxation studies show that the cleavage site of the lipid-signaling protein intermediate bears rigid α-helical topology. An evolutionary conserved motif, by contrast, interrupts the secondary structure ∼9-10 residues C-terminal of the scissile bond and acts as an inducer of conformational flexibility within the carboxyl-terminal transmembrane region. These results are consistent with molecular dynamics simulations. Topology, stability, and site-resolved dynamics data suggest that the cleavage of the α-helical substrate in the case of RIP may be associated with a hinge motion triggered by the molecular environment.

Dynamic Descriptions of Highly Flexible Molecules from NMR Dipolar Couplings: Physical Basis and Limitations.

Biomolecules that control physiological function by changing their conformation play key roles in biology and remain poorly characterized. NMR dipolar couplings (DCs) depend intrinsically on both molecular shape and structural fluctuations, thereby providing the enticing prospect of tracking these conformational changes at atomic detail. Although this dual dependence has until now severely complicated analysis of DCs from highly dynamic systems, general approaches have recently been proposed that simplify interpretation of experimental DCs, by entirely eliminating molecular alignment from the analysis. Using simple and intuitive simulation of target ensembles, we investigate the impact of such approaches on the resulting descriptions of the conformational energy landscape. We find that ensemble descriptions of highly flexible systems derived from DCs without explicit consideration of the alignment properties of the constituent conformations can be compromised and inaccurate, despite exhibiting high correlation with experimental measurement.

Structure of the eukaryotic translation initiation factor eIF4E in complex with 4EGI-1 reveals an allosteric mechanism for dissociating eIF4G.

The interaction of the eukaryotic translation initiation factor eIF4E with the initiation factor eIF4G recruits the 40S ribosomal particle to the 5' end of mRNAs, facilitates scanning to the AUG start codon, and is crucial for eukaryotic translation of nearly all genes. Efficient recruitment of the 40S particle is particularly important for translation of mRNAs encoding oncoproteins and growth-promoting factors, which often harbor complex 5' UTRs and require efficient initiation. Thus, inhibiting the eIF4E/eIF4G interaction has emerged as a previously unpursued route for developing anticancer agents. Indeed, we discovered small-molecule inhibitors of this eIF4E/eIF4G interaction (4EGIs) that inhibit translation initiation both in vitro and in vivo and were used successfully in numerous cancer-biology and neurobiology studies. However, their detailed molecular mechanism of action has remained elusive. Here, we show that the eIF4E/eIF4G inhibitor 4EGI-1 acts allosterically by binding to a site on eIF4E distant from the eIF4G binding epitope. Data from NMR mapping and high-resolution crystal structures are congruent with this mechanism, where 4EGI-1 attaches to a hydrophobic pocket of eIF4E between ß-sheet2 (L60-T68) and α-helix1 (E69-N77), causing localized conformational changes mainly in the H78-L85 region. It acts by unfolding a short 310-helix (S82-L85) while extending α-helix1 by one turn (H78-S82). This unusual helix rearrangement has not been seen in any previous eIF4E structure and reveals elements of an allosteric inhibition mechanism leading to the dislocation of eIF4G from eIF4E.

Analytical Description of NMR Relaxation Highlights Correlated Dynamics in Intrinsically Disordered Proteins.

The dynamic fluctuations of intrinsically disordered proteins (IDPs) define their function. Although experimental nuclear magnetic resonance (NMR) relaxation reveals the motional complexity of these highly flexible proteins, the absence of physical models describing IDP dynamics hinders their mechanistic interpretation. Combining molecular dynamics simulation and NMR, we introduce a framework in which distinct motions are attributed to local libration, backbone dihedral angle dynamics and longer-range tumbling of one or more peptide planes. This model provides unique insight into segmental organization of dynamics in IDPs and allows us to investigate the presence and extent of the correlated motions that are essential for function.

Identification of Dynamic Modes in an Intrinsically Disordered Protein Using Temperature-Dependent NMR Relaxation.

The dynamic modes and time scales sampled by intrinsically disordered proteins (IDPs) define their function. Nuclear magnetic resonance (NMR) spin relaxation is probably the most powerful tool for investigating these motions delivering site-specific descriptions of conformational fluctuations from throughout the molecule. Despite the abundance of experimental measurement of relaxation in IDPs, the physical origin of the measured relaxation rates remains poorly understood. Here we measure an extensive range of auto- and cross-correlated spin relaxation rates at multiple magnetic field strengths on the C-terminal domain of the nucleoprotein of Sendai virus, over a large range of temperatures (268-298 K), and combine these data to describe the dynamic behavior of this archetypal IDP. An Arrhenius-type relationship is used to simultaneously analyze up to 61 relaxation rates per amino acid over the entire temperature range, allowing the measurement of local activation energies along the chain, and the assignment of physically distinct dynamic modes. Fast (τ ≤ 50 ps) components report on librational motions, a dominant mode occurs on time scales around 1 ns, apparently reporting on backbone sampling within Ramachandran substates, while a slower component (5-25 ns) reports on segmental dynamics dominated by the chain-like nature of the protein. Extending the study to three protein constructs of different lengths (59, 81, and 124 amino acids) substantiates the assignment of these contributions. The analysis is shown to be remarkably robust, accurately predicting a broad range of relaxation data measured at different magnetic field strengths and temperatures. The ability to delineate intrinsic modes and time scales from NMR spin relaxation will improve our understanding of the behavior and function of IDPs, adding a new and essential dimension to the description of this biologically important and ubiquitous class of proteins.

The Role of Dynamics and Allostery in the Inhibition of the eIF4E/eIF4G Translation Initiation Factor Complex.

Lack of regulation of the interaction between the eIF4E/eIF4G subunits of the translation initiation factor complex eIF4F is a hallmark of cancer. The inhibitor 4EGI-1 binds to eIF4E, thereby preventing association with eIF4G through an allosteric mechanism. NMR spectroscopy and MD simulations were used to obtain a mechanistic description of the role of correlated dynamics in this allosteric regulation. We show that binding of 4EGI-1 perturbs native correlated motions and increases correlated fluctuations in part of the eIF4G binding site.

Distribution of Pico- and Nanosecond Motions in Disordered Proteins from Nuclear Spin Relaxation.

Intrinsically disordered proteins and intrinsically disordered regions (IDRs) are ubiquitous in the eukaryotic proteome. The description and understanding of their conformational properties require the development of new experimental, computational, and theoretical approaches. Here, we use nuclear spin relaxation to investigate the distribution of timescales of motions in an IDR from picoseconds to nanoseconds. Nitrogen-15 relaxation rates have been measured at five magnetic fields, ranging from 9.4 to 23.5 T (400-1000 MHz for protons). This exceptional wealth of data allowed us to map the spectral density function for the motions of backbone NH pairs in the partially disordered transcription factor Engrailed at 11 different frequencies. We introduce an approach called interpretation of motions by a projection onto an array of correlation times (IMPACT), which focuses on an array of six correlation times with intervals that are equidistant on a logarithmic scale between 21 ps and 21 ns. The distribution of motions in Engrailed varies smoothly along the protein sequence and is multimodal for most residues, with a prevalence of motions around 1 ns in the IDR. We show that IMPACT often provides better quantitative agreement with experimental data than conventional model-free or extended model-free analyses with two or three correlation times. We introduce a graphical representation that offers a convenient platform for a qualitative discussion of dynamics. Even when relaxation data are only acquired at three magnetic fields that are readily accessible, the IMPACT analysis gives a satisfactory characterization of spectral density functions, thus opening the way to a broad use of this approach.

Large-Scale Conformational Dynamics Control H5N1 Influenza Polymerase PB2 Binding to Importin α.

Influenza A RNA polymerase complex is formed from three components, PA, PB1, and PB2. PB2 is independently imported into the nucleus prior to polymerase reconstitution. All crystallographic structures of the PB2 C-terminus (residues 536-759) reveal two globular domains, 627 and NLS, that form a tightly packed heterodimer. The molecular basis of the affinity of 627-NLS for importins remained unclear from these structures, apparently requiring large-scale conformational changes prior to importin binding. Using a combination of solution-state NMR, small-angle neutron scattering, small-angle X-ray scattering (SAXS), and Förster resonance energy transfer (FRET), we show that 627-NLS populates a temperature-dependent dynamic equilibrium between closed and open states. The closed state is stabilized by a tripartite salt bridge involving the 627-NLS interface and the linker, that becomes flexible in the open state, with 627 and NLS dislocating into a highly dynamic ensemble. Activation enthalpies and entropies associated with the rupture of this interface were derived from simultaneous analysis of temperature-dependent chemical exchange saturation transfer measurements, revealing a strong temperature dependence of both open-state population and exchange rate. Single-molecule FRET and SAXS demonstrate that only the open-form is capable of binding to importin α and that, upon binding, the 627 domain samples a dynamic conformational equilibrium in the vicinity of the C-terminus of importin α. This intrinsic large-scale conformational flexibility therefore enables 627-NLS to bind importin through conformational selection from a temperature-dependent equilibrium comprising both functional forms of the protein.

Cross-correlated relaxation measurements under adiabatic sweeps: determination of local order in proteins.

Adiabatically swept pulses were originally designed for the purpose of broadband spin inversion. Later, unexpected advantages of their utilization were also found in other applications, such as refocusing to excite spin echoes, studies of chemical exchange or fragment-based drug design. Here, we present new experiments to characterize fast (ps-ns) protein dynamics, which benefit from little-known properties of adiabatic pulses. We developed a strategy for measuring cross-correlated cross-relaxation (CCCR) rates during adiabatic pulses. This experiment provides a linear combination of longitudinal and transverse CCCR rates, which is offset-independent across a typical amide (15)N spectrum. The pulse sequence can be recast to provide accurate transverse CCCR rates weighted by the populations of exchanging states. Sensitivity can be improved in systems in slow exchange. Finally, the experiments can be easily modified to yield residue-specific correlation times. The average correlation time of motions can be determined with a single experiment while at least two different experiments had to be recorded until now.

Challenges in preparing, preserving and detecting para-water in bulk: overcoming proton exchange and other hurdles.

Para-water is an analogue of para-hydrogen, where the two proton spins are in a quantum state that is antisymmetric under permutation, also known as singlet state. The populations of the nuclear spin states in para-water are believed to have long lifetimes just like other Long-Lived States (LLSs). This hypothesis can be verified by measuring the relaxation of an excess or a deficiency of para-water, also known as a "Triplet-Singlet Imbalance" (TSI), i.e., a difference between the average population of the three triplet states T (that are symmetric under permutation) and the population of the singlet state S. In analogy with our recent findings on ethanol and fumarate, we propose to adapt the procedure for Dissolution Dynamic Nuclear Polarization (D-DNP) to prepare such a TSI in frozen water at very low temperatures in the vicinity of 1.2 K. After rapid heating and dissolution using an aprotic solvent, the TSI should be largely preserved. To assess this hypothesis, we studied the lifetime of water as a molecular entity when diluted in various solvents. In neat liquid H2O, proton exchange rates have been characterized by spin-echo experiments on oxygen-17 in natural abundance, with and without proton decoupling. One-dimensional exchange spectroscopy (EXSY) has been used to study proton exchange rates in H2O, HDO and D2O mixtures diluted in various aprotic solvents. In the case of 50 mM H2O in dioxane-d8, the proton exchange lifetime is about 20 s. After dissolving, one can observe this TSI by monitoring intensities in oxygen-17 spectra of H2O (if necessary using isotopically enriched samples) where the AX2 system comprising a "spy" oxygen A and two protons X2 gives rise to binomial multiplets only if the TSI vanishes. Alternatively, fast chemical addition to a suitable substrate (such as an activated aldehyde or ketone) can provide AX2 systems where a carbon-13 acts as a spy nucleus. Proton signals that relax to equilibrium with two distinct time constants can be considered as a hallmark of a TSI. We optimized several experimental procedures designed to preserve and reveal dilute para-water in bulk.

A specific phosphorylation-dependent conformational switch in SARS-CoV-2 nucleocapsid protein inhibits RNA binding.

The nucleocapsid protein of severe acute respiratory syndrome coronavirus 2 encapsidates the viral genome and is essential for viral function. The central disordered domain comprises a serine-arginine-rich (SR) region that is hyperphosphorylated in infected cells. This modification regulates function, although mechanistic details remain unknown. We use nuclear magnetic resonance to follow structural changes occurring during hyperphosphorylation by serine arginine protein kinase 1, glycogen synthase kinase 3, and casein kinase 1, that abolishes interaction with RNA. When eight approximately uniformly distributed sites have been phosphorylated, the SR domain binds the same interface as single-stranded RNA, resulting in complete inhibition of RNA binding. Phosphorylation by protein kinase A does not prevent RNA binding, indicating that the pattern resulting from physiologically relevant kinases is specific for inhibition. Long-range contacts between the RNA binding, linker, and dimerization domains are abrogated, phenomena possibly related to genome packaging and unpackaging. This study provides insight into the recruitment of specific host kinases to regulate viral function.

Time scales of slow motions in ubiquitin explored by heteronuclear double resonance.

Understanding how proteins function at the atomic level relies in part on a detailed characterization of their dynamics. Ubiquitin, a small single-domain protein, displays rich dynamic properties over a wide range of time scales. In particular, several regions of ubiquitin show the signature of chemical exchange, including the hydrophobic patch and the ß4-α2 loop, which are both involved in many interactions. Here, we use multiple-quantum relaxation techniques to identify the extent of chemical exchange in ubiquitin. We employ our recently developed heteronuclear double resonance method to determine the time scales of motions that give rise to chemical exchange. Dispersion profiles are obtained for the backbone NH(N) pairs of several residues in the hydrophobic patch and the ß4-α2 loop, as well as the C-terminus of helix α1. We show that a single time scale (ca. 50 µs) can be used to fit the data for most residues. Potential mechanisms for the propagation of motions and the possible extent of correlation of these motions are discussed.

Boosting the sensitivity of ligand-protein screening by NMR of long-lived states.

A new NMR method for the study of ligand-protein interactions exploits the unusual lifetimes of long-lived states (LLSs). The new method provides better contrast between bound and free ligands and requires a protein-ligand ratio ca. 25 times lower than for established T(1ρ) methods, thus saving on costly proteins. The new LLS method was applied to the screening of inhibitors of urokinase-type plasminogen activator (uPA), which is a prototypical target of cancer research. With only 10 µM protein, a dissociation constant (K(D)) of 180 ± 20 nM was determined for the strong ligand (inhibitor) UK-18, which can be compared with K(D) = 157 ± 39 nM determined by the established surface plasmon resonance method.