Hydrogen skeleton, mobility and protein architecture.

The mobility of the proton-proton radial vectors is introduced as a quantitative measure for the structural dynamics of organic materials, especially protein molecules. As defined for the entire molecule, the hydrogen mobility (HM) is proposed as an "order parameter," which describes the effect of motional narrowing on inter-proton dipole-dipole interactions. HM satisfies all requirements of an order parameter in the Landau molecular field theory of phase transitions. The wide-line NMR second moments needed to obtain HM are exactly defined and measurable physical quantities, which are not produced by mathematical fitting and do not carry the limitations and restrictions of any model (theoretical formalism). We first demonstrate the usefulness of HM on small organic molecules with data taken form the literature. We outline its link with structural and functional characteristics on a range of proteins: HM provides a model-free parameter based on first principles that can clearly distinguish between globular and intrinsically disordered proteins, and can also provide insight into the behavior of disease-related mutants.

Introduction

The importance of the dynamic nature of protein molecules cannot be overestimated: understanding their function in its entirety cannot be achieved by considering them as rigid molecules. Their dynamic state is affected by their thermal and chemical surroundings. The limitations of information inherent in the applied experimental method(s) also cannot be neglected. Beyond the general and conceptual question of motion (to the effect as it is detailed in Appendix A), the identification and proportion of mobile parts (residues), the physical characteristics of their motion and the effect of mobility on their reactivity/function all require special attention.

Addressing mobility is motivated by works approaching this subject from different angles. Halle criticized experimental approaches and results on aqueous protein solutions, by stating that “the progress is less and erratic and the results given by different experimental methods are contradictory and more or less model dependent in the interpretation”. Saito and Kobayashi examined carefully the physical basis of protein architecture. Wlodawer and colleagues discussed the limitations of structural information obtained by X-ray crystallography. They emphasized the lack of direct determination of hydrogen positions, which suggests the role of proton magnetic resonance experiments to complement such measurements very effectively. The regular and generally used NMR procedures were summarized e.g., by Dyson and Wright, with particular focus on intrinsically disordered proteins, by Duer in the case of studying molecular motions in solids, by Antzutkin in molecular structure determination in biology, and by Smith in the case of studying protein structures. The change in protein structural paradigm from order to disorder (from globular to intrinsically disordered (ID) proteins) was discussed in reference 8. Realizing the limitations of all approaches applied to date, we would like to add a novel wide-line NMR approach to complement the usual set of methods; the experimental analyses (including wide-line NMR) of the IDPs were introduced earlier in.,

Our experimental approach relies on wide-line NMR spectrometry, as we seek global, non-selective information. The elementary probe we use is the atomic nucleus of hydrogen (the proton magnetic dipolar moment) residing in the protein molecules. The relevant variants of wide-line NMR spectrometry are detailed in the book chapters 12 in reference 9, and 13 in reference 10 The wide-line spectrum and its even moments, especially the second moment (a main indicator of internal motion in solids; see Appendix B) are used here as a first step. In prior works, we outlined the physical basis of the applied method and detailed the evaluation and interpretation protocols. The key factor in our approach is the <span class="Chemical">hydrogen mobility (i.e., the time dependence of proton-proton radial vectors), which is indicative of NMR-visible motions and not the high-frequency lattice vibrations. This hydrogen mobility introduced depends on the molecular structure, temperature, and also on the chemical environment in the case of solutions.

In the paper, we describe in detail the theoretical basis of the second moment of the NMR spectrum. We first consider the case of a rigid lattice; the motional narrowing is treated subsequently. We give the definition of the novel order parameter we term hydrogen mobility, and put it into the context of order parameters in general. We then review the experimental second moment data found in the literature for several organic compounds and we calculate the corresponding <span class="Chemical">hydrogen mobility values. We test and give reasons for the application of the hydrogen mobility. Based on experimental data, we calculate the hydrogen mobility factors for some selected proteins, and discuss the importance of hydrogen mobility as a global, model-free order parameter in describing function-related features of flexibility of proteins.

NMR Spectrum, Even Moments, and Multi-Component Wide-Line Spectrum of “Solids”

Rigid lattice

This section is prefaced with the remark that there is no analytical description of the sp.ectra of multi-sp.in systems (e.g., molecules made of numerous atomic nuclei of various types). The series expansion technique can be used, to obtain a moments expansion of the time-domain NMR-signal (FID or solid echo).

