PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
J Phys Chem B. Author manuscript; available in PMC 2017 March 24.
Published in final edited form as:
PMCID: PMC4837759
NIHMSID: NIHMS768045

Local Ordering at Mobile Sites in Proteins from Nuclear Magnetic Resonance Relaxation: the Role of Site Symmetry

Abstract

Restricted motions in proteins (e.g., N–H bonds dynamics) are studied effectively with NMR. By analogy with restricted motions in liquid crystals (LC), the local ordering has in the past been primarily represented by potentials comprising the L = 2, |K| = 0, 2 spherical harmonics. However, probes dissolved in LC’s experience non-polar ordering, often referred to as alignment, while protein-anchored probes experience polar ordering, often referred to as orientation. In this study we investigate the role of local (site) symmetry in the context of the polarity of the local ordering. We find that potentials comprising the L = 1, |K| = 0, 1 spherical harmonics represent adequately polar ordering. It is useful to characterize potential symmetry in terms of the irreducible representations of D2h point group, which is already implicit in the definition of the rotational diffusion tensor. Thus, the relevant rhombic L = 1 potentials have B1u and B3u symmetry whereas the relevant rhombic L = 2 potentials have Ag symmetry. A comprehensive scheme where local potentials and corresponding probability density functions (PDFs) are represented in Cartesian and spherical coordinates clarifies how they are affected by polar and non-polar ordering. The Cartesian coordinates are chosen so that the principal axis of polar axial PDF is pointing along the z-axis, whereas the principal axis of the non-polar axial PDF is pointing along ±z. Two-term axial potentials with 1 ≤ L ≤ 3 exhibit substantial diversity; they are expected to be useful in NMR-relaxation-data-fitting. It is shown how potential coefficients are reflected in the experimental order parameters. The comprehensive scheme representing local potentials and PDFs is exemplified for the L = 2 case using experimental data from 15N-labeled plexin-B1 and thioredoxin, 2H-, and 13C-labeled benzenehexa-n-alkanoates, and nitroxide-labeled T4 lysozyme. Future prospects for improved ordering analysis based on combined atomistic and mesoscopic approaches are delineated.

Keywords: restricted motions, symmetry of local ordering, alignment versus orientation, spherical harmonics

TOC Figure

An external file that holds a picture, illustration, etc.
Object name is nihms768045u1.jpg

1. Introduction

Molecules may be characterized by 3D structures, and dynamics which are often associated with motional rates. We consider, for example, a rigid molecule reorienting in an isotropic medium.13 Such motion may be represented by a 3D diffusion tensor,3 whose principal values are a function of molecular shape. However, when the molecule is reorienting in an anisotropic medium, or parts of a flexible molecule are reorienting in anisotropic local surroundings, the description of the motion also has to include proper representation of the local spatial constraints.47 The latter can be expressed in terms of orienting potentials,47 which represent the local structures. Orienting potentials are associated with the preferential arrangement of the molecule in the surrounding local space, and thus reflect the local geometric arrangement. In other words, structure-related and geometry-related kinetic features can be elucidated by studying motions that take place in anisotropic environments.419

A class of restricted motions of particular interest are intra-molecular reorientation processes of segments of proteins (more generally, bio-macromolecules).1019 Examples include small moieties such as amide or methyl groups1319 that execute localized motions, as well as large loops or domains that execute collective motions.1719 These processes are associated with the internal mobility/flexibility of the protein. In this study we focus on structure-related and geometry-related aspects of the local ordering in proteins in terms of their respective symmetries.

Restricted motions are often treated within the scope of atomistic molecular dynamics (MD) simulations. Usually this is done by taking into account sets of relatively slow coordinates.2025 Using effective sampling and appropriate analysis of MD trajectories, probability density functions (PDFs) yielding approximate potentials of mean force (PMFs) have been determined (e.g., ref 25). To our knowledge, the more mesoscopic (collective) local features, such as explicit forms of local anisotropic potentials and local molecular geometry associated with these potentials, as well as effective diffusion tensors, have not been provided by MD methods.

On the other hand, mesoscopic models for treating restricted motions are convenient for interpreting magnetic resonance experiments.419 They encompass order parameters and PDFs in a straightforward manner.6,7 These stochastic approaches utilize a tensorial representation of diffusion and ordering. Using NMR or ESR spectroscopy as reporters of the structural dynamics, one is able to determine the structure-related principal values and the geometry-related orientations of the ordering and diffusion tensors.

The orienting potential, U(Ω), is typically given in mesoscopic stochastic models by the expansion in the Wigner rotation matrix elements, DMKL(Ω):

u(Ω)U(Ω)/kBT=L=1M=L+LK=L+LcMKLDMKL(Ω)
(1)

with u(Ω) and the coefficients cM,KL being dimensionless. The Euler angles, Ω, describe the molecular orientation relative to the director frame associated with the orienting surroundings. This general expansion needs to be truncated at some value L = Lmax before attempting to reproduce experimental data, thus leaving just a few leading terms. However, the “quantum numbers” M and K are just restricted to |M| ≤ L and |K| ≤ L, which still yields too many cMKL fitting parameters for the typically limited experimental data available. We thus need to make sensible choices whereby a sufficiently reduced number of parameters are used in the fitting process.

Let us first take a hint from the symmetry of the rotational diffusion tensor, which is a second rank symmetric tensor consistent with D2h point group symmetry. That is, however complex the shape of a molecule, we are approximating its diffusive properties in terms of D2h point group symmetry. Thus, it seems reasonable that we impose similar effective point group symmetry on the orientational potential of eq 1. Let us examine what has been done in this context in the past.

Stochastic models for studying restricted motions of rigid molecules dissolved in liquid crystalline (LC) media by ESR and NMR have been developed in early work.46 The usual assumption is the existence of a uniaxial LC director. Consequently the “quantum number”, M, in eq 1 is zero, yielding:47

u(θ,φ)=L=1K=L+LcKLD0KL(0,θ,φ)
(2)

In this case only the two polar angles, θ and [var phi], are needed to describe the molecular orientation relative to the director frame. This is a typical simplification, in that it requires just the D0KL(0,θ,φ) functions, which are proportional to the corresponding spherical harmonics, YLK(θ,φ) (e.g., see ref 26).

In a simple sense, the ordering of a large solute probe reorienting in a macroscopically aligned liquid crystal (LC) solvent which exhibits inversion symmetry with respect to the origin of the director frame may be taken as reflecting this inversion symmetry. Such a scenario represents “non-polar” ordering. Potentials describing non-polar ordering comprise only the even-L terms of eq 2. This may be thought of as preferential arrangement in space of double-headed-arrows, a setup often referred to as alignment27,28 The leading terms here are the functions D0K2 with |K| = 0, 1, 2 (more generally, L = 2, 4, 6…).

Envision a smaller probe residing in one leaflet of a lipid bilayer. The local environment is most certainly polar, requiring odd-L terms to be retained in eq 2. A similar requirement would apply to a probe dissolved in a thermotropic liquid crystal formed of polar molecules (unless the probe is much larger than the liquid crystal molecule). One may think of polar ordering as preferential arrangement in space of single-headed arrows, which is often referred to as orientation.27,28

