Nanosecond-Timescale Dynamics and Conformational Heterogeneity in Human GCK Regulation and Disease

Human glucokinase (GCK) is the prototypic example of an emerging class of proteins with allosteric-like behavior that originates from intrinsic polypeptide dynamics. High-resolution NMR investigations of GCK have elucidated millisecond-timescale dynamics underlying allostery. In contrast, faster motions have remained underexplored, hindering the development of a comprehensive model of cooperativity. Here, we map nanosecond-timescale dynamics and structural heterogeneity in GCK using a combination of unnatural amino acid incorporation, time-resolved fluorescence, and ¹⁹F nuclear magnetic resonance spectroscopy. We find that a probe inserted within the enzyme's intrinsically disordered loop samples multiple conformations in the unliganded state. Glucose binding and disease-associated mutations that suppress cooperativity alter the number and/or relative population of these states. Together, the nanosecond kinetics characterized here and the millisecond motions known to be essential for cooperativity provide a dynamical framework with which we address the origins of cooperativity and the mechanism of activated, hyperinsulinemia-associated, noncooperative variants.

Describing Complex Structure-Function Relationships in Biomolecules at Equilibrium

One of the great ambitions of structural biology is to describe structure-function relationships quantitatively. Statistical thermodynamics is a powerful, general tool for computing the behavior of biological macromolecules at equilibrium because it establishes a direct link between structure and function. Complex behavior emerges as equilibria of multiple reactions are coupled. Analytical treatment of linked equilibria scales poorly with increasing numbers of reactions and states as the algebraic constructs rapidly become unwieldy. We therefore developed a generalizable, but straightforward computational method to handle arbitrarily complex systems. To demonstrate this approach, we collected a multidimensional fluorescence landscape of an engineered fluorescent glucose biosensor and showed that its features could be modeled with ten intricately linked ligand-binding and conformational exchange reactions. This protein represents a minimalist model of sufficient complexity to encompass fundamental biomolecular structure-function relationships: two-state and multistate conformational ensembles, conformational hierarchies, osmolytes, coupling between different binding sites and coupling between ligand binding and conformational change. The successful fit of this complex, multifaceted system demonstrates generality of the method.

Large-Scale Biomolecular Conformational Transitions Explored by a Combined Elastic Network Model and Enhanced Sampling Molecular Dynamics

Biomolecules often undergo large-scale conformational transitions when carrying out their functions. However, it is still challenging for conventional molecular dynamics simulations to provide adequate structural dynamics information to interpret associated mechanisms. Here, we present a combined elastic network model and enhanced sampling-based strategy (iterANM-IaMD) by adopting iterANM to construct initial conformation space and enhanced sampling IaMD to explore the free energy landscape along specific large-scale conformational transitions. We applied this strategy to three functionally and structurally distinct proteins (adenylate kinase, calmodulin, and p38α kinase), which undergo striking conformational change upon ligand binding. The simulation results for both free and ligand-bound proteins show qualitative and quantitative agreement with existing studies, suggesting iterANM-IaMD as an accurate and efficient tool to investigate structural dynamics involved in complicated biological processes. Our work also provides insights into the relationship between the dynamics and functionality of biomolecules.

Enhancing NMR derived ensembles with kinetics on multiple timescales

Nuclear magnetic resonance (NMR) has the unique advantage of elucidating the structure and dynamics of biomolecules in solution at physiological temperatures, where they are in constant movement on timescales from picoseconds to milliseconds. Such motions have been shown to be critical for enzyme catalysis, allosteric regulation, and molecular recognition. With NMR being particularly sensitive to these timescales, detailed information about the kinetics can be acquired. However, nearly all methods of NMR-based biomolecular structure determination neglect kinetics, which introduces a large approximation to the underlying physics, limiting both structural resolution and the ability to accurately determine molecular flexibility. Here we present the Kinetic Ensemble approach that uses a hierarchy of interconversion rates between a set of ensemble members to rigorously calculate Nuclear Overhauser Effect (NOE) intensities. It can be used to simultaneously refine both temporal and structural coordinates. By generalizing ideas from the extended model free approach, the method can analyze the amplitudes and kinetics of motions anywhere along the backbone or side chains. Furthermore, analysis of a large set of crystal structures suggests that NOE data contains a surprising amount of high-resolution information that is better modeled using our approach. The Kinetic Ensemble approach provides the means to unify numerous types of experiments under a single quantitative framework and more fully characterize and exploit kinetically distinct protein states. While we apply the approach here to the protein ubiquitin and cross validate it with previously derived datasets, the approach can be applied to any protein for which NOE data is available.

Conformational sampling and kinetics changes across a non-Arrhenius break point in the enzyme thermolysin

Numerous studies have suggested a significant role that protein dynamics play in optimizing enzyme catalysis, and changes in conformational sampling offer a window to explore this role. Thermolysin from Bacillus thermoproteolyticus rokko, which is a heat-stable zinc metalloproteinase, serves here as a model system to study changes of protein function and conformational sampling across a temperature range of 16-36 °C. The temperature dependence of kinetics of thermolysin showed a biphasic transition at 26 °C that points to potential conformational and dynamic differences across this temperature. The non-Arrhenius behavior observed resembled results from previous studies of a thermophilic alcohol dehydrogenase enzyme, which also indicated a biphasic transition at ambient temperatures. To explore the non-Arrhenius behavior of thermolysin, room temperature crystallography was applied to characterize structural changes in a temperature range across the biphasic transition temperature. The alternate conformation of side chain fitting to electron density of a group of residues showed a higher variability in the temperature range from 26 to 29 °C, which indicated a change in conformational sampling that correlated with the non-Arrhenius break point.

Patterns of Dynamics Comprise a Conserved Evolutionary Trait

