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.     

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.     

Characterization of the mobility of various chemical groups in the purple membrane of halobacterium halobium by 13C, 31P and 2H solid state N M R

Lyophilized purple membrane sheets have been investigated by C-13- and P-31-cross polarization/magic angle spinning N M R spectroscopy. The high-resolution C-13 spectrum and its non-quaternary suppression version indicate fast protein side-chain motions but a rigid backbone structure on a time scale of roughly < 0.001 to 0.01 s. Three components of exchangeable hydrogen have been detected by deuterium N M R. The mean exchange time of the peptide hydrogens must be longer than 1 μs. The medium component is attributed to mobile side-chains. In addition a narrow line has been observed which is assigned to the residual hydration water.     

Analysis of 31P nuclear magnetic resonance lineshapes and transversal relaxation of bacteriophage M13 and tobacco mosaic virus.

The experimentally observed 31P lineshapes and transversal relaxation of 15% (wt/wt) M13, 30% M13, and 30% tobacco mosaic virus (TMV) are compared with lineshapes and relaxation curves that are simulated for various types of rotational diffusion using the models discussed previously (Magusin, P. C. M. M., and M. A. Hemminga. 1993. Biophys. J. 64:1851-1860). It is found that isotropic diffusion cannot explain the observed lineshape effects. A rigid rod diffusion model is only successful in describing the experimental data obtained for 15% M13. For 30% M13 the experimental lineshape and relaxation curve cannot be interpreted consistently and the TMV lineshape cannot even be simulated alone, indicating that the rigid rod diffusion model does not generally apply. A combined diffusion model with fast isolated motions of the encapsulated nucleic acid dominating the lineshape and a slow overall rotation of the virion as a whole, which mainly is reflected in the transversal relaxation, is able to provide a consistent picture for the 15 and 30% M13 samples, but not for TMV. Strongly improved lineshape fits for TMV are obtained assuming that there are three binding sites with different mobilities. The presence of three binding sites is consistent with previous models of TMV. The best lineshapes are simulated for a combination of one mobile and two static sites. Although less markedly, the assumption that two fractions of DNA with different mobilities exist within M13 also improves the simulated lineshapes. The possible existence of two 31P fractions in M13 sheds new light on the nonintegral ratio 2.4:1 between the number of nucleotides and protein coat subunits in the phage: 83% of the viral DNA is less mobile, suggesting that the binding of the DNA molecule to the protein coat actually occurs at the integral ratio of two nucleotides per protein subunit.     

New approach to study fast and slow motions in lipid bilayers: application to dimyristoylphosphatidylcholine-cholesterol interactions.

Natural abundance 13C solid-state nuclear magnetic resonance spectroscopy was used to investigate the effect of the incorporation of cholesterol on the dynamics of dimyristoylphosphatidylcholine (DMPC) bilayers in the liquid-crystalline phase. In particular, the use of a combination of the cross-polarization and magic angle spinning techniques allows one to obtain very high resolution spectra from which can be distinguished several resonances attributed to the polar head group, the glycerol backbone, and the acyl chains of the lipid molecule. To examine both the fast and slow motions of the lipid bilayers, 1H spin-lattice relaxation times as well as proton and carbon spin-lattice relaxation times in the rotating frame were measured for each resolved resonance of DMPC. The use of the newly developed ramped-amplitude cross-polarization technique results in a significant increase in the stability of the cross-polarization conditions, especially for molecular groups undergoing rapid motions. The combination of T1 and T1 rho measurements indicates that the presence of cholesterol significantly decreases the rate and/or amplitude of both the high and low frequency motions in the DMPC bilayers. This effect is particularly important for the lipid acyl chains and the glycerol backbone region.     

Secondary structure and backbone dynamics of Escherichia coli diacylglycerol kinase,as revealed by site-directed solid-state 13C NMR