The magnetic dipolar interaction plays a key role in the wide-line spectroscopy of solids. The magnitudes (or lengths) of the vectors connecting the interacting nuclei and their directions with respect to the external magnetic field directly determine the architecture (geometrical arrangement, structure or topology) of the nuclear spin system. This is the reason why the direct measurement of dipolar interactions is essential for describing the structure of proteins. The basic quantum theory of the dipolar interaction and properties of the NMR spectrum was given by Van Vleck in his Nobel laureate work. Determining the second moments is viewed as a first step, used subsequently to characterize the rigid lattice and the internal motions of molecules.

To understand the NMR characteristics, the wide-line NMR spectrum and the even moments, and to lay the groundwork for planning and interpreting the experiments, we invoke short parts from Abragam’s (chapters IV and X, ref. 14) and Slichter’s (chapter 13, ref. 15) basic monographs.

On the basis of the Van Vleck theory, the second moment coming from dipole-dipole interaction between identical spins (II) is (p. 79, ref. 15)

where is the gyromagnetic ratio, is the magnetic constant, N is the number of nuclei, r is the magnitude of the vector connecting the ith and jth nuclei, θ is the angle between the internuclear vector r and the time independent magnetic induction vector B0, and the double summation means averaging over all the spins and all the neighbors. If all the resonant spins are located in equivalent positions, the double summation is reduced to a single one, being independent from one of the indices. There are then N equivalent sums, one for each value of j. The second moment is

in that case. (We have to emphasize that it is not the case with protons in proteins, in spite of the general claim in some papers.) Each term in Eqs. 1 and 2 is clearly of the order of, where is the contribution of the k spin to the local field at sp.in j. Equations 1 and 2 precisely define the local field, which enables one to compare an exactly defined theoretical quantity with the measured (experimental) value. The summation has a short range character because of the strong distance dependence of the term.

It is also necessary to address what “equivalent positions” mean in general and especially in a protein molecule? “The equivalent arrangements of the non-zero nuclear magnets around the resonant nuclei” gives the answer. In our case, not all the protons are equivalent, and consequently, a few terms from one part of the double summation are to be used in the second moment calculation. How many? Decomposition of the components in the wide-line spectrum helps to estimate the number.

If, besides protons, there are also other nuclear species, they also contribute to the local fields and to the total second moment (the contribution is given in ref. 15 not necessary for us at this moment). In the case of proton resonance, the dipolar contributions of other nuclei (e.g., 13C, 15N, 17O) are zero or negligible; this is not true for other sp.ecies (e.g., 23Na or 35Cl).

As a consequence, the hydrogen nuclei (the protons) give the topological map of the immobile protein molecule and the rigid protein-<span class="Chemical">water system (which we may term the hydrogen skeleton). The use of Eq. 1 or Eq. 2 gives the theoretical basis for controlling the topological construction of the hydrogen skeleton of our molecules. It is to be mentioned, again, that the arrangement of the hydrogen atoms is missing from the X-ray maps.

Both the length and direction of the r proton-proton vectors are given as structural elements in single crystals. Only the length of the proton-proton vector exists in a powder (polycrystalline or lyophilized) sample because of “space-averaging.” The second moment then simply depends on the length (magnitude) of r. For interacting sp.in-pairs, the contribution to the second moment of the rigid lattice is

The fourth moment is only used here as a control for the analytical form of the line shape in Chapter 4 (wide-line NMR spectra and second moments of proteins), i.e., the details are not cited, only the references 10, 12, and 13 are mentioned again. The important point for our treatment is that the ratio of the fourth moment to the second moment-square is 3 for a Gaussian line.

Motional narrowing

Observing the temperature dependence of the systems studied, the position of protons and the dipole-dipole interaction between the neighboring pairs will be time dependent. The second moments (Eqs. 1 and 2) will be time dependent by the r and θ quantities as a consequence of atomic motions. The phenomenon is known as motional narrowing, and the results are the narrowing of the NMR sp.ectrum and reduced moments as a consequence of time averaging. For a spin pair, Abragam and Slichter following Andrew and Eades gave the description of time averaging by rotation around a given axis