The importance of protein dynamics in function may suggest an evolutionary selection on large-scale protein motions. Here we systematically studied the dynamic characteristics in 2221 protein domains (58477 sequences) of the Pfam database. We defined the patterns of dynamics (PODs) based on the estimated NMR order parameters and the predicted degree of disorder, and found a significant correlation between them in families of both structured and disordered protein domains. We demonstrate that conservation of dynamic patterns frequently exceeds conservation of sequence and is comparable to the patterns of hydropathy and nonspecific interaction potential. Similarity of dynamic patterns is weakly correlated to structure similarity and to the degree of disorder. We illustrate that POD alignments could be applied to sequentially divergent or intrinsically disordered regions. We propose that patterns of dynamics comprise a conserved evolutionary trait, which could be used to infer evolutionary relationships as an alternative to sequence and structure.

Role of hydration water in the onset of protein structural dynamics

Proteins are the molecular workhorses in a living organism. Their 3D structures are animated by a multitude of equilibrium fluctuations and specific out-of-equilibrium motions that are required for proteins to be biologically active. When studied as a function of temperature, functionally relevant dynamics are observed at and above the so-called protein dynamical transition (~240 K) in hydrated, but not in dry proteins. In this review we present and discuss the main experimental and computational results that provided evidence for the dynamical transition, with a focus on the role of hydration water dynamics in sustaining functional protein dynamics. The coupling and mutual influence of hydration water dynamics and protein dynamics are discussed and the hypotheses illustrated that have been put forward to explain the physical origin of their onsets.

Unsaturated lipid bilayers at cryogenic temperature: librational dynamics of chain-labeled lipids from pulsed and CW-EPR

Fully hydrated bilayers of monounsaturated palmitoyloleoylphosphatidylcholine (POPC) and diunsaturated dioleoylphosphatidylcholine (DOPC) lipids have low main phase transition temperatures (271 K for POPC and 253 K for DOPC). Two-pulse echo detected spectra, combined with continuous wave electron paramagnetic resonance spectroscopy, are employed to study the low-temperature lamellar phases of the POPC and DOPC unsaturated bilayers that are usually studied in the fluid state. Phosphatidylcholine spin-labeled at C-5 and C-16 carbon atom positions along the acyl chain were used and the temperature varied over the range 77-270 K. Segmental chain librational oscillations of small amplitude and with correlation time in the subnanosecond to nanosecond range are found in both membranes. The mean-square angular amplitude, α2, of librations increases with temperature, is larger close to the bilayer midplane than close to the first acyl chain segments, and is larger in diunsaturated than in monounsaturated bilayers. In the inner hydrocarbon region of both lipid matrices, α2 increases first slowly and linearly with temperature and then more rapidly, and a dynamical transition is detected in the range 190-210 K. Compared to dipalmitoylphosphatidylcholine bilayers of fully saturated symmetric chain lipids, the presence of double bonds in the acyl chain enhances the intensity of librational motion which is characterized by larger angular variations at the terminal methyl ends. These findings highlight biophysical properties of unsaturated bilayers in the frozen state, including a detailed characterization of segmental chain dynamics and the evidence of a dynamical transition that appears to be a generic feature in hydrated macromolecular systems. These results can also be relevant in regulating membrane physical properties and function at higher physiological temperatures.

Localized and Collective Motions in HET-s(218-289) Fibrils from Combined NMR Relaxation and MD Simulation

Nuclear magnetic resonance (NMR) relaxation data and molecular dynamics (MD) simulations are combined to characterize the dynamics of the fungal prion HET-s(218-289) in its amyloid form. NMR data is analyzed with the dynamics detector method, which yields timescale-specific information. An analogous analysis is performed on MD trajectories. Because specific MD predictions can be verified as agreeing with the NMR data, MD was used for further interpretation of NMR results: for the different timescales, cross-correlation coefficients were derived to quantify the correlation of the motion between different residues. Short timescales are the result of very local motions, while longer timescales are found for longer-range correlated motion. Similar trends on ns- and μs-timescales suggest that μs motion in fibrils is the result of motion correlated over many fibril layers.

Control of Solvent Dynamics around the B12-Dependent Ethanolamine Ammonia-Lyase Enzyme in Frozen Aqueous Solution by Using Dimethyl Sulfoxide Modulation of Mesodomain Volume

The temperature-dependent structure and dynamics of two concentric solvent phases, the protein-associated domain (PAD) and the mesodomain, that surround the ethanolamine ammonia-lyase (EAL) protein from Salmonella typhimurium in frozen polycrystalline aqueous solution are addressed by using electron paramagnetic resonance spectroscopy of the paramagnetic nitroxide spin probe, TEMPOL, over the temperature ( T) range 190-265 K. Dimethyl sulfoxide (DMSO), added at 0.5, 2.0, and 4.0% v/v and present at the maximum freeze concentration at T ≤ 245 K, varies the volume of the interstitial aqueous DMSO mesodomain ( V_meso) relative to a fixed PAD volume ( V_PAD). The increase in V_meso/ V_PAD from 0.8 to 6.0 is quantified by the partitioning of TEMPOL between the two phases. As V_meso/ V_PAD is increased, the Arrhenius parameters for activated TEMPOL rotational motion in the mesodomain remain uniform, whereas the parameters for TEMPOL in the PAD show a progressive transformation toward the mesodomain values (higher mobility). An order-disorder transition (ODT) in the PAD is detected by the exclusion of TEMPOL from the PAD into the mesodomain. The ODT T value is systematically lowered by increased V_meso/ V_PAD (from 215 to 200 K), and PAD ordering kinks the mesodomain Arrhenius dependence. Thus there is reciprocity in PAD-mesodomain solvent coupling. The results are interpreted as a dominant influence of ice-boundary confinement on the PAD solvent structure and dynamics, which is transmitted through the mesodomain and which decreases with mesodomain volume at increased added DMSO. The systematic tuning of PAD and mesodomain solvent dynamics by the variation of added DMSO is an incisive approach for the resolution of contributions of protein-solvent dynamical coupling to EAL catalysis.