Now consider a probe reorienting in the anisotropic internal space of a protein while anchored at it. The fact that the probe is bonded to the protein implies a single-headed-arrow-type, i.e., polar, ordering. We represent this scenario by taking the local director to be uniaxial and polar. In this case eq 2, with odd-L terms retained, should be used. The leading terms here are D0K1 with |K| = 0, 1 (more generally, L = 1, 3, 5…).

Comparison of the polar and non-polar scenarios then clarifies the distinction between orientation (visualized as a vector) and alignment (visualized as a symmetric tensor), as we have defined them above. The respective potentials may then be classified in terms of the irreducible representations of D2h point group.

The comparison between the polar and non-polar cases and their admixture is carried out within the scope of a comprehensive scheme which includes potential shapes, and forms of the associated PDFs, represented in both Cartesian-coordinate and spherical-coordinate systems. For axial probes (K = 0), the terms with L = 3 and L = 4 are also considered.

Unlike the LC director which is fixed in space (although fluctuations do occur), a local director such as prevailing in internally mobile proteins is itself tumbling in solution, being defined in terms of the protein axes. This represents a two-body (protein and probe) mesoscopic coupled-rotator problem. Freed et al. solved this problem within the scope of the slowly relaxing local structure (SRLS) approach.8 SRLS, and its limit where the protein motion is frozen, the microscopic-order-macroscopic-disorder (MOMD) approach,29 have been applied extensively within the scope of ESR spectroscopy.1012 Applications to model and biological membranes, complex fluids, nitroxide-labeled proteins and DNA fragments, etc., have been reported.30,31 We have also applied SRLS extensively to NMR relaxation in proteins dissolved in aqueous solution.1319 Recently we have applied MOMD to the analysis of 2H NMR lineshapes from polycrystalline proteins experiencing internal mobility.32,33 All of these approaches use potentials expanded in the basis set of the D0KL(0,θ,φ) functions. Thus, the ordering-related analysis in this work is relevant to SRLS, MOMD, and any other stochastic model for analyzing NMR and ESR lineshapes and/or relaxation parameters in proteins utilizing the D0KL(0,θ,φ) functions (i.e., the spherical harmonics) as basis-sets for the expansion of the local potential. We also describe ties to atomistic MD methods, as described below.

In previous SRLS-based and MOMD-based studies the potential typically included the leading even-L (L = 2) terms.9,10,1419,29,32,33 In a few cases the L = 4, K = 0, 2 terms were also included in the analysis.9 A potential given by the L = 1, K = 0 term has been considered by Polimeno and Freed in ref 8. It was found that using potentials given by the L = 1, K = 0 term or the L = 2, K = 0 term yields results for motional correlation functions that are similar in nature (although the potential coefficients determined have different magnitudes) when, as in ESR or NMR, the effect of the motion consists of averaging second-rank magnetic tensors.8 In this study we also take into account the leading odd-L (L = 1) terms; this constitutes a substantial enhancement (currently within the scope of designing the local potential) of the SRLS and MOMD potentials.

Given that order parameters rather than potential coefficients are typically determined experimentally, we show how the latter are reflected by the former. The comprehensive potential/PDF/order parameter description delineated above is exemplified for L = 2 potentials using previously-determined order parameters or potential coefficients from 15N-labeled plexin-B1 in solution,25,34 T4 lysozyme in solution,35 polycrystalline thioredoxin,36 and discotic and solid phases of benzenehexa-n-alkanoates (with n = 6 and 7).37

Our results are expected to be relevant to a broad range of experiments. Liquid-crystal phases of phospholipids require odd-L potentials for proper analysis of 31P spin relaxation.38 Rhombic polar order-parameters have been invoked in the context of twist-bend nematic phases.39 However, they have been associated with the conical angle and the pitch of the medium, rather than potentials formed of odd-L terms of eq 2. Potentials given by the L = 1, K = 0 term have been used to analyze dielectric relaxation of ferroelectric systems,40 to treat electric-field-implied polar order,41 and interpret the orientation of asymmetric paramagnetic molecules in strong magnetic fields. To our knowledge, applications where the L = 1, |K| = 1 terms, or both L = 1 and L = 2 terms, were also included in the expression for the potential, as required for nonaxial cases, have not been reported.

2. Theoretical Background

As pointed out above, we take the local director to be uniaxial. Next we have to ensure that the potential is real; this requires that cKL=()K(cKL) (ref 9). The potential which is real is given by:

u(θ,φ)=L{c0LD00L(0,θ,0)+K>0LcK±L[D0KL(0,θ,φ)+()jD0KL(0,θ,φ)]}
(3)

where j = 1 for odd-L and j = 0 for even-L, and cK±LcKL±cKL.

Let us focus on the L = 2 potential and let us assume for convenience, as is usually done, that the local ordering tensor defined in terms of this potential is diagonal in the same molecular frame as the diffusion tensor. The coefficients c02,c2±2c22±c22 are themselves irreducible tensor components. The Cartesian tensor with components cij2=cji2 is diagonal (i.e., Tr(cii2)=0) in the same frame as the local ordering tensor,42,43 with complete specification given by c02 and c2+2 (for simplicity, the plus sign in c2+2 will be omitted below). The corresponding D002 and (D022+D022) terms in eq 3 have symmetry Ag in D2h point group.

Nest we focus on the L = 1 potential. Three real potentials, given by the spherical harmonic functions Y10, 12(Y11Y11) and i12(Y11+Y11), may be defined. The respective symmetries in the D2h point group are B1u, B3u, and B2u. These functions are similar in character to the atomic orbitals pz, px, and py, respectively.26 When the primary polar axis is taken as the z-axis of the diffusion tensor only the term Y10 survives in eq 3. When Y10 is linearly combined with 12(Y11Y11)(i12(Y11+Y11)) with positive coefficients, the primary polar axis is tilted in the zx (zy) plane (negative coefficients change the sign/signs of the x, y and/or z coordinates). We choose Y10 linearly combined with 12(Y11Y11) as polar orientational potential (i.e., we assume that in the rhombic case the principal polar axis is tilted in the |zx| plane rather than the |zy| plane). The plus sign in the designation of the coefficient c1+1 of 12(Y11Y11) will be omitted below.

The relevant order parameters are thus:

S0L=D00L(0,θ,0),L=1,2,3,4
(4a)

S11=(D011(0,θ,φ)D011(0,θ,φ))
(4b)

S22=(D022(0,θ,φ)+D022(0,θ,φ))
(4c)

The normalized probability density function is given by:

D0KL(0,θ,φ)eu(θ,φ)sinθdθdφeu(θ,φ)sinθdθdφ
(5)

The part of the Character Table of D2h point group relevant to this study is given in Table S1 of the Supporting Information.26 It depicts the symmetry operations that characterize the various potential functions considered.