where γ is the angle between the radius vector r and the rotation axis. (If γ = 90°, the time-averaged reduced moment is the rigid lattice value reduced by the factor of 4.) Equation 4 forms the basis of the magic-angle spinning (MAS) method, in which the sample is spun physically at the angle γ = 54°44′ (cosγ = 1/3) with respect to the direction of the magnetic field. This way, the reduced second moment takes zero value for each spin pair. When there are internal molecular motions, the averaging is done also over the possible γ angles, which does not result in zero value. The very simple reduction factor in Equation 4 is only valid for the intramolecular contribution to the second moment. The intermolecular contribution that results from interactions between sp.ins that belong to different molecules is affected by the rotation in a more complicated way since both the distances between spins and the orientations of the spin-spin vectors change.

The 2 types of averaging, space-average and time-average, can be done for every topological and motional model. The summations (Eqs. 1–2) and several possible motions should be considered, not forgetting that the summation means that the individual contributions are added up. This procedure involves numerous assumptions about the topology and the motional states and, consequently, the results are model dependent. Instead, we chose to determine the second moments experimentally: both the reduced second moment and the rigid-lattice value can be measured. In the typical temperature-dependence of the NMR second moment, i.e., on the frequency of the internal motion, the decrease of the second moment occurs in a temperature range, the position and width of which depends on the material studied.

Order parameter

An order parameter, which is characteristic of the motional state of the proton-proton pair, can be theoretically defined on the basis of Equations 1–4. This Hydrogen-The motional order parameter we intend to introduce, can be calculated by introducing the notation as

where M2(T) is the second moment measured at an arbitrary temperature T, M2(RL) is the second moment measured in the rigid lattice state at a sufficiently low temperature, and HM refers to the average mobility of the proton-proton vectors according to Equations 1–4. If the spectrum can be decomposed into e.g., 2 components as found in several cases, then 2 order parameters, i.e., 2 mobility values, can be measured. Equation 5 then takes the form

where means the “total,” “broad,” and “narrow” spectral components, respectively. HM is zero for an immobile (rigid/ordered) system and it is one for highly mobile (fluid-like) systems. The former state is probably realized at the temperature of liquid helium, whereas the latter one is only achieved above the melting point. There are 2 important points to be addressed: (1) What is the quantity to normalize with, i.e., what is the second moment of the rigid lattice? (2) Which spectral component can be assigned to which part of the molecule?

“Order” is connected with the time dependence or mobility of the proton spin (nuclear magnetic moment) as HM = 0 for a rigid system and HM = 1 for the liquid state; both these states and intermediate states can also be measured. In Appendix C, some considerations are given for the “order parameters” as defined through the Lipari-Szabo formalism of protein-solvent systems.

About order parameters in general

We have learned from the thermal properties of solids, e.g., in reference 19 and especially in the Landau mean-field theory of phase transition that a large variety of systems can be described by a single, temperature-dependent order parameter. The system might be the magnetization in ferromagnets, the dielectric polarization in a ferroelectric system, the fraction of superconducting electron-pairs in a superconductor, the director distribution in nematic liquid crystals, or the fraction of neighbor A-B bonds of the total bonds in an alloy AB or A3B. The physical quantity, which the order parameter refers to, is clearly given and its value must be between 0 and 1. The proposed HM(T) satisfies these criteria.

Experimental Second Moment Data in the Literature and Estimated Hydrogen Mobility

To appreciate the insight and relevance of HM as an order parameter in describing proteins, first we look into relevant data in the literature on small organic compounds. These results are related to motions of small molecules and parts of molecules, which helps us familiarize ourselves with the relevant orders of magnitude. The HM(T) values were calculated for cases where the second moment values were given at 2 temperatures at least.

The second moment measured for samples containing water of crystallization is M2 ~28·10−8 T2. The expected second moment for groups is only M2 ~3·10−8 T2.,

