A suite of Mathematica notebooks for the analysis of protein main chain 15N NMR relaxation data

A suite of Mathematica notebooks has been designed to ease the analysis of protein main chain ¹⁵N NMR relaxation data collected at a single magnetic field strength. Individual notebooks were developed to perform the following tasks: nonlinear fitting of ¹⁵N-T ₁ and -T ₂ relaxation decays to a two parameter exponential decay, calculation of the principal components of the inertia tensor from protein structural coordinates, nonlinear optimization of the principal components and orientation of the axially symmetric rotational diffusion tensor, model-free analysis of ¹⁵N-T ₁, -T ₂, and {¹H}–¹⁵N NOE data, and reduced spectral density analysis of the relaxation data. The principle features of the notebooks include use of a minimal number of input files, integrated notebook data management, ease of use, cross-platform compatibility, automatic visualization of results and generation of high-quality graphics, and output of analyses in text format.L. Spyracopoulos is an AHFMR Medical Research Senior Scholar     

Refinement of protein structure against non-redundant carbonyl <Superscript>13</Superscript>C NMR relaxation

Carbonyl ¹³C′ relaxation is dominated by the contribution from the ¹³C′ chemical shift anisotropy (CSA). The relaxation rates provide useful and non-redundant structural information in addition to dynamic parameters. It is straightforward to acquire, and offers complimentary structural information to the ¹⁵N relaxation data. Furthermore, the non-axial nature of the ¹³C′ CSA tensor results in a T₁/T₂ value that depends on an additional angular variable even when the diffusion tensor of the protein molecule is axially symmetric. This dependence on an extra degree of freedom provides new geometrical information that is not available from the NH dipolar relaxation. A protocol that incorporates such structural restraints into NMR structure calculation was developed within the program Xplor-NIH. Its application was illustrated with the yeast Fis1 NMR structure. Refinement against the ¹³C′ T₁/T₂ improved the overall quality of the structure, as evaluated by cross-validation against the residual dipolar coupling as well as the ¹⁵N relaxation data. In addition, possible variations of the CSA tensor were addressed. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Cross-correlation suppressed T1 and NOE experiments for protein side-chain 13CH2 groups

Relaxation measurements of side-chain ¹³CH₂-groups of uniformly ¹³C labeled human ubiquitin were performed at 600 MHz and 800 MHz magnetic field strength at 30°C. Dipole-dipole cross-correlated relaxation effects in T₁ experiments were suppressed by the combination of radio-frequency pulses and pulsed field gradients during the relaxation delay leading to monoexponential relaxation decays that allow a more accurate extraction of the ¹³C T₁ relaxation times. Heteronuclear ¹H-¹³C NOEs obtained by using different proton saturation schemes indicate that the influence of cross-correlation is small. The experimental T₁ and NOE data were interpreted in a model-free way in terms of a generalized order parameter and an internal correlation time.     

High-resolution heteronuclear NMR of human ubiquitin in an aqueous liquid crystalline medium

A mixture of dihexanoyl phosphatidylcholine and dimyristoylphosphatidylcholine in water forms disc-shaped particles, often referred toas bicelles [Sanders and Schwonek (1992) Biochemistry, 31, 8898–8905].These adopt an ordered, liquid crystalline phase, which can be maintained atvery low concentrations of the bicelles (down to 3% w/v). At thisconcentration the spacing between individual bicelles, on average, exceeds300 Å. The bicelles are shown to have a negligible effect on therotational diffusion of ubiquitin as judged by the ¹⁵NT₁ values of the backbone amides relative to those inisotropic aqueous solution. The protein exhibits a residual degree ofalignment which is proportional to the bicelle concentration, andapproximately collinear with ubiquitin's rotational diffusion tensor. Thedegree of alignment obtained offers unique opportunities for studying theprotein's structure and dynamics.     

Domain orientation in β-cyclodextrin-loaded maltose binding protein: Diffusion anisotropy measurements confirm the results of a dipolar coupling study

Maltose binding protein (MBP) is a 370-residue two-domain molecule involved in bacterial chemotaxis and sugar uptake. Rotational diffusion tensors were calculated for a complex between MBP and -cyclodextrin using backbone ¹⁵N T₁ and T₁ relaxation times and steady state ¹H-¹⁵N NOE values. The tensors obtained for each of the two domains in the protein were subsequently used to determine the relative domain orientation in the molecule. The average domain orientation determined using this approach agrees well with results from dipolar coupling data, but differs significantly from the domain orientation deduced from X-ray studies of the complex.     