As pointed out above, for probes with axial symmetry the functions D001(0,θ,0), D002(0,θ,0), D003(0,θ,0), and D004(0,θ,0), proportional to Y10, Y20, Y30, Y40, respectively, are also considered; they are given in Table S2 of the Supporting Information. The odd-L functions have B1u symmetry whereas the even-L functions have Ag symmetry (Tables S1 and S2 of the Supporting Information).

The real spherical harmonics Y10, 12(Y11Y11), i12(Y11+Y11), Y20 and 12(Y22+Y22) are given in Table S3 of the Supporting Information. Note that Y20 and 12(Y22+Y22) are similar in character to the atomic orbitals dz2 and dx2y2, respectively.26 These symmetries are also helpful in determining the non-zero matrix elements of the potential terms in the Stochastic Liouville Equation (SLE) developed within the scope of standard theoretical methods.

The polar and non-polar scenarios are compared in the following manner. 3D potential surfaces are generated as a function of the spherical coordinates, θ and [var phi]. In this representation the value of the potential, u(θ, [var phi]), is depicted on the z-axis, and color-coding is used to illustrate variations. 3D surfaces of u are also shown in Cartesian coordinates, x = sin θ cos [var phi]|u|, y = sin θ sin [var phi] |u| and z = cos θ |u|. In this representation the value of the potential is given by the distance from the origin of the coordinate frame to any point on the 3D surface. Positive and negative values are colored red and blue, respectively, by analogy with positive and negative phases of the p-orbitals for the L = 1 case, and the dz2 and dx2y2 orbitals for the L = 2 case. Note, however, that “blue” potential values are populated preferentially. The utility of the Cartesian representation of the potential will be realized in cases of compounded potential forms.

PDFs are shown in both Cartesian44 (x = sin θ cos [var phi] exp(−u), y = sin θ sin [var phi] exp (−u) and z = cos θ exp (−u)) and spherical coordinates (θ and [var phi]). Because of the PDF exponential functional dependence on u, large positive potential values hardly contribute to the visual display of the PDF. It is useful to examine simultaneously the Cartesian representations of corresponding potential and PDF pairs. The information provided by the four representations is thus complementary.

3. Results and Discussion

3.1. Axial Ordering (K = 0)

Order Parameters

Let us consider axially symmetric molecules described by potentials comprising a single term. The L = 1 and K = 0 potential is u(θ)=c01cosθ and the corresponding order parameter, S01, is given by eq 4a. The L = 2 and K = 0 potential is u(θ)=c02[32cos2θ12] and the order parameter, S02, is also given by eq 4a, etc.

Thus, in these simple cases the order parameter, S0L, is obtained from the corresponding potential coefficient, c0L, in a straightforward manner by solving eq 4a. However, in many cases order parameters rather than potential coefficients are obtained experimentally.36,37 It is of interest to derive c0L from S0L, given that potential forms are substantially more sensitive to details of the local ordering than order parameters (see below).

This can be accomplished by pre-calculating graphs of S0L as a function of c0L, and using them subsequently to derive c0L from known S0L values. Figure 1 shows such graphs for L = 1, 2, 3, 4, with c0L ranging from −30 to 30. The trigonometric expressions underlying the potentials for L = 1–4 are given in Table S2 of the Supporting Information.

Figure 1
Order parameters, S0L, shown as a function of potential coefficients, c0L, for L = 1, 2, 3, 4 and K = 0, calculated according to eqs 3 and 4. The respective functional forms are given in Table S2 of the Supporting Information.

For any given positive value of c0L the extent of ordering increases along the series S04, S03, S02 and S01, with S01S02 somewhat larger than S03S04. All of the order parameters converge to 1 for large and positive values of c0L. For large and negative values of c0L the order parameters S01 and S03 converge to −1, S02 converges to −1/2, and S04 converges to −3/7. For any given c0L<0 the extent of ordering increases along the series |S04|, |S02|, |S03| and |S01|. In the region 15c0L0 all four order parameters differ. For c0L15 (very strong perpendicular ordering), all four order parameters converge to their respective lower limits.

Figure 1 is useful in cases where one seeks ascribing appropriate single-term axial potentials to specific systems.

Potentials

The potential forms corresponding to L = 1, 2, 3 are shown in Figures 2A–C for c0L=10. Interestingly, the potential with L = 3 and K = 0 exhibits two minima: one for θ = 0, and the other for θ=cos1(15) ([congruent with] 116.6°) (Figure 2C).

Figure 2
Axial potentials, u(θ), shown as a function of θ (we show 3D rather than 2D plots for easy comparison with rhombic potentials, u(θ,[var phi])). The various potentials are defined by the coefficients c01=10 (A), c02=10 (B), ...

Figures 2D–F show examples of potentials comprising the L = 1 and L = 2 (part D), L = 1 and L = 3 (part E), and L = 1, L = 2 and L = 3 terms (part F) spherical harmonics. All of the coefficients are equal to 10, representing thus strong local ordering.1319 Figure 2 indicates that even within the scope of axial symmetry diverse potential forms can be devised, although the corresponding order parameters are quite similar (Figure 1).

Let us focus on case E, where the potential is given by the L = 1 and L = 3 spherical harmonics. The angle where the minimum with θ ≠ 0 occurs (call it θmin) is given by the ratio Rc01/c03 according to the expression:

θmin=12cos1[4R+915].
(6)

Thus, one can control θmin by varying the ratio R. The requirement that the absolute value of the cosine function be smaller than, or equal to, 1, implies that −6 ≤ R ≤ 1.5; consequently 0° < θmin < 60°. We found empirically that the range of 60° < θmin ≤ 70° (70° < θmin ≤ 90°) can be accessed by combining the axial L = 1, L = 2 and L = 3 (L = 2 and L = 3) spherical harmonics.

Given the diversity in their forms, axial potentials of the kind examined above are expected to be useful for analyzing (with SRLS or other stochastic models) NMR relaxation in proteins. We found previously that one has to account for structural asymmetry in order to obtain physically relevant result.1416 This can be accomplished in a straightforward manner by allowing for rhombic potentials.1416 However, the respective data-fitting calculations are demanding. The computation will become substantially more effective if structural asymmetry is accounted for within the scope of axial potentials by allowing for separate local ordering and local diffusion frames.4547 This strategy, although less straightforward than allowing for rhombic potentials, was found to be appropriate. In refs 45–47 the potential comprises a single term with L = 2, K = 0. The axial potentials considered above comprise at least two terms. Using them in data-fitting might obviate the need to separate the local ordering and local diffusion frames, while maintaining the axial symmetry of the potential. Such a scenario is both effective and economical.

PDFs

To gain further insight into multi-term axial potentials, u(θ), we examine this scenario within the scope of the associated PDFs given in Cartesian coordinates. The underlying potentials are specified by the coefficients, c0L, shown in rows A–F of Table 1. All of these potentials are strong.

Table 1
Potential Coefficients, c0L, Associated with the PDFs of Figures 3A3F, and angle θmaxa.