Works on the molecular motions detected in solid hydrocarbons- indicate the possible form of molecular motions, which range from torsional and rotational motions to translational diffusion. The second moment values measured for long-chain paraffins (n-CH2+2) at T = 82–99 K can be considered as rigid-lattice values with M2 values given in Table 1; these relevant spectra showed no fine structure. If the structure of a crystal is known and the lattice can be assumed rigid, the second moment of the proton resonance line can be calculated by applying the Van Vleck formula (Eq. 1). The crystal structure of n-C29H60 was given among the results of X-ray investigations of a series of normal paraffins (n-CH2+2; i = 6 – 44).- Although X-ray diffraction cannot give the locations of the hydrogen atoms, it accurately measures the C-C bond lengths. The geometric positions of the hydrogen atoms (protons) were then determined by assuming tetrahedral symmetry for the C-H bonds and by using the spectroscopic C-H bond lengths. The long-chain paraffins n-octadecane (n-C18H36), n-octacosane (n-C28H58) and dicetyl (n-C32H66) were investigated by continuous-wave NMR-sp.ectroscopy. These molecules consist of methylene and methyl groups. The rigid-lattice second moments measured at 82–99 K are 23.9·10−8 T2 (i = 18), 26.6·10−8 T2 (i = 28), and 27.1·10−8 T2 (i = 32). Fine structure in the sp.ectra at that temperature was not found. On the basis of the Van Vleck theory and Equation 1, the theoretical M2 values obtained are 27.3·10−8 T2 (i = 18), 27.0·10−8 T2 (i = 28), and 26.9·10−8 T2 (i = 32). The good agreement between the values measured and the ones calculated for n-octacosane and dicetyl proves the rigidity of the H skeleton, whereas the difference in the case of n-octadecane points to a small hydrogen mobility. The deviation from the two other paraffins is not surprising, because the lower melting point of n-octadecane also shows a looser crystal structure. The value of the HM mobility “order parameter” introduced here is 0.12.

Table 1. Second moments and hydrogen mobility for benzene and its deuterated versions

Temperature	Motional state	M2 [10−8 T2]	HM
C6H6
T < 68 K	Rigid lattice	9.72	0
145 K < T < 230 K	Rotation around C axis	1.61	0.83
C6H5D
T < 68 K	Rigid lattice	7.94	0
145 K < T < 230 K	Rotation around C axis	1.86	0.85
1,3,5-C6H3D3
T < 68 K	Rigid lattice	3.57	0
145 K < T < 230 K	Rotation around C axis	0.44	0.88

Systematic investigations were done on benzene, by replacing H atoms partly with D (<span class="Chemical">deuteron). The H-D change dilutes the magnetization of the H lattice, because of the much smaller nuclear moment of D, due to which its contribution to the second moment is approximately 2% of the proton magnetization. The great HM value (Table 1) can be explained by the rotation around the hexad axis of the full proton system (H-skeleton), that is, by the rotation of the benzene molecule. HM is practically independent of the H/D ratio.

In the case of N,N-dimethylanilin, the connection of the –N(CH3)2 group to the <span class="Chemical">benzene ring gives an estimated rigid lattice value of M2 = 18.4·10−8 T2. The measured second moment is M2 = 9.4…8.5·10−8 T2 in the temperature range of -190°C to 0°C, and it slowly decreases with increasing temperature, because of the lattice dilatation. The HM = 0.51 value is the consequence of the decrease of symmetry as compared that of the benzene. The actual value can be interpreted by the rotation of –CH3 groups and by the immobility of benzene ring.

For urea, the measured second moment is M2 = 20.8·10−8 T2 at ~195 K and it is M2 = 6.9·10−8 T2 at room temperature, and the <span class="Chemical">hydrogen mobility estimated by us is HM = 0.87. The high value exhibits a molecular motion of high degree of symmetry, higher than that of the –NH2 groups. Xylenes and mesitylene give measured second moments of 9.8·10−8 T2 to 9.9·10−8 T2 at 95 K, which decrease upon heating to 205 K due to the lattice dilatation. The interpretation of the results is similar to that of the N,N-dimethylanilin, that is, rotation of methyl groups and the immobile benzene ring.

The experimental work on hexamethylbenzene is outstanding regarding that several NMR characteristics were measured between 2 K and 450 K. The solid-echo radiofrequency pulse combination was used for the second moment measurements. The X-ray structure is known and the calculated rigid lattice second moment is M2 = 32.7·10−8 T2. The measured second moments are 20·10−8 T2 to 14.5·10−8 T2 in the 2 K to 50 K range, they take the values of M2 = 13.0·10−8 T2 at 90 K, and M2 = 2.5·10−8 T2 above 210 K. The corresponding HM values are 0.60 and 0.98, respectively. The whole molecule rotates above 210 K. One of the surprising conclusions is that methyl groups are in a motional state at 2 K.

For the compounds 1,4-dicyclohexylcyclohexane, <span class="Chemical">cyclohexane and adamantane (Bokor, M. et al., to be published), the results are summarized in Table 2. Isotropic (not uniaxial) rotation exists in the high symmetry cyclohexane and adamantane molecules in the solid phase.

Table 2. Second moments and hydrogen mobility for 1,4-dicyclohexylcyclohexane, cyclohexane and adamantane

