Molecular motion of spin labeled side chains in the C-terminal domain of RGL2 protein: a SDSL-EPR and MD study

Five singly spin labeled side chains at surface sites in the C-terminal domain of RGL2 protein have been analyzed to investigate the general relationship between nitroxide side chain mobility and protein structure. At these sites, the structural perturbation produced by replacement of a native residue with a nitroxide side chain appears to be very slight at the level of the backbone fold. The primary determinants of the nitroxide side chain mobility are backbone dynamics and tertiary interactions. On the exposed surfaces of alpha-helices, the side chain mobility is not restricted by tertiary interactions but appears to be determined by backbone dynamics, while in loop sites, the side chain mobility is even higher. For a better understanding of the changes in the EPR spectral line shape, molecular dynamics simulations were performed and found in agreement with EPR spectral data.     

Motion of spin-labeled side chains in T4 lysozyme: effect of side chain structure

Previous studies have shown that the mobility of nitroxide side chains in a protein, inferred from the electron paramagnetic resonance (EPR) spectra, can be used to classify particular sites as helix surface sites, tertiary contact sites, buried sites, or loop sites. In addition, the sequence dependence of mobility can identify regular secondary structure. However, in the most widely used side chain, an apparent interaction of the nitroxide ring with the protein at some helix surface sites gives rise to EPR spectra degenerate with those at tertiary contact sites. In the present study, we use selected sites in T4 lysozyme to evaluate novel nitroxide side chains designed to resolve this degeneracy. The results indicate that the reagent 3-(methanesulfonylthiomethyl)-2,2, 5,5-tetramethylpyrrolidin-1-yloxy reacts with cysteine to give a nitroxide side chain that has a high contrast in mobility between helix surface and tertiary contact sites, effectively resolving the degeneracy. The reagent 3-(iodomercuriomethyl)-2,2,5,5-tetramethyl-2, 5-dihydro-1H-pyrrol-1-yloxy reacts with cysteine to provide a mercury-linked nitroxide that also shows reduced interaction with the protein at most helix surface sites. Thus, these new side chains may be the preferred choices for structure determination using site-directed spin labeling.     

Characterization of threonine side chain dynamics in an antifreeze protein using natural abundance 13C NMR spectroscopy

The dynamics of threonine side chains of the Tenebrio molitor antifreeze protein (TmAFP) were investigated using natural abundance (13)C NMR. In TmAFP, the array of threonine residues on one face of the protein is responsible for conferring its ability to bind crystalline ice and inhibit its growth. Heteronuclear longitudinal and transverse relaxation rates and the [(1)H]-(13)C NOE were determined in this study. The C alpha H relaxation measurements were compared to the previously measured (15)N backbone parameters and these are found to be in agreement. For the analysis of the threonine side chain motions, the model of restricted rotational diffusion about the chi(1) dihedral angle was employed [London and Avitabile (1978) J. Am. Chem. Soc., 100, 7159-7165]. We demonstrate that the motion experienced by the ice binding threonine side chains is highly restricted, with an approximate upper limit of less than +/-25 degrees.     

Main chain and side chain dynamics of oxidized flavodoxin from Cyanobacterium anabaena.

Oxidized flavodoxin from Cyanobacterium anabaena PCC 7119 is used as a model system to investigate the fast internal dynamics of a flavin-bearing protein. Virtually complete backbone and side chain resonance NMR assignments of an oxidized flavodoxin point mutant (C55A) have been determined. Backbone and side chain dynamics in flavodoxin (C55A) were investigated using (15)N amide and deuterium methyl NMR relaxation methods. The squared generalized order parameters (S(NH)(2)) for backbone amide N-H bonds are found to be uniformly high ( approximately 0.923 over 109 residues in regular secondary structure), indicating considerable restriction of motion in the backbone of the protein. In contrast, methyl-bearing side chains are considerably heterogeneous in their amplitude of motion, as indicated by obtained symmetry axis squared generalized order parameters (S(axis)(2)). However, in comparison to nonprosthetic group-bearing proteins studied with these NMR relaxation methods, the side chains of oxidized flavodoxin are unusually rigid.     

Temperature dependence of the internal dynamics of a calmodulin-peptide complex