High resolution proton magnetic resonance of sonicated phospholipids

We have recorded high resolution proton magnetic resonance spectra of sonicated phospholipid vesicles. The following lipids were used in separate experiments: phosphatidylglycerol, phosphatidylserine, phosphatidylethanolamine, and phosphatidylcholine from egg yolk as well as dimyristoyl phosphatidylcholine. Mixed lipid vesicles were also investigated. Assignments of the peaks associated with the various protons of the different lipids are presented. It is shown that in favorable cases, it is possible to resolve the different phospholipid head groups of mixed lipid samples. Spin lattice relaxation times (T₁) of each peak were collected at 500 MHz and 90 MHz. The influence of the addition of a small concentration of spin labeled phospholipid on i) the linewidths ii) the spin lattice relaxation times, was determined. It is shown that nitroxide radicals selectively broaden the peaks associated with the protons localized at a comparable depth of the bilayer. On the other hand, T₁ are less selectively perturbed. Potential applicability of ¹H-NMR for the investigation of lipid-proton specificity in membranes is discussed.     

Determination of 13Cα relaxation times in uniformly 13C/15N-enriched proteins

Summary Relaxation times of ¹³C carbons of uniformly ¹³C/¹⁵N-enriched probes have been investigated. The relaxation behaviour was analyzed in terms of a multispin system. Pulse sequences for the determination of T₁, T₂ and the heteronuclear NOE of ¹³C in uniformly ¹³C/¹⁵N-enriched ribonuclease T1 are presented. The experiments performed in order to obtain T₁ and the heteronuclear NOE were similar to those of the corresponding ¹⁵N experiments published previously. The determination of T₂ for the C-carbon in a completely labeled protein is more complicated, since the magnetization transfer during the T₂ evolution period owing to the scalar coupling of C–C must be suppressed. Various different pulse sequences for the T₂ evolution period were simulated in order to optimize the bandwidth for which reliable T₂ relaxation times can be obtained. A proof for the quality of these pulse sequences is given by fitting the intensity decay of individual ¹H–¹³C cross peaks, in a series of (¹H, ¹³C)-ct-HSQC spectra with a modified CPMG sequence as well as a T_1p sequence for the transverse relaxation time, to a single exponential using a simplex algorithm.     

Backbone dynamics of free barnase and its complex with barstar determined by 15N NMR relaxation study

Backbone dynamics of uniformly ¹⁵N-labeled free barnase and its complex with unlabelled barstar have been studied at 40°C, pH 6.6, using ¹⁵N relaxation data obtained from proton-detected 2D {¹H}-¹⁵N NMR spectroscopy. ¹⁵N spin-lattice relaxation rate constants (R₁), spin-spin relaxation rate constants (R₂), and steady-state heteronuclear {¹H}-¹⁵N NOEs have been measured at a magnetic field strength of 14.1 Tesla for 91 residues of free barnase and for 90 residues out of a total of 106 in the complex (excluding three prolines and the N-terminal residue) backbone amide ¹⁵N sites of barnase. The primary relaxation data for both the cases have been analyzed in the framework of the model-free formalism using both isotropic and axially symmetric models of the rotational diffusion tensor. As per the latter, the overall rotational correlation times (_m) are 5.0 and 9.5 ns for the free and complexed barnase, respectively. The average order parameter is found to be 0.80 for free barnase and 0.86 for the complex. However, the changes are not uniform along the backbone and for about 5 residues near the binding interface there is actually a significant decrease in the order parameters on complex formation. These residues are not involved in the actual binding. For the residues where the order parameter increases, the magnitudes vary significantly. It is observed that the complex has much less internal mobility, compared to free barnase. From the changes in the order parameters, the entropic contribution of NH bond vector motion to the free energy of complex formation has been calculated. It is apparent that these motions cause significant unfavorable contributions and therefore must be compensated by many other favorable contributions to effect tight complex formation. The observed variations in the motion and their different locations with regard to the binding interface may have important implications for remote effects and regulation of the enzyme action.     

Transverse relaxation optimized triple-resonance NMR experiments for nucleic acids