To gain insight into secondary structure and backbone dynamics, we have recorded ¹³C NMR spectra of [3-¹³C]Ala-, [1-¹³C]Val-labeled Escherichia coli diacylglycerol kinase (DGK), using cross-polarization-magic angle spinning (CP-MAS) and single-pulse excitation with dipolar decoupled-magic angle spinning (DD-MAS) methods. DGK was either solubilized in n-decyl-β-maltoside (DM) micelle or integrated into 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) or 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) bilayers. Surprisingly, the ¹³C NMR spectra were broadened to yield rather featureless peaks at physiological temperatures, both in DM solution or lipid bilayers at liquid crystalline phase, due to interference of motional frequencies of DGK with frequencies of magic angle spinning (MAS) or proton decoupling (10⁴ or 10⁵ Hz, respectively). In gel phase lipids, however, up to six distinct ¹³C NMR peaks were well-resolved due to lowered fluctuation frequencies (<10⁵ Hz) for the transmembrane region, the amphipathic α-helices and loops. While DGK can be tightly packed in gel phase lipids, DGK is less tightly packed at physiological temperatures, where it becomes more mobile. The fact that the enzymatic activity is low under conditions where motion is restricted and high when conformational fluctuations can occur suggests that acquisition of low frequency backbone motions, on the microsecond to millisecond time scale, may facilitate the efficient enzymatic activity of DGK.     

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.     

Protein dynamics detected in a membrane-embedded potassium channel using two-dimensional solid-state NMR spectroscopy

We report longitudinal ¹⁵N relaxation rates derived from two-dimensional (¹⁵N, ¹³C) chemical shift correlation experiments obtained under magic angle spinning for the potassium channel KcsA-Kv1.3 reconstituted in multilamellar vesicles. Thus, we demonstrate that solid-state NMR can be used to probe residue-specific backbone dynamics in a membrane-embedded protein. Enhanced backbone mobility was detected for two glycine residues within the selectivity filter that are highly conserved in potassium channels and that are of core relevance to the filter structure and ion selectivity.     

Irreversible conformational change of bacterio-opsin induced by binding of retinal during its reconstitution to bacteriorhodopsin, as studied by (13)C NMR

We compared (13)C NMR spectra of [3-(13)C]Ala- and [1-(13)C]Val-labeled bacterio-opsin (bO), produced either by bleaching bR with hydroxylamine or from a retinal-deficient strain, with those of bacteriorhodopsin (bR), in order to gain insight into the conformational changes of the protein backbone that lead to correct folding after retinal is added to bO. The observed (13)C NMR spectrum of bO produced by bleaching is not greatly different from that of bR, except for the presence of suppressed or decreased peak-intensities. From careful evaluation of the intensity differences between cross polarization magic angle spinning (CP-MAS) and dipolar decoupled-magic angle spinning (DD-MAS) spectra, it appears that the reduced peak-intensities arise from reduced efficiency of cross polarization or interference of internal motions with proton decoupling frequencies. In particular, the E-F and F-G loops and some transmembrane helices of the bleached bO have acquired internal motions whose frequencies interfere with proton decoupling frequencies. In contrast, the protein backbone of the bO from the retinal-negative cells is incompletely folded. Although it contains mainly a-helices, its very broad (13)C NMR signals indicate that its tertiary structure is different from bR. Importantly, this changed structure is identical in form to that of bleached bO from wild-type bR after it was regenerated with retinal in vitro, and bleached with hydroxylamine. We conclude that the binding of retinal is essential for the correct folding of bR after it is inserted in vitro into the lipid bilayer, and the final folded state does not revert to the partially folded form upon removal of the retinal.     

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.     

Internal motional amplitudes and correlated bond rotations in an alpha-helical peptide derived from 13C and 15N NMR relaxation