Temperature	Motional state	M2 [10−8 T2]	HM
1,4-dicyclohexylcyclohexane
T < 310 K	Rigid lattice	24.9	0
335 K < T < 344 K	Torsional oscillation	15.9	0.36
T > 344 K	Rotation around long axis	6.8	0.73
cyclohexane
T < 150 K	Rigid lattice	26.5	0
156 K < T < 186 K	Rotation around C axis	6.4	0.76
T > 186 K	Isotropic rotation	1.4	0.95
T > ~200 K	Isotropic rotation + translational diffusion	~0	~1
adamantane
T ~4 K	Rigid lattice	24.8	0
room temperature	Isotropic rotation	0.9	0.95

The results for a few amino acids are also shown in Table 3. Amino acids are in zwitterionic form in the crystalline state. Rapid random reorientation of the –NH3 and the –CH3 groups in the solid amino acids were found above 150 K., There are 2 motional processes manifesting themselves in the narrowing of the NMR spectra with increasing temperature. The second moments fall in 2 successive steps with temperature as the frequency of motion of each process becomes comparable with the spectral width.

Table 3. Second moments and hydrogen mobility for selected amino acids

	M2 [10−8 T2]	M2 [10−8 T2]	M2 [10−8 T2]	HM
	T = 4.2 K	T = 77 K	room temp.	( Eq. 5 )
glycine		30.3	12.7	0.58
α-alanine		31.8	9.2	0.71
β-alanine		34.9	15.0	0.57
isoleucine	30		11	0.66
leucine	26		8	0.74
methionine	21		nd	nd
threonine	27.5		11.5	0.59
valine	26		11	0.67

A few generalizations are already apparent upon analyzing and comparing data on these organic compounds. Each molecule shows extensive internal mobility in the solid phase at room temperature. In many instances, the mobility of the H-H vectors can be as high as resulting in HM = 0.95. The mobile state exists also at lower temperatures, even at T = 4 K in some cases. The HM order parameter for the hydrogen mobility is a quantitative measure of the dynamic state, and it requires no model to be applied, in contrast with the second moment measured at a given temperature.

Wide-Line NMR Spectra and Second Moments of Proteins

Studies on small organic molecules and amino acids suggest that the order parameter introduced adequately captures the extent of hydrogen mobility, without invoking a model for interpreting experimental observations. Diakova et al. studied the wide-line 1H-NMR FID signals of dry and hydrated lysozyme powders. The FID of the hydrated protein powder was fitted to a sum of Gaussian functions (plus a constant at the highest water contents). It was found that the residual water at levels below 5 wt% is not rotationally mobile and the water-proton NMR signal is indistinguishable from the FID of the solid protein protons. The fast component of the FID was isolated and analyzed, which corresponds to the majority of the protein protons that do not experience significant dynamical averaging. These experiments infer that our approach might also be used for characterizing the general structural behavior of proteins.

To this end, we selected proteins that fall into 2 broad structural classes, ordered (folded) and intrinsically disordered (ID or ID protein, <span class="Gene">IDP), anticipating that IDPs show much more hydrogen mobility and much less order than folded ones. As a representative of ordered proteins, we chose lysozyme, the enzyme that has been amply studied for its structure, function and disease-causing mutations. We then extended our studies on several IDPs which have already been characterized in much detail. Thymosin β4 is a small actin-regulatory protein, which has been shown experimentally to be highly disordered (Á. Tantos, et al., to be published). α-synuclein is also a fully disordered protein,, involved in Parkinson disease where it undergoes a transition to a highly structured amyloid state. Familial mutations (A30P, E46K and A53T) promote this transition, possibly via altering the conformational equilibrium and structural flexibility of the protein., ERD14 is a plant IDP of stress-related functions,, also shown to be highly disordered. Our goal here is to characterize these proteins in detail to determine their HM and correlate their behavior with their physiological function and relationship with disease. It is to be noted that some of our conclusions rely on preliminary results, and further experimental results we intend to publish in later detailed publications.

Methods Applied

The applied NMR measurement and evaluation methods were summarized in chapter 13 of. In addition to the generally used FID signals, solid echo signals were also detected. Results for the entire temperature range covered (extended down to 4.2 K) are discussed in separate publications. The NMR signals detected in the time domain were analyzed and transformed to spectra similarly as it was described in reference 44.