Triple resonance HCN and HCNCH experiments are reliable methods of establishing sugar-to-base connectivity in the NMR spectra of isotopicaly labeled oligonucleotides. However, with larger molecules the sensitivity of the experiments is drastically reduced due to relaxation processes. Since the polarization transfer between ¹³C and ¹⁵N nuclei relies on rather small heteronuclear coupling constants (11–12 Hz), the long evolution periods (up to 30–40 ms) in the pulse sequences cannot be avoided. Therefore any effort to enhance sensitivity has to concentrate on manipulating the spin system in such a way that the spin–spin relaxation rates would be minimized. In the present paper we analyze the efficiency of the two known approaches of relaxation rate control, namely the use of multiple-quantum coherence (MQ) and of the relaxation interference between chemical shift anisotropy and dipolar relaxation – TROSY. Both theoretical calculations and experimental results suggest that for the sugar moiety (H1-C1-N1/9) the MQ approach is clearly preferable. For the base moiety (H6/8-C6/8-N1/9), however, the TROSY shows results superior to the MQ suppression of the dipole–dipole relaxation at moderate magnetic fields (500 MHz) and the sensitivity improvement becomes dramatically more pronounced at very high fields (800 MHz). The pulse schemes of the triple-resonance HCN experiments with sensitivity optimized performance for unambiguous assignments of intra-residual sugar-to-base connectivities combining both approaches are presented.     

Selectively 13C-enriched DNA: Evidence from 13C1′ relaxation rate measurements of an internal dynamics sequence effect in the lac operator

Summary In order to study some internal dynamic processes of the lac operator sequence, the ¹³C-labeled duplex 5d(C₀G₁C₂T₃C₄A₅C₆A₇A₈T₉T₁₀) · d(A₁₀A₉T₈T₇G₆T₅G₄A₃G₂C₁G₀)3 was used. The spreading of both the H1 and C1 resonances brought about an excellent dispersion of the ¹H1-¹³C1 correlations. The spinlattice relaxation parameters R(C_z), R(C_x,y) and R(H_zC_z) were measured for each residue of the two complementary strands, except for the 3-terminal residues which were not labeled. Variation of the relaxation rates was found along the sequence. These data were analyzed in the context of the model-free formalism proposed by Lipari and Szabo [(1982) J. Am. Chem. Soc., 104, 4546–4570] and extended to three parameters by Clore et al. [(1990) Biochemistry, 29, 7387–7401; and (1990) J. Am. Chem. Soc., 112, 4989–4991]. A careful analysis using a least-squares program showed that our data must be interpreted in terms of a three-parameter spectral density function. With this approach, the global correlation time was found to be the same for each residue. All the C1-H1 fragments exhibited both slow (_s = 1.5) and fast (_f = 20 ps) restricted libration motions (S _inf2 ^sups =0.74 to 1.0 and S _inf2 ^supf =0.52 to 0.96). Relaxation processes were described as governed by the motion of the sugar relative to the base and in terms of bending of the whole duplex. The possible role played by the special structure of the AATT sequence is discussed. No evident correlation was found between the amplitude motions of the complementary residues. The 5-terminal residues showed large internal motions (S²=0.5), which describe the fraying of the double helix. Global examination of the microdynamical parameters S _inf2 ^supf and S _inf2 ^sups along the nucleotide sequence showed that the adenine residues exhibit more restricted fast internal motions (S _inf2 ^supf =0.88 to 0.96) than the others, whereas the measured relaxation rates of the four nucleosides in solution were mainly of dipolar origin. Moreover, the fit of both R(C_z) and R(HzC_z) experimental relaxation rates using an only global correlation time for all the residues, gave evidence of a supplementary relaxation pathway affecting R(C_x,y) for the purine residues in the (53) G₄A₃ and A₁₀A₉T₈T₇ sequences. This relaxation process was analyzed in terms of exchange stemming from motions of the sugar around the glycosidic bond on the millisecond time scale. It should be pointed out that these residues gave evidence of close contacts with the protein in the complex with the lac operator [Boelens et al. (1987) J. Mol. Biol., 193, 213–216] and that these motions could be implied in the lac-operator-lac-repressor recognition process.     

Addressing the overlap problem in the quantitative analysis of two dimensional NMR spectra: Application to <Superscript>15</Superscript>N relaxation measurements

A quantitative analysis of 2D ¹H-¹⁵N spectra is often complicated by resonance overlap. Here a simple method is presented for resolving overlapped correlations by recording 2D projection planes from HNCO data sets. Applications are presented involving the measurement of ¹⁵N T₁ relaxation rates in a high molecular weight protein, malate synthase G, and in a system that exchanges between folded and unfolded states, the drkN SH3 domain. By supplementing relaxation data recorded in the conventional way as a series of 2D ¹H-¹⁵N data sets with a series of a pair of projection planes the number of dynamics probes is increased significantly for both systems studied.     

