Dynamics of biomolecules from picoseconds to seconds at atomic resolution

Studying micro to millisecond protein dynamics using simple amide 15N CEST experiments supplemented with major-state R2 and visible peak-position constraints

Over the last decade amide 15N CEST experiments have emerged as a popular tool to study protein dynamics that involves exchange between a 'visible' major state and sparsely populated 'invisible' minor states. Although initially introduced to study exchange between states that are in slow exchange with each other (typical exchange rates of, 10 to 400 s-1), they are now used to study interconversion between states on the intermediate to fast exchange timescale while still using low to moderate (5 to 350 Hz) 'saturating' B1 fields. The 15N CEST experiment is very sensitive to exchange as the exchange delay TEX can be quite long (~0.5 s) allowing for a large number of exchange events to occur making it a very powerful tool to detect minor sates populated ([Formula: see text]) to as low as 1%. When systems are in fast exchange and the 15N CEST data has to be described using a model that contains exchange, the exchange parameters are often poorly defined because the [Formula: see text] versus [Formula: see text] and [Formula: see text] versus exchange rate ([Formula: see text]) plots can be quite flat with shallow or no minima and the analysis of such 15N CEST data can lead to wrong estimates of the exchange parameters due to the presence of 'spurious' minima. Here we show that the inclusion of experimentally derived constraints on the intrinsic transverse relaxation rates and the inclusion of visible state peak-positions during the analysis of amide 15N CEST data acquired with moderate B1 values (~50 to ~350 Hz) results in convincing minima in the [Formula: see text] versus [Formula: see text] and the [Formula: see text] versus [Formula: see text] plots even when exchange occurs on the 100 μs timescale. The utility of this strategy is demonstrated on the fast-folding Bacillus stearothermophilus peripheral subunit binding domain that folds with a rate constant ~104 s-1. Here the analysis of 15N CEST data alone results in [Formula: see text] versus [Formula: see text] and [Formula: see text] versus [Formula: see text] plots that contain shallow minima, but the inclusion of visible-state peak positions and restraints on the intrinsic transverse relaxation rates of both states during the analysis of the 15N CEST data results in pronounced minima in the [Formula: see text] versus [Formula: see text] and [Formula: see text] versus [Formula: see text] plots and precise exchange parameters even in the fast exchange regime ([Formula: see text]~5). Using this strategy we find that the folding rate constant of PSBD is invariant (~10,500 s-1) from 33.2 to 42.9 °C while the unfolding rates (~70 to ~500 s-1) and unfolded state populations (~0.7 to ~4.3%) increase with temperature. The results presented here show that protein dynamics occurring on the 10 to 104 s-1 timescale can be studied using amide 15N CEST experiments.

SRLS Analysis of 15N-1H NMR Relaxation from the Protein S100A1: Dynamic Structure, Calcium Binding, and Related Changes in Conformational Entropy

We report on amide (N-H) NMR relaxation from the protein S100A1 analyzed with the slowly relaxing local structure (SRLS) approach. S100A1 comprises two calcium-binding "EF-hands" (helix-loop-helix motifs) connected by a linker. The dynamic structure of this protein, in both calcium-free and calcium-bound form, is described as the restricted local N-H motion coupled to isotropic protein tumbling. The restrictions are given by a rhombic potential, u (∼10 kT), the local motion by a diffusion tensor with rate constant D₂ (∼10⁹ s^-1), and principal axis tilted from the N-H bond at angle β (10-20°). This parameter combination provides a physically insightful picture of the dynamic structure of S100A1 from the N-H bond perspective. Calcium binding primarily affects the C-terminal EF-hand, among others slowing down the motion of helices III and IV approximately 10-fold. Overall, it brings about significant changes in the shape of the local potential, u, and the orientation of the local diffusion axis, β. Conformational entropy derived from u makes an unfavorable entropic contribution to the free energy of calcium binding estimated at 8.6 ± 0.5 kJ/mol. The N-terminal EF-hand undergoes moderate changes. These findings provide new insights into the calcium-binding process. The same data were analyzed previously with the extended model-free (EMF) method, which is a simple limit of SRLS. In that interpretation, the protein tumbles anisotropically. Locally, calcium binding increases ordering in the loops of S100A1 and conformational exchange (R_ex) in the helices of its N-terminal EF-hand. These are very unusual features. We show that they most likely stem from problematic data-fitting, oversimplifications inherent in EMF, and experimental imperfections. R_ex is shown to be mainly a fit parameter. By reanalyzing the experimental data with SRLS, which is largely free of these deficiencies, we obtain-as delineated above-physically-relevant structural, kinetic, geometric, and binding information.