1H NMR measurements and data acquisition were accomplished by Bruker AVANCE III and Bruker SXP 4–100 NMR pulse spectrometers at ω0/2π = 82.4 MHz with a stability better than ± 10−6, and with 2 ppm inhomogeneity of the magnet. The temperature was controlled by an open-cycle Janis cryostat with a stability of ± 0.1°C, the uncertainty of the temperature scale was ± 1°C. The data points in the figures are based on spectra recorded by averaging signals to reach the signal/noise ratio 50. The number of averaged NMR signals was varied to achieve the desired signal quality for each sample.

Results and Discussion

The sNMR-signals are limited to the lowest temperature of 4.2 K and to room temperature. The proteins lysozyme (T. Verebélyi, et al., to be published) and thymosin β4 were chosen as examples representing the expected lowest and highest HM values. The spectra of the FID and the solid echo signals differ only slightly from each other. The spectrum of the FID is somewhat wider as a consequence of the local magnetic fields coming from inhomogeneous proton-proton contributions. Unlike to the data reported earlier, the shape of the spectrum was not Gaussian either at T = 4.2 K or at room temperature (Figs. 1, 2, 3, and 4). The spectra could be decomposed into at least 2 components at both low and high temperatures. The presence of different spectral components indicates the heterogeneity of the proton spin-systems. The relative weights of the components vary with the protein types and the temperatures. The room-temperature narrow-spectrum component was wider than a signal coming from “free” water alone. The other proteins listed in Table 4 show similar behaviors.

Figure 1.1H NMR FID (left) and spectrum (right) of lyophilized lysozyme measured at T = 4.2 K, ν0 = 82.4 MHz. The sp.ectrum is 41.9 ± 0.2 kHz wide. The second moment (M2) of the spectrum is (16.1 ± 0.1)·10−8 T2, the fourth moment is M4 = (87 ± 8)·10−17 T. The shape factor of the spectrum is M22/M4 = 3.3 ± 0.3 (it is 3 for a Gaussian). Estimated errors are given.

Figure 2.1H NMR FID (left) and spectrum (right) of lyophilized lysozyme measured at T = 295.1 K, ν0 = 82.4 MHz. The spectrum was decomposed into a broad rigid lattice line and a motionally narrowed line. The FID shows fast rigid lattice and slow motionally narrowed relaxation accordingly. The second moment (M2) of the spectrum is (10.4 ± 0.1)·10−8 T2, the broad line is of M2 = (11.0 ± 0.1)·10−8 T2 while the narrow line has M2 = (0.9 ± 0.1)·10−8 T2. The spectral widths are 29.1 ± 0.1 kHz, 35.3 ± 0.2 kHz, and 7.1 ± 0.1 kHz, respectively. (The corresponding shape parameters are 3.50 ± 0.3, 3.32 ± 0.3, and 3.51 ± 0.3, respectively.) Estimated errors are given.

Figure 3.1H NMR FID (left) and sp.ectrum (right) of lyophilized Tβ4 measured at T = 7.4 K, ν0 = 82.4 MHz. The spectrum is 49.5 ± 0.2 kHz wide. The second moment (M2) of the spectrum is (20.6 ± 0.1)·10−8 T2, the fourth moment is M4 = (117 ± 10)·10−17 T. The shape factor of the spectrum is M22/M4 = 2.7 ± 0.3 (it is 3 for a Gaussian). Estimated errors are given.

Figure 4.1H NMR FID (left) and spectrum (right) of lyophilized Tβ4 measured at T = 300 K, ν0 = 82.4 MHz. The spectrum was decomposed into a broad rigid lattice line and a narrow motionally narrowed line. The FID shows fast rigid lattice and slow motionally narrowed relaxation accordingly. The second moment (M2) of the spectrum is (6.2 ± 0.1)·10−8 T2, the broad line is of M2 = (7.0 ± 0.1)·10−8 T2 while the narrow line has M2 = (0.5 ± 0.1)·10−8 T2. The spectral widths are 3.8 ± 0.1 kHz, 27.9 ± 0.1 kHz, and 0.5 ± 0.05 kHz, respectively. (The corresponding shape parameters are 3.8 ± 0.3, 3.0 ± 0.3, and 11.2 ± 0.3, respectively.) Estimated errors are given.

Table 4. The results on a polycrystalline globular and a few lyophilized powdered ID proteins

	Lysozyme	α-synuclein variants				ERD14	Tβ4
		A30P	E46K	WT	A53T