The temperature dependence of the fast internal dynamics of calcium-saturated calmodulin in complex with a peptide corresponding to the calmodulin-binding domain of the smooth muscle myosin light chain kinase is examined using 15N and 2H NMR relaxation methods. NMR relaxation studies of the complex were carried out at 13 temperatures that span 288-346 K. The dynamics of the backbone and over four dozen methyl-bearing side chains, distributed throughout the calmodulin molecule, were probed. The side chains show a much more variable and often considerably larger response to temperature than the backbone. A significant variation in the temperature dependence of the amplitude of motion of individual side chains is seen. The amplitude of motion of some side chains is essentially temperature-independent while many show a simple roughly linear temperature dependence. In a few cases, angular order increases with temperature, which is interpreted as arising from interactions with neighboring residues. In addition, a number of side chains display a nonlinear temperature dependence. The significance of these and other results is illuminated by several simple interpretative models. Importantly, analysis of these models indicates that changes in generalized order parameters can be robustly related to corresponding changes in residual entropy. A simple cluster model that incorporates features of cooperative or conditional motion reproduces many of the unusual features of the experimentally observed temperature dependence and illustrates that side chain interactions result in a dynamically changing environment that significantly influences the motion of internal side chains. This model also suggests that the intrinsic entropy of interacting clusters of side chains is only modestly reduced from that of independent side chain motion. Finally, estimates of protein heat capacity support the view that the major contribution to the heat capacity of protein solutions largely arises from local bond vibrations and solvent interactions and not from torsional oscillations of side chains.     

A theoretical study of rotational diffusion models for rod-shaped viruses. The influence of motion on 31P nuclear magnetic resonance lineshapes and transversal relaxation.

Information about the interaction between nucleic acids and coat proteins in intact virus particles may be obtained by studying the restricted backbone dynamics of the incapsulated nucleic acids using 31P nuclear magnetic resonance (NMR) spectroscopy. In this article, simulations are carried out to investigate how reorientation of a rod-shaped virus particle as a whole and isolated nucleic acid motions within the virion influence the 31P NMR lineshape and transversal relaxation dominated by the phosphorus chemical shift anisotropy. Two opposite cases are considered on a theoretical level. First, isotropic rotational diffusion is used as a model for mobile nucleic acids that are loosely or partially bound to the protein coat. The effect of this type of diffusion on lineshape and transversal relaxation is calculated by solving the stochastic Liouville equation by an expansion in spherical functions. Next, uniaxial rotational diffusion is assumed to represent the mobility of phosphorus in a virion that rotates as a rigid rod about its length axis. This type of diffusion is approximated by an exchange process among discrete sites. As turns out from these simulations, the amplitude and the frequency of the motion can only be unequivocally determined from experimental data by a combined analysis of the lineshape and the transversal relaxation. In the fast motional region both the isotropic and the uniaxial diffusion model predict the same transversal relaxation as the Redfield theory. For very slow motion, transversal relaxation resembles the nonexponential relaxation as observed for water molecules undergoing translational diffusion in a magnetic field gradient. In this frequency region T2e is inversely proportional to the cube root of the diffusion coefficient. In addition to the isotropic and uniaxial diffusion models, a third model is presented, in which fast restricted nucleic acid backbone motions dominating the lineshape are superimposed on a slow rotation of the virion about its length axis, dominating transversal relaxation. In an accompanying article the models are applied to the 31P NMR results obtained for bacteriophage M13 and tobacco mosaic virus.     

Structure and dynamics of a helical hairpin and loop region in annexin 12: a site-directed spin labeling study

A 30-residue nitroxide scan encompassing a helical hairpin and an extended loop in soluble annexin 12 (helices D and E in repeat 2; residues 134-163) has been analyzed in terms of nitroxide side chain mobility and accessibility to collision with paramagnetic reagents (Pi). Values of Pi for both O(2) and a Ni(II) metal complex (NiEDDA) are remarkably well correlated with the fractional solvent accessibility of the native side chains at the corresponding positions computed from the known crystal structure. This result demonstrates the utility of Pi as an experimental measure of side chain accessibility in solution, as well as the lack of structural perturbation due to the presence of the nitroxide side chain. The pattern of side chain mobility is also in excellent agreement with predictions from the crystal structure. The results presented here extend the correlations between mobility and structure described in earlier work on other helical proteins, and suggest their generality. The periodic dependence of Pi and mobility along the sequence of annexin 12 reveals the helical segments and their orientation in the fold, as expected for a nonperturbing nitroxide side chain. However, these data do not distinguish the helix-loop-helix motif from a continuous helix, because immobilized side chains in the short loop sequence maintain the periodicity. As shown here, the ratio of Pi values for O(2) and NiEDDA clearly delineates the loop region, due to size exclusion effects between the two reagents. A new feature evident in a nitroxide scan through multiple secondary elements is a modulation of the basic Pi and mobility patterns along the sequence, apparently due to differences in helix packing and backbone motion. Thus, in the short helix D, residues are consistently more mobile and accessible throughout the sequence compared to the residues in the longer, less-solvated and more ordered helix E.     