Figure 3A shows the PDF corresponding to the strong potential u(θ)=−10 cos θ. For the implied small values of θ one may approximate u(θ) by its Taylor expansion around θ = 0. This yields u(θ) [congruent with] −10×(1 − θ2/2), which has been used in the past to describe wobble-in-a-cone.48 Thus, Figure 3A may be considered as representing the wobble-in-a-cone PDF. Figure 3B shows the PDF corresponding to the potential given by c01=30 and c03=5.1. Based on its shape, it may be considered as representing a generalized version of the wobble-in-a-cone PDF.

Figure 3
Axial probability densities, exp(−u), given in Cartesian coordinates (x = sinθ cos[var phi] |exp(−u)|, y = sinθ sin[var phi] |exp(−u)|, z = cosθ |exp(−u)|). The potentials, u, are defined by the ...

The PDFs shown in Figures 3C3E (corresponding to the potentials C–E in Table 1) are cone-shaped, with various vertex angles and thicknesses. The vertex angle is given by θmax, which is the same as θmin of the corresponding potential. θmax is approximately 21, 57 and 72° for the PDFs of Figures 3C, ,3D3D and and3E,3E, respectively (Table 1). A simple model for describing internal motions is Woessner’s model, where the probe is confined to reside on the infinitely thin surface of a cone with variable vertex angle.49 The potentials underlying the PDFs shown in Figures 3C3E represent generalizations of Woessner’s model. The PDF shown in Figure 3B corresponds to θmax [congruent with] 7°; one may consider the underlying potential as representing a generalization of either wobble-in-a-cone48 or Woessner’s model49.

Finally, the PDF of Figure 3F may be considered a generalization of ideal perpendicular ordering, where the probe is confined to lie in the xy plane of the local director frame.

As pointed out above, the entire range of 0 ≤ θ 90° can be accessed with appropriate combinations of the L = 1–3, K = 0 axial spherical harmonics.

Orientation versus Alignment

Figure 4 shows results for polar and non-polar scenarios, with large positive potential coefficients representing strong parallel ordering, and large negative potential coefficients representing strong perpendicular ordering. The polar case with strong parallel ordering is shown in Figure 4A. It is represented by the potential u(θ)= −10 cos θ Part a shows the Cartesian representation of u which (as pointed out above) resembles the atomic orbital pz. Part c shows the spherical-coordinate-representation of u with minimum at θ = 0 maximum at θ = π. The PDF is shown in Cartesian coordinates in part b. It is directed along +z, as implied by c01 being positive. In addition, it is lobe-shaped as compared to the sphere-shaped potential of Figure 4Aa, because (as pointed out above) for large c01>0 potential values that are positive virtually do not contribute to |exp(−u)|, while they contribute to |u|. The spherical-coordinate-representation of the PDF appears in part d; it shows the function values on the z-axis, and illustrates their variations.

Figure 4
4-panel representations of the axial scenarios L = 1, K = 0 and c01=10 (part A, panels a–d); L = 2, K = 0 and c02=10 (part B, panels a–d); L = 1, K = 0 and c01=10 (part C, panels a–d); and L = 2, K = 0 and c02= ...

The non-polar case with strong parallel ordering is shown in Figure 4B. It is represented by the potential u(θ)=10[32cos2θ12] Part a shows its Cartesian-coordinate-representation which resembles the atomic orbital dz2.26 Part c shows the spherical-coordinate-representation of the potential with minima at θ = 0 and π, and maximum at θ=π2. The PDF is shown in Cartesian coordinates in part b. It comprises two lobes disposed symmetrically along +z and −z. They are thinner than the PDF lobe corresponding to L = 1 (Figure 4Ab), indicating reduced spatial distribution, as implied by the functional form of non-polar Y20 as compared to polar Y10. The spherical-coordinate-representation of the PDF appears in Figure 4Bd, illustrating the variations in PDF magnitude; comparison with Figures 4Ad shows larger spatial distribution in the L= 1 case.

The distinction between orientation, implied by L = 1 with symmetry B1u, and alignment, implied by L = 2 with symmetry Ag, is illustrated with particular clarity by the Cartesian PDFs (Figures 4Ab and 4Bb): for L = 1 one has directionality (only positive z-values are featured) whereas for L = 2 there is no directionality (both positive and negative z-values are featured).

Figure 4C shows the results obtained for L = 1 and c01=10. Rendering c01 negative implies conversion of z into −z for both the Cartesian-, and the spherical-coordinate-representations.

Figure 4D shows the results obtained for L = 2 and c02=10. Rendering c02 negative implies conversion of z into −z for the spherical-coordinate-representations. For the Cartesian-coordinate-representations this change consists of u being converted into −u.

3.2. Rhombic Ordering

Order Parameters

Axial orienting potentials, u(θ)=c0LD00L(0,θ,0), are defined by a single potential coefficient, c0L. The derivation of c0L from S0L, based on pre-calculated profiles (Figure 1), is unambiguous. The rhombic orienting potentials and the two-term axial orienting potentials considered in this study are defined by two potential coefficients, from which two order parameters are derived. It is not obvious that two order parameters, e.g., S02 and S22, correspond to a unique potential form, i.e., a unique pair of potential coefficients c02 and c22.

We examine this matter below for the following combinations: c01 and c11 derived from S01 and S11; c02 and c22 derived from S02 and S22; and c01 and c03 derived from S01 and S03. The coefficients c01 and c11 have been varied from −20 to 20 in steps of 0.5. For each point with coordinates (c01,c11) we calculated S01 and S11 according to eq 3 with j = 1, and eqs 4a, 4b and 5. This yielded the 3D surfaces S01(c01,c11) and S11(c01,c11). Figure 5A displays the contours of these surfaces in the c01/c11 plane. This procedure was also carried out for c02 and c22 (using eq 3 with j = 0, and eqs 4a, 4c and 5) to obtain the surfaces S02(c02,c22) and S22(c02,c22) with contours shown in Figure 5B. Finally, c01 and c03 yielded the surfaces S01(c01,c03) and S03(c01,c03), with contours shown in Figure 5C.

Figure 5
Contour plots of the surfaces S01(c01,c11) and S11(c01,c11) (part A), S02(c02,c22) and S22(c02,c22) (part B), and S01(c01,c03) and S03(c01,c03) (part C), as a function of the respective potential coefficients. The designations a–d represent ...

It can be seen that if two “isolines” intersect they will do so only once. Thus, if a solution exists, it is unique. The “isolines” are very dense for c01 and c11 smaller than approximately 5 in absolute value. This implies c01 and c11 being determined with good certainty in this parameter range. The isolines become much sparser for c01 and c11 of approximately |10|, and very sparse for c01 and c11 larger than |10|. This implies c01 and c11 being determined with reduced certainty when their absolute values are large. Note that Figure 5A is symmetric with respect to both the horizontal and the vertical lines through the point (c01,c11)=(0,0). The designations a–d refer to the extreme values of |S01| and |S11| in the (−20, 20) parameter range for c01 and c11. The extreme limits of S01 are given in section 3.1. above; the extreme limit of |S11| is 2.