M2(He)	16.1	22.2	21.7	17.9 (22.0)	17.2 (22.0)	20.1	22.0
M2t(RT)	10.4	8.2	7.8	6.7 (8.6)	8.6 (8.6)	7.6	6.0
M2b(RT)	11.0	10.7	9.4	8.9 (11.0)	10.1 (10.1)	10.4	6.8
M2n(RT)	0.9	1.7	0.2	0.4	0.4	0.50	0.9
HMt	0.35	0.63	0.64	0.53 (0.72)	0.50 (0.72)	0.62	0.73
HMb	0.32	0.52	0.57	0.50 (0.50)	0.41 (0.54)	0.48	0.69
HMn	0.94	0.92	0.99	0.98	0.98	0.98	0.96

The parameter values in parentheses are given for a hypothetical rigid lattice state. The second moment (M2) values are given in 10−8 T2 = 1 G2 units, for the HM parameters are defined by Equation 6. The designations He and RT mean T ~4 K (liquid helium temperature) and room temperature, respectively. The index t refers to the whole spectrum, b stands for the broad while and n for the narrow spectrum component.

The parameter values in parentheses are given for a hypothetical rigid lattice state. The second moment (M2) values are given in 10−8 T2 = 1 G2 units, for the HM parameters are defined by Equation 6. The designations He and RT mean T ~4 K (liquid <span class="Chemical">helium temperature) and room temperature, respectively. The index t refers to the whole spectrum, b stands for the broad while and n for the narrow spectrum component.

The second moments calculated from the spectra and the HM hydrogen mobility parameters are summarized in Table 4; one should remember that the t-indexed parameters apply to whole molecules and indices b and n refer to the 2-component resolutions. The 2-component spectra and the relevant moments measured at room temperature indicate 2 types of hydrogen-hydrogen radial vector mobility, characteristic of 2 individual populations of residues of the molecule. The difference between the very-low-temperature second moment and the room-temperature second moment of the broad component shows that the molecule has no part made up entirely of static hydrogen atoms (immobile proton-proton radial vectors) at elevated temperatures. Even in the solid phase, the mobility HMt(T) characteristic of the whole molecule is considerably greater for the IDP molecules than for globular lysozyme. In accord, the value of HMt(T) is presumably in close association with structural disorder of proteins. It is of note that these types of molecular motions do not become frozen in aqueous solutions either. The second moments measured for the protein lysozyme at low and high temperatures ( and , respectively) give the smallest HMt value, which we use as a reference point for globular (structured) proteins. Thymosin β4 represents the other extremity, which shows the highest proton-proton radial vector mobility. All the experimental data measured by different methods proves the substantial disorder of the molecule. ERD14 also has a high mobility, which is in line with its largely disordered character as a molecular chaperone.

The results on α-synuclein mutants are important from another point of view. HM provides a qualitative measure of distinction, it is model-free and quantitative, which can be related to the function of IDPs. Whereas it is premature to draw general conclusions, this perfectly fits into the current trend of linking quantitative description structural features of IDPs with their function (unstructure-function relationship), it is of further note that proton-proton mobility of α-synuclein point mutants already shows significant differences at the temperature of liquid helium in the case of WT and A53T, in contrast with E46K and A30P. The smaller values for the α-synuclein variants wild-type (WT) and A53T (compared with the 2 other mutants) can only be explained by assuming that the molecule is not totally rigid even at such low temperatures and they show a mobility of HM ~0.20. It was presumed that the rigid-lattice second moments of the WT and the other α-synuclein mutants are the same with the value M2(He) = 22·10−8 T2 (numbers in parentheses, Table 4). This presumption is reasonable as these α-synuclein variants differ from each other in only one amino acid of 141. By the criterion of HM, wild-type and A53T are the most disordered of the four α-synuclein variants, i.e., they have the most dynamic molecular groups (hydrogen atoms). It is of note that earlier NMR data showed a similar behavior on their frozen aqueous solutions, and these differences in behavior seem to have bearing on the effect of familial mutations on local ordering being conducive to their transition to the disease-related amyloid state.

Conclusions and Outlook

The hydrogen mobility order parameter HM(T) proposed here represents the relative missing part of the time-independent (rigid lattice) <span class="Chemical">1H NMR second moment coming from the proton-proton dipole-dipole interaction in hydrogen rich molecules, including organic compounds, amino acids and proteins. HM can take values between 0 and 1, characteristic of a rigid lattice and liquid state, respectively, and the actual values are temperature dependent quantities.

