Dynamics of Hydrophobic Core Phenylalanine Residues Probed by Solid-State Deuteron NMR

Persistence of Methionine Side Chain Mobility at Low Temperatures in a Nine-Residue Low Complexity Peptide, as Probed by 2 H Solid-State NMR

Methionine side chains are flexible entities which play important roles in defining hydrophobic interfaces. We utilize deuterium static solid-state NMR to assess rotameric inter-conversions and other dynamic modes of the methionine in the context of a nine-residue random-coil peptide (RC9) with the low-complexity sequence GGKGMGFGL. The measurements in the temperature range of 313 to 161 K demonstrate that the rotameric interconversions in the hydrated solid powder state persist to temperatures below 200 K. Removal of solvation significantly reduces the rate of the rotameric motions. We employed 2 H NMR line shape analysis, longitudinal and rotation frame relaxation, and chemical exchange saturation transfer methods and found that the combination of multiple techniques creates a significantly more refined model in comparison with a single technique. Further, we compare the most essential features of the dynamics in RC9 to two different methionine-containing systems, characterized previously. Namely, the M35 of hydrated amyloid-β1-40 in the three-fold symmetric polymorph as well as Fluorenylmethyloxycarbonyl (FMOC)-methionine amino acid with the bulky hydrophobic group. The comparison suggests that the driving force for the enhanced methionine side chain mobility in RC9 is the thermodynamic factor stemming from distributions of rotameric populations, rather than the increase in the rate constant.

Rigidifying of the internal dynamics of amyloid-beta fibrils generated in the presence of synaptic plasma vesicles

We investigated the changes in internal flexibility of amyloid-β1-40 (Aβ) fibrils grown in the presence of rat synaptic plasma vesicles. The fibrils are produced using a modified seeded growth protocol, in which the Aβ concentration is progressively increased at the expense of the decreased lipid to protein ratio. The morphologies of each generation are carefully assessed at several fibrils' growth time points using transmission electron microscopy. The side-chain dynamics in the fibrils is investigated using deuterium solid-state NMR measurements, with techniques spanning line shapes analysis and several NMR relaxation rates measurements. The dynamics is probed in the site-specific fashion in the hydrophobic C-terminal domain and the disordered N-terminal domain. An overall strong rigidifying effect is observed in comparison with the wild-type fibrils generated in the absence of the membranes. In particular, the overall large-scale fluctuations of the N-terminal domain are significantly reduced, and the activation energies of rotameric inter-conversion in methyl-bearing side-chains of the core (L17, L34, M35, V36), as well as the ring-flipping motions of F19 are increased, indicating a restricted core environment. Membrane-induced flexibility changes in Aβ aggregates can be important for the re-alignment of protein aggregates within the membrane, which in turn would act as a disruption pathway of the bilayers' integrity.

Effect of Cross-Seeding of Wild-Type Amyloid-β1-40 Peptides with Post-translationally Modified Fibrils on Internal Dynamics of the Fibrils Using Deuterium Solid-State NMR

Post-translationally modified (PTM) amyloid-β (Aβ) species can play an important role in modulating Alzheimer's disease pathology. These relatively less populated modifications can cross-seed the wild-type Aβ peptides to produce fibrils that retain many structural and functional features of the original PTM variants. We focus on studies of internal flexibility in the cross-seeded Aβ1-40 fibrils originating from seeding with two PTM variants with modifications in the disordered N-terminal domain: ΔE3 truncation and S8-phosphorylation. We employ an array of 2H solid-state NMR techniques, including line shape analysis over a broad temperature range, longitudinal relaxation, and quadrupolar CPMG, to assess the dynamics of the cross-seeded fibrils. The focus is placed on selected side-chain sites in the disordered N-terminal domain (G9 and V12) and hydrophobic core methyl and aromatic groups (L17, L34, M35, V36, and F19). We find that many of the essential features of the dynamics present in the original PTM seeds persist in the cross-seeded fibrils, and several of the characteristic features are even enhanced. This is particularly true for the activation energies of the rotameric motions and large-scale rearrangements of the N-terminal domain. Thus, our results on the dynamics complement prior structural and cell toxicity studies, suggesting that many PTM Aβ species can aggressively cross-seed the wild-type peptide in a manner that propagates the PTM's signature.

NMR Studies of Aromatic Ring Flips to Probe Conformational Fluctuations in Proteins

