A study of phenylalanine side-chain dynamics in surface-adsorbed peptides using solid-state deuterium NMR and rotamer library statistics

Carbon-detected deuterium solid-state NMR rotating frame relaxation measurements for protein methyl groups under magic angle spinning

Deuterium rotating frame solid-state NMR relaxation measurements (2H R1ρ) are important tools in quantitative studies of molecular dynamics. We demonstrate how 2H to 13C cross-polarization (CP) approaches under 10-40 kHz magic angle spinning rates can be combined with the 2H R1ρ blocks to allow for extension of deuterium rotating frame relaxation studies to methyl groups in biomolecules. This extension permits detection on the 13C nuclei and, hence, for the achievement of site-specific resolution. The measurements are demonstrated using a nine-residue low complexity peptide with the sequence GGKGMGFGL, in which a single selective -13CD3 label is placed at the methionine residue. Carbon-detected measurements are compared with the deuterium direct-detection results, which allows for fine-tuning of experimental approaches. In particular, we show how the adiabatic respiration CP scheme and the double adiabatic sweep on the 2H and 13C channels can be combined with the 2H R1ρ relaxation rates measurement. Off-resonance 2H R1ρ measurements are investigated in addition to the on-resonance condition, as they extent the range of effective spin-locking field.

Effect of Post-Translational Modifications and Mutations on Amyloid-β Fibrils Dynamics at N Terminus

We investigate the variability in the dynamics of the disordered N-terminal domain of amyloid-β fibrils (Aβ), comprising residues 1-16 of Aβ1-40, due to post-translational modifications and mutations in the β-bend regions known to modulate aggregation properties. Using 2H static solid-state NMR approaches, we compare the dynamics in the wild-type Aβ fibrils in the threefold symmetric polymorph with the fibrils from three post-translational modification sequences: isoaspartate-D7, the phosphorylation of S8, and an N-terminal truncation ΔE3. Additional comparisons are made with the mutants in the β-bend region (residues 21-23) corresponding to the familial Osaka E22Δ deletion and D23N Iowa mutation. We also include the aggregates induced by Zn2+ ions. The dynamics are probed at the F4 and G9 positions. The main motional model involves two free states undergoing diffusion and conformational exchanges with the bound state in which the diffusion is quenched because of transient interactions involving fibril core and other intrastrand contacts. The fraction of the bound state increases in a sigmoidal fashion with a decrease in temperature. There is clear variability in the dynamics: the phosphorylation of S8 variant is the most rigid at the G9 site in line with structural studies, the ΔE3 fibrils are more flexible at the G9 site in line with the morphological fragmentation pattern, the Zn-induced aggregates are the most mobile, and the two β-bend mutants have the strongest changes at the F4 site toward higher rigidity. Overall, the changes underlie the potential role of conformational ensembles in setting the stage for aggregation-prone states.

Solid-State NMR and MD Study of the Structure of the Statherin Mutant SNa15 on Mineral Surfaces

Elucidation of the structure and interactions of proteins at native mineral interfaces is key to understanding how biological systems regulate the formation of hard tissue structures. In addition, understanding how these same proteins interact with non-native mineral surfaces has important implications for the design of medical and dental implants, chromatographic supports, diagnostic tools, and a host of other applications. Here, we combine solid-state NMR spectroscopy, isotherm measurements, and molecular dynamics simulations to study how SNa15, a peptide derived from the hydroxyapatite (HAP) recognition domain of the biomineralization protein statherin, interacts with HAP, silica (SiO2), and titania (TiO2) mineral surfaces. Adsorption isotherms are used to characterize the binding affinity of SNa15 to HAP, SiO2, and TiO2. We also apply 1D 13C CP MAS, 1D 15N CP MAS, and 2D 13C-13C DARR experiments to SNa15 samples with uniformly 13C- and 15N-enriched residues to determine backbone and side-chain chemical shifts. Different computational tools, namely TALOS-N and molecular dynamics simulations, are used to deduce secondary structure from backbone and side-chain chemical shift data. Our results show that SNa15 adopts an α-helical conformation when adsorbed to HAP and TiO2, but the helix largely unravels upon adsorption to SiO2. Interactions with HAP are mediated in general by acidic and some basic amino acids, although the specific amino acids involved in direct surface interaction vary with surface. The integrated experimental and computational approach used in this study is able to provide high-resolution insights into adsorption of proteins on interfaces.