Peptide GFSKAELAKARAAKRGGY folds in an alpha-helical conformation that is stabilized by formation of a hydrophobic staple motif and an N-terminal capping box (Munoz V. Blanco FJ, Serrano L, 1995, Struct Biol 2:380-385). To investigate backbone and side-chain internal motions within the helix and hydrophobic staple, residues F2, A5, L7, A8, and A10 were selectively 13C- and 15N-enriched and NMR relaxation experiments were performed in water and in water/trifluoroethanol (TFE) solution at four Larmor frequencies (62.5, 125, 150, and 200 MHz for 13C). Relaxation data were analyzed using the model free approach and an anisotropic diffusion model. In water, angular variances of motional vectors range from 10 to 20 degrees and backbone phi,psi bond rotations for helix residues A5, L7, A8, and A10 are correlated indicating the presence of Calpha-H, Calpha-Cbeta, and N-H rocking-type motions along the helix dipole axis. L7 side-chain CbetaH2 and CgammaH motions are also correlated and as motionally restricted as backbone CalphaH, suggesting considerable steric hindrance with neighboring groups. In TFE which stabilizes the fold, internal motional amplitudes are attenuated and rotational correlations are increased. For the side chain of hydrophobic staple residue F2, wobbling-in-a-cone type motions dominate in water, whereas in TFE, the Cbeta-Cgamma bond and phenyl ring fluctuate more simply about the Calpha-Cbeta bond. These data support the Daragan-Mayo model of correlated bond rotations (Daragan VA, Mayo KH, 1996, J Phys Chem 100:8378-8388) and contribute to a general understanding of internal motions in peptides and proteins.     

NMR of bicelles: orientation and mosaic spread of the liquid-crystal director under sample rotation

Model-membrane systems composed of liquid-crystalline bicellar phases can be uniaxially oriented with respect to a magnetic field, thereby facilitating structural and dynamics studies of membrane-associated proteins. Here we quantitatively characterize a method that allows the manipulation of the direction of this uniaxial orientation. Bicelles formed from DMPC/DHPC are examined by ³¹P NMR under variable-angle sample-spinning (VAS) conditions, confirming that the orientation of the liquid-crystalline director can be influenced by sample spinning. The director is perpendicular to the rotation axis when (the angle between the sample-spinning axis and the magnetic field direction) is smaller than the magic angle, and is parallel to the rotation axis when is larger than the magic angle. The new ³¹P NMR VAS data presented are considerably more sensitive to the orientation of the bicelle than earlier ²H studies and the analysis of the sideband pattern allows the determination of the orientation of the liquid-crystal director and its variation over the sample, i.e., the mosaic spread. Under VAS, the mosaic spread is small if deviates significantly from the magic angle but becomes very large at the magic angle.     

Phosphorus-31 Transverse Relaxation Rate Measurements by NMR Spectroscopy: Insight into Conformational Exchange Along the Nucleic Acid Backbone

Phosphorus ((31)P) NMR spectroscopy can provide important information about the dynamics of nucleic acids. In this communication, we propose an inversely detected (31)P transverse relaxation rate ( R (2)) measurement experiment. This experiment enables fast measurement of accurate (31)P transverse relaxation rates and provides the possibility to detect slow motions mapped by the phosphorus nuclei along the nucleic acid backbone. Dispersion curves show some (31)P nuclei experiencing chemical exchange in the millisecond time scale. Under the assumption of a two-state exchange process, the reduced lifetimes of the exchanging sites (tau(ex)) obtained are in accordance with base pair lifetime estimates deduced from imino proton exchange measurements.     

Molecular organization and motions of crystalline monoacylglycerols and diacylglycerols: a C-13 MASNMR study.