Aromatic Ring Dynamics, Thermal Activation, and Transient Conformations of a 468 kDa Enzyme by Specific 1H-13C Labeling and Fast Magic-Angle Spinning NMR

Aromatic residues are located at structurally important sites of many proteins. Probing their interactions and dynamics can provide important functional insight but is challenging in large proteins. Here, we introduce approaches to characterize the dynamics of phenylalanine residues using ¹H-detected fast magic-angle spinning (MAS) NMR combined with a tailored isotope-labeling scheme. Our approach yields isolated two-spin systems that are ideally suited for artifact-free dynamics measurements, and allows probing motions effectively without molecular weight limitations. The application to the TET2 enzyme assembly of ∼0.5 MDa size, the currently largest protein assigned by MAS NMR, provides insights into motions occurring on a wide range of time scales (picoseconds to milliseconds). We quantitatively probe ring-flip motions and show the temperature dependence by MAS NMR measurements down to 100 K. Interestingly, favorable line widths are observed down to 100 K, with potential implications for DNP NMR. Furthermore, we report the first ¹³C R_1ρ MAS NMR relaxation-dispersion measurements and detect structural excursions occurring on a microsecond time scale in the entry pore to the catalytic chamber and at a trimer interface that was proposed as the exit pore. We show that the labeling scheme with deuteration at ca. 50 kHz MAS provides superior resolution compared to 100 kHz MAS experiments with protonated, uniformly ¹³C-labeled samples.

Resolution enhancement and proton proximity probed by 3D TQ/DQ/SQ proton NMR spectroscopy under ultrafast magic-angle-spinning beyond 70 kHz

Proton nuclear magnetic resonance (NMR) in solid state has gained significant attention in recent years due to the remarkable resolution and sensitivity enhancement afforded by ultrafast magic-angle-spinning (MAS). In spite of the substantial suppression of ¹H-¹H dipolar couplings, the proton spectral resolution is still poor compared to that of ¹³C or ¹⁵N NMR, rendering it challenging for the structural and conformational analysis of complex chemicals or biological solids. Herein, by utilizing the benefits of double-quantum (DQ) and triple-quantum (TQ) coherences, we propose a 3D single-channel pulse sequence that correlates proton triple-quantum/double-quantum/single-quantum (TQ/DQ/SQ) chemical shifts. In addition to the two-spin proximity information, this 3D TQ/DQ/SQ pulse sequence enables more reliable extraction of three-spin proximity information compared to the regular 2D TQ/SQ correlation experiment, which could aid in revealing the proton network in solids. Furthermore, the TQ/DQ slice taken at a specific SQ chemical shift only reveals the local correlations to the corresponding SQ chemical shift, and thus it enables accurate assignments of the proton peaks along the TQ and DQ dimensions and simplifies the interpretation of proton spectra especially for dense proton networks. The high performance of this 3D pulse sequence is well demonstrated on small compounds, L-alanine and a tripeptide, N-formyl-L-methionyl-L-leucyl-L-phenylalanine (MLF). We expect that this new methodology can inspire the development of multidimensional solid-state NMR pulse sequences using the merits of TQ and DQ coherences and enable high-throughput investigations of complex solids using abundant protons.

Dynamic Nuclear Polarization Magic-Angle Spinning Nuclear Magnetic Resonance Combined with Molecular Dynamics Simulations Permits Detection of Order and Disorder in Viral Assemblies

We report dynamic nuclear polarization (DNP)-enhanced magic-angle spinning (MAS) NMR spectroscopy in viral capsids from HIV-1 and bacteriophage AP205. Viruses regulate their life cycles and infectivity through modulation of their structures and dynamics. While static structures of capsids from several viruses are now accessible with near-atomic-level resolution, atomic-level understanding of functionally important motions in assembled capsids is lacking. We observed up to 64-fold signal enhancements by DNP, which permitted in-depth analysis of these assemblies. For the HIV-1 CA assemblies, a remarkably high spectral resolution in the 3D and 2D heteronuclear data sets permitted the assignment of a significant fraction of backbone and side-chain resonances. Using an integrated DNP MAS NMR and molecular dynamics (MD) simulation approach, the conformational space sampled by the assembled capsid at cryogenic temperatures was mapped. Qualitatively, a remarkable agreement was observed for the experimental ¹³C/¹⁵N chemical shift distributions and those calculated from substructures along the MD trajectory. Residues that are mobile at physiological temperatures are frozen out in multiple conformers at cryogenic conditions, resulting in broad experimental and calculated chemical shift distributions. Overall, our results suggest that DNP MAS NMR measurements in combination with MD simulations facilitate a thorough understanding of the dynamic signatures of viral capsids.

Solvent-dependent segmental dynamics in intrinsically disordered proteins

Protein and water dynamics have a synergistic relationship, which is particularly important for intrinsically disordered proteins (IDPs), although the details of this coupling remain poorly understood. Here, we combine temperature-dependent molecular dynamics simulations using different water models with extensive nuclear magnetic resonance (NMR) relaxation to examine the importance of distinct modes of solvent and solute motion for the accurate reproduction of site-specific dynamics in IDPs. We find that water dynamics play a key role in motional processes internal to "segments" of IDPs, stretches of primary sequence that share dynamic properties and behave as discrete dynamic units. We identify a relationship between the time scales of intrasegment dynamics and the lifetime of hydrogen bonds in bulk water. Correct description of these motions is essential for accurate reproduction of protein relaxation. Our findings open important perspectives for understanding the role of hydration water on the behavior and function of IDPs in solution.