Figure 5B shows the “isolines” of S02 and S22 as a function of c02 and c22. The parameter range where potential coefficients can be determined with good certainty include relatively small values of |c02 | and |c22|, as in the L = 1 case. In addition, it also includes large values of |c02| and |c22|, provided they are similar in magnitude. Figure 5B is only symmetrical with respect to the vertical line through the point (c02,c22)=(0,0). The designations a–d refer to the extreme values of |S02| and |S22| in the (−20, 20) parameter range for c02 and c22. The extreme limits of S02 are given in section 3.1 above; the extreme limit of |S22| is 32.

Figure 5C shows “isolines” of S01 and S03 as a function of c01 and c03. Here the parameter ranges where potential coefficients can be determined with good certainty include 5c015 and 20c0320. This figure is asymmetric. The designations a–d refer to the extreme values of |S01| and |S03| in the −20, 20 parameter range for c01 and c03. The extreme limits of |S01| and |S03| are given in section 3.1. above.

The information inherent in Figure 5 is important for predicting the certainty of the results to be obtained from scenarios with different symmetries.

Orientation versus Alignment

Figure 6 shows results for rhombic polar and non-polar scenarios represented by potential coefficients c01=c11=10 (part A) and c02=c22=10 (part B). An additional example, represented by polar-potential coefficients c01=10,c11=10 (part C) and non-polar-potential coefficients c02=10,c22=10 (part D), is provided.

Figure 6
4-panel representations of the rhombic scenarios L = 1, K = 0, 1 and c01=c11=10 (part A, panels a–d); L = 2, K = 0, 2 and c02=c22=10 (part B, panels a–d); L = 1, K = 0, 1 and c02=c22=10 (part C, panels a–d); and L = ...

Let us focus on Figure 6A, associated with the linear combination of Y10 and 12(Y11Y11) with coefficients c01=c11=10; the respective potential is given by u(θ,φ)=10cosθ102sinθcosφ (Table S3 of the Supporting Information). Figure 6B is associated with the linear combination of Y20 and 12(Y22Y22) with coefficients c02=c22=10; the respective potential is given by u(θ,φ)=10(32cos2θ12)1032sin2θcos(2φ) (Table S3 of the Supporting Information). The 4-panel representations illustrating these scenarios are shown in Figures 6A and and6B,6B, respectively.

The Cartesian L = 1 potential of Figure 6Aa has negative values in the xyz quadrant and positive values in the (−x)(−y)(−z) quadrant. Both this function and the Cartesian PDF (Figure 6Ab) have the primary polar axis in the zx plane. The spherical-coordinate representation of the potential has minimum, whereas the spherical-coordinate representation of the PDF has maximum, for θ=tan120.955rad(54.7°) and [var phi]= 0 (Figures 6Ac and 6Ad).

The Cartesian L = 2 potential has negative values in a relatively thin region centered at the zx plane (the coordinates have both positive and negative values) with maximum |z| exceeding maximum |x|, and positive values in lobe-shaped regions along the y-axis (the coordinates have both positive and negative values) (Figure 6Ba). The Cartesian PDF is aligned preferentially along ±z and distributed in a relatively thin regions centered at the zx plane (Figure 6Bb). The spherical-coordinate-representation of the potential has minimum value for θ = 0 and any [var phi], as well as θ = π and any [var phi]; it has maximum value for θ = [var phi] = π/2; the spherical-coordinate-representation of the PDF has reversed extrema (Figures 6Bc and 6Bd).

The difference between orientation and alignment is demonstrated with particular clarity by the Cartesian representation of the PDF (Figures 6Ab and 6Bb). Figure 6A illustrates comprehensively a polar ordering scenario based on one of the two lowest odd-L rhombic potentials. Figure 6B illustrates comprehensively a non-polar ordering scenario based on the lowest even-L rhombic potential.

We examine now the rhombic ordering scenario represented by the polar potential given by c01=10 and c11=10, shown in Figure 6C. The Cartesian potential has negative values in the (−x)yz quadrant and positive values in the x(−y)(−z) quadrant (Figure 6Aa). Both this function and the Cartesian PDF (Figure 6Cb) have the primary polar axis in the z(−x) plane. Figure 6Ca can be obtained from Figure 6Aa, and Figure 6Cb from Figure 6Ab, by inversion with respect to the origin followed by reflection with respect to the xy plane. The spherical-coordinate representations of the potential and the PDF are shown in Figures 6Cc and 6Cd, respectively.

The rhombic ordering scenario represented by the non-polar potential given by c02=10 and c22=10 is shown in Figure 6D. The Cartesian potential and PDF (Figures 6Da and 6Db) are obtained from their counterparts in Figures 6Ba and 6Bb by a π/2 clockwise rotation around z. The spherical-coordinate representations of the potential and the PDF are shown in Figures 6Dc and 6Dd, respectively.

Figure S1 of the Supporting Information is analogous with Figure 6 except that c01=10 and c11=10 in Figure S1A; c02=10 and c22=10 in Figure S1B; c01=10 and c11=10 in Figure S1C; c02=10 and c22=10 in Figure S1D.