In the context of the present study, the term motion is needed for the internal motions of a molecule (not for the translational diffusion of molecules in a liquid), and the time scale is set by the motions visible to NMR. The above results together with the introduction of the order parameter represent a theoretical and practical framework that can guide further investigations and assist in classifying proteins.

It is worthy to address why we chose the order parameter and not the second moment to characterize mobility? To calculate the second moment of a rigid lattice, precise atomic coordinates are needed, which can be obtained only by applying models. In the case of active molecular motions, the details of these motions should also be modeled. It is fair to ask how many parameters should be introduced to interpret a single measured quantity. Molecular motions result in almost twice as many parameters than a rigid system. The hydrogen mobility introduced here provides an overall dynamic attribute, i.e., a model-independent quantitative value, which characterizes the general internal dynamics of the molecule at the actual temperature. It is hard to overemphasize that only 2 measured quantities are used in the hydrogen mobility factor, namely the reduced second moment measured at the actual temperature and a reference second moment value measured for the rigid lattice state. It is not at all intended to exclude models as tools to help the interpretation of the measured second moments. For example, the 4 variants of α-synuclein have almost identical primary structures, yet their NMR spectra and second moments show explicit differences even at T < 10 K.

The knowledge of the measured M2(RL) gives significant help in explaining the NMR relaxation times. The coupling constant of the relaxation formalisms (see e.g., Chapter 12 in ref. 9, and Chapter 13 in refs.10,15 and 45) can be determined directly from the second moment of the rigid lattice. Therefore, not all 3 of the quantities—coupling constant, activation energy, and correlation time— only the latter two, should be determined from the relaxation time vs. temperature curves. We previously found that the three-parameter fitting of the relaxation time model produces absurd results in the case of the IDPs. The spatial and the temporal averages are immiscible categories for heterogeneous systems such as the protein molecules.

Our results prove that the proposed order parameter (HM) for hydrogen mobility is a well-defined physical quantity, which can be experimentally determined in a model-free way for small molecules and proteins alike. Furthermore, it does not involve any speculation or reference to “fitting parameters.” HM provides a theoretical and practical framework for the investigation of the correlations between global structural order and internal molecular mobility in solids. It provides a simple yet powerful method for the fast and quantitative distinction between globular and ID proteins. Whereas order parameters of individual resonances have the power of resolving detailed structural features of proteins (globular and ID alike), there is a strong trend in the <span class="Gene">IDP literature that uses global parameters of structure, disorder, and flexibility for approaching their function. In this sense, the novel parameter fits into the general characterization of IDPs. It even has the power to show quantitative differences between the dynamics of proteins, the sequence of which differ only in one amino acid. It provides a bridge connecting the dynamic properties measured for aqueous solutions and lyophilized proteins. And finally, HM calls attention to the importance of precise measurements and the adequate selection of reference temperatures.

Appendices

Appendix A

Thermal energy in a solid can produce motions ranging from lattice vibrations to restricted rotations of parts or the whole of molecules and even translational diffusion. The effect they produce depends on their frequency spectrum, spatial amplitude and the relationship of the motional axes with respect to those defining nuclear sp.in interactions. When motional averaging is present, the effective interaction which contributes to the spectral shape and width is the average over this motion.

Appendix B

Moment is a quantity that expresses the average or expected value of the first, second, third, or fourth power of the deviation of each component of a frequency distribution from some given value, typically mean or zero. The first moment is the mean, the second moment the variance, the third moment the skew, and the fourth moment the kurtosis. The dipolar interaction between protons gives symmetrical wide-line NMR spectra and consequently zero odd and non-zero even moments. The detailed information theory of the NMR line shape and moments are given in the work of Powles and Carazza.

Appendix C

For fast internal motions of proteins, the (generalized, or effective, or rigidity, or mobility) order parameter S, is used in the literature (see also e.g., refs. 50 and 51), which is the measure of the degree of spatial restriction of the motion. These order parameters meet a persistent and common claim and have the following common properties. (1) They are used for the interpretation of the relaxation characteristics (R1, R2, NMRD, correlation function) measured for aqueous protein solutions and not for solid state proteins as in the examples analyzed here by us. (2) The mathematical relations, which describe the above relaxation characteristics, were determined relying upon multi-parameter fittings. (3) That is, we think that it is not expedient to make a comparison between the order parameters S and HM, in the absence of a critical review on the whole field (and before additional measurements of HM for several more proteins).