Side chain mobility and ligand interactions of the G strand of tear lipocalins by site-directed spin labeling.

Side chain mobility, accessibility, and backbone motion were studied by site-directed spin labeling of sequential cysteine mutants of the G strand in tear lipocalins (TL). A nitroxide scan between residues 98 and 105 revealed the alternating periodicity of mobility and accessibility to NiEDDA and oxygen, characteristic of a beta-strand. Residue 99 was the most inaccessible to NiEDDA and oxygen. EPR spectra with the fast relaxing agent, K(3)Fe(CN)(6), exhibited two nitroxide populations for most residues. The motionally constrained population was relatively less accessible to K(3)Fe(CN)(6) because of dynamic tertiary contact, probably with side chain residues of adjacent strands. With increasing concentrations of sucrose, the spectral contribution of the immobile component was greater, indicating a larger population with tertiary contact. Increased concentrations of sucrose also resulted in a restriction of mobility of spin-labeled fatty acids which were bound within the TL cavity. The data suggest that sucrose enhanced ligand affinity by slowing the backbone motion of the lipocalin. The correlation time of an MTSL derivative (I) attached to F99C resulted in the lack of side chain motion and therefore reflects the overall rotation of the TL complex. The correlation time of F99C in tears (13.5 ns) was the same as that in buffer and indicates that TL exists as a dimer under native conditions. TL-spin-labeled ligand complexes have a shorter correlation time than the protein alone, indicating that the fatty acids are not rigidly anchored in the cavity, but move within the pocket. This segmental motion of the ligand was modulated by protein backbone fluctuations. Accessibility studies with oxygen and NiEDDA were performed to determine the orientation and depth of a series of fatty acid derivatives in the cavity of TL. Fatty acids are oriented with the hydrocarbon tail buried in the cavity and the carboxyl group oriented toward the mouth. In general, the mobility of the nitroxide varied according to position such that nitroxides near the mouth had greater mobility than those located deep in the cavity. Nitroxides positioned up to 16 carbon units from the hydrocarbon tail of the ligand are motionally restricted and inaccessible, indicating the cavity extends to at least this depth. EPR spectra obtained with and without sucrose showed that the intracavitary position of lauric acid in TL is similar to that in beta-lactoglobulin. However, unlike beta-lactoglobulin, TL binds 16-doxyl stearic acid, suggesting less steric hindrance and greater promiscuity for TL.     

Thermodynamic interpretation of protein dynamics from NMR relaxation measurements

Protein dynamics and thermodynamics can be characterized through measurements of relaxation rates of side chain (2)H and (13)C, and backbone (15)N nuclei using NMR spectroscopy. The rates reflect protein motions on timescales from picoseconds to milliseconds. Backbone and methyl side chain NMR relaxation measurements for several proteins are beginning to reveal the role of protein dynamics in protein stability and ligand binding.     

Structural origins of nitroxide side chain dynamics on membrane protein α-helical sites

Understanding the structure and dynamics of membrane proteins in their native, hydrophobic environment is important to understanding how these proteins function. EPR spectroscopy in combination with site-directed spin labeling (SDSL) can measure dynamics and structure of membrane proteins in their native lipid environment; however, until now the dynamics measured have been qualitative due to limited knowledge of the nitroxide spin label's intramolecular motion in the hydrophobic environment. Although several studies have elucidated the structural origins of EPR line shapes of water-soluble proteins, EPR spectra of nitroxide spin-labeled proteins in detergents or lipids have characteristic differences from their water-soluble counterparts, suggesting significant differences in the underlying molecular motion of the spin label between the two environments. To elucidate these differences, membrane-exposed α-helical sites of the leucine transporter, LeuT, from Aquifex aeolicus, were investigated using X-ray crystallography, mutational analysis, nitroxide side chain derivatives, and spectral simulations in order to obtain a motional model of the nitroxide. For each crystal structure, the nitroxide ring of a disulfide-linked spin label side chain (R1) is resolved and makes contacts with hydrophobic residues on the protein surface. The spin label at site I204 on LeuT makes a nontraditional hydrogen bond with the ortho-hydrogen on its nearest neighbor F208, whereas the spin label at site F177 makes multiple van der Waals contacts with a hydrophobic pocket formed with an adjacent helix. These results coupled with the spectral effect of mutating the i ± 3, 4 residues suggest that the spin label has a greater affinity for its local protein environment in the low dielectric than on a water-soluble protein surface. The simulations of the EPR spectra presented here suggest the spin label oscillates about the terminal bond nearest the ring while maintaining weak contact with the protein surface. Combined, the results provide a starting point for determining a motional model for R1 on membrane proteins, allowing quantification of nitroxide dynamics in the aliphatic environment of detergent and lipids. In addition, initial contributions to a rotamer library of R1 on membrane proteins are provided, which will assist in reliably modeling the R1 conformational space for pulsed dipolar EPR and NMR paramagnetic relaxation enhancement distance determination.     