A CEST NMR experiment to obtain glycine 1Hα chemical shifts in 'invisible' minor states of proteins

Chemical exchange saturation transfer (CEST) experiments are routinely used to study protein conformational exchange between a 'visible' major state and 'invisible' minor states because they can detect minor states with lifetimes varying from ~ 3 to ~ 100 ms populated to just ~ 0.5%. Consequently several ¹H, ¹⁵N and ¹³C CEST experiments have been developed to study exchange and obtain minor state chemical shifts at almost all backbone and sidechain sites in proteins. Conspicuously missing from this extensive set of CEST experiments is a ¹H CEST experiment to study exchange at glycine (Gly) ¹H^α sites as the existing ¹H CEST experiments that have been designed to study dynamics in amide ¹H-¹⁵N spin systems and methyl ¹³CH₃ groups with three equivalent protons while suppressing ¹H-¹H NOE induced dips are not suitable for studying exchange in methylene ¹³CH₂ groups with inequivalent protons. Here a Gly ¹H^α CEST experiment to obtain the minor state Gly ¹H^α chemical shifts is presented. The utility of this experiment is demonstrated on the L99A cavity mutant of T4 Lysozyme (T4L L99A) that undergoes conformational exchange between two compact conformers. The CEST derived minor state Gly ¹H^α chemical shifts of T4L L99A are in agreement with those obtained previously using CPMG techniques. The experimental strategy presented here can also be used to obtain methylene proton minor state chemical shifts from protein sidechain and nucleic acid backbone sites.

Tuning a timing device that regulates lateral root development in rice

Peptidyl Prolyl Isomerases (PPIases) accelerate cis-trans isomerization of prolyl peptide bonds. In rice, the PPIase LRT2 is essential for lateral root initiation. LRT2 displays in vitro isomerization of a highly conserved W-P peptide bond (¹⁰⁴W-P¹⁰⁵) in the natural substrate OsIAA11. OsIAA11 is a transcription repressor that, in response to the plant hormone auxin, is targeted to ubiquitin-mediated proteasomal degradation via specific recognition of the cis isomer of its ¹⁰⁴W-P¹⁰⁵ peptide bond. OsIAA11 controls transcription of specific genes, including its own, that are required for lateral root development. This auxin-responsive negative feedback circuit governs patterning and development of lateral roots along the primary root. The ability to tune LRT2 activity via mutagenesis is crucial for understanding and modeling the role of this bimodal switch in the auxin circuit and lateral root development. We present characterization of the thermal stability and isomerization rates of several LRT2 mutants acting on the OsIAA11 substrate. The thermally stable mutants display activities lower than that of wild-type (WT) LRT2. These include binding diminished but catalytically active P125K, binding incompetent W128A, and binding capable but catalytically incompetent H133Q mutations. Additionally, LRT2 homologs hCypA from human, TaCypA from Triticum aestivum (wheat) and PPIB from E. coli were shown to have 110, 50 and 60% of WT LRT2 activity on the OsIAA11 substrate. These studies identify several thermally stable LRT2 mutants with altered activities that will be useful for establishing relationships between cis-trans isomerization, auxin circuit dynamics, and lateral root development in rice.

Methyl-Based NMR Spectroscopy Methods for Uncovering Structural Dynamics in Large Proteins and Protein Complexes