Solid-state NMR reveals a comprehensive view of the dynamics of the flexible, disordered N-terminal domain of amyloid-β fibrils

Amyloid fibril deposits observed in Alzheimer's disease comprise amyloid-β (Aβ) protein possessing a structured hydrophobic core and a disordered N-terminal domain (residues 1-16). The internal flexibility of the disordered domain is likely essential for Aβ aggregation. Here, we used ²H static solid-state NMR methods to probe the dynamics of selected side chains of the N-terminal domain of Aβ_1-40 fibrils. Line shape and relaxation data suggested a two-state model in which the domain's free state undergoes a diffusive motion that is quenched in the bound state, likely because of transient interactions with the structured C-terminal domain. At 37 °C, we observed freezing of the dynamics progressively along the Aβ sequence, with the fraction of the bound state increasing and the rate of diffusion decreasing. We also found that without solvation, the diffusive motion is quenched. The solvent acted as a plasticizer reminiscent of its role in the onset of global dynamics in globular proteins. As the temperature was lowered, the fraction of the bound state exhibited sigmoidal behavior. The midpoint of the freezing curve coincided with the bulk solvent freezing for the N-terminal residues and increased further along the sequence. Using ²H R _1ρ measurements, we determined the conformational exchange rate constant between the free and bound states under physiological conditions. Zinc-induced aggregation leads to the enhancement of the dynamics, manifested by the faster conformational exchange, faster diffusion, and lower freezing-curve midpoints.

Conformational Excitation and Nonequilibrium Transition Facilitate Enzymatic Reactions: Application to Pin1 Peptidyl-Prolyl Isomerase

Conformational flexibility of protein is essential for enzyme catalysis. Yet, how protein's conformational rearrangements and dynamics contribute to catalysis remains highly controversial. To unravel protein's role in catalysis, it is inevitable to understand the static and dynamic mechanisms simultaneously. To this end, here the Pin1-catalyzed isomerization reaction is studied from the two perspectives. The static view indicates that the hydrogen bonds involving Pin1 rearrange in a tightly coupled manner with isomerization. In sharp contrast, the isomerization dynamics are found to be very rapid; protein's slow conformational rearrangements thus cannot occur simultaneously with isomerization, and the reaction proceeds in a nonequilibrium manner. The distinctive protein conformations necessary to stabilize the transition state are prepared a priori, i.e., as conformational excited states. The present result suggests that enzymatic reaction is not a simple thermal activation from equilibrium directly to the transition state, thus adding a novel perspective to Pauling's view.

Recent advances in solid-state nuclear magnetic resonance spectroscopy of exotic nuclei

We present a review of recent advances in solid-state nuclear magnetic resonance (SSNMR) studies of exotic nuclei. Exotic nuclei may be spin-1/2 or quadrupolar, and typically have low gyromagnetic ratios, low natural abundances, large quadrupole moments (when I > 1/2), or some combination of these properties, generally resulting in low receptivities and/or prohibitively broad line widths. Some nuclides are little studied for other reasons, also rendering them somewhat exotic. We first discuss some of the recent progress in pulse sequences and hardware development which continues to enable researchers to study new kinds of materials as well as previously unfeasible nuclei. This is followed by a survey of applications to a wide range of exotic nuclei (including e.g., ⁹Be, ²⁵Mg, ³³S, ³⁹K, ⁴³Ca, ^47/49Ti, ⁵³Cr, ⁵⁹Co, ⁶¹Ni, ⁶⁷Zn, ⁷³Ge, ⁷⁵As, ⁸⁷Sr, ¹¹⁵In, ¹¹⁹Sn, ^121/123Sb, ^135/137Ba, ^185/187Re, ²⁰⁹Bi), most of them quadrupolar. The scope of the review is the past ten years, i.e., 2007-2017.

Towards a Stochastic Paradigm: From Fuzzy Ensembles to Cellular Functions

The deterministic sequence → structure → function relationship is not applicable to describe how proteins dynamically adapt to different cellular conditions. A stochastic model is required to capture functional promiscuity, redundant sequence motifs, dynamic interactions, or conformational heterogeneity, which facilitate the decision-making in regulatory processes, ranging from enzymes to membraneless cellular compartments. The fuzzy set theory offers a quantitative framework to address these problems. The fuzzy formalism allows the simultaneous involvement of proteins in multiple activities, the degree of which is given by the corresponding memberships. Adaptation is described via a fuzzy inference system, which relates heterogeneous conformational ensembles to different biological activities. Sequence redundancies (e.g., tandem motifs) can also be treated by fuzzy sets to characterize structural transitions affecting the heterogeneous interaction patterns (e.g., pathological fibrillization of stress granules). The proposed framework can provide quantitative protein models, under stochastic cellular conditions.

Molecular Dynamics Simulations of Macromolecular Crystals

The structures of biological macromolecules would not be known to their present extent without X-ray crystallography. Most simulations of globular proteins in solution begin by surrounding the crystal structure of the monomer in a bath of water molecules, but the standard simulation employing periodic boundary conditions is already close to a crystal lattice environment. With simple protocols, the same software and molecular models can perform simulations of the crystal lattice, including all asymmetric units and solvent to fill the box. Throughout the history of molecular dynamics, studies of crystal lattices have served to investigate the quality of the underlying force fields, correlate the simulated ensembles to experimental structure factors, and extrapolate the behavior in lattices to behavior in solution. Powerful new computers are enabling molecular simulations with greater realism and statistical convergence. Meanwhile, the advent of exciting new methods in crystallography, including femtosecond free-electron lasers and image reconstruction for time-resolved crystallography on slurries of small crystals, is expanding the range of structures accessible to X-ray diffraction. We review past fusions of simulations and crystallography, then look ahead to the ways that simulations of crystal structures will enhance structural biology in the future.