Hydration dependence of backbone and side chain polylysine dynamics: a 13C solid-state NMR and IR spectroscopy study

The molecular dynamics of solid poly-L-lysine has been studied by the following natural abundance (13)C-NMR relaxation methods: measurements of the relaxation times T(1) at two resonance frequencies, off-resonance T(1rho) at two spin-lock frequencies, and proton-decoupled T(1rho). Experiments were performed at different temperatures and hydration levels (up to 17% H(2)O by weight). The natural abundance (13)C-CPMAS spectrum of polylysine provides spectral resolution of all types of backbone and side chain carbons and thus, dynamic parameters could be determined separately for each of them. At the same time, the conformational properties of polylysine were investigated by Fourier transform infrared spectroscopy. The data obtained from the different NMR experiments were simultaneously analyzed using the correlation function formalism and model-free approach. The results indicate that in dry polylysine both backbone and side chains take part in two low amplitude motions with correlation times of the order of 10(-4) s and 10(-9) s. Upon hydration, the dynamic parameters of the backbone remain almost constant except for the amplitude of the slower process that increases moderately. The side chain dynamics reveals a much stronger hydration response: the amplitudes of both slow and fast motions increase significantly and the correlation time of the slow motion shortens by about five orders of magnitude, and at hydration levels of more than 10% H(2)O fast and slow side chain motions are experimentally indistinguishable. These changes in the molecular dynamics cannot be ascribed to any hydration-dependent conformational transitions of polylysine because IR spectra reveal almost no hydration dependence in either backbone or side chain absorption domains. The physical nature of the fast and slow motions, their correlation time distributions, and hydration dependence of microdynamic parameters are discussed.     

Dynamics of a de novo designed three-helix bundle protein studied by 15N, 13C, and 2H NMR relaxation methods

Understanding how the amino acid sequence of a polypeptide chain specifies a unique, functional three-dimensional structure remains an important goal, especially in the context of the emerging discipline of de novo protein design. Alpha3D is a single chain protein of 73 amino acids resulting from a de novo design effort. Previous solution nuclear magnetic resonance studies of alpha3D confirm that the protein adopts the designed structure of a three-helix bundle. Furthermore, alpha3D has been previously shown to possess all of the major thermodynamic and structural characteristics of natural proteins, though it shares no sequence homology to any protein sequence in the database. In this work, the backbone and side-chain dynamics of alpha3D were investigated using 15N, 13C, and 2H nuclear magnetic resonance relaxation methods with the aim of assessing the character of the internal motions of this native-like protein of de novo design. At the backbone level, both 15N and 13C(alpha) relaxation studies indicate highly restrictive motion on the picosecond to nanosecond time scale in the alpha-helical regions of alpha3D, with increasing mobility at the ends of the alpha-helices and in the two loop regions. This is largely consistent with what is seen in proteins of natural origin. Overall, the view provided by both 2H and 13C methyl relaxation methods suggest that the side chains of alpha3D are more dynamic compared to natural proteins. Regions of relative flexibility bound clusters of rigid methyl-bearing side-chain groups that are interspersed with aromatic and beta-branched amino acids. The time scale of motions associated with methyl-bearing side chains of alpha3D are significantly longer than that seen in natural proteins. These results indicate that the strategies underlying the design of alpha3D have largely, but not completely, captured both the structural and dynamic character of natural proteins.     

The analysis of NMR relaxation data in terms of multiple internal motions

A novel formalism for estimating the complex motions of proteins and other flexible macromolecules from NMR relaxation measurements is applied to ¹³C NMR relaxation data on the Bovine Pancreatic Trypsin Inhibitor (M. W. 6,500). Six experimental parameters measured at two field strengths are accounted for by a minimum of three motions at each carbon group. Low frequency components make small but finite contribution to the relaxation of all resonances, suggesting a general low frequency distortion of the backbone. Rotational diffusion of the protein makes a relatively minor contribution to the relaxation process. For aliphatic groups, rotation of side chains dominates relaxation.     