Protein dynamics in the solid-state from 2H NMR lineshape analysis. III. MOMD in the presence of Magic Angle Spinning

We report on a new approach to the analysis of dynamic NMR lineshapes from polycrystalline (i.e., macroscopically disordered) samples in the presence of Magic Angle Spinning (MAS). This is an application of the Stochastic Liouville Equation developed by Freed and co-workers for treating restricted (i.e., microscopically ordered) motions. The ²H nucleus in an internally-mobile C-CD₃ moiety serves as a prototype probe. The acronym is ²H/MOMD/MAS, where MOMD stands for "microscopic-order-macroscopic-disorder." The key elements describing internal motions - their type, the local spatial restrictions, and related features of local geometry - are treated in MOMD generally, within their rigorous three-dimensional tensorial requirements. Based on this representation a single physically well-defined model of local motion has the capability of reproducing experimental spectra. There exist other methods for analyzing dynamic ²H/MAS spectra which advocate simple motional modes. Yet, to reproduce satisfactorily the experimental lineshapes, one has either to use unusual parameter values, or combine several simple motional modes. The multi-simple-mode reasoning assumes independence of the constituent modes, features ambiguity as different simple modes may be used, renders inter-system comparison difficult as the overall models differ, and makes possible model-improvement only by adding yet another simple mode, i.e., changing the overall model. ²H/MOMD/MAS is free of such limitations and inherently provides a clear physical interpretation. These features are illustrated. The advantage of ²H/MOMD/MAS in dealing with sensitive but hardly investigated slow-motional lineshapes is demonstrated by applying it to actual experimental data. The results differ from those obtained previously with a two-site exchange scheme that yielded unusual parameters.

Solid state deuterium NMR study of LKα14 peptide aggregation in biosilica

In nature, organisms including diatoms, radiolaria, and marine sponges use proteins, long chain polyamines, and other organic molecules to regulate the assembly of complex silica-based structures. Here, the authors investigate structural features of small peptides, designed to mimic the silicifying activities of larger proteins found in natural systems. LKα14 (Ac-LKKLLKLLKKLLKL-C), an amphiphilic lysine/leucine repeat peptide with an α-helical secondary structure at polar/apolar interfaces, coprecipitates with silica to form nanospheres. Previous 13C magic angle spinning studies suggest that the tetrameric peptide bundles that LKα14 is known to form in solution may persist in the silica-complexed form, and may also function as catalysts and templates for silica formation. To further investigate LKα14 aggregation in silica, deuterium solid-state nuclear magnetic resonance (2H ssNMR) was used to establish how leucine side-chain dynamics differ in solid LKα14 peptides isolated from aqueous solution, from phosphate-buffered solution, and in the silica-precipitated states. Modeling the 2H ssNMR line shapes probed the mechanisms of peptide preaggregation and silica coprecipitation. The resulting NMR data indicates that the peptide bundles in silica preserve the hydrophobic interior that they display in the hydrated solid state. However, NMR data also indicate free motion of the leucine residues in silica, a condition that may result from structural deformation of the aggregates arising from interactions between the surface lysine side chains and the surrounding silica matrix.

Prediction and clarification of structures of (bio)molecules on surfaces

The design of future materials for biotechnological applications via deposition of molecules on surfaces will require not only exquisite control of the deposition procedure, but of equal importance will be our ability to predict the shapes and stability of individual molecules on various surfaces. Furthermore, one will need to be able to predict the structure patterns generated during the self-organization of whole layers of (bio)molecules on the surface. In this review, we present an overview over the current state of the art regarding the prediction and clarification of structures of biomolecules on surfaces using theoretical and computational methods.

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