Six saturated acylglycerols (1-myristoyl-sn-glycerol, 1-palmitoyl-sn-glycerol, 1,2-dimyristoyl-sn-glycerol, 1,2-dipalmitoyl-sn-glycerol, 1,2-dipalmitoyl-rac-glycerol, and 1,3-dimyristoylglycerol) were studied in their various polymorphic forms (sub-alpha, alpha, beta') by natural abundance C-13 nuclear magnetic resonance (NMR) with magic angle spinning (MASNMR). C-13 MASNMR does not require single crystals and can observe relatively disordered crystals, distinct advantages over crystallographic diffraction methods. Well resolved spectra were obtained for each acylglycerol, and the chemical shifts of corresponding carbons were different for each crystalline phase and the isotropic liquid phase; moreover, in the case of monoacylglycerols, the symmetrically nonequivalent molecules in the same crystalline structure gave distinct C-13 resonances for the same carbon. The C-13 chemical shifts corresponding to each polymorphic phase were interpreted in terms of differences in intramolecular bond distances, intermolecular interactions (such as H bonding), and molecular motions. Mobilities of the glycerol backbone and acyl chains were assessed by the C-13 linewidths and the C-H dipolar relaxation rates. The chemical shift anisotropy(ies) (delta sigma) of the carbonyl group(s) of each acylglycerol was determined from slow-spinning MAS spectra, and was discussed in terms of the conformational and/or motional changes for the carbonyl carbon(s).     

Conformation and dynamics of melittin bound to magnetically oriented lipid bilayers by solid-state (31)P and (13)C NMR spectroscopy

The conformation and dynamics of melittin bound to the dimyristoylphosphatidylcholine (DMPC) bilayer and the magnetic orientation in the lipid bilayer systems were investigated by solid-state (31)P and (13)C NMR spectroscopy. Using (31)P NMR, it was found that melittin-lipid bilayers form magnetically oriented elongated vesicles with the long axis parallel to the magnetic field above the liquid crystalline-gel phase transition temperature (T(m) = 24 degrees C). The conformation, orientation, and dynamics of melittin bound to the membrane were further determined by using this magnetically oriented lipid bilayer system. For this purpose, the (13)C NMR spectra of site-specifically (13)C-labeled melittin bound to the membrane in the static, fast magic angle spinning (MAS) and slow MAS conditions were measured. Subsequently, we analyzed the (13)C chemical shift tensors of carbonyl carbons in the peptide backbone under the conditions where they form an alpha-helix and reorient rapidly about the average helical axis. Finally, it was found that melittin adopts a transmembrane alpha-helix whose average axis is parallel to the bilayer normal. The kink angle between the N- and C-terminal helical rods of melittin in the lipid bilayer is approximately 140 degrees or approximately 160 degrees, which is larger than the value of 120 degrees determined by x-ray diffraction studies. Pore formation was clearly observed below the T(m) in the initial stage of lysis by microscope. This is considered to be caused by the association of melittin molecules in the lipid bilayer.     

Alteration of conformation and dynamics of bacteriorhodopsin induced by protonation of Asp 85 and deprotonation of Schiff base as studied by 13C NMR

According to previous X-ray diffraction studies, the D85N mutant of bacteriorhodopsin (bR) with unprotonated Schiff base assumes a protein conformation similar to that in the M photointermediate. We recorded (13)C NMR spectra of [3-(13)C]Ala- and [1-(13)C]Val-labeled D85N and D85N/D96N mutants at ambient temperature to examine how conformation and dynamics of the protein backbone are altered when the Schiff base is protonated (at pH 7) and unprotonated (at pH 10). Most notably, we found that the peak intensities of three to four [3-(13)C]Ala-labeled residues from the transmembrane alpha-helices, including Ala 39, 51, and 53 (helix B) and 215 (helix G), were suppressed in D85N and D85N/D96N both from CP-MAS (cross polarization-magic angle spinning) and DD-MAS (dipolar decoupled-magic angle spinning) spectra, irrespective of the pH. This is due to conformational change and subsequent acquisition of intermediate time-range motions, with correlation times in the order of 10(-)(5) or 10(-)(4) s, which interferes with proton decoupling frequency or frequency of magic angle spinning, respectively, essential for an attempted peak-narrowing to achieve high-resolution NMR signals. Greater changes were achieved, however, at pH 10, which indicate large-amplitude motions of transmembrane helices upon deprotonation of Schiff base and the formation of the M-like state in the absence of illumination. The spectra detected more rapid motions in the extracellular and/or cytoplasmic loops, with correlation times increasing from 10(-)(4) to 10(-)(5) s. Conformational changes in the transmembrane helices were located at helices B, G, and D as viewed from the above-mentioned spectral changes, as well as at 1-(13)C-labeled Val 49 (helix B), 69 (B-C loop), and [3-(13)C]Ala-labeled Ala 126 (D-helix) signals, in addition to the cytoplasmic and extracellular loops. Further, we found that in the M-like state the charged state of Asp 96 at the cytoplasmic side substantially modulated the conformation and dynamics of the extracellular region through long-distance interaction.     

Identifying anisotropic constraints in multiply labeled bacteriorhodopsin by 15N MAOSS NMR: a general approach to structural studies of membrane proteins

Structural models of membrane proteins can be refined with sets of multiple orientation constraints derived from structural NMR studies of specifically labeled amino acids. The magic angle oriented sample spinning (MAOSS) NMR approach was used to determine a set of orientational constraints in bacteriorhodopsin (bR) in the purple membrane (PM). This method combines the benefits of magic angle spinning (MAS), i.e., improved sensitivity and resolution, with the ability to measure the orientation of anisotropic interactions, which provide important structural information. The nine methionine residues in bacteriorhodopsin were isotopically (15)N labeled and spectra simplified by deuterium exchange before cross-polarization magic angle spinning (CPMAS) experiments. The orientation of the principal axes of the (15)N chemical shift anisotropy (CSA) tensors was determined with respect to the membrane normal for five of six residual resonances by analysis of relative spinning sideband intensities. The applicability of this approach to large proteins embedded in a membrane environment is discussed in light of these results.     

The chemical and dynamical influence of the anti-viral drug amantadine on the M2 proton channel transmembrane domain

The M(2) proton channel plays a vital role in the life cycle of the influenza A virus. His(37), the key residue in the M(2) transmembrane domain (M(2)-TMD), plays a central role in the proton conductance mechanism. The anti-influenza drug, amantadine, inhibits the channel activity through binding to the pore of the M(2) channel. The nuclear spin relaxation data and polarization inversion spin exchange at the magic angle spectra show that both the polypeptide backbone and His(37) side chain are more constrained in the presence of amantadine. Using (15)N cross polarization magic-angle spinning NMR spectroscopy, the protonation of His(37) of M(2)-TMD in lipid bilayers was monitored in the absence and presence of amantadine as a function of pH. Binding amantadine lowers the His(37) pK(a) values by approximately three orders of magnitude compared with the first pK(a) of histidine in amantadine-free M(2)-TMD. Amantadine's influence on the His(37) chemical properties suggests a novel mechanism by which amantadine may inhibit proton conductance.     

Repeat motions and backbone flexibility in designed proteins with different numbers of identical consensus tetratricopeptide repeats

The tetratricopeptide repeat (TPR) is a 34-residue helix-turn-helix motif that occurs as three or more tandem repeats in a wide variety of proteins. We have determined the repeat motions and backbone fluctuations of proteins containing two or three consensus TPR repeats (CTPR2 and CPTR3, respectively) using 15N NMR relaxation measurements. Rotational diffusion tensors calculated from these data for each repeat within each TPR protein indicate that there is a high degree of motional correlation between different repeats in the same protein. This is consistent with the prevailing view that repeat proteins, such as CTPR2 and CTPR3, behave as single cooperatively folded domains. The internal motions of backbone NH groups were determined using the Lipari-Szabo model-free formalism. For most residues, there was a clear separation between the influence of internal motion and the influence of global rotational tumbling on the observed magnetic relaxation. The local internal motions are highly restricted in most of the helical elements, with slightly greater flexibility in the linker elements. Comparisons between CTPR2 and CTPR3 indicate that an addition of a TPR repeat to the C-terminus (before the solvation helix) of CTPR2 slightly reduces the flexibility of the preceding helix.