Aromatic residues form a significant part of the protein core, where they make tight interactions with multiple surrounding side chains. Despite the dense packing of internal side chains, the aromatic rings of phenylalanine and tyrosine residues undergo 180° rotations, or flips, which are mediated by transient and large-scale "breathing" motions that generate sufficient void volume around the aromatic ring. Forty years after the seminal work by Wagner and Wüthrich, NMR studies of aromatic ring flips are now undergoing a renaissance as a powerful means of probing fundamental dynamic properties of proteins. Recent developments of improved NMR methods and isotope labeling schemes have enabled a number of advances in addressing the mechanisms and energetics of aromatic ring flips. The nature of the transition states associated with ring flips can be described by thermodynamic activation parameters, including the activation enthalpy, activation entropy, activation volume, and also the isothermal volume compressibility of activation. Consequently, it is of great interest to study how ring flip rate constants and activation parameters might vary with protein structure and external conditions like temperature and pressure. The field is beginning to gather such data for aromatic residues in a variety of environments, ranging from surface exposed to buried. In the future, the combination of solution and solid-state NMR spectroscopy together with molecular dynamics simulations and other computational approaches is likely to provide detailed information about the coupled dynamics of aromatic rings and neighboring residues. In this Perspective, we highlight recent developments and provide an outlook toward the future.

Aromatic ring flips in differently packed ubiquitin protein crystals from MAS NMR and MD

Probing the dynamics of aromatic side chains provides important insights into the behavior of a protein because flips of aromatic rings in a protein's hydrophobic core report on breathing motion involving a large part of the protein. Inherently invisible to crystallography, aromatic motions have been primarily studied by solution NMR. The question how packing of proteins in crystals affects ring flips has, thus, remained largely unexplored. Here we apply magic-angle spinning NMR, advanced phenylalanine ¹H-¹³C/²H isotope labeling and MD simulation to a protein in three different crystal packing environments to shed light onto possible impact of packing on ring flips. The flips of the two Phe residues in ubiquitin, both surface exposed, appear remarkably conserved in the different crystal forms, even though the intermolecular packing is quite different: Phe4 flips on a ca. 10-20 ns time scale, and Phe45 are broadened in all crystals, presumably due to µs motion. Our findings suggest that intramolecular influences are more important for ring flips than intermolecular (packing) effects.

Characterizing Fast Conformational Exchange of Aromatic Rings Using Residual Dipolar Couplings: Distinguishing Jumplike Flips from Other Exchange Mechanisms

Aromatic ring flips are a hallmark of protein dynamics. They are experimentally studied by NMR spectroscopy, where recent advances have led to improved characterization across a wide range of time scales. Results on different proteins have been interpreted as continuous diffusive ring rotations or jumplike flips, leading to diverging views of the protein interior as being fluidlike or solidlike, respectively. It is challenging to distinguish between these mechanisms and other types of conformational exchange because chemical-shift-mediated line broadening provides only conclusive evidence for ring flips only if the system can be moved from the slow- to intermediate/fast-exchange regime. Moreover, whenever the chemical shift difference between the two symmetry-related sites is close to zero, it is not generally possible to determine the exchange time scale. Here we resolve these issues by measuring residual dipolar coupling (RDC)-mediated exchange contributions using NMR relaxation dispersion experiments on proteins dissolved in dilute liquid crystalline media. Excellent agreement is found between the experimental difference in RDC between the two symmetry-related sites and the value calculated from high-resolution X-ray structures, demonstrating that dynamics measured for F52 in the B1 domain of protein G reports on distinct, jumplike flips rather than other types of conformational exchange.

Comparative Hydrophobic Core Dynamics Between Wild-Type Amyloid-β Fibrils, Glutamate-3 Truncation, and Serine-8 Phosphorylation

Post-translational modifications (PTMs) of amyloid-β (Aβ) species are implicated in the modulation of overall toxicities and aggregation propensities. We investigated the internal dynamics in the hydrophobic core of the truncated ΔE3 mutant fibrils of Aβ_1-40 and compared them with prior and new data for wild-type fibrils as well as with phosphorylated S8 fibrils. Deuteron static solid-state NMR techniques, spanning line-shape analysis, longitudinal relaxation, and chemical exchange saturation transfer methods, were employed to assess the rotameric jumps of several methyl-bearing and aromatic groups in the core of the fibrils. Taken together, the results indicate the rather significant influence of the PTMs on the hydrophobic core dynamics, which propagates far beyond the local site of the chemical modification. The phosphorylated S8 fibrils display an overall rigidifying of the core based on the higher activation barriers of motions than the wild-type fibrils, whereas the ΔE3 fibrils induce a broader variety of changes, some of which are thermodynamic in nature rather than the kinetic ones.