NMR spectroscopy is particularly adept at site-specifically monitoring dynamic processes in proteins, such as protein folding, domain movements, ligand binding, and side-chain rotations. By coupling the favorable spectroscopic properties of highly dynamic side-chain methyl groups with transverse-relaxation-optimized spectroscopy (TROSY), it is now possible to routinely study such dynamic processes in high-molecular-weight proteins and complexes approaching 1 MDa. In this Perspective, we describe many elegant methyl-based NMR experiments that probe slow (second) to fast (picosecond) dynamics in large systems. To demonstrate the power of these methods, we also provide interesting examples of studies that utilized each methyl-based NMR technique to uncover functionally important dynamics. In many cases, the NMR experiments are paired with site-directed mutagenesis and/or other biochemical assays to put the dynamics and function into context. Our vision of the future of structural biology involves pairing methyl-based NMR spectroscopy with biochemical studies to advance our knowledge of the motions large proteins and macromolecular complexes use to choreograph complex functions. Such studies will be essential in elucidating the critical structural dynamics that underlie function and characterizing alterations in these processes that can lead to human disease.

Simple and Robust Study of Backbone Dynamics of Crystalline Proteins Employing 1H-15N Dipolar Coupling Dispersion

We report a new solid-state multidimensional NMR approach based on the cross-polarization with variable-contact pulse sequence [ Paluch , P. ; Pawlak , T. ; Amoureux , J.-P. ; Potrzebowski , M. J. J. Magn. Reson. 233 , 2013 , 56 ], with ¹H inverse detection and very fast magic angle spinning (ν_R = 60 kHz), dedicated to the measurement of local molecular motions of ¹H-¹⁵N vectors. The introduced three-dimensional experiments, ¹H-¹⁵N-¹H and hCA(N)H, are particularly useful for the study of molecular dynamics of proteins and other complex structures. The applicability and power of this methodology have been revealed by employing as a model sample the GB-1 small protein doped with Na₂CuEDTA. The results clearly prove that the dispersion of ¹H-¹⁵N dipolar coupling constants well correlates with higher order structure of the protein. Our approach complements the conventional studies and offers a fast and reasonably simple method.

Conformational switch of the bacterial adhesin FimH in the absence of the regulatory domain: Engineering a minimalistic allosteric system

For many biological processes such as ligand binding, enzymatic catalysis, or protein folding, allosteric regulation of protein conformation and dynamics is fundamentally important. One example is the bacterial adhesin FimH, where the C-terminal pilin domain exerts negative allosteric control over binding of the N-terminal lectin domain to mannosylated ligands on host cells. When the lectin and pilin domains are separated under shear stress, the FimH-ligand interaction switches in a so-called catch-bond mechanism from the low- to high-affinity state. So far, it has been assumed that the pilin domain is essential for the allosteric propagation within the lectin domain that would otherwise be conformationally rigid. To test this hypothesis, we generated mutants of the isolated FimH lectin domain and characterized their thermodynamic, kinetic, and structural properties using isothermal titration calorimetry, surface plasmon resonance, nuclear magnetic resonance, and X-ray techniques. Intriguingly, some of the mutants mimicked the conformational and kinetic behaviors of the full-length protein and, even in absence of the pilin domain, conducted the cross-talk between allosteric sites and the mannoside-binding pocket. Thus, these mutants represent a minimalistic allosteric system of FimH, useful for further mechanistic studies and antagonist design.

A J-modulated protonless NMR experiment characterizes the conformational ensemble of the intrinsically disordered protein WIP

Intrinsically disordered proteins (IDPs) are multi-conformational polypeptides that lack a single stable three-dimensional structure. It has become increasingly clear that the versatile IDPs play key roles in a multitude of biological processes, and, given their flexible nature, NMR is a leading method to investigate IDP behavior on the molecular level. Here we present an IDP-tailored J-modulated experiment designed to monitor changes in the conformational ensemble characteristic of IDPs by accurately measuring backbone one- and two-bond J(¹⁵N,¹³Cα) couplings. This concept was realized using a unidirectional (H)NCO ¹³C-detected experiment suitable for poor spectral dispersion and optimized for maximum coverage of amino acid types. To demonstrate the utility of this approach we applied it to the disordered actin-binding N-terminal domain of WASp interacting protein (WIP), a ubiquitous key modulator of cytoskeletal changes in a range of biological systems. One- and two-bond J(¹⁵N,¹³Cα) couplings were acquired for WIP residues 2-65 at various temperatures, and in denaturing and crowding environments. Under native conditions fitted J-couplings identified in the WIP conformational ensemble a propensity for extended conformation at residues 16-23 and 45-60, and a helical tendency at residues 28-42. These findings are consistent with a previous study of the based upon chemical shift and RDC data and confirm that the WIP^2-65 conformational ensemble is biased towards the structure assumed by this fragment in its actin-bound form. The effects of environmental changes upon this ensemble were readily apparent in the J-coupling data, which reflected a significant decrease in structural propensity at higher temperatures, in the presence of 8 M urea, and under the influence of a bacterial cell lysate. The latter suggests that crowding can cause protein unfolding through protein-protein interactions that stabilize the unfolded state. We conclude that J-couplings are a useful measureable in characterizing structural ensembles in IDPs, and that the proposed experiment provides a practical method for accurately performing such measurements, once again emphasizing the power of NMR in studying IDP behavior.

