Slowly Relaxing Local Structure Research Articles

Slowly Relaxing Local Structure Analysis of 15N Relaxation from the Proteins p50 and Human Neutrophil Gelatinase-Associated Lipocalin: New Insights into the Dynamic Structure of β-Barrel Proteins

Nuclear magnetic resonance relaxation analysis is a powerful method for studying the internal mobility of proteins. We have developed for analysis the slowly relaxing local structure (SRLS) approach. SRLS is general in its nature in several respects, including the tensorial representation of the physical quantities comprising the dynamic model. By controlling tensor symmetry, a broad range of systems can be treated with physical relevance, typically with data-fitting techniques. In simple limits, SRLS yields the traditional model-free (MF) method. In the present context, MF simplicity means featuring the highest possible tensor symmetry. This renders MF-based data-fitting susceptible to the usage of fit parameters, yielding physically ill-defined results. A typical candidate is the Rex term, devised to represent ms-μs motions but often invoked by the fitting scheme just to improve the statistics. Here, we consider two such cases using the N-H bond as probe and the proteins p50 and human neutrophil gelatinase-associated lipocalin as paradigm systems. We illustrate the harm caused by the physically unjustified involvement of Rex in MF-based 15N relaxation analysis. Then, we show that forgoing the usage of Rex, SRLS analysis of the very same experimental data provides interesting new information.

Structural Dynamics from NMR Relaxation by SRLS Analysis: Local Geometry, Potential Energy Landscapes, and Spectral Densities.