Dynamics of Serine-8 Side-Chain in Amyloid-β Fibrils and Fluorenylmethyloxycarbonyl Serine Amino Acid, Investigated by Solid-State Deuteron NMR

Serine side-chains are strategic sites of post-translational modifications, and it is important to establish benchmarks of their internal dynamics. In this work, we compare the dynamics of serine side-chains in several biologically important systems: serine-8 in the disordered domain of Aβ_1-40 fibrils in the hydrated and dry states and fluorenylmethyloxycarbonyl (Fmoc) serine with the bulky group that mimics the hydrophobicity of the fibril contacts yet lacks the complexity of the protein system. Using deuterium solid-state NMR static line shape and longitudinal relaxation techniques in the 310 to 180 K temperature range, we compare the main features of the dynamics in these systems. The main motional modes in the fibrils are large-scale fluctuations in the hydrated state of the fibrils as well as local motions such as 3-site jumps of the C^β deuterons at high temperatures and small-angle fluctuations of the C^α-C^β axis at low temperatures. In the hydrated fibrils, two distinct states are present with vastly different extents of large-scale diffusive motions and 3-site-jump rate constants. The hydrated state at the physiological conditions is dominated by the "free" state undergoing large-scale diffusive motions and very fast local 3-site jumps, while in the "bound" state, these large-scale motions are quenched due to transient inter- and intramolecular interactions. Additionally, in the bound state, the 3-site-jump motions are orders of magnitude slower. Details of the dynamics in the serine side-chain are dependent on fine structural features and hydration levels of the systems.

Combining Molecular and Spin Dynamics Simulations with Solid-State NMR: A Case Study of Amphiphilic Lysine-Leucine Repeat Peptide Aggregates

Interpreting dynamics in solid-state molecular systems requires characterization of the potentially heterogeneous environmental contexts of molecules. In particular, the analysis of solid-state nuclear magnetic resonance (ssNMR) data to elucidate molecular dynamics (MD) involves modeling the restriction to overall tumbling by neighbors, as well as the concentrations of water and buffer. In this exploration of the factors that influence motion, we utilize atomistic MD trajectories of peptide aggregates with varying hydration to mimic an amorphous solid-state environment and predict ssNMR relaxation rates. We also account for spin diffusion in multiply spin-labeled (up to 19 nuclei) residues, with several models of dipolar-coupling networks. The framework serves as a general approach to determine essential spin couplings affecting relaxation, benchmark MD force fields, and reveal the hydration dependence of dynamics in a crowded environment. We demonstrate the methodology on a previously characterized amphiphilic 14-residue lysine-leucine repeat peptide, LKα14 (Ac-LKKLLKLLKKLLKL-c), which has an α-helical secondary structure and putatively forms leucine-burying tetramers in the solid state. We measure the R₁ relaxation rates of uniformly ¹³C-labeled and site-specific ²H-labeled leucines in the hydrophobic core of LKα14 at multiple hydration levels. Studies of 9 and 18 tetramer bundles reveal the following: (a) for the incoherent component of ¹³C relaxation, the nearest-neighbor spin interactions dominate, while the ¹H-¹H interactions have minimal impact; (b) the AMBER ff14SB dihedral barriers for the leucine C_γ-C_δ bond ("methyl rotation barriers") must be lowered by a factor of 0.7 to better match the ²H data; (c) proton-driven spin diffusion explains some of the discrepancy between experimental and simulated rates for the C_β and C_α nuclei; and (d) ¹³C relaxation rates are mostly underestimated in the MD simulations at all hydrations, and the discrepancies identify likely motions missing in the 50 ns MD trajectories.

Local ordering and dynamics in anisotropic media by magnetic resonance: from liquid crystals to proteins