Refocusing CSA during magic angle spinning rotating-frame relaxation experiments

We examine coherent evolution of spin-locked magnetization during magic-angle spinning (MAS), in the context of relaxation experiments designed to probe chemical exchange (rotating-frame relaxation (R_1ρ)). Coherent evolution is expected in MAS based rotating-frame relaxation decay experiments if matching conditions are met (such as, ω₁ = nω_r) and if the chemical shielding anisotropy (CSA) is substantial. We show here using numerical simulations and experiments that even when such matching requirements are avoided (e.g., ω₁ < 0.5ω_r, ∼1.5ω_r, >2.5ω_r), coherent evolution of spin-locked magnetization with large CSA is still considerable. The coherent evolution has important consequences on the analysis of relaxation decay and the ability to extract accurate information of interest about dynamics. We present a pulse sequence that employs rotary echoes and refocuses CSA contributions, allowing for more sensitive measurement of rotating-frame relaxation with less interference from coherent evolution. In practice, the proposed pulse sequence, REfocused CSA Rotating-frame Relaxation (RECRR) is robust to carrier frequency offset, B₁-field inhomogeneity, and slight miscalibrations of the refocusing pulses.

Direct Experimental Characterization of Contributions from Self-Motion of Hydrogen and from Interatomic Motion of Heavy Atoms to Protein Anharmonicity

One fundamental challenge in biophysics is to understand the connection between protein dynamics and its function. Part of the difficulty arises from the fact that proteins often present local atomic motions and collective dynamics on the same time scales, and challenge the experimental identification and quantification of different dynamic modes. Here, by taking lyophilized proteins as the example, we combined deuteration technique and neutron scattering to separate and characterize the self-motion of hydrogen and the collective interatomic motion of heavy atoms (C, O, N) in proteins on the pico-to-nanosecond time scales. We found that hydrogen atoms present an instrument-resolution-dependent onset for anharmonic motions, which can be ascribed to the thermal activation of local side-group motions. However, the protein heavy atoms exhibit an instrument-resolution-independent anharmonicity around 200 K, which results from unfreezing of the relaxation of the protein structures on the laboratory equilibrium time (100-1000 s), softening of the entire bio-macromolecules.

Lipid- and Cholesterol-Mediated Time-Scale-Specific Modulation of the Outer Membrane Protein X Dynamics in Lipid Bilayers

Membrane protein function fundamentally depends on lipid-bilayer fluidity and the composition of the biological membrane. Although dynamic interdependencies between membrane proteins and the surrounding lipids are suspected, a detailed description is still missing. To uncover lipid-modulated membrane protein backbone dynamics, time-scale-specific NMR relaxation experiments with residue-resolution were recorded. The data revealed that lipid order, modified either biochemically or biophysically, changes the dynamics of the immersed membrane protein in a specific and time-scale-dependent manner. A temperature-dependent dynamics analysis furthermore suggests a direct coupling between lipid and protein dynamics in the picosecond-nanosecond, microsecond, and millisecond time scales, caused by the lipid's trans-gauche isomerization, the segmental and rotational motion of lipids, and the fluidity of the lipid phase, respectively. These observations provide evidence of a direct modulatory capability of the membrane to regulate protein function through lipid dynamics ranging from picoseconds to milliseconds.

Protein-Protein Interfaces Probed by Methyl Labeling and Proton-Detected Solid-State NMR Spectroscopy

Proton detection and fast magic-angle spinning have advanced biological solid-state NMR, allowing for the backbone assignment of complex protein assemblies with high sensitivity and resolution. However, so far no method has been proposed to detect intermolecular interfaces in these assemblies by proton detection. Herein, we introduce a concept based on methyl labeling that allows for the assignment of these moieties and for the study of protein-protein interfaces at atomic resolution.

Rotational Dynamics of Proteins from Spin Relaxation Times and Molecular Dynamics Simulations

Conformational fluctuations and rotational tumbling of proteins can be experimentally accessed with nuclear spin relaxation experiments. However, interpretation of molecular dynamics from the experimental data is often complicated, especially for molecules with anisotropic shape. Here, we apply classical molecular dynamics simulations to interpret the conformational fluctuations and rotational tumbling of proteins with arbitrarily anisotropic shape. The direct calculation of spin relaxation times from simulation data did not reproduce the experimental data. This was successfully corrected by scaling the overall rotational diffusion coefficients around the protein inertia axes with a constant factor. The achieved good agreement with experiments allowed the interpretation of the internal and overall dynamics of proteins with significantly anisotropic shape. The overall rotational diffusion was found to be Brownian, having only a short subdiffusive region below 0.12 ns. The presented methodology can be applied to interpret rotational dynamics and conformation fluctuations of proteins with arbitrary anisotropic shape. However, a water model with more realistic dynamical properties is probably required for intrinsically disordered proteins.

Slow Dynamics around a Protein and Its Coupling to Solvent

Solvent is essential for protein dynamics and function, but its role in regulating the dynamics remains debated. Here, we employ saturation transfer electron spin resonance (ST-ESR) to explore the issue and characterize the dynamics on a longer (from μs to s) time scale than has been extensively studied. We first demonstrate the reliability of ST-ESR by showing that the dynamical changeovers revealed in the spectra agree to liquid-liquid transition (LLT) in the state diagram of the glycerol/water system. Then, we utilize ST-ESR with four different probes to systematically map out the variation in local (site-specific) dynamics around a protein surface at subfreezing temperatures (180-240 K) in 10 mol % glycerol/water mixtures. At highly exposed sites, protein and solvent dynamics are coupled, whereas they deviate from each other when temperature is greater than LLT temperature (∼190 K) of the solvent. At less exposed sites, protein however exhibits a dynamic, which is distinct from the bulk solvent, throughout the temperature range studied. Dominant dynamic components are thus revealed, showing that (from low to high temperatures) the overall structural fluctuation, rotamer dynamics, and internal side-chain dynamics, in turn, dominate the temperature dependence of spin-label motions. The structural fluctuation component is relatively slow, collective, and independent of protein structural segments, which is thus inferred to a fundamental dynamic component intrinsic to protein. This study corroborates that bulk solvent plasticizes protein and facilitates rather than slaves protein dynamics.