We have developed the two-body coupled-rotator slowly relaxing local structure (SRLS) approach for elucidating protein dynamics by nuclear magnetic resonance (NMR) relaxation. The rotators are represented by diffusion tensors D1 for overall protein tumbling and D2 for locally ordered probe motion. D1 and D2 are coupled dynamically by a potential, u, typically given by linear combinations of the Wigner functions D002 and (D022 + D0-22). Until now, our SRLS analyses provided the tensors, D1 and D2, the potential, u, and the geometric link between SRLS and NMR. Here we enhance this description by also examining the SRLS spectral densities obtained by solving the SRLS Smoluchowski equation. In addition, we show that the form of u specified above complies with two NMR-detected potential energy landscapes representing preferential ordering along N-H or Cα-Cα. Pictorial illustrations thereof are provided. The extended SRLS analysis is applied to 15N-H relaxation from the carbohydrate recognition domain of galectin-3 (Gal3C) in complex with two diastereomeric ligands, S and R. We find that D2 is isotropic with a principal value, D2, of 1010 s-1 on average, and it is faster in the strands β3, β5, and β8. The potential, u, is strong (∼20 kT); it is slightly rhombic when N-H is the main ordering axis and highly rhombic when Cα-Cα is the main ordering axis. Gal3C-S exhibits primarily preferential ordering along Cα-Cα; Gal3C-R exhibits both types of ordering. The binding-associated polypeptide chain segment of Gal3C-S is homogeneous, whereas that of Gal3C-R is diversified, with regard to D2 and ordering preference. We associate these features with the previously determined diminished binding constant of Gal3C-R in comparison with Gal3C-S. Thus, the present study enhances the SRLS analysis, in general, and provides new insights into the dynamic structure and binding properties of Gal3C-S and Gal3C-R, in particular.

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 D2 (∼109 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 (Rex) 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. Rex 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.

Conformational Entropy from Mobile Bond Vectors in Proteins: A Viewpoint that Unifies NMR Relaxation Theory and Molecular Dynamics Simulation Approaches.

A new method for determining conformational entropy in proteins is reported. Proteins prevail as conformational ensembles, p ∝ exp(-u). By selecting a bond vector (e.g., N-H) as a conformation representative, molecular dynamics simulations can provide (relative to a reference structure) p as approximate Boltzmann probability density and u as N-H potential of mean force (POMF). The latter is as accurate as implied by the force field but statistical in character; this limits the insights it can provide and its utilization. Conformational entropy is given exclusively by u. Deriving it from POMFs renders it accurate but statistical in character. Previously, we devised explicit (i.e., analytical but not exact) potentials made of Wigner functions, D0KL, with L ≤ 4, which closely resemble the corresponding POMFs in form; hence, they also approach the latter in accuracy. Such potentials can be beneficially characterized/compared in terms of composition, symmetry, and associated order parameters. In this study, we develop a method for deriving conformational entropy from them, which also features the benefits specified above. The method developed is applied to the dimerization of the Rho GTPase-binding domain of plexin-B1. Insights into local ordering, entropy compensation, and features of allostery are gained. In previous work, we developed the slowly relaxing local structure (SRLS) approach for the analysis of NMR relaxation from restricted bond vector motion in proteins. SRLS comprises explicit (restricting) potentials of the kind developed here. It also comprises diffusion tensors describing the local motion and related features of local geometry. The complete model fits experimental data. In future work, the explicit potentials developed here will be inserted unchanged in SRLS-based data fitting, thereby improving the picture of structural dynamics. Given that SRLS is unique in featuring potentials that can closely approach the corresponding POMFs in accuracy, the present study is an important step toward generally improving protein dynamics by NMR relaxation.

Local Ordering at Mobile Sites in Proteins: Combining Perspectives from NMR Relaxation and Molecular Dynamics.

We report on progress toward improving NMR relaxation analysis in proteins in terms of the slowly relaxing local structure (SRLS) approach by developing a method that combines SRLS with molecular dynamics (MD) simulations. 15N-H bonds from the Rho GTPase binding domain of plexin-B1 are used as test case. We focus on the locally restricting/ordering potential of mean force (POMF), u(θ,φ), at the N-H site (θ and φ specify the orientation of the N-H bond in the protein). In SRLS, u(θ,φ) is expanded in the basis set of the real linear combinations of the Wigner rotation matrix elements with M = 0, D L,| K|(θ,φ). Because of limited data sensitivity, only the lowest ( L = 2) terms are preserved; this potential function is denoted by u(SRLS). In MD, the force-field-parametrized POMF is the potential, u(MD), defined in terms of the probability distribution, Peq(MD) ∝ exp(- u(MD)). Peq(MD), and subsequently u(MD), can be derived from the MD trajectory as histograms. One might contemplate utilizing u(MD) instead of u(SRLS); however, histograms cannot be used in SRLS analyses. Here, we approximate u(θ,φ) in terms of linear combinations of the D L,| K| functions with L = 1-4 and appropriate symmetry, denoted by u(DLK), and optimize the latter (via Peq) against u(MD). This yields for every N-H bond an analytical ordering potential, u(DLK-BEST), which exceeds u(SRLS) considerably in accuracy. u(DLK-BEST) can be used fixed in SRLS data fitting, thereby enabling the determination of additional parameters. This yields a substantially improved picture of structural dynamics, which is a significant benefit. The primary achievement of this work is to have employed for the first time MD data to derive a suitable (in terms of composition and symmetry) approximation to the SRLS POMF.

Conformational Entropy from Slowly Relaxing Local Structure Analysis of 15N-H Relaxation in Proteins: Application to Pheromone Binding to MUP-I in the 283-308 K Temperature Range.

The slowly relaxing local structure (SRLS) approach is applied to 15N-H relaxation from the major urinary protein I (MUP-I), and its complex with pheromone 2-sec-butyl-4,5-dihydrothiazol. The objective is to elucidate dynamics, and binding-induced changes in conformational entropy. Experimental data acquired previously in the 283-308 K temperature range are used. The N-H bond is found to reorient globally with correlation time, τ1,0, and locally with correlation time, τ2,0, where τ1,0 ≫ τ2,0. The local motion is restricted by the potential u = -c02D002, where D002 is the Wigner rotation matrix element for L = 2, K = 0, and c02 evaluates the strength of the potential. u yields straightforwardly the order parameter, ⟨D002⟩, and the conformational entropy, Sk, both given by Peq = exp(-u). The deviation of the local ordering/local diffusion axis from the N-H bond, given by the angle β, is also determined. We find that c02 ≅ 18 ± 4 and τ2,0 = 0-170 ps for ligand-free MUP-I, whereas c02 ≅ 15 ± 4 and τ2,0 = 20-270 ps for ligand-bound MUP-I. β is in the 0-10° range. c02 and τ2,0 decrease, whereas β increases, when the temperature is increased from 283 to 308 K. Thus, SRLS provides physically well-defined structure-related (c02 and ⟨D002⟩), motion-related (τ2,0), geometry-related (β), and binding-related (Sk) local parameters, and their temperature-dependences. Intriguingly, upon pheromone binding the conformational entropy of MUP-I decreases at high temperature and increases at low temperature. The very same experimental data were analyzed previously with the model-free (MF) method which yielded "global" (in this context, "relating to the entire 283-308 K range") amplitude (S2) and rate (τe) of the local motion, and a phenomenological exchange term (Rex). S2 is found to decrease (implying implicitly "global" increase in Sk) upon pheromone binding.

15N-H-Related Conformational Entropy Changes Entailed By Plexin-B1 RBD Dimerization: Combined Molecular Dynamics/NMR Relaxation Approach.

We report on a new method for determining function-related conformational entropy changes in proteins. Plexin-B1 RBD dimerization serves as example, and internally mobile N-H bonds serve as probes. Sk (entropy in units of kBT) is given by -∫(PeqlnPeq)dΩ, where Peq = exp(-u) is the probability density for probe orientation, and u the local potential. Previous slowly relaxing local structure (SRLS) analyses of 15N-H relaxation in proteins determined linear combinations of D002(Ω) and (D022(Ω) + D0-22(Ω)) (D0KL(Ω) represents a Wigner rotation matrix element in uniaxial local medium) as "best-fit" form of u. SRLS also determined the "best-fit" orientation of the related ordering tensor. On the basis of this information the coefficients (in the linear combination) of the terms specified above are determined with molecular dynamics (MD) simulations. With the explicit expression for u thus in hand, Sk is calculated. We find that in general Sk decreases, i.e., the local order increases, upon plexin-B1 RBD dimerization. The largest decrease in Sk occurs in the helices α1 and α2, followed by the α2/β6 turn. Only the relatively small peripheral β2 strand, β2/α1 turn, and L3 loop become more disordered. That α-helices dominate ΔSk = Sk(dimer) - Sk(monomer), a few peripheral outliers partly counterbalance the overall decrease in Sk, and the probability density function, Peq, has rhombic symmetry given that the underlying potential function, u, has rhombic symmetry, are interesting features. We also derive S2 (the proxy of u in the simple "model-free (MF)" limit of SRLS) with MD. Its conversion into a potential requires assumptions and adopting a simple axial form of u. Ensuing ΔSk(MF) profiles are u-dependent and differ from ΔSk(SRLS). A method that provides consistent, general, and accurate Sk, atomistic/mesoscopic in nature, has been developed. Its ability to provide new insights in protein research has been illustrated.

Conformational Entropy from NMR Relaxation in Proteins: The SRLS Perspective.

Conformational entropy changes associated with bond-vector motions in proteins contribute to the free energy of ligand-binding. To derivesuch contributions, we apply the slowly relaxing local structure (SRLS) approachtoNMRrelaxation from 15N-H bonds or C-CDH2 moieties of several proteinsinfree and ligand-bound form. The spatial restraints on probe motion,whichdetermine the extent of local order, are expressed in SRLS byawell-definedpotential, u(θ). The latter yields the orientational probabilitydensity,Peq=exp(-u(θ)), and hence the related conformational entropy, Ŝ=-∫Peq(θ) ln[Peq(θ)] sin θ dθ (Ŝ is "entropy" in units of kBT, and θ represents the bond-vector orientation in the protein). SRLS is applied to 4-oxalocrotonate tautomerase (4-OT), the acyl-coenzyme A binding protein (ACBP), the C-terminal SH2 domain of phospholipase Cγ1 (PLCγ1C SH2), the construct dihydrofolate reductase-E:folate (DHFR-E:folate), and their complexes with appropriate ligands, to determine ΔŜ. Eglin C and its V18A and V34A mutants are also studied. Finally, SRLS is applied to the structurally homologous proteins TNfn3 and FNfn10 to characterize within its scope the unusual "dynamics" of the TNfn3 core. Upon ligand-binding, the backbones of 4-OT, ACBP, and PLCγ1C SH2 show limited, increased, and decreased order, respectively; the cores of DHFR-E:folate and PLCγ1C SH2 become more ordered. The V18A (V34A) mutation increases (decreases) the order within the eglin C core. The core of TNfn3 is less ordered structurally and more mobile kinetically. Secondary structure versus loops, surface-binding versus core insertion, and ligand size emerged as being important in rationalizing ΔŜ. The consistent and general tool developed herein is expected to provide further insights in future work.

An SRLS Study of 2H Methyl-Moiety Relaxation and Related Conformational Entropy in Free and Peptide-Bound PLCγ1C SH2.

The two-body (protein and probe) coupled-rotator slowly relaxing local structure (SRLS) approach for NMR relaxation in proteins is extended to derive conformational entropy, Ŝ. This version of SRLS is applied to deuterium relaxation from the C-CDH2 bonds of free and peptide-bound PLCγ1C SH2. Local C-CDH2 motion is described by a correlation time for local diffusion, τ2, and a Maier-Saupe potential, u. On average, τ2, which largely fulfills τ2 ≪ τ1 (τ1 - correlation time for global tumbling), is 270 ± 41 ps and u is 2 ± 0.1 kBT. The PLCγ1C SH2 data were analyzed previously with the model-free (MF) method. SRLS is a generalization of MF, assumed so far to yield the latter for τ2 ≪ τ1 and simple local geometry. Despite these conditions being fulfilled, we find here that τ2 and u differ substantially from their MF counterparts. This is shown to stem from MF (a) disregarding mode-coupling of the first type (see below) and (b) parametrizing the methyl-moiety-related spectral density function (SDF). Our main interest lies in ΔŜ, the conformational entropy difference between the peptide-bound and free PLCγ1C SH2 forms. We find that ΔŜ is rendered inaccurate in MF because factors a and b above impair the accuracy of Saxis, the parameter on which the calculation of ΔŜ is based. Conformational entropy was obtained previously using various simple system-specific models. SRLS is unique in obtaining this important thermodynamic quantity based on a general physically well-defined local potential. It is also unique in its ability to extract the information inherent in 2H relaxation parameters from methyl moieties in protein with accuracy commensurate with data sensitivity.

Polar Versus Non-polar Local Ordering at Mobile Sites in Proteins: Slowly Relaxing Local Structure Analysis of (15)N Relaxation in the Third Immunoglobulin-Binding Domain of Streptococcal Protein G.

We developed recently the slowly relaxing local structure (SRLS) approach for studying restricted motions in proteins by NMR. The spatial restrictions have been described by potentials comprising the traditional L = 2, K = 0, 2 spherical harmonics. However, the latter are associated with non-polar ordering whereas protein-anchored probes experience polar ordering, described by odd-L spherical harmonics. Here we extend the SRLS potential to include the L = 1, K = 0, 1 spherical harmonics and analyze (15)N-(1)H relaxation from the third immunoglobulin-binding domain of streptococcal protein G (GB3) with the polar L = 1 potential (coefficients c0(1) and c1(1)) or the non-polar L = 2 potential (coefficients c0(2) and c2(2)). Strong potentials, with ⟨c0(1)⟩ ∼ 60 for L = 1 and ⟨c0(2)⟩ ∼ 20 for L = 2 (in units of kBT), are detected. In the α-helix of GB3 the coefficients of the rhombic terms are c1(1) ∼ c2(2) ∼ 0; in the preceding (following) chain segment they are ⟨c1(1)⟩ ∼ 6 for L = 1 and ⟨c2(2)⟩ ∼ 14 for L = 2 (⟨c1(1)⟩ ∼ 3 for L = 1 and ⟨c2(2)⟩ ∼ 7 for L = 2). The local diffusion rate, D2, lies in the 5 × 10(9)-1 × 10(11) s(-1) range; it is generally larger for L = 1. The main ordering axis deviates moderately from the N-H bond. Corresponding L = 1 and L = 2 potentials and probability density functions are illustrated for residues A26 of the α-helix, Y3 of the β1-strand, and L12 of the β1/β2 loop; they differ considerably. Polar/orientational ordering is shown to be associated with GB3 binding to its cognate Fab fragment. The polarity of the local ordering is clearly an important factor.

Structural Dynamics of the Potassium Channel Blocker ShK: SRLS Analysis of (15)N Relaxation.

The 35-residue ShK peptide binds with high affinity to voltage-gated potassium channels. The dynamics of the binding surface was studied recently with (microsecond to millisecond) (15)N relaxation dispersion and (picosecond to nanosecond) (15)N spin relaxation of the N-H bonds. Relaxation dispersion revealed microsecond conformational-exchange-mediated exposure of the functionally important Y23 side chain to the peptide surface. The spin relaxation parameters acquired at 14.1 and 16.45 T have been subjected to model-free (MF) analysis, which yielded a squared generalized order parameter, S(2), of approximately 0.85 for virtually all of the N-H bonds. Only a "rigid backbone" evaluation could be inferred. We ascribe this limited information to the simplicity of MF in the context of challenging data. To improve the analysis, we apply the slowly relaxing local structure (SRLS) approach, which is a generalization of MF. SRLS describes N-H bond dynamics in ShK in terms of a local potential, u, ranging from 10 to 18.5 kBT, and a local diffusion rate, D2, ranging from 4.2 × 10(8) to 2.4 × 10(10) s(-1). This analysis shows that u is outstandingly strong for Y23 and relatively weak for K22, whereas D2 is slow for Y23 and fast for K22. These observations are relevant functionally because of the key role of the K22-Y23 dyad in ShK binding to potassium channels. The disulfide-bond network exhibits a medium-strength potential and an alternating wave-like D2 pattern. This is indicative of moderate structural restraints and motional plasticity, in support of, although not directly correlated with, the microsecond binding-related conformational exchange process detected previously. Thus, new information on functionally important residues in ShK and its overall conformational stability emerged from the SRLS analysis, as compared with the previous MF-based estimate of backbone dynamics as backbone rigidity.

NMR spin relaxation in proteins: The patterns of motion that dissipate power to the bath

We developed in recent years the two-body coupled-rotator slowly relaxing local structure (SRLS) approach for the analysis of NMR relaxation in proteins. The two bodies/rotators are the protein (diffusion tensor D1) and the spin-bearing probe, e.g., the 15N−1H bond (diffusion tensor, D2), coupled by a local potential (u). A Smoluchowski equation is solved to yield the generic time correlation functions (TCFs), which are sums of weighted exponentials (eigenmodes). By Fourier transformation one obtains the generic spectral density functions (SDFs) which underlie the experimental relaxation parameters. The typical paradigm is to characterize structural dynamics in terms of the best-fit values of D1, D2, and u. Additional approaches we pursued employ the SRLS TCFs, SDFs, or eigenmodes as descriptors. In this study we develop yet another perspective. We consider the SDF as function of the angular velocity associated with the fluctuating fields underlying NMR relaxation. A parameter called j-fraction, which represents the relative contribution of eigenmode, i, to a given value of the SDF function at a specific frequency, ω, is defined. j-fraction profiles of the dominant eigenmodes are derived for 0 ≤ ω ≤ 1012 rad/s. They reveal which patterns of motion actuate power dissipation at given ω-values, what are their rates, and what is their relative contribution. Simulations are carried out to determine the effect of timescale separation, D1/D2, axial potential strength, and local diffusion axiality. For D1/D2 ≤ 0.01 and strong local potential of 15 kBT, power is dissipated by global diffusion, renormalized (by the strong potential) local diffusion, and probe diffusion on the surface of a cone (to be called cone diffusion). For D1/D2 = 0.1, power is dissipated by mixed eigenmodes largely of a global-diffusion-type or cone-diffusion-type, and a nearly bare renormalized-local-diffusion eigenmode. For D1/D2 > 0.1, most eigenmodes are of a mixed type. The analysis is affected substantially by reducing the potential strength from 15 to 5 kBT, and/or allowing for axial D2 with D2,∥/D2,⊥ = 10. The scheme developed is applied to 15N−1H relaxation from the β-sheet residue K19 and the α-helix residue A34 of the third immunoglobulin-binding domain of streptococcal protein G. Previous studies revealed rhombic local potentials with different rhombicity around ${\rm C}_{i - 1}^\alpha {\rm - C}_i^\alpha $Ci−1α−Ciα, and different timescale separation (0.047 for K19 and 0.102 for A34). Here, we find that K19 and A34 dissipate power to the bath through global diffusion, mixed cone-diffusion-related and mixed renormalized-local-diffusion-related motions. At small ω-values, A34 is more effective than K19 in dissipating power. In general, it executes faster cone-diffusion-type, and slower renormalized-local-diffusion-type and local-probe-fluctuation-type motions. K19 experiences faster N−H fluctuations than A34. Eigenmode clustering, experienced by K19 to a larger extent, is observed in the fast-probe-fluctuation regime. New information on the effect of the structural context on N−H bond dynamics has been obtained. The patterns of motion that dissipate NMR-relaxation-related power illuminate protein dynamics from a new perspective. They constitute yet another qualifier of N−H bond dynamics. This study sets the stage for developing ways for enhancing the contribution of desired pathways for power dissipation at selected angular velocities.

The Slowly Relaxing Local Structure Perspective of Protein Dynamics by NMR Relaxation

AbstractNMR relaxation is a powerful method for elucidating structural dynamics. Standard stochastic dynamic models generate time correlation functions (TCFs) that feature physically well‐defined parameters. We developed such a model, called the slowly relaxing local structure (SRLS) approach, for proteins. SRLS is a two‐body (protein and probe) coupled‐rotator approach. Given that the protein (featuring diffusion tensor, D1) restricts the probe (featuring diffusion tensor, D2), the two “bodies” are inherently coupled dynamically. This is substantiated by a local potential, u, associated with a local ordering tensor, S. SRLS allows for general tensorial properties of D1, D2, S and the magnetic NMR tensors, and a general form of u. The TCFs are multi‐exponential, in accordance with the degree of generality of the various tensors. The traditional model‐free (MF) method is based on a different conceptualization. According to it a mode‐decoupling bi‐exponential (one term for each rotator) TCF captures adequately the detectable features of structural dynamics. Hence, stochastic approaches are unnecessary. Here, we show that this (amply proven) oversimplification leads to physically vague constructs/composites as descriptors of structural dynamics. We illustrate misleading results obtained with MF when mode coupling, or S tensor asymmetry, dominate the analysis. Finally, we delineate the substantial advantage in using SRLS TCF as quantity to be compared with its atomistic molecular dynamicsbased counterpart.

The eigenmode perspective of NMR spin relaxation in proteins

We developed in recent years the two-body (protein and probe) coupled-rotator slowly relaxing local structure (SRLS) approach for elucidating protein dynamics from NMR spin relaxation. So far we used as descriptors the set of physical parameters that enter the SRLS model. They include the global (protein-related) diffusion tensor, D1, the local (probe-related) diffusion tensor, D2, and the local coupling∕ordering potential, u. As common in analyzes based on mesoscopic dynamic models, these parameters have been determined with data-fitting techniques. In this study, we describe structural dynamics in terms of the eigenmodes comprising the SRLS time correlation functions (TCFs) generated by using the best-fit parameters as input to the Smoluchowski equation. An eigenmode is a weighted exponential with decay constant given by an eigenvalue of the Smoluchowski operator, and weighting factor determined by the corresponding eigenvector. Obviously, both quantities depend on the SRLS parameters as determined by the SRLS model. Unlike the set of best-fit parameters, the eigenmodes represent patterns of motion of the probe-protein system. The following new information is obtained for the typical probe, the (15)N-(1)H bond. Two eigenmodes, associated with the protein and the probe, dominate when the time scale separation is large (i.e., D2 >> D1), the tensorial properties are simple, and the local potential is either very strong or very weak. When the potential exceeds these limits while the remaining conditions are preserved, new eigenmodes arise. The multi-exponentiality of the TCFs is associated in this case with the restricted nature of the local motion. When the time scale separation is no longer large, the rotational degrees of freedom of the protein and the probe become statistically dependent (coupled dynamically). The multi-exponentiality of the TCFs is associated in this case with the restricted nature of both the local and the global motion. The effects of local diffusion axiality, potential strength, and extent of mode-coupling on the eigenmode setup are investigated. We detect largely global motional or largely local motional eigenmodes. In addition, we detect mixed eigenmodes associated with correlated∕prograde or anti-correlated∕retrograde rotations of the global (D1) and local (D2) motional modes. The eigenmode paradigm is applied to N-H bond dynamics in the β-sheet residue K19, and the α-helix residue A34, of the third immunoglobulin-binding domain of streptococcal protein G. The largest contribution to the SRLS TCFs is made by mixed anti-correlated D1 and D2 eigenmodes. The next largest contribution is made by D1-dominated eigenmodes. Eigenmodes dominated by the local motion contribute appreciably to A34 and marginally to K19. Correlated D1 and D2 eigenmodes contribute exclusively to K19 and do not contribute above 1% to A34. The differences between K19 and A34 are delineated and rationalized in terms of the best-fit SRLS parameters and mode-mixing. It may be concluded that eigenmode analysis is complementary and supplementary to data-fitting-based analysis.

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.

SRLS Analysis of 15N Relaxation from Bacteriophage T4 Lysozyme: a Tensorial Perspective That Features Domain Motion

Bacteriophage T4L lysozyme (T4L) comprises two domains connected by a helical linker. Several methods detected ns domain motion associated with the binding of the peptidoglycan substrate. An ESR study of nitroxide-labeled T4L, based on the slowly relaxing local structure (SRLS) approach, detected ns local motion involving the nitroxide and the helix housing it. (15)N−H spin relaxation data from T4L acquired at magnetic fields of 11.7 and 18.8 T, and 298 K, were analyzed previously with the model-free (MF) method. The results did not detect domain motion. SRLS is the generalization of MF. Here, we apply it to the same data analyzed previously with MF. The restricted local N−H motion is described in terms of tilted axial local ordering (S) and local diffusion (D(2)) tensors; dynamical coupling to the global tumbling is accounted for. We find that D(2,⊥) is 1.62 × 10(7) (1.56 × 10(7)) s(−1) for the N-terminal (C-terminal) domain. This dynamic mode represents domain motion. For the linker D(2,⊥) is the same as the rate of global tumbling, given by (1.46 ± 0.04) × 10(7) s(−1). D(2,∥) is 1.3 × 10(9), 1.8 × 10(9) and 5.3 × 10(9) s(−1) for the N-terminal domain, the C-terminal domain, and the linker, respectively. This dynamic mode represents N−H bond vector fluctuations. The principal axis of D(2) is virtually parallel to the N−H bond. The order parameter, S(0)(2), is 0.910 ± 0.046 for most N−H bonds. The principal axis of S is tilted from the C(i−1)(α) −C(i)(α) axis by −2° to 6° for the N-, and C-terminal domains, and by 2.5° for the linker. The tensorial-perspective-based and mode-coupling-based SRLS picture provides new insights into the structural dynamics of bacteriophage T4 lysozyme.

Slowly Relaxing Local Structure (SRLS) Analysis of 15N–H Relaxation from the Prototypical Small Proteins GB1 and GB3

15N-H relaxation parameters from the first (GB1) and third (GB3) immunoglobulin-binding domains of streptococcal protein G were analyzed previously with the traditional model-free (MF) method. These proteins comprise an α-helix and a four-stranded β-sheet. An extensive study of GB1 (GB3) used combined three-field (five-field) data acquired in the 278-323 K range (at 297 K). For successful analysis of the GB3 data, it was necessary to allow for variations in the 15N chemical shift anisotropy (CSA) tensor and virtually eliminate the local motion. In the case of GB1, the spectral density was parametrized. Here, we analyze these data with the slowly relaxing local structure (SRLS) approach, which is the generalization of MF in allowing for general tensorial properties, and accounting for mode-coupling. A standard (featuring constant magnetic tensors) SRLS fitting scheme is used. This analysis accounts for the important asymmetry of the local spatial restrictions; it provides physical order parameters, local diffusion rates, related activation energies, and key features of local geometry. Using data from GB3 we show that the main local ordering axis is C(i-1)(α) - C(i)(α), and the average axial (rhombic) order parameter is -0.457 ± 0.017 (1.156 ± 0.015) for the α-helix and -0.484 ± 0.002 (1.10 ± 0.04) for the rest of the polypeptide chain. The N-H bonds within (outside of) the α-helix reorient locally with an average correlation time, (τ), of 310 (130) ps, as compared to 3.33 ns for the global tumbling. Several N-H bonds in the loops β1/β2, β2/α-helix, and α-helix/β3 have (τ) of 380, 320, and 750 ps, respectively. The distinctive experimental data of the α-helix are due to relatively weak and substantially rhombic local ordering and slow local motion. For GB1, we derive activation energies from local diffusion rates. They are 43.3 ± 7.1 kJ/mol for the β-strands, 24.7 ± 3.9 kJ/mol for the α-helix (and approximately for the loop β3/β4), and 18.9 ± 1.8 kJ/mol for the other loops. The physical SRLS description provides new insights into the backbone dynamics of GB1 and GB3 in particular, and proteins in general.

SRLS Analysis of 15N Spin Relaxation from E. coli Ribonuclease HI: The Tensorial Perspective

15N–H relaxation parameters from ribonuclease HI (RNase H), acquired in previous work at magnetic fields of 14.1 and 18.8 T, and at 300 K, are analyzed with the mode-coupling slowly relaxing local structure (SRLS) approach. In accordance with standard theoretical treatments of restricted motions, SRLS approaches N-H bond dynamics from a tensorial perspective. As shown previously, a physically adequate description of this phenomenon has to account for the asymmetry of the local spatial restrictions. So far, we used rhombic local ordering tensors; this is straightforward but computationally demanding. Here, we propose substantiating the asymmetry of the local spatial restrictions in terms of tilted axial local ordering (S) and local diffusion (D2) tensors. Although less straightforward, this description provides physically sound structural and dynamic information and is efficient computationally. We find that the local order parameter, S(0)2, is on average 0.89 (0.84, and may be as small as 0.6) for the secondary structure elements (loops). The main local ordering axis deviates from the C(i-1)α-C(i)α axis by less than 6°. At 300 K, D(2,perpendicular) is virtually the same as the global diffusion rate, D1 = 1.8 × 10(7) s(-1). The correlation time 1/6D(2,parallel) ranges from 3-125 (208-344) ps for the secondary structure elements (loops) and is on average 125 ps for the C-terminal segment. The main local diffusion axis deviates from the N-H bond by less than 2° (10°) for the secondary structure elements (loops). An effective data-fitting protocol, which leads in most cases to unambiguous results with limited uncertainty, has been devised. A physically sound and computationally effective methodology for analyzing 15N relaxation in proteins, that provides a new picture of N–H bond structural dynamics in proteins, has been set forth.

Methyl dynamics of a Ca2+-calmodulin-peptide complex from NMR/SRLS.

We developed the slowly relaxing local structure (SRLS) approach for analyzing NMR spin relaxation in proteins. SRLS accounts for dynamical coupling between the tumbling of the protein and the local motion of the probe and for general tensorial properties. It is the generalization of the traditional model-free (MF) method, which does not account for mode-coupling and treats only simple tensorial properties. SRLS is applied herein to ²H relaxation of ¹³CDH₂ groups in the complex of Ca(2+)-calmodulin with the peptide smMLCKp. Literature data comprising ²H T₁ and T₂ acquired at 14.1 and 17.6 T, and 288, 295, 308, and 320 K, are used. We find that mode-coupling is a small effect for methyl dynamics. On the other hand, general tensorial properties are important. In particular, it is important to allow for the asymmetry of the local spatial restrictions, which can be represented in SRLS by a rhombic local ordering tensor with components S(0)(2) and S(2)(2). The principal axes frame of this tensor is obviously different from the axial frames of the magnetic tensors. Here, we find that -0.2 ≤ S(0)(2) ≤ 0.5 and -0.4 ≤ S(2)(2) ≤ 0. MF features a single "generalized" order parameter, S, confined to the 0-0.316 range; the local geometry is inherently simple. The parameter S is inaccurate, having absorbed unaccounted for effects, notably S(2)(2) ≠ 0. We find that the methionine methyls (the other methyl types) reorient with rates of 8.6 × 10⁹ to 21.4 × 10⁹ (0.67 × 10⁹ to 6.5 × 10⁹) 1/s. The corresponding activation energies are 10 (10-27) kJ/mol. By contrast, MF yields inaccurate effective local motional correlation times, τ(e), with nonphysical temperature dependence. Thus, the problematic S- and τ(e)-based MF picture of methyl dynamics has been replaced with an insightful physical picture based on a local ordering tensor related to structural features, and a local diffusion tensor that yields accurate activation energies.

Backbone dynamics of deoxy and carbonmonoxy hemoglobin by NMR/SRLS.

The slowly relaxing local structure (SRLS) approach, developed for NMR spin relaxation analysis in proteins, is applied herein to amide ¹⁵N relaxation in deoxy and carbonmonoxy hemoglobin. Experimental data including ¹⁵N T₁, T₂ and ¹⁵N-{¹H} NOE, acquired at 11.7 and 14.1 T, and 29 and 34 °C, are analyzed. The restricted local motion of the N-H bond is described in terms of the principal value (S(0)(2)) and orientation (β(D)) of an axial local ordering tensor, S, and the principal values (R(||)(L) and R(⊥)(L)) and orientation (β(O)) of an axial local diffusion tensor, R(L). The parameters c₀² (the potential coefficient in terms of which S(0)(2) is defined), R(||)(L), β(D), and β(O) are determined by data fitting; R(⊥)(L) is set equal to the global motional rate, R(C), found previously to be (5.2-5.8) × 10⁶ 1/s in the temperature range investigated. The principal axis of S is (nearly) parallel to the C(i-1)(α)-C(i)(α) axis; when the two axes are parallel, β(D) = -101.3° (in the frame used). The principal axis of R(L) is (nearly) parallel to the N-H bond; when the two axes are parallel, β(O) = -101.3°. For "rigid" N-H bonds located in secondary structure elements the best-fit parameters are S(0)(2) = 0.88-0.95 (corresponding to local potentials of 8.6-19.9 k(B)T), R(||)(L) = 10⁹-10¹⁰ 1/s, β(D) = -101.3° ± 2.0°, and β(O) = -101.3° ± 4°. For flexible N-H bonds located in loops the best-fit values are S(0)(2) = 0.75-0.80 (corresponding to local potentials of 4.5-5.5 k(B)T), R(||)(L) = (1.0-6.3) × 10⁸ 1/s, β(D) = -101.3° ± 4.0°, and β(O) = -101.3° ± 10°. These results are important in view of their physical clarity, inherent potential for further interpretation, consistency, and new qualitative insights provided (vide infra).