Direct Investigation of Slow Correlated Dynamics in Proteins via Dipolar Interactions

The synchronization of native state motions as they transition between microstates influences catalysis kinetics, mediates allosteric interactions, and reduces the conformational entropy of proteins. However, it has proven difficult to describe native microstates because they are usually minimally frustrated and may interconvert on the micro- to millisecond time scale. Direct observation of concerted equilibrium fluctuations would therefore be an important tool for describing protein native states. Here we propose a strategy that relates NMR cross-correlated relaxation (CCR) rates between dipolar interactions to residual dipolar couplings (RDCs) of individual consecutive H(N)-N and H(α)-C(α) bonds, which act as a proxy for the peptide planes and the side chains, respectively. Using Xplor-NIH ensemble structure calculations restrained with the RDC and CCR data, we observe collective motions on time scales slower than nanoseconds in the backbone for GB3. To directly access the correlations from CCR, we develop a structure-free data analysis. The resulting dynamic correlation map is consistent with the ensemble-restrained simulations and reveals a complex network. In general, we find that the bond motions are on average slightly correlated and that the local environment dominates many observations. Despite this, some patterns are typical over entire secondary structure elements. In the β-sheet, nearly all bonds are weakly correlated, and there is an approximately binary alternation in correlation intensity corresponding to the solvent exposure/shielding alternation of the side chains. For α-helices, there is also a weak correlation in the H(N)-N bonds. The degree of correlation involving H(α)-C(α) bonds is directly affected by side-chain fluctuations, whereas loops show complex and nonuniform behavior.

Investigating protein conformational energy landscapes and atomic resolution dynamics from NMR dipolar couplings: a review

Nuclear magnetic resonance spectroscopy is exquisitely sensitive to protein dynamics. In particular inter-nuclear dipolar couplings, that become measurable in solution when the protein is dissolved in a dilute liquid crystalline solution, report on all conformations sampled up to millisecond timescales. As such they provide the opportunity to describe the Boltzmann distribution present in solution at atomic resolution, and thereby to map the conformational energy landscape in unprecedented detail. The development of analytical methods and approaches based on numerical simulation and their application to numerous biologically important systems is presented.

Analysis of local molecular motions of aromatic sidechains in proteins by 2D and 3D fast MAS NMR spectroscopy and quantum mechanical calculations

We report a new multidimensional magic angle spinning NMR methodology, which provides an accurate and detailed probe of molecular motions occurring on timescales of nano- to microseconds, in sidechains of proteins. The approach is based on a 3D CPVC-RFDR correlation experiment recorded under fast MAS conditions (ν(R) = 62 kHz), where (13)C-(1)H CPVC dipolar lineshapes are recorded in a chemical shift resolved manner. The power of the technique is demonstrated in model tripeptide Tyr-(d)Ala-Phe and two nanocrystalline proteins, GB1 and LC8. We demonstrate that, through numerical simulations of dipolar lineshapes of aromatic sidechains, their detailed dynamic profile, i.e., the motional modes, is obtained. In GB1 and LC8 the results unequivocally indicate that a number of aromatic residues are dynamic, and using quantum mechanical calculations, we correlate the molecular motions of aromatic groups to their local environment in the crystal lattice. The approach presented here is general and can be readily extended to other biological systems.

An optimized method for (15)N R(1) relaxation rate measurements in non-deuterated proteins