Effect of glycerol on photobleaching of cytochrome Raman lines in frozen yeast cells

We applied a Raman spectroscopy approach to investigate the effect of a cryoprotectant on the redox state of cytochromes on freezing yeast cells. The redox activity of cytochromes was studied using time-resolved photobleaching of the resonance Raman lines. It is found that ice formation causes a drastic change in the redox state of cytochromes in cells frozen without cryoprotectant, whereas in the presence of glycerol the effects of ice formation are more gradual. The photobleaching rate of cells frozen in glycerol solution shows a gradual slowing with temperature decrease and an abrupt slowdown below - 48 °C. This abrupt decrease was interpreted as originating from changes in protein conformational dynamics. Our findings provide important new insights into the transition from active to inactive cytochrome states as cells undergo freezing in the presence and absence of cryoprotectant.

Hidden motions and motion-induced invisibility: Dynamics-based spectral editing in solid-state NMR

Solid-state nuclear magnetic resonance (ssNMR) spectroscopy enables the structural characterization of a diverse array of biological assemblies that include amyloid fibrils, non-amyloid aggregates, membrane-associated proteins and viral capsids. Such biological samples feature functionally relevant molecular dynamics, which often affect different parts of the sample in different ways. Solid-state NMR experiments' sensitivity to dynamics represents a double-edged sword. On the one hand, it offers a chance to measure dynamics in great detail. On the other hand, certain types of motion lead to signal loss and experimental inefficiencies that at first glance interfere with the application of ssNMR to overly dynamic proteins. Dynamics-based spectral editing (DYSE) ssNMR methods leverage motion-dependent signal losses to simplify spectra and enable the study of sub-structures with particular motional properties.

DNP enhanced NMR with flip-back recovery

DNP methods can provide significant sensitivity enhancements in magic angle spinning solid-state NMR, but in systems with long polarization build up times long recycling periods are required to optimize sensitivity. We show how the sensitivity of such experiments can be improved by the classic flip-back method to recover bulk proton magnetization following continuous wave proton heteronuclear decoupling. Experiments were performed on formulations with characteristic build-up times spanning two orders of magnitude: a bulk BDPA radical doped o-terphenyl glass and microcrystalline samples of theophylline, l-histidine monohydrochloride monohydrate, and salicylic acid impregnated by incipient wetness. For these systems, addition of flip-back is simple, improves the sensitivity beyond that provided by modern heteronuclear decoupling methods such as SPINAL-64, and provides optimal sensitivity at shorter recycle delays. We show how to acquire DNP enhanced 2D refocused CP-INADEQUATE spectra with flip-back recovery, and demonstrate that the flip-back recovery method is particularly useful in rapid recycling regimes. We also report Overhauser effect DNP enhancements of over 70 at 592.6 GHz/900 MHz.

New applications of solid-state NMR in structural biology

Various recent developments in solid-state nuclear magnetic resonance (ssNMR) spectroscopy have enabled an array of new insights regarding the structure, dynamics, and interactions of biomolecules. In the ever more integrated world of structural biology, ssNMR studies provide structural and dynamic information that is complementary to the data accessible by other means. ssNMR enables the study of samples lacking a crystalline lattice, featuring static as well as dynamic disorder, and does so independent of higher-order symmetry. The present study surveys recent applications of biomolecular ssNMR and examines how this technique is increasingly integrated with other structural biology techniques, such as (cryo) electron microscopy, solution-state NMR, and X-ray crystallography. Traditional ssNMR targets include lipid bilayer membranes and membrane proteins in a lipid bilayer environment. Another classic application has been in the area of protein misfolding and aggregation disorders, where ssNMR has provided essential structural data on oligomers and amyloid fibril aggregates. More recently, the application of ssNMR has expanded to a growing array of biological assemblies, ranging from non-amyloid protein aggregates, protein–protein complexes, viral capsids, and many others. Across these areas, multidimensional magic angle spinning (MAS) ssNMR has, in the last decade, revealed three-dimensional structures, including many that had been inaccessible by other structural biology techniques. Equally important insights in structural and molecular biology derive from the ability of MAS ssNMR to probe information beyond comprehensive protein structures, such as dynamics, solvent exposure, protein–protein interfaces, and substrate–enzyme interactions.

Mesodomain and Protein-Associated Solvent Phases with Temperature-Tunable (200-265 K) Dynamics Surround Ethanolamine Ammonia-Lyase in Globally Polycrystalline Aqueous Solution Containing Dimethyl Sulfoxide

Electron paramagnetic resonance spectroscopy of the spin probe, TEMPOL, is used to resolve solvent phases that surround the ethanolamine ammonia-lyase (EAL) protein from Salmonella typhimurium at low temperature (T) in frozen, globally polycrystalline aqueous solution and to report on the T dependence of their detectably rigid and fluid states. EAL plays a role in human gut microbiome-based disease conditions, and physicochemical studies provide insight into protein structure and mechanism, toward potential therapeutics. Temperature dependences of the rotational correlation times (τ_c; detection range, 10^-11 ≤ τ_c ≤ 10^-7 s) and the corresponding weights of TEMPOL tumbling components from 200 to 265 K in the presence of EAL are measured in two frozen systems: (1) water-only and (2) 1% v/v dimethyl sulfoxide (DMSO). In the water-only condition, a protein-vicinal solvent component detectably fluidizes at 230 K and melts the surrounding ice-crystalline region with increasing T, creating a bounded, relatively high-viscosity aqueous solvent domain, up to 265 K. In the EAL, 1% v/v DMSO condition, two distinct concentric solvent phases are resolved around EAL: protein-associated domain (PAD) and mesodomain. The DMSO aqueous mesodomain fluidizes at 200 K, followed by PAD fluidization at 210 K. The interphase dynamical coupling is consistent with the spatial arrangement and significant contact areas of the phases, indicated by the experimentally determined mean volume ratio, V(mesodomain)/V(PAD)/V(protein) = 0.5:0.3:1.0. The results provide a rationale for native chemical reactions of EAL at T < 250 K and an advance toward precise control of solvent dynamics as a tunable parameter for quantifying the coupling between solvent and protein fluctuations and chemical reaction steps in EAL and other enzymes.