Perspective: revisiting the field dependence of TROSY sensitivity

The discovery of the TROSY effect (Pervushin et al. in Proc Natl Acad Sci USA 94:12366–12371, 1997) for reducing transverse relaxation and line sharpening through selecting pathways in which dipole–dipole and CSA Hamiltonians partially cancel each other had a tremendous impact on solution NMR studies of macromolecules. Together with the methyl TROSY (Tugarinov and Kay in J Biomol NMR 28:165–172, 2004) it enabled structural and functional studies of significantly larger systems. The optimal field strengths for TROSY have been estimated to be on spectrometers operating around 900 MHz (21.14 T) for the ¹H_N TROSY (Pervushin et al. in Proc Natl Acad Sci USA 94:12366–12371, 1997) while the aromatic ¹³C (¹³C_aro) TROSY is posited to be optimal at around 600 MHz (14.09 T) (Pervushin et al. in J Am Chem Soc 120:6394–6400, 1998b; Pervushin in Q Rev Biophys 33:161–197, 2000). The initial rational was based on the consideration of where the quadratic B₀ field dependences of the TROSY relaxation rates reach a minimum. For sensitivity consideration, however, it is interesting to estimate which field strengths yield the tallest peaks. Recent studies of ¹⁵N-detected TROSYs suggested that maximal peak heights are expected at 1.15 GHz (27.01 T) although the slowest relaxation rates or longest transverse relaxation times T₂ are indeed expected around 900 MHz (21.14 T) (Takeuchi in J Biomol NMR 63:323–331, 2015; Takeuchi et al. in J Biomol NMR 64:143–151, 2016). This was based on the fact that the heights of Lorentzian lines are proportional to B _o ^3/2 * T₂ (B_o). Thus, multiplying the parabolic T₂(B_o) dependence with the increasing function of B _o ^3/2 shifts the maxima of peak-height field dependence from the T₂ maximum at 900 MHz to higher fields. Moreover, besides shifting the peak height maximum for ¹⁵N TROSY, this analysis yields estimates for optimal peak heights for ¹H_N detected TROSY to 1.5 GHz, and to 900 MHz for ¹³C-detected ¹³C_aroTROSY as is detailed below. To our knowledge, this aspect of field dependence of TROSY sensitivity has not been in the attention of the NMR community but may affect perspectives of NMR at ultra-high fields.     

The relative orientation of the fibronectin 6F11F2 module pair: A 15N NMR relaxation study

The structure of a pair of modules (⁶F1¹F2), that forms part of the collagen-binding region of fibronectin, is refined using heteronuclear relaxation data. A structure of the pair was previously derived from ¹H-¹H NOE and ³ J _HHN data [Bocquier et al. (1999) Structure, 7, 1451–1460] and a weak module–module interface, comprising Leu19 and Leu28, in ⁶F1, and Tyr68 in ²F1, was identified. In this study, the definition of the average relative orientation of the two modules is improved using the dependence of ¹⁵N relaxation on rotational diffusion anisotropy. This structure refinement is based on the selection of a subset of structures from sets calculated with NOE and ³ J _HHN data alone, using the quality of the fits to the relaxation data as the selection criterion. This simple approach is compared to a refinement strategy where ¹⁵N relaxation data are included in the force field as additional restraints [Tjandra et al. (1997) Nat. Struct. Biol., 4, 443–449].     

Characterization of the overall and local dynamics of a protein with intermediate rotational anisotropy: Differentiating between conformational exchange and anisotropic diffusion in the B3 domain of protein G

Because the overall tumbling provides a major contribution to protein spectral densities measured in solution, the choice of a proper model for this motion is critical for accurate analysis of protein dynamics. Here we study the overall and backbone dynamics of the B3 domain of protein G using ¹⁵N relaxation measurements and show that the picture of local motions is markedly dependent on the model of overall tumbling. The main difference is in the interpretation of the elevated R ₂ values in the -helix: the isotropic model results in conformational exchange throughout the entire helix, whereas no exchange is predicted by anisotropic models that place the longitudinal axis of diffusion tensor almost parallel to the helix axis. Due to small size (fast tumbling) of the protein, the T ₁ values have low sensitivity to NH bond orientation. The diffusion tensor derived from orientation dependence of R ₂/R ₁ is anisotropic (D _par/D _perp=1.4), with a small rhombic component. In order to distinguish the correct picture of motion, we apply model-independent methods that are sensitive to conformational exchange and do not require knowledge of protein structure or assumptions about its dynamics. A comparison of the CSA/dipolar cross-correlation rate constants with ¹⁵N relaxation rates and the estimation of R _ex terms from relaxation data at 9.4 and 14.1 T indicate no conformational exchange in the helix, in support of the anisotropic models. The experimentally derived diffusion tensor is in excellent agreement with theoretical predictions from hydrodynamic calculations; a detailed comparison with various hydrodynamic models revealed optimal parameters for hydrodynamic calculations.     