(15)N longitudinal relaxation rates are extensively used for the characterization of protein dynamics; however, their accurate measurement is hindered by systematic errors. (15)N CSA/(1)H-(15)N dipolar cross-correlated relaxation (CC) and amide proton exchange saturation transfer from water protons are the two main sources of systematic errors in the determination of (15)N R1 rates through (1)H-(15)N HSQC-based experiments. CC is usually suppressed through a train of 180° proton pulses applied during the variable (15)N relaxation period (T), which can perturb water magnetization. Thus CC cancellation is required in such a way as to minimize water saturation effects. Here we examined the level of water saturation during the T period caused by various types of inversion proton pulses to suppress CC: (I) amide-selective IBURP-2; (II) cosine-modulated IBURP-2; (III) Watergate-like blocks; and (IV) non-selective hard. We additionally demonstrate the effect of uncontrolled saturation of aliphatic protons on (15)N R1 rates. In this paper we present an optimized pulse sequence that takes into account the crucial effect of controlling also the saturation of the aliphatic protons during (15)N R1 measurements in non-deuterated proteins. We show that using cosine-modulated IBURP-2 pulses spaced 40 ms to cancel CC in this optimized pulse program is the method of choice to minimize systematic errors coming from water and aliphatic protons saturation effects.

NMR studies of dynamic biomolecular conformational ensembles

Multidimensional heteronuclear NMR approaches can provide nearly complete sequential signal assignments of isotopically enriched biomolecules. The availability of assignments together with measurements of spin relaxation rates, residual spin interactions, J-couplings and chemical shifts provides information at atomic resolution about internal dynamics on timescales ranging from ps to ms, both in solution and in the solid state. However, due to the complexity of biomolecules, it is not possible to extract a unique atomic-resolution description of biomolecular motions even from extensive NMR data when many conformations are sampled on multiple timescales. For this reason, powerful computational approaches are increasingly applied to large NMR data sets to elucidate conformational ensembles sampled by biomolecules. In the past decade, considerable attention has been directed at an important class of biomolecules that function by binding to a wide variety of target molecules. Questions of current interest are: "Does the free biomolecule sample a conformational ensemble that encompasses the conformations found when it binds to various targets; and if so, on what time scale is the ensemble sampled?" This article reviews recent efforts to answer these questions, with a focus on comparing ensembles obtained for the same biomolecules by different investigators. A detailed comparison of results obtained is provided for three biomolecules: ubiquitin, calmodulin and the HIV-1 trans-activation response RNA.

Chimeric Avidin--NMR structure and dynamics of a 56 kDa homotetrameric thermostable protein

Chimeric avidin (ChiAVD) is a product of rational protein engineering remarkably resistant to heat and harsh conditions. In quest of the fundamentals behind factors affecting stability we have elucidated the solution NMR spectroscopic structure of the biotin–bound form of ChiAVD and characterized the protein dynamics through ¹⁵N relaxation and hydrogen/deuterium (H/D) exchange of this and the biotin–free form. To surmount the challenges arising from the very large size of the protein for NMR spectroscopy, we took advantage of its high thermostability. Conventional triple resonance experiments for fully protonated proteins combined with methyl–detection optimized experiments acquired at 58°C were adequate for the structure determination of this 56 kDa protein. The model–free parameters derived from the ¹⁵N relaxation data reveal a remarkably rigid protein at 58°C in both the biotin–bound and the free forms. The H/D exchange experiments indicate a notable increase in hydrogen protection upon biotin binding.

Structural basis of conformational transitions in the active site and 80's loop in the FK506-binding protein FKBP12