Magnetic resonance methods have been used extensively for over 50 years to elucidate molecular structure and dynamics of liquid crystals (LCs), providing information quite unique in its rigour and extent. The ESR- or NMR-active probe is often a solute molecule reporting on characteristics associated with the surrounding (LC) medium, which exerts the spatial restrictions on the probe. The theoretical approaches developed for LCs are applicable to anisotropic media in general. Of particular interest is the interior space of a globular protein labelled, e.g. with a nitroxide moiety or a ¹⁵N-¹H bond. The ESR or NMR label plays the role of the probe and the internal protein surroundings the role of the anisotropic medium. A general feature of the restricted motions is the local ordering, i.e. the nature, magnitude and symmetry of the spatial restraints exerted at the site of the moving probe. This property is the main theme of the present review article. We outline its treatment in our work from both the theoretical and the experimental points of view, highlighting the new physical insights gained. Our illustrations include studies on thermotropic (nematic and smectic) and lyotropic liquid crystals formed by phospholipids, in addition to studies of proteins.

Phenyl-Ring Dynamics in Amyloid Fibrils and Proteins: The Microscopic-Order-Macroscopic-Disorder Perspective

We have developed the microscopic-order-macroscopic-disorder (MOMD) approach for studying internal mobility in polycrystalline proteins with ²H lineshape analysis. The motion itself is expressed by a diffusion tensor, R, the local spatial restraints by a potential, u, and the "local geometry" by the relative orientation of the model-related and nuclear magnetic resonance-related tensors. Here, we apply MOMD to phenyl-ring dynamics in several Αβ₄₀-amyloid-fibrils, and the villin headpiece subdomain (HP36). Because the available data are limited in extent and sensitivity, we adjust u and R in the relevant parameter ranges, fixing the "local geometry" in accordance with standard stereochemistry. This yields a physically well-defined and consistent picture of phenyl-ring dynamics, enabling comparison between different systems. In the temperature range of 278-308 K, u has a strength of (1.7-1.8) kT and a rhombicity of (2.4-2.6) kT, and R has components of 5.0 × 10² ≤ R_⊥ ≤ 2.0 × 10³ s^-1 and 6.3 × 10⁵ ≤ R_∥ ≤ 2.0 × 10⁶ s^-1. At 278 K, fibril hydration increases the axiality of both u and R; HP36 hydration has a similar effect at 295 K, reducing R_⊥ considerably. The D23N mutation slows down the motion of the probe; Aβ₄₀ polymorphism affects both this motion and the related local potential. The present study identifies the impact of various factors on phenyl-ring mobility in amyloid fibrils and globular proteins; the difference between the two protein forms is considerable. The distinctive impact of hydration on phenyl-ring motion and previously studied methyl-group motion is also examined. The ²H lineshapes considered here were analyzed previously with various multi-simple-mode (MSM) models, where several simple motional modes are combined. The MOMD and MSM interpretations differ in essence.

Basic experiments in 2H static NMR for the characterization of protein side-chain dynamics

The focus of this review is the basic methodology for applications of static deuteron NMR for studies of dynamics in the side chains of proteins. We review experimental approaches for the measurements of static line shapes and relaxation rates as well as signal enhancement strategies using the multiple echo acquisition scheme. Further, we describe computational strategies for modeling jump and diffusive motions underlying experimental data. Applications are chosen from studies of amyloid fibrils comprising the amyloid-β protein.

Acetone mobility in zeolite cages with new features in the deuteron NMR spectra and relaxation

We studied deuteron NMR spectra and spin-lattice relaxation of deuterated acetone-d₆, adsorbed into zeolites NaX (1.3) and NaY(2.4) at 100% coverage of sodium cations. At temperatures roughly below 160 K the deuterons are localized and their NMR characteristics are determined by CD₃ rotation and rotational oscillations of acetone molecules. In NaX the CD₃ rotation and rotational oscillations about the twofold axis of acetone dominate the spectra below 100 K, while above it oscillations also about other axes become important. In NaY dominant features are related to methyl tunnelling and to a smaller extent to rigid acetones, before the rotational oscillations about twofold axis start to prevail above 40 K. The analysis of the strongly non-exponential magnetization recovery was done by applying the recently introduced method (Ylinen et al., 2015 [12]), improved here to take into account the limited fast recovery at the level crossings, 10% at ω_t=ω₀ and 28% at ω_t=2ω₀. At first the experimental recovery is fitted by three exponentials with adjustable weights and decay rates. Then these quantities are calculated from activation energy distributions and known expressions for the deuteron relaxation rate. In NaY two distinctly separate activation energy distributions were needed, the dominant one being very broad. The use of three distributions, two of them covering practically the same energies as the broad one, lead to a somewhat better agreement with experiment. In general the theoretical results agree with experiment within experimental scatter. As the final result the mean activation energies and widths are obtained for activation energy distributions.