Developmental changes in the calcium currents in embryonic chick ventricular myocytes

Summary Using the patch-clamp technique, we recorded whole-cell calcium current from isolated cardiac myocytes dissociated from the apical ventricles of 7-day and 14-day chick embryos. In 70% of 14-day cells after 24 hr in culture, two component currents could be separated from totalI _Ca activated from a holding potential (V _h) of –80 mV. L-type current (I _L) was activated by depolarizing steps fromV _h –30 or –40 mV. The difference current (I _T) was obtained by subtractingI _L, fromI _Ca.I _T could also be distinguished pharmacologically fromI _L in these cells.I _T was selectively blocked by 40–160 m Ni²⁺, whereasI _L was suppressed by 1 m D600 or 2 m nifedipine. The Ni²⁺-resistant and D600-resistant currents had activation thresholds and peak voltages that were near those ofI _T andI _L defined by voltage threshold, and resembled those in adult mammalian heart. In 7-day cells,I _T andI _L could be distinguished by voltage threshold in 45% (S cells), while an additional 45% of 7-day cells were nonseparable (NS) by activation voltage threshold. Nonetheless, in mostNS cells,I _Ca was partly blocked by Ni²⁺ and by D600 given separately, and the effects were additive when these agents were given together. Differences among the cells in the ability to separateI _T andI _L by voltage threshold resulted largely from differences in the position of the steady-state inactivation and activation curves along the voltage axis. In all cells at both ages in which the steady-state inactivation relation was determined with a double-pulse protocol, the half-inactivation potential (V _1/2) of the Ni²⁺-resistant currentI _L averaged –18 mV. In contrast,V _1/2 of the Ni²⁺-sensitiveI _T was –60 mV in 14-day cells, –52 mV in 7-dayS cells, and –43 mV in 7-day NS cells. The half-activation potential was near –2 mV forI _L at both ages, but that ofI _T was –38 mV in 14-day and –29 mV in 7-day cells. Maximal current density was highly variable from cell to cell, but showed no systematic differences between 7-day and 14-day cells. These results indicate that the main developmental change that occurs in the components ofI _Ca is a negative shift with, embryonic age in the activation and inactivation relationships ofI _T along the voltage axis.     

[31P]-Nuclear magnetic resonance spin lattice relaxation in lecithin reverse micelles

[³¹P] -Nuclear magnetic resonance (NMR) spin lattice relaxation times (T₁) have been measured for lecithin-nonpolar solvent-water as a function of added water for three solvents, namely, benzene, carbon tetrachloride and cyclohexane. In benzene and carbon tetrachloride systems, where spherical reverse micelles are formed, [³¹P]-NMR T₁, values increase linearly with added water. However, in cyclohexane, the trends in the [³¹P]-T₁ values indicate very different micellisation processes. Even at the lowest concentration of added water, the [³¹P]-T₁ values in this solvent are substantially larger than the corresponding values in benzene and carbon tetrachloride, which is attributed to the intramolecular chlorinephosphate interaction being the weakest in cyclohexane. At a higher water content of six mols of water per mol of lecithin in cyclohexane solvent, the [³¹P]-T₁ values show a sharp decrease indicating a sudden change in the dynamics of the phosphate group, and this confirms the on set of ‘reverse micelle-to-liquid crystalline’ phase transition observed in this system by other spectroscopic and physical techniques.     

A C-detection NMR approach for large glycoproteins