The extensive set of NMR doublings exhibited by the immunophilin FKBP12 (FK506-binding protein 12) arose from a slow transition to the cis-peptide configuration at Gly⁸⁹ near the tip of the 80′s loop, the site for numerous protein-recognition interactions for both FKBP12 and other FKBP domain proteins. The 80′s loop also exhibited linebroadening, indicative of microsecond to millisecond conformational dynamics, but only in the trans-peptide state. The G89A variant shifted the trans–cis peptide equilibrium from 88:12 to 33:67, whereas a proline residue substitution induced fully the cis-peptide configuration. The 80′s loop conformation in the G89P crystal structure at 1.50 Å resolution differed from wild-type FKBP12 primarily at residues 88, 89 and 90, and it closely resembled that reported for FKBP52. Structure-based chemical-shift predictions indicated that the microsecond to millisecond dynamics in the 80′s loop probably arose from a concerted main chain (ψ₈₈ and ϕ₈₉) torsion angle transition. The indole side chain of Trp⁵⁹ at the base of the active-site cleft was reoriented ~90^o and the adjacent backbone was shifted in the G89P crystal structure. NOE analysis of wild-type FKBP12 demonstrated that this indole populates the perpendicular orientation at 20%. The ¹⁵N relaxation analysis was consistent with the indole reorientation occurring in the nanosecond timeframe. Recollection of the G89P crystal data at 1.20 Å resolution revealed a weaker wild-type-like orientation for the indole ring. Differences in the residues that underlie the Trp⁵⁹ indole ring and altered interactions linking the 50′s loop to the active site suggested that reorientation of this ring may be disfavoured in the other six members of the FKBP domain family that bear this active-site tryptophan residue.

Extensive resonance doubling arises from a cis–trans peptide transition at Gly⁸⁹, whereas linebroadening appears due to a concerted shift in the neighbouring torsion angles. The active site Trp⁵⁹ ring adopts a perpendicular orientation at a population of 20%.

NMR determines transient structure and dynamics in the disordered C-terminal domain of WASp interacting protein

WASp-interacting protein (WIP) is a 503-residue proline-rich polypeptide expressed in human T cells. The WIP C-terminal domain binds to Wiskott-Aldrich syndrome protein (WASp) and regulates its activation and degradation, and the WIP-WASp interaction has been shown to be critical for actin polymerization and implicated in the onset of WAS and X-linked thrombocytopenia. WIP is predicted to be an intrinsically disordered protein, a class of polypeptides that are of great interest because they violate the traditional structure-function paradigm. In this first (to our knowledge) study of WIP in its unbound state, we used NMR to investigate the biophysical behavior of WIP(C), a C-terminal domain fragment of WIP that includes residues 407-503 and contains the WASp-binding site. In light of the poor spectral dispersion exhibited by WIP(C) and the high occurrence (25%) of proline residues, we employed 5D-NMR(13)C-detected NMR experiments with nonuniform sampling to accomplish full resonance assignment. Secondary chemical-shift analysis, (15)N relaxation rates, and protection from solvent exchange all concurred in detecting transient structure located in motifs that span the WASp-binding site. Residues 446-456 exhibited a propensity for helical conformation, and an extended conformation followed by a short, capped helix was observed for residues 468-478. The (13)C-detected approach allows chemical-shift assignment in the WIP(C) polyproline stretches and thus sheds light on their conformation and dynamics. The effects of temperature on chemical shifts referenced to a denatured sample of the polypeptide demonstrate that heating reduces the structural character of WIP(C). Thus, we conclude that the disordered WIP(C) fragment is comprised of regions with latent structure connected by flexible loops, an architecture with implications for binding affinity and function.

The time correlation function perspective of NMR relaxation in proteins