Atomic resolution conformational dynamics of intrinsically disordered proteins from NMR spin relaxation

Nuclear magnetic resonance (NMR) spectroscopy is one of the most powerful experimental approaches for investigating the conformational behaviour of intrinsically disordered proteins (IDPs). IDPs represent a significant fraction of all proteomes, and, despite their importance for understanding fundamental biological processes, the molecular basis of their activity still remains largely unknown. The functional mechanisms exploited by IDPs in their interactions with other biomolecules are defined by their intrinsic dynamic modes and associated timescales, justifying the considerable interest over recent years in the development of technologies adapted to measure and describe this behaviour. NMR spin relaxation delivers information-rich, site-specific data reporting on conformational fluctuations occurring throughout the molecule. Here we review recent progress in the use of ¹⁵N relaxation to identify local backbone dynamics and long-range chain-like motions in unfolded proteins.

Low-Temperature Dynamics of Chain-Labeled Lipids in Ester- and Ether-Linked Phosphatidylcholine Membranes

Continuous wave electron paramagnetic resonance spectroscopy and two-pulse echo detected spectra of chain-labeled lipids are used to study the dynamics of frozen lipid membranes over the temperature range 77-260 K. Bilayers of ester-linked dihexadecanoylphosphatidylcholine (DPPC) with noninterdigitated chains and ether-linked dihexadecyl phosphatidylcholine (DHPC) with interdigitated chains are considered. Rapid stochastic librations of small angular amplitude are found in both lipid matrices. In noninterdigitated DPPC bilayers, the mean-square angular amplitude, [Formula: see text], of the motion increases with temperature and it is larger close to the chain termini than close to the polar/apolar interface. In contrast, in interdigitated DHPC lamellae, [Formula: see text] is small and temperature and label-position independent at low temperature and increases steeply at high temperature. The rotational correlation time, τ_c, of librations lies in the subnanosecond range for DPPC and in the nanosecond range for DHPC. In all membrane samples, the temperature dependence of [Formula: see text] resembles that of the mean-square atomic displacement revealed by neutron scattering and a dynamical transition is detected in the range 210-240 K. The results highlight the librational oscillations and the glass-like behavior in bilayer and interdigitated lipid membranes.

The dynamics of peptide-water interactions in dialanine: An ultrafast amide I 2D IR and computational spectroscopy study

We present a joint experimental and computational study of the dynamic interactions of dialanine (Ala-Ala) with water, comparing the results of ultrafast 2D IR and infrared transient absorption spectroscopy of its amide I vibration with spectra modeled from molecular dynamics (MD) simulations. The experimental data are analyzed to describe vibrational frequency fluctuations, vibrational energy relaxation, and chemical exchange processes. The origin of these processes in the same underlying fluctuating forces allows a common description in terms of the fluctuations and conformational dynamics of the peptide and associated solvent. By comparing computational spectroscopy from MD simulations with multiple force fields and water models, we describe how the dynamics of water hydrogen bond fluctuations and switching processes act as a source of friction that governs the dephasing and vibrational relaxation, and provide a description of coupled water and peptide motions that give rise to spectroscopic exchange processes.

Dynamics of Hydration Water Plays a Key Role in Determining the Binding Thermodynamics of Protein Complexes

Interfacial waters are considered to play a crucial role in protein–protein interactions, but in what sense and why are they important? Here, using molecular dynamics simulations and statistical thermodynamic analyses, we demonstrate distinctive dynamic characteristics of the interfacial water and investigate their implications for the binding thermodynamics. We identify the presence of extraordinarily slow (~1,000 times slower than in bulk water) hydrogen-bond rearrangements in interfacial water. We rationalize the slow rearrangements by introducing the “trapping” free energies, characterizing how strongly individual hydration waters are captured by the biomolecular surface, whose magnitude is then traced back to the number of water–protein hydrogen bonds and the strong electrostatic field produced at the binding interface. We also discuss the impact of the slow interfacial waters on the binding thermodynamics. We find that, as expected from their slow dynamics, the conventional approach to the water-mediated interaction, which assumes rapid equilibration of the waters’ degrees of freedom, is inadequate. We show instead that an explicit treatment of the extremely slow interfacial waters is critical. Our results shed new light on the role of water in protein–protein interactions, highlighting the need to consider its dynamics to improve our understanding of biomolecular bindings.

Temperature-Independent Nuclear Quantum Effects on the Structure of Water