NMR spectroscopy has great potential to provide us with information on structure and dynamics at atomic resolution of glycoproteins in solution. In larger glycoproteins, however, the detrimental fast ¹H transverse relaxation renders the conventional ¹H-detected NMR experiments difficult. ¹³C direct detection potentially offers a valuable alternative to ¹H detection to overcome the fast T₂ relaxation. Here, we applied ¹³C-detected NMR methods to observe the NMR signals of ¹³C-labeled glycans attached to the Fc fragment of immunoglobulin G with a molecular mass of 56 kDa. Spectral analysis revealed that a ¹³C-detected ¹³C-¹³C NOESY experiment is highly useful for spectral assignments of the glycans of large glycoproteins. This approach would be, in part, complementary to ¹³C-¹³C TOCSY and ¹H-detection experiments.     

The role of ion antiporters in the maintenance of intracellular pH in rat vascular smooth muscle cells

Vascular smooth muscle intracellular pH is maintained by the Na⁺/H⁺ and Cl^–/HCO ₃ ^– antiporters. The Na⁺/H⁺ exchanger is a major route of H⁺ extrusion in most eukaryotic cells and is present in vascular smooth muscle cells in a similar capacity. It extrudes H^– into the extracellular space in exchange for Na⁺. The Cl^–/HCO ₃ ^– exchanger plays an analogous role to lower the pH of vascular smooth muscle cells when increases in intracellular pH occur. Its activity has also been demonstrated in A7r5 and A10 vascular smooth muscle cells. The Na⁺/H⁺ exchanger is regulated by a number of agents which act through inositol trisphosphate/diacylglycerol, to stimulate the antiporter. Calcium-calmodulin dependent protein kinase may also activate the antiporter in vivo. Phosphorylation of the Cl^–/HCO ₃ ^– exchanger has also been observed but its physiological role is not known. Both these antiporters exist in the plasma membrane as integral proteins with free acidic cytoplasmic termini. These regions may be important in sensing changes in intracellular pH, to which these antiporters respond.Abbreviations CaM Calmodulin - DCCD Dicylohexyl-Carbodiimide - DG Diacylglycerol - DIDS-4 4-Diisthiocyanostilbene-2,2-Disulfonic Acid - IP₃ Inositol Trisphosphate - PKC protein Kinase C - SITS-4 4-Acetamido-4-Isothiocyanstilbene-2,2-Disulfonate - VSMC Vascular Smooth Muscle Cell     

Backbone dynamics of an alamethicin in methanol and aqueous detergent solution determined by heteronuclear 1H−15N NMR spectroscopy

Summary The ¹⁵N relaxation rates of the -aminoisobutyric acid (Aib)-rich peptide alamethicin dissolved in methanol at 27°C and 5°C, and dissolved in aqueous sodium dodecylsulfate (SDS) at 27°C, were measured using inverse-detected one-and two-dimensional ¹H–¹⁵N NMR spectroscopy. Measurements of ¹⁵N longitudinal (R_N(N_z)) and transverse (R_N(N_x,y)) relaxation rates and the {¹H} ¹⁵N nuclear Overhauser enhancement (NOE) at 11.7 Tesla were used to calculate (quasi-) spectral density values at 0, 50, and 450 MHz for the peptide in methanol and in SDS. Spectral density mapping at 0, 50, 450, 500, and 550 MHz was done using additional measurements of the ¹H–¹⁵N lingitudinal two-spin order, R_NH(2H _infZ ^supN N_Z), two-spin antiphase coherence, R_NH(2H _infN ^supZ N_x,y), and the proton longitudinal relaxation rate, R_H(H _infN ^supZ ), for the peptide dissolved in methanol only. The spectral density of motions was also modeled using the three-parameter Lipari-Szabo function. The overall rotational correlation times were determined to be 1.1, 2.5, and 5.7 ns for alamethicin in methanol at 27°C and 5°C, and in SDS at 27°C, respectively. From the rotational correlation time determined in SDS the number of detergent molecules associated with the peptide was estimated to be about 40. The average order parameter was about 0.7 and the internal correlation times were about 70 ps for the majority of backbone amide ¹⁵N sites of alamethicin in methanol and in SDS. The relaxation data, spectral densities, and order parameters suggest that the peptide N-H vectors of alamethicin are not as highly constrained as the core regions of folded globular proteins. However, the peptide backbone is clearly not as mobile as the most unconstrained regions of folded proteins, such as those found in the frayed C-and N-termini of some proteins, or in randomcoil peptides. The data also suggest significant mobility at both ends of the peptide dissolved in methanol. In SDS the mobility in the middle and at the ends of the peptide is reduced. The implications of the results with respect to the sterically hindered Aib residues and the biological activities of the peptide are discussed.To whom correspondence should be addressed.