Mapping Membrane Protein Backbone Dynamics: A Comparison of Site-Directed Spin Labeling with NMR 15N-Relaxation Measurements

The ability to detect nanosecond backbone dynamics with site-directed spin labeling (SDSL) in soluble proteins has been well established. However, for membrane proteins, the nitroxide appears to have more interactions with the protein surface, potentially hindering the sensitivity to backbone motions. To determine whether membrane protein backbone dynamics could be mapped with SDSL, a nitroxide was introduced at 55 independent sites in a model polytopic membrane protein, TM0026. Electron paramagnetic resonance spectral parameters were compared with NMR ¹⁵N-relaxation data. Sequential scans revealed backbone dynamics with the same trends observed for the R₁ relaxation rate, suggesting that nitroxide dynamics remain coupled to the backbone on membrane proteins.     

Mapping Membrane Protein Backbone Dynamics: A Comparison of Site-Directed Spin Labeling with NMR N-Relaxation Measurements

The ability to detect nanosecond backbone dynamics with site-directed spin labeling (SDSL) in soluble proteins has been well established. However, for membrane proteins, the nitroxide appears to have more interactions with the protein surface, potentially hindering the sensitivity to backbone motions. To determine whether membrane protein backbone dynamics could be mapped with SDSL, a nitroxide was introduced at 55 independent sites in a model polytopic membrane protein, TM0026. Electron paramagnetic resonance spectral parameters were compared with NMR ¹⁵N-relaxation data. Sequential scans revealed backbone dynamics with the same trends observed for the R₁ relaxation rate, suggesting that nitroxide dynamics remain coupled to the backbone on membrane proteins.     

Temperature dependence of fast dynamics in proteins

The temperature dependence of the internal dynamics of recombinant human ubiquitin has been measured using solution NMR relaxation techniques. Nitrogen-15 relaxation has been employed to obtain a measure of the amplitude of subnanosecond motion at amide N-H sites in the protein. Deuterium relaxation has been used to obtain a measure of the amplitude of motion of methyl-groups in amino-acid side chains. Data was obtained between 5 and 55 degrees C. The majority of amide N-H and methyl groups show a roughly linear (R(2)>0.75) temperature dependence of the associated Lipari-Szabo model-free squared generalized-order parameter (O(2)) describing the amplitude of motion. Interestingly, for those sites showing a linear response, the temperature dependence of the backbone is distinct from that of the methyl-bearing side chains with the former being characterized by a significantly larger Lambda-value, where Lambda is defined as d ln(1 - O)/d lnT. These results are comparable to the sole previous such study of the temperature dependence of protein motion obtained for a calmodulin-peptide complex. This suggests that the distinction between the main chain and methyl-bearing side chains may be general. Insight into the temperature dependence is gathered from a simple two-state step potential model.     

Solid-state NMR investigation of the dynamics of the soluble and membrane-bound colicin Ia channel-forming domain.

Solid-state NMR spectroscopy was employed to study the molecular dynamics of the colicin Ia channel domain in the soluble and membrane-bound states. In the soluble state, the protein executes small-amplitude librations (with root-mean-square angular fluctuations of 0-10 degrees ) in the backbone and larger-amplitude motions (16-17 degrees ) in the side chains. Upon membrane binding, the motional amplitudes increase significantly for both the backbone (12-16 degrees ) and side chains (23-29 degrees ), as manifested by the reduction in the C-H and H-H dipolar couplings and (15)N chemical shift anisotropy. These motions occur not only on the pico- to nanosecond time scales, but also on the microsecond time scale, as revealed by the (1)H rotating-frame spin-lattice relaxation times. Average motional correlation times of 0.8 and 1.2 micros were extracted for the soluble and membrane-bound states, respectively. In comparison, both forms of the colicin Ia channel domain are completely immobile on the millisecond scale. These results indicate that the colicin Ia channel domain has enhanced conformational mobility in the lipid bilayer compared to the soluble state. This membrane-induced mobility increase is consistent with the loss of tertiary structure of the protein in the membrane, which was previously suggested by the extended helical array model [Zakharov et al. (1998) Proc. Natl. Acad. Sci. U.S.A. 95, 4282-4287]. An extended structure would also facilitate protein interactions with the mobile lipids and thus increase the protein internal motions. We speculate that the large mobility of the membrane-bound colicin Ia channel domain is a prerequisite for channel opening in the presence of a voltage gradient.