Nuclear quantum effects (NQEs) have a significant influence on the hydrogen bonds in water and aqueous solutions and have thus been the topic of extensive studies. However, the microscopic origin and the corresponding temperature dependence of NQEs have been elusive and still remain the subject of ongoing discussion. Previous x-ray scattering investigations indicate that NQEs on the structure of water exhibit significant temperature dependence [Phys. Rev. Lett. 94, 047801 (2005)PRLTAO0031-900710.1103/PhysRevLett.94.047801]. Here, by performing wide-angle x-ray scattering of H_{2}O and D_{2}O droplets at temperatures from 275 K down to 240 K, we determine the temperature dependence of NQEs on the structure of water down to the deeply supercooled regime. The data reveal that the magnitude of NQEs on the structure of water is temperature independent, as the structure factor of D_{2}O is similar to H_{2}O if the temperature is shifted by a constant 5 K, valid from ambient conditions to the deeply supercooled regime. Analysis of the accelerated growth of tetrahedral structures in supercooled H_{2}O and D_{2}O also shows similar behavior with a clear 5 K shift. The results indicate a constant compensation between NQEs delocalizing the proton in the librational motion away from the bond and in the OH stretch vibrational modes along the bond. This is consistent with the fact that only the vibrational ground state is populated at ambient and supercooled conditions.

Solvent-Driven Dynamical Crossover in the Phenylalanine Side-Chain from the Hydrophobic Core of Amyloid Fibrils Detected by 2H NMR Relaxation

Aromatic residues are important markers of dynamical changes in proteins' hydrophobic cores. In this work we investigated the dynamics of the F19 side-chain in the core of amyloid fibrils across a wide temperature range of 300 to 140 K. We utilized solid-state 2H NMR relaxation to demonstrate the presence of a solvent-driven dynamical crossover between different motional regimes, often also referred to as the dynamical transition. In particular, the dynamics are dominated by small-angle fluctuations at low temperatures and by π-flips of the aromatic ring at high temperatures. The crossover temperature is more than 43 degrees lower for the hydrated state of the fibrils compared to the dry state, indicating that interactions with water facilitate π-flips. Further, crossover temperatures are shown to be very sensitive to polymorphic states of the fibrils, such as the 2-fold and 3-fold symmetric morphologies of the wild-type protein as well as D23N mutant protofibrils. We speculate that these differences can be attributed, at least partially, to enhanced interactions with water in the 3-fold polymorph, which has been shown to have a water-accessible cavity. Combined with previous studies of methyl group dynamics, the results highlight the presence of multiple dynamics modes in the core of the fibrils, which was originally believed to be quite rigid.

Transmembrane allosteric energetics characterization for strong coupling between proton and potassium ion binding in the KcsA channel

The slow spontaneous inactivation of potassium channels exhibits classic signatures of transmembrane allostery. A variety of data support a model in which the loss of K⁺ ions from the selectivity filter is a major factor in promoting inactivation, which defeats transmission, and is allosterically coupled to protonation of key channel activation residues, more than 30 Å from the K⁺ ion binding site. We show that proton binding at the intracellular pH sensor perturbs the potassium affinity at the extracellular selectivity filter by more than three orders of magnitude for the full-length wild-type KcsA, a pH-gated bacterial channel, in membrane bilayers. Studies of F103 in the hinge of the inner helix suggest an important role for its bulky sidechain in the allosteric mechanism; we show that the energetic strength of coupling of the gates is strongly altered when this residue is mutated to alanine. These results provide quantitative site-specific measurements of allostery in a bilayer environment, and highlight the power of describing ion channel gating through the lens of allosteric coupling.

Bacteriophage Tail-Tube Assembly Studied by Proton-Detected 4D Solid-State NMR

Obtaining unambiguous resonance assignments remains a major bottleneck in solid-state NMR studies of protein structure and dynamics. Particularly for supramolecular assemblies with large subunits (>150 residues), the analysis of crowded spectral data presents a challenge, even if three-dimensional (3D) spectra are used. Here, we present a proton-detected 4D solid-state NMR assignment procedure that is tailored for large assemblies. The key to recording 4D spectra with three indirect carbon or nitrogen dimensions with their inherently large chemical shift dispersion lies in the use of sparse non-uniform sampling (as low as 2 %). As a proof of principle, we acquired 4D (H)COCANH, (H)CACONH, and (H)CBCANH spectra of the 20 kDa bacteriophage tail-tube protein gp17.1 in a total time of two and a half weeks. These spectra were sufficient to obtain complete resonance assignments in a straightforward manner without use of previous solution NMR data.

Solid-State NMR Provides Evidence for Small-Amplitude Slow Domain Motions in a Multispanning Transmembrane α-Helical Protein

Proteins are dynamic entities and populate ensembles of conformations. Transitions between states within a conformational ensemble occur over a broad spectrum of amplitude and time scales, and are often related to biological function. Whereas solid-state NMR (SSNMR) spectroscopy has recently been used to characterize conformational ensembles of proteins in the microcrystalline states, its applications to membrane proteins remain limited. Here we use SSNMR to study conformational dynamics of a seven-helical transmembrane (TM) protein, Anabaena Sensory Rhodopsin (ASR) reconstituted in lipids. We report on site-specific measurements of the 15N longitudinal R1 and rotating frame R1ρ relaxation rates at two fields of 600 and 800 MHz and at two temperatures of 7 and 30 °C. Quantitative analysis of the R1 and R1ρ values and of their field and temperature dependencies provides evidence of motions on at least two time scales. We modeled these motions as fast local motions and slower collective motions of TM helices and of structured loops, and used the simple model-free and extended model-free analyses to fit the data and estimate the amplitudes, time scales and activation energies. Faster picosecond (tens to hundreds of picoseconds) local motions occur throughout the protein and are dominant in the middle portions of the TM helices. In contrast, the amplitudes of the slower collective motions occurring on the nanosecond (tens to hundreds of nanoseconds) time scales, are smaller in the central parts of helices, but increase toward their cytoplasmic sides as well as in the interhelical loops. ASR interacts with a soluble transducer protein on its cytoplasmic surface, and its binding affinity is modulated by light. The larger amplitude of motions on the cytoplasmic side of the TM helices correlates with the ability of ASR to undergo large conformational changes in the process of binding/unbinding the transducer.