We applied over a decade ago the two-body coupled-rotator slowly relaxing local structure (SRLS) approach to NMR relaxation in proteins. One rotator is the globally moving protein and the other rotator is the locally moving probe (spin-bearing moiety, typically the (15)N-(1)H bond). So far we applied SRLS to (15)N-H relaxation from seven different proteins within the scope of the commonly used data-fitting paradigm. Here, we solve the SRLS Smoluchowski equation using typical best-fit parameters as input, to obtain the corresponding generic time correlation functions (TCFs). The following new information is obtained. For actual rhombic local ordering and main ordering axis pointing along C(i-1)(α)-C(i)(α), the measurable TCF is dominated by the (K,K') = (-2,2), (2,2), and (0,2) components (K is the order of the rank 2 local ordering tensor), determined largely by the local motion. Global diffusion axiality affects the analysis significantly when the ratio between the parallel and perpendicular components exceeds approximately 1.5. Local diffusion axiality has a large and intricate effect on the analysis. Mode-coupling becomes important when the ratio between the global and local motional rates falls below 0.01. The traditional method of analysis--model-free (MF)--represents a simple limit of SRLS. The conditions under which the MF and SRLS TCFs are the same are specified. The validity ranges of wobble-in-a-cone and rotation on the surface of a cone as local motions are determined. The evolution of the intricate Smoluchowski operator from the simple diffusion operator for a sphere reorienting in isotropic medium is delineated. This highlights the fact that SRLS is an extension of the established stochastic theories for treating restricted motions. This study lays the groundwork for TCF-based comparison between mesoscopic SRLS and atomistic molecular dynamics.

Atomic-Resolution Structural Dynamics in Crystalline Proteins from NMR and Molecular Simulation

Solid-state NMR can provide atomic-resolution information about protein motions occurring on a vast range of time scales under similar conditions to those of X-ray diffraction studies and therefore offers a highly complementary approach to characterizing the dynamic fluctuations occurring in the crystal. We compare experimentally determined dynamic parameters, spin relaxation, chemical shifts, and dipolar couplings, to values calculated from a 200 ns MD simulation of protein GB1 in its crystalline form, providing insight into the nature of structural dynamics occurring within the crystalline lattice. This simulation allows us to test the accuracy of commonly applied procedures for the interpretation of experimental solid-state relaxation data in terms of dynamic modes and time scales. We discover that the potential complexity of relaxation-active motion can lead to significant under- or overestimation of dynamic amplitudes if different components are not taken into consideration.

Force-field dependence of chignolin folding and misfolding: comparison with experiment and redesign

We study the folding of the designed hairpin chignolin, using simulations with four different force fields. Interestingly, we find a misfolded, out-of-register, structure comprising 20-50% of the ordered structures with three force fields, but not with a fourth. A defining feature of the misfold is that Gly-7 adopts a β(PR) conformation rather than α(L). By reweighting, we show that differences between the force fields can mostly be attributed to differences in glycine properties. Benchmarking against NMR data suggests that the preference for β(PR) is not a force-field artifact. For chignolin, we show that including the misfold in the ensemble results in back-recalculated NMR observables in slightly better agreement with experiment than parameters calculated from a folded ensemble only. For comparison, we show by NMR and circular dichroism spectroscopy that the G7K mutant of chignolin, in which formation of this misfold is impossible, is well folded with stability similar to the wild-type and does not populate the misfolded state in simulation. Our results highlight the complexity of interpreting NMR data for small, weakly structured, peptides in solution, as well as the importance of accurate glycine parameters in force fields, for a correct description of turn structures.

Site-resolved measurement of microsecond-to-millisecond conformational-exchange processes in proteins by solid-state NMR spectroscopy

We demonstrate that conformational exchange processes in proteins on microsecond-to-millisecond time scales can be detected and quantified by solid-state NMR spectroscopy. We show two independent approaches that measure the effect of conformational exchange on transverse relaxation parameters, namely Carr-Purcell-Meiboom-Gill relaxation-dispersion experiments and measurement of differential multiple-quantum coherence decay. Long coherence lifetimes, as required for these experiments, are achieved by the use of highly deuterated samples and fast magic-angle spinning. The usefulness of the approaches is demonstrated by application to microcrystalline ubiquitin. We detect a conformational exchange process in a region of the protein for which dynamics have also been observed in solution. Interestingly, quantitative analysis of the data reveals that the exchange process is more than 1 order of magnitude slower than in solution, and this points to the impact of the crystalline environment on free energy barriers.

Tyrosine residues mediate fibril formation in a dynamic light chain dimer interface

Light chain amyloidosis is an incurable protein misfolding disease where monoclonal immunoglobulin light chains misfold and deposit as amyloid fibrils, causing organ failure and death. Previously, we determined that amyloidogenic light chains AL-09 and AL-103 do not form fibrils at pH 10 (tyrosine pK(a)). There are three tyrosine residues (32, 91, and 96) clustered in the dimer interface, interacting differently in the two light chain proteins due to their two different dimer conformations. These tyrosines may be ionized at pH 10, causing repulsion and inhibiting fibril formation. Here, we characterize single and double Tyr-to-Phe mutations in AL-09 and AL-103. All AL-09 Tyr-to-Phe mutants form fibrils at pH 10, whereas none of the AL-103 mutants form fibrils at pH 10. NMR studies suggest that although both AL-09 and AL-103 present conformational heterogeneity, only AL-09 favors dimer conformations where tyrosine residues mediate crucial interactions for amyloid formation.