We conducted a detailed investigation of the dynamics of two phenylalanine side chains in the hydrophobic core of the villin headpiece subdomain protein (HP36) in the hydrated powder state over the 298-80 K temperature range. Our main tools were static deuteron NMR measurements of longitudinal relaxation and line shapes supplemented with computational modeling. The temperature dependence of the relaxation times reveals the presence of two main mechanisms that can be attributed to the ring-flips, dominating at high temperatures, and small-angle fluctuations, dominating at low temperatures. The relaxation is nonexponential at all temperatures with the extent of nonexponentiality increasing from higher to lower temperatures. This behavior suggests a distribution of conformers with unique values of activation energies. The central values of the activation energies for the ring-flipping motions are among the smallest reported for aromatic residues in peptides and proteins and point to a very mobile hydrophobic core. The analysis of the widths of the distributions, in combination with the earlier results on the dynamics of flanking methyl groups (Vugmeyster et al. J. Phys. Chem. B 2013, 117, 6129-6137), suggests that the hydrophobic core undergoes slow concerted fluctuations. There is a pronounced effect of dehydration on the ring-flipping motions, which shifts the distribution toward more rigid conformers. The crossover temperature between the regions of dominance of the small-angle fluctuations and ring-flips shifts from 195 K in the hydrated protein to 278 K in the dry one. This result points to the role of solvent in softening the core and highlights aromatic residues as markers of the protein dynamical transitions.

Protein dynamics in the solid state from 2H NMR line shape analysis: a consistent perspective

Deuterium line shape analysis of CD3 groups has emerged as a particularly useful tool for studying microsecond-millisecond protein motions in the solid state. The models devised so far consist of several independently conceived simple jump-type motions. They are comprised of physical quantities encoded in their simplest form; improvements are only possible by adding yet another simple motion, thereby changing the model. The various treatments developed are case-specific; hence comparison among the different systems is not possible. Here we develop a new methodology for (2)H NMR line shape analysis free of these limitations. It is based on the microscopic-order-macroscopic-disorder (MOMD) approach. In MOMD motions are described by diffusion tensors, spatial restrictions by potentials/ordering tensors, and geometric features by relative tensor orientations. Jump-type motions are recovered in the limit of large orientational potentials. Model improvement is accomplished by monitoring the magnitude, symmetry, and orientation of the various tensors. The generality of MOMD makes possible comparison among different scenarios. CD3 line shapes from the Chicken Villin Headpiece Subdomain and the Streptomyces Subtilisin Inhibitor are used as experimental examples. All of these spectra are reproduced by using rhombic local potentials constrained for simplicity to be given by the L = 2 spherical harmonics, and by axial diffusion tensors. Potential strength and rhombicity are found to be ca. 2-3 k(B)T. The diffusion tensor is tilted at 120° from the C-CD3 axis. The perpendicular (parallel) correlation times for local motion are 0.1-1.0 ms (3.3-30 μs). Activation energies in the 1.1-8.0 kcal/mol range are estimated. Future prospects include extension to the (2)H relaxation limit, application to the (15)N and (13)C NMR nuclei, and accounting for collective motions and anisotropic media.

Internal dynamics of dynactin CAP-Gly is regulated by microtubules and plus end tracking protein EB1

CAP-Gly domain of dynactin, a microtubule-associated activator of dynein motor, participates in multiple cellular processes, and its point mutations are associated with neurodegenerative diseases. Recently, we have demonstrated that conformational plasticity is an intrinsic property of CAP-Gly. To understand its origin, we addressed internal dynamics of CAP-Gly assembled on polymeric microtubules, bound to end-binding protein EB1 and free, by magic angle spinning NMR and molecular dynamics simulations. The analysis of residue-specific dynamics of CAP-Gly on time scales spanning nano- through milliseconds reveals its unusually high mobility, both free and assembled on polymeric microtubules. On the contrary, CAP-Gly bound to EB1 is significantly more rigid. Molecular dynamics simulations indicate that these motions are strongly temperature-dependent, and loop regions are surprisingly mobile. These findings establish the connection between conformational plasticity and internal dynamics in CAP-Gly, which is essential for the biological functions of CAP-Gly and its ability to bind to polymeric microtubules and multiple binding partners. In this work, we establish an approach, for the first time, to probe atomic resolution dynamic profiles of a microtubule-associated protein assembled on polymeric microtubules. More broadly, the methodology established here can be applied for atomic resolution analysis of dynamics in other microtubule-associated protein assemblies, including but not limited to dynactin, dynein, and kinesin motors assembled on microtubules.