Let us assume that a given protein site with largely known site symmetry is subject to investigation. As shown above, well-defined symmetry operations convert parts ad of Figure 6A ((6B)6B) into the corresponding parts of Figure 6C, S1A and S1C ((6D,6D, S1B and S1D). Thus, one can devise a potential form that largely matches the known site symmetry. Further refinement can then be carried out with analysis of NMR relaxation data.

3.3. Applications

To illustrate the generality of our approach we select four essentially different examples. They include 15N relaxation from a protein dissolved in aqueous solution, 15N relaxation from a polycrystalline protein, ESR lineshape analysis from a nitroxide-labeled protein dissolved in aqueous solution, and 2H and 13C relaxation in one of the phases of a medium-size molecule exhibiting polymorphism in the solid-state.

15N relaxation in proteins in the liquid-state

We studied previously with SRLS 15N spin relaxation of N–H bonds in plexin-B1 in aqueous solution.25 The order parameters S02 and S22 have been determined as the plateau values of time correlation functions calculated from a 110 ns MD trajectory, after eliminating the global motion. This is a commonly used strategy, based on the assumption that large timescale separation prevails between the global and local motions. Let us denote these parameters as S02(MD) and S22(MD). The potential coefficients c02 and c22 were obtained from S02(MD) and S22(MD) by solving the set of non-linear equations:25

S02(c02,c22)S02(MD)=0
(7a)

S22(c02,c22)S22(MD)=0
(7b)

Values of c02=3.56 and c22=3.47 were obtained in this manner for the local potential of the N–H bond of Gln56. We show the 4-panel representation for this residue (Figure 7A); PDF contours instead of 3D surfaces are shown.

Figure 7
4-panel representation associated with the potential given by c02=3.56 and c22=3.47, obtained with SRLS analysis of 15N relaxation from the N–H bond of Gln56 of plexin-B1 in solution25 (part A). 4-panel representation associated ...

PDF contours we call “atomistic” have also been calculated in ref 25 for Gln56 directly from the MD trajectory. Let us compare them with the “mesoscopic” PDF contours associated with the analytical potential with coefficients c02=3.56 and c22=3.47. Figure 8(c) shows the "atomistic" PDF contours;25 Figure 8(d) (which is the same as Figure 7Ad) shows the “mesoscopic” PDF counterparts. Figure 8(c) is considerably more asymmetric than Figure 8(d). We show in Figure 7B the 4-panel representation associated with the “enhanced mesoscopic” potential given by c02=3.56, c22=3.47, c03=3.00, and c13=5.00. The “enhanced mesoscopic” PDF contours (Figure 7Bd) are considerably more asymmetric than the “mesoscopic” PDF contours (Figure 7Ad), and quite similar to the “atomistic” PDF contours (Figure 8(c)).

Figure 8
PDF contour plots (in spherical coordinates, θ and [var phi]) for residue Gln56 of plexin-B1 calculated directly from the MD trajectory25 (c). PDF contour plots for residue Gln56 of plexin-B1 corresponding to the potential with coefficients ...

Thus, one can improve mesoscopic PDF contours by adjusting their forms so as to match the corresponding atomistic PDF contours. This is accomplished by linearly combining different spherical harmonics. The best-fit mesoscopic potential underlying the best-fit mesoscopic PDF may be considered a good approximation to the atomistic potential implicit in the atomistic PDF. Enhanced mesoscopic potentials obtained in this manner may be compared among sites, providing new insights. Importantly, they may be used without further adjustment in data-fitting schemes, reducing thereby the number of parameters allowed to vary.

The following comments are in order. (1) One of the future prospects of this work is to develop computational methods50 for extracting the best-fit mesoscopic PMF directly from the MD-based PDF. In this process the form of the PMF will be enhanced by adding spherical-Harmonics-based terms systematically; error estimations on the potential coefficients will emerge automatically. (2) In the context of protein dynamics by NMR relaxation there are previous reports on MD-derived PMFs used to interpret order parameters. For example, it was shown that a parabolic PMF (given by c02Y20 in the limit of strong axial potential) applies to rapidly fluctuating N–H bonds.51 However, this PMF would not apply to more general cases, such as illustrated by Figure 8(c). On the other hand, in our approach the PMF may comprise even-L and/or odd-L spherical harmonics of axial and/or rhombic symmetry; clearly, this is a general approach.

15N relaxation in proteins in the solid-state

The authors of ref 36 studied 15N spin relaxation in polycrystalline thioredoxin. Assuming that the local motions are in the extreme motional narrowing limit, the experimental powder patterns have been interpreted in terms of effective partially-reduced rhombic 15N chemical shift anisotropy (CSA) tensors. We calculated from the tensor components of Gly21 the potential coefficients c02=1.71 and c22=4.21 using eqs 7a,b, and employed these values to generate the 4-panel representation shown in Figure 9A.

Figure 9
4-panel representation based on the potential given by c02=1.71 and c22=4.21, obtained from the order parameters of the N–H bond of Gly21 of polycrystalline thioredoxin (part A). 4-panel representation based on the potential given by c02=3.31 ...

It can be seen that the coefficients of the local potentials at the N–H site of Gln56 of plexin-B1 in the liquid state ( c02=3.56 and c22=3.47), and the N–H site of Gly21 of thioredoxin in the solid-state ( c02=1.71 and c22=4.21), differ moderately in their absolute values, substantially in the ratios c22/c02, and have opposite signs. Thus, the main difference lies in the symmetry of the potential.

Nitroxide relaxation in T4L Lysozyme in solution

In Figure 9B we show the 4-panel representation for residue 72 of T4 Lysozyme which has been mutated into cysteine and labeled with a nitroxide.35 Liquid-state ESR lineshapes were analyzed with SRLS. Two conformers have been detected at this site; we show conformer 2, for which a local potential with c02=3.31 and c22=1.67 was determined at 25 °C.

This case illustrates yet another type of local potential with magnitude similar to, but symmetry different from, the local potentials at the N–H sites of Gln56 of plexin-B1 in solution and Gly21 of thioredoxin in the solid-state.

NMR relaxation in benzenehexa-n-alkanoates

Our last example is related to restricted motion in the discotic mesophase of benzenehexa-7-alkanoate, and in phase I of benzenehexa-6-alkanoate. Having been determined to be in their extreme motional narrowing limit, these motions give rise to partially averaged rhombic powder spectra.37 We converted the axial and rhombic components of these tensors into the order parameters S02 and S22, calculated c02 and c22 using eqs 7a,b,25 and generated the 4-panel representations shown in Figure 10.

Figure 10
4-panel representation associated with the potential given by c02=1.49 and c22=1.81 obtained from the order parameters of 2H-labeled β-hydrogen of BHA7 in its discotic mesophase (part A). 4-panel representation associated with the potential ...

In part A we show results obtained for the potential with c02=1.49 and c22=1.81, associated with 2H relaxation from the β-hydrogen of benzenehexa-7-alkanoate (BHA7) in its discotic mesophase. In part B we show results obtained for the potential with c02=1.28 and c02=2.68, associated with 13C relaxation from the 13C-labeled carbonyl of benzenehexa-7-alkanoate (BHA7) in phase I. Once again, we detect local potentials which differ moderately in magnitude, but quite a bit in symmetry, from the other local potentials discussed above.

Both Gln56 of plexin-B1 and Gly21 of thioredoxin are rather flexible residues. The nitroxide of T4 lysozyme is attached to the protein through a long flexible tether. The side-chains of the benzenehexa-n-alkanoates, which have been labeled with 2H or 13C, are flexible moieties. In accordance, our examples are associated with medium-strength local potentials. N–H sites in proteins residing in well-structured regions of the polypeptide were found to be associated with strong local potentials.1416 On the other hand, side-chain methyl sites in proteins were found to be associated with weak local potentials.1416 Although the strength of the local potentials in our examples might not be typical of many protein N–H and methyl sites, the diversity in symmetry is likely to be general in nature. Further studies in this context are required.

We show in Table 2 the experimental order parameters S02 and S22 of the five examples considered, and the corresponding potential coefficients. We also show the Saupe order parameters Sxx, Syy and Szz, related to S02 and S22 as:

Sxx=+(3/2S22S02)/2
(8a)

Syy=(3/2S22+S02)/2
(8b)

Szz=S02
(8c)
Table 2
Potential Coefficients and Corresponding Order Parameters in Irreducible Tensor Notation ( S02 and S22) and in Saupe Tensor Notation (Sxx, Syy, Szz,) (Eqs 3, 4, 5 and 8).a

Table 2, together with the 4-panel representations illustrated in Figures 7, ,99 and and10,10, constitute an insightful and comprehensive description of the local spatial restrictions at mobile sites of molecules, in particular proteins.

Let us focus on the Cartesian potential and PDF representations shown in Figures 7Aa, 9Aa, 9Ba, 10Aa and 10Ba. The potentials and PDFs shown in Figures 9Aa and 10Ba are similar given that both coefficients are positive and the ratio c22/c02 is comparable: 2.46 in Figure 9Aa and 2.09 in Figure 10Ba. The potential in Figure 10Aa also has positive coefficients but c22/c02=1.21; hence its shape is different. This illustrates the sensitivity of potential and PDF shapes to the relative magnitudes of c02 and c22.

The signs of c02and c22 are also important factors in determining potential and PDF shapes. Let us focus on the spherical-coordinate-representations of the potentials (Figures with the designation “c”). Figure 7Ac, where both coefficients negative, shows two potential wells. Figure 9Ac, where c02 is positive and c22 is negative, shows two potential barriers. Figure 9Bc, where both coefficients are positive, shows three potential barriers.

4. Conclusions

Potentials consisting of linear combinations of the L = 1 spherical harmonics Y10 and 12(Y11Y11) provide a useful description of polar local ordering at mobile sites in proteins. Comparison was made with corresponding non-polar L = 2 potentials consisting of linear combinations of the spherical harmonics Y20 and (Y22+Y22).

The shapes of the L = 1 and L = 2 potentials and associated PDFs differ substantially. The orientational polarity of L = 1 potentials is reflected particularly well by the Cartesian-coordinate-representations of the PDF. Thus, for positive (negative) potential coefficients the primary polar axis points in a direction in the zx (z(−x)) plane determined by their magnitudes. By contrast, the L = 2 PDFs have their main symmetry axis pointing symmetrically along +z and −z.

A representation where potentials and probability density functions are shown in both Cartesian and spherical coordinate systems provides an extensive characterization of the local ordering.

Substantial diversity was detected among potentials obtained as linear combinations of L = 1, L = 2 and/or L = 3, and K = 0, spherical harmonics. Being multi-term axial functions, these potentials are expected to help in devising economical and computationally effective data-fitting schemes for NMR relaxation analysis in proteins.

The representation specified above was applied to potential coefficients c02and c22 from fits to several different experimental systems. For the examples considered we found similarity in potential strength and diversity in potential symmetry.

A method whereby experimental mesoscopic PDFs can be improved by using their atomistic counterparts as a benchmark is proposed. The underlying best-fit mesoscopic potentials may then be considered as reasonable approximations to the atomistic potentials implicit in the atomistic PDFs.

Supplementary Material

S1

Acknowledgments

This work was supported by the Israel–U.S.A. Binational Science Foundation (Grant No. 2010185 to E.M. and J.H.F.), and the Israel Science Foundation (Grants No. 437/11 and No. 369/15 to E.M.). This work was also supported by NIH/NIGMS grant P41GM103521 to J.H.F.

Footnotes

Supporting Information Available: Table S1, showing the relevant Part of the D2h Point Group Character Table, and Table S2, showing the functions D00L(0,θ,0) for L = 1−4, are provided. In addition figure S1, analogous with Figure 6 except that for c01=10 and c11=10 in Figure S1A; c22=10 and c22=10 in Figure S1B; c01=10 and c11=10 in Figure S1C; c02=10 and c22=10 in Figure S1D, is provided. This material is available free of charge via the internet at http://pubs.acs.org.

References

1. Debye P. Polar Molecules. Dover Publications, Inc; New York: 1929.
2. Perrin F. Mouvement Brownien d’un Ellipsoide (I). Dispersion Dielectrique pour des Molecules Ellipsoidales. Le J de Phys ét le Radium. 1934;5(10):497–511.
3. Freed JH. Anisotropic Rotational Diffusion and Electron Spin Resonance Linewidths. J Chem Phys. 1964;41:2077–2083.
4. Nordio PL, Busolin P. Electron Spin Resonance Line Shapes in Partially Oriented Systems. J Chem Phys. 1971;55:5485.
5. Polnaszek CF, Bruno GV, Freed JH. ESR Lineshapes in the Slow-Motional Region: Anisotropic Liquids. J Chem Phys. 1973;58:3185–3199.
6. Emsley JW, editor. NMR of Liquid Crystals. Reidel; Dordrecht: 1983.
7. Freed JH, Nayeem A, Rananavare SB. In: The Molecular Dynamics of Liquid Crystals. Luckhurst GR, Veracini CA, editors. Kluwer Academic Publishers; The Netherlands: 1994. p. 271. Chapter 12.
8. Polimeno A, Freed JH. A Many-Body Stochastic Approach to Rotational Motions in Liquids. Adv Chem Phys. 1993;83:89–204.
9. Polimeno A, Freed JH. Slow-Motional ESR in Complex Fluids: the Slowly Relaxing Local Structure Model of Solvent Cage Effects. J Phys Chem. 1995;99:10995–11006.
10. Liang Z, Freed JH. An Assessment of the Applicability of Multifrequency ESR to Study the Complex Dynamics of Biomolecules. J Phys Chem B. 1999;103:6384–6396.
11. Liang Z, Freed JH, Keyes RS, Bobst AM. An Electron Spin Resonance Study of DNA Dynamics Using the Slowly Relaxing Local Structure Model. J Phys Chem B. 2000;104:5372–5381.
12. Liang Z, Lou Y, Freed JH, Columbus L, Hubbell WL. A Multifrequency Electron Spin Resonance Study of T4 Lysozyme Using the Slowly Relaxing Local Structure Model. J Phys Chem B. 2004;108:17649–1765.
13. Tugarinov V, Liang Z, Shapiro Yu E, Freed JH, Meirovitch E. A Structural Mode-Coupling Approach to 15N NMR Relaxation in Proteins. J Am Chem Soc. 2001;123:3055–3063. [PubMed]
14. Meirovitch E, Shapiro Yu E, Polimeno A, Freed JH. Protein Dynamics from NMR: the SRLS Analysis Compared with MF Analysis. J Phys Chem A. 2006;110:8366–8396. [PMC free article] [PubMed]
15. Meirovitch E, Shapiro Yu E, Polimeno A, Freed JH. Structural Dynamics of Biomolecules by NMR: the Slowly Relaxing Local Structure Approach. Prog NMR Spectrosc. 2010;56:360–405. [PMC free article] [PubMed]
16. Meirovitch E, Polimeno A, Freed JH. Protein Dynamics by NMR Spin Relaxation: The Slowly Relaxing Local Structure Perspective. In: Harris RK, Wasylishen RE, editors. Encyclopedia of Magnetic Resonance. John Wiley; Chichester, U.K.: 2012. p. 3591.
17. Tugarinov V, Shapiro Yu E, Liang Z, Freed JH, Meirovitch E. A novel view of domain flexibility in E. coli adenylate kinase based on structural mode-coupling. J Mol Biol. 2002;315:171–186. [PubMed]
18. Shapiro Yu E, Kahana E, Tugarinov V, Liang Z, Freed JH, Meirovitch E. Domain flexibility in ligand-free and inhibitor-bound Escherichia coli adenylate kinase based on a mode-coupling analysis of 15N spin relaxation. Biochemistry. 2002;41:6271–6281. [PubMed]
19. Shapiro Yu E, Kahana E, Meirovitch E. Domain Mobility in Proteins from NMR/SRLS. J Phys Chem B. 2009;113:12050–12060. [PubMed]
20. Laio A, Parrinello M. Escaping free-energy minima. Proc Nat Acad Sci. 2002;20:12562–12566. [PubMed]
21. Marini F, Camilloni C, Provasi D, Broglia RA, Tiana G. Metadynamic sampling of free-energy landscapes of proteins coupled with a Monte Carlo algorithm. Gene. 2008;422:37–40. [PubMed]
22. Darve E, Rodríguez-Gómez D, Pohorille A. Adaptive biasing force method for scalar and vector free energy calculations. J Chem Phys. 2008;128:144120–144132. [PubMed]
23. Minoukadeh K, Chipot C, Leliévre T. Potential of mean force calculations: a multiple-walker adaptive biasing force approach. J Chem Theory Comput. 2010;6:1008–1017.
24. Su L, Cukier RI. Hamiltonian replica exchange method study of Escherichia coli and Yersinia pestis HPPK. J Phys Chem B. 2009;113:16197–16208. [PubMed]
25. Zerbetto M, Buck M, Meirovitch E, Polimeno A. Integrated computational approach to the analysis of NMR relaxation in proteins: application to ps-ns main-chain 15N-1H and global dynamics of the Rho GTPase binding domain of plexin-B1. J Phys Chem B. 2011;115:376–388. [PMC free article] [PubMed]
26. e-Chemical Portal, D2h point group. http://www.webqc.org/printable-symmetrypointgroup-d2h.html.
27. de Gennes PG, Prost J. The Physics of Liquid Crystals. Oxford Universiy Press; 1993.
28. Zare RN. Angular Momentum. Wiley; 1988. p. 226. Ch. 5, Application 13.
29. Meirovitch E, Nayeem A, Freed JH. Analysis of Protein-Lipid Interactions Based pn Model Simulations of Electron Spin Resonance Spectra. J Phys Chem. 1984;88:3454–3465.
30. Misra SK, Freed JH. In: Molecular Motions in Multifrequency Electron Paramagnetic Resonance. Misra SK, editor. Vol. 11. Wiley-VCH; New-York: 2011. pp. 497–544.
31. Dzikowski B, Freed JH. Encyclopedia of Biophysics. In: Roberts GCK, editor. Membrane Fluidity. Springer-Verlag; Berlin, Heidelberg: 2012.
32. Meirovitch E, Liang Z, Freed JH. Protein Dynamics in the Solid State from 2H NMR Lineshape Analysis: a Consistent perspective. J Phys Chem B. 2015;119:2857–2868. [PMC free article] [PubMed]
33. Meirovitch E, Liang Z, Freed JH. Protein Dynamics in the Solid State from 2H NMR Lineshape Analysis. II MOMD Appplied to C–D and C–CD3 Probes. J Phys Chem B. 2015;119:14022–14032. [PMC free article] [PubMed]
34. Zerbetto M, Anderson R, Rech M, Bouguet-Bonnet S, Liqun Zhang L, Meirovitch E, Polimeno A, Matthias Buck M. Analysis of NMR relaxation in proteins by the Modelfree formalism and by a combined experimental and computational approach: Application to ps-ns mainchain 15N-1H and global dynamics of the Rho GTPase binding domain of plexin-B1 in the dimeric state. J Phys Chem B. 2013;117:174–184. [PMC free article] [PubMed]
35. Zhang Z, Fleissner MR, Tipikin DS, Liang Z, Moscicki JK, Earle KA, Hubbell WL, Freed JH. Multifrequency Electron Spin Resonance Study of the Dynamics of Spin-Labeled T4 Lysozyme. J Phys Chem. 2010;114:5503–5521. [PMC free article] [PubMed]
36. Yang Y, Tasayco ML, Polenova T. Dynamics of Reasssembled Thioredoxin Studies with Magic Angle NMR: Snapshots from different Timescales. J Am Chem Soc. 2009;131:13690–13702. [PMC free article] [PubMed]
37. Maliniak A, Greenbaum S, Poupko R, Zimmermann H, Luz Z. Deuterium and Carbon-13 NMR of the Solid and Discotic Phases of Three Benzenehexa-n-alkanoates. J Phys Chem. 1993;97:4832–4840.
38. Ferrarini A, Nordio PL, Moro GJ. In: The Molecular Dynamics of Liquid Crystals. Luckhurst GR, Veracini CA, editors. Kluwer Academic Publishers; The Netherlands: 1994. Chapter 3.
39. Greco C, Marini A, Frezza E, Ferrarini A. From the Molecular Structure of Spectroscopic and Material Properties: Computational Investigation of a Bent-Core Nematic Liquid Crystal. ChemPhysChem. 2014;15:1336–1344. [PubMed]
40. Biscarini F, Chiccoli C, Parsini P, Zannoni C. On a Simple Model of Ferroelectricity in Liquid Crystals. J of Non-Crystalline Solids. 1991;131–133:1190–1193.
41. van der Vorst CPJM, Picken SJ. Electric field poling of acceptor-donor molecules. J Opt Soc Am B. 1990;7:320–325.
42. Polnaszek CF, Freed JH. Electron Spin Resonance Studies of Anisotropic Ordering, Spin Relaxation, and Slow Tumbling in Liquid Crystalline Solvents. J Phys Chem. 1975;79:2283–2306.
43. Meirovitch E, Igner D, Igner E, Moro G, Freed JH. Electron-Spin Relaxation and Ordering in Smectic and Supercooled Nematic Liquid Crystals. J Chem Phys. 1982;77:3915–3938.
44. Barnes JP, Freed JH. Dynamics and Ordering in Mixed Model Membranes of Dimyristoyl-phosphatadyl-choline and Di-myristoyl-phosphatadyl-serine: a 250 GHz Electron Spin Resonance Study Using Cholestane. Biophys J. 1998;75:2532–2546. [PubMed]
45. Meirovitch E, Zerbetto M, Polimeno A, Freed JH. Backbone dynamics of deoxy and carbonmonoxy hemoglobin by NMR/SRLS. J Phys Chem B. 2011;115:143–157. [PMC free article] [PubMed]
46. Meirovitch E, Shapiro Yu E, Zerbetto M, Polimeno A. SRLS analysis of the 15N spin relaxation from E. coli Ribonuclease HI: the tensorial perspective. J Phys Chem B. 2012;116:886–894. [PubMed]
47. Meirovitch E. SRLS analysis of 15N relaxation from bacteriophage T4 lysozyme: a tensorial perspective that features domain motion. J Phys Chem B. 2012;116:6118–6127. [PubMed]
48. Petersen NO, Chan SI. More on the Motional State of Lipid Bilayer Membranes: Interpretation of Order Parameters Obtained from Nuclear Magnetic Resonance Experiments. Biochemistry. 1977;16:2657–2667. [PubMed]
49. Woessner DE. Spin Relaxation Processes in a Two-Spin System Undergoing Anisotropic Reorientation. J Chem Phys. 1962;36:1–4.
50. Budil DE, Sale KL, Khairy KA, Fajer PG. Calculating Slow-Motional Electron Paramagnetic Resonance Spectra from Molecular Dynamics Using a Diffusion Operator Approach. J Phys Chem A. 2006;110:3703–3713. [PubMed]
51. Massi F, Palmer AG., III Temperature Dependence of NMR Order Parameters and Protein Dynamics. J Am Chem Soc. 2003;125:11158–11159. [PubMed]