Comparative Dynamics of Methionine Side-Chain in FMOC-Methionine and in Amyloid Fibrils

We compared the dynamics of key methionine methyl groups in the water-accessible hydrophobic cavity of amyloid fibrils and Fluorenylmethyloxycarbonyl-Methionine (FMOC-Met), which renders general hydrophobicity to the environment without the complexity of the protein. Met35 in the hydrated cavity was recently found to undergo a dynamical cross-over from the dominance of methyl rotations at low temperatures to the dominance of diffusive motion of methyl axis at high temperatures. Current results indicate that in FMOC-Met this cross-over is suppressed, similar to what was observed for the dry fibrils, indicating that hydration of the cavity is driving the onset of the dynamical transition.

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.

Static solid-state 2H NMR methods in studies of protein side-chain dynamics

In this review, we discuss the experimental static deuteron NMR techniques and computational approaches most useful for the investigation of side-chain dynamics in protein systems. Focus is placed on the interpretation of line shape and relaxation data within the framework of motional modeling. We consider both jump and diffusion models and apply them to uncover glassy behaviors, conformational exchange and dynamical transitions in proteins. Applications are chosen from globular and membrane proteins, amyloid fibrils, peptide adsorbed on surfaces and proteins specific to connective tissues.

Ubiquitin immobilized on mesoporous MCM41 silica surfaces - Analysis by solid-state NMR with biophysical and surface characterization

Deriving the conformation of adsorbed proteins is important in the assessment of their functional activity when immobilized. This has particularly important bearings on the design of contemporary and new encapsulated enzyme-based drugs, biosensors, and other bioanalytical devices. Solid-state nuclear magnetic resonance (NMR) measurements can expand our molecular view of proteins in this state and of the molecular interactions governing protein immobilization on popular biocompatible surfaces such as silica. Here, the authors study the immobilization of ubiquitin on the mesoporous silica MCM41 by NMR and other techniques. Protein molecules are shown to bind efficiently at pH 5 through electrostatic interactions to individual MCM41 particles, causing their agglutination. The strong attraction of ubiquitin to MCM41 surface is given molecular context through evidence of proximity of basic, carbonyl and polar groups on the protein to groups on the silica surface using NMR measurements. The immobilized protein exhibits broad peaks in two-dimensional 13C dipolar-assisted rotational resonance spectra, an indication of structural multiplicity. At the same time, cross-peaks related to Tyr and Phe sidechains are missing due to motional averaging. Overall, the favorable adsorption of ubiquitin to MCM41 is accompanied by conformational heterogeneity and by a major loss of motional degrees of freedom as inferred from the marked entropy decrease. Nevertheless, local motions of the aromatic rings are retained in the immobilized state.

LK peptide side chain dynamics at interfaces are independent of secondary structure

Real-time observation of the ultrafast motions of leucine side chains within model peptides at the water–air interface with representative folds – α-helix, 3₁₀-helix, β-strand – show that interfacial dynamics are mostly determined by surface interactions.

Fast Motions of Key Methyl Groups in Amyloid-β Fibrils

Amyloid-β (Aβ) peptide is the major component of plaques found in Alzheimer's disease patients. Using solid-state 2H NMR relaxation performed on selectively deuterated methyl groups, we probed the dynamics in the threefold symmetric and twofold symmetric polymorphs of native Aβ as well as the protofibrils of the D23N mutant. Specifically, we investigated the methyl groups of two leucine residues that belong to the hydrophobic core (L17 and L34) as well as M35 residues belonging to the hydrophobic interface between the cross-β subunits, which has been previously found to be water-accessible. Relaxation measurements performed over 310-140 K and two magnetic field strengths provide insights into conformational variability within and between polymorphs. Core packing variations within a single polymorph are similar to what is observed for globular proteins for the core residues, whereas M35 exhibits a larger degree of variability. M35 site is also shown to undergo a solvent-dependent dynamical transition in which slower amplitude motions of methyl axes are activated at high temperature. The motions, modeled as a diffusion of methyl axis, have activation energy by a factor of 2.7 larger in the twofold compared with the threefold polymorph, whereas D23N protofibrils display a value similar to the threefold polymorph. This suggests enhanced flexibility of the hydrophobic interface in the threefold polymorph. This difference is only observed in the hydrated state and is absent in the dry fibrils, highlighting the role of solvent at the cavity. In contrast, the dynamic behavior of the core is hydration-independent.