Deuterium solid-state nuclear magnetic resonance studies of methyl group dynamics in bacteriorhodopsin and retinal model compounds: evidence for a 6-s-trans chromophore in the protein

Solid-state NMR spectroscopy

Solid-state nuclear magnetic resonance (NMR) spectroscopy is an atomic-level method used to determine the chemical structure, three-dimensional structure, and dynamics of solids and semi-solids. This Primer summarizes the basic principles of NMR as applied to the wide range of solid systems. The fundamental nuclear spin interactions and the effects of magnetic fields and radiofrequency pulses on nuclear spins are the same as in liquid-state NMR. However, because of the anisotropy of the interactions in the solid state, the majority of high-resolution solid-state NMR spectra is measured under magic-angle spinning (MAS), which has profound effects on the types of radiofrequency pulse sequences required to extract structural and dynamical information. We describe the most common MAS NMR experiments and data analysis approaches for investigating biological macromolecules, organic materials, and inorganic solids. Continuing development of sensitivity-enhancement approaches, including 1H-detected fast MAS experiments, dynamic nuclear polarization, and experiments tailored to ultrahigh magnetic fields, is described. We highlight recent applications of solid-state NMR to biological and materials chemistry. The Primer ends with a discussion of current limitations of NMR to study solids, and points to future avenues of development to further enhance the capabilities of this sophisticated spectroscopy for new applications.

Dynamic Nuclear Polarization of Biomembrane Assemblies

While atomic scale structural and dynamic information are hallmarks of nuclear magnetic resonance (NMR) methodologies, sensitivity is a fundamental limitation in NMR studies. Fully exploiting NMR capabilities to study membrane proteins is further hampered by their dilution within biological membranes. Recent developments in dynamic nuclear polarization (DNP), which can transfer the relatively high polarization of unpaired electrons to nuclear spins, show promise for overcoming the sensitivity bottleneck and enabling NMR characterization of membrane proteins under native-like conditions. Here we discuss fundamental aspects of DNP-enhanced solid-state NMR spectroscopy, experimental details relevant to the study of lipid assemblies and incorporated proteins, and sensitivity gains which can be realized in biomembrane-based samples. We also present unique insights which can be gained from DNP measurements and prospects for further development of the technique for elucidating structures and orientations of membrane proteins in native lipid environments.

Concepts and Methods of Solid-State NMR Spectroscopy Applied to Biomembranes

Concepts of solid-state NMR spectroscopy and applications to fluid membranes are reviewed in this paper. Membrane lipids with ²H-labeled acyl chains or polar head groups are studied using ²H NMR to yield knowledge of their atomistic structures in relation to equilibrium properties. This review demonstrates the principles and applications of solid-state NMR by unifying dipolar and quadrupolar interactions and highlights the unique features offered by solid-state ²H NMR with experimental illustrations. For randomly oriented multilamellar lipids or aligned membranes, solid-state ²H NMR enables direct measurement of residual quadrupolar couplings (RQCs) due to individual C-²H-labeled segments. The distribution of RQC values gives nearly complete profiles of the segmental order parameters S_CD⁽ⁱ⁾ as a function of acyl segment position (i). Alternatively, one can measure residual dipolar couplings (RDCs) for natural abundance lipid samples to obtain segmental S_CH order parameters. A theoretical mean-torque model provides acyl-packing profiles representing the cumulative chain extension along the normal to the aqueous interface. Equilibrium structural properties of fluid bilayers and various thermodynamic quantities can then be calculated, which describe the interactions with cholesterol, detergents, peptides, and integral membrane proteins and formation of lipid rafts. One can also obtain direct information for membrane-bound peptides or proteins by measuring RDCs using magic-angle spinning (MAS) in combination with dipolar recoupling methods. Solid-state NMR methods have been extensively applied to characterize model membranes and membrane-bound peptides and proteins, giving unique information on their conformations, orientations, and interactions in the natural liquid-crystalline state.

Protein dynamics in the solid state from 2H NMR line shape analysis: a consistent perspective

Deuterium line shape analysis of CD3 groups has emerged as a particularly useful tool for studying microsecond-millisecond protein motions in the solid state. The models devised so far consist of several independently conceived simple jump-type motions. They are comprised of physical quantities encoded in their simplest form; improvements are only possible by adding yet another simple motion, thereby changing the model. The various treatments developed are case-specific; hence comparison among the different systems is not possible. Here we develop a new methodology for (2)H NMR line shape analysis free of these limitations. It is based on the microscopic-order-macroscopic-disorder (MOMD) approach. In MOMD motions are described by diffusion tensors, spatial restrictions by potentials/ordering tensors, and geometric features by relative tensor orientations. Jump-type motions are recovered in the limit of large orientational potentials. Model improvement is accomplished by monitoring the magnitude, symmetry, and orientation of the various tensors. The generality of MOMD makes possible comparison among different scenarios. CD3 line shapes from the Chicken Villin Headpiece Subdomain and the Streptomyces Subtilisin Inhibitor are used as experimental examples. All of these spectra are reproduced by using rhombic local potentials constrained for simplicity to be given by the L = 2 spherical harmonics, and by axial diffusion tensors. Potential strength and rhombicity are found to be ca. 2-3 k(B)T. The diffusion tensor is tilted at 120° from the C-CD3 axis. The perpendicular (parallel) correlation times for local motion are 0.1-1.0 ms (3.3-30 μs). Activation energies in the 1.1-8.0 kcal/mol range are estimated. Future prospects include extension to the (2)H relaxation limit, application to the (15)N and (13)C NMR nuclei, and accounting for collective motions and anisotropic media.

Temperature dependence of internal motions of protein side-chain NH3(+) groups: insight into energy barriers for transient breakage of hydrogen bonds

Although charged side chains play important roles in protein function, their dynamic properties are not well understood. Nuclear magnetic resonance methods for investigating the dynamics of lysine side-chain NH3(+) groups were established recently. Using this methodology, we have studied the temperature dependence of the internal motions of the lysine side-chain NH3(+) groups that form ion pairs with DNA phosphate groups in the HoxD9 homeodomain-DNA complex. For these NH3(+) groups, we determined order parameters and correlation times for bond rotations and reorientations at 15, 22, 28, and 35 °C. The order parameters were found to be virtually constant in this temperature range. In contrast, the bond-rotation correlation times of the NH3(+) groups were found to depend strongly on temperature. On the basis of transition state theory, the energy barriers for NH3(+) rotations were analyzed and compared to those for CH3 rotations. Enthalpies of activation for NH3(+) rotations were found to be significantly higher than those for CH3 rotations, which can be attributed to the requirement of hydrogen bond breakage. However, entropies of activation substantially reduce the overall free energies of activation for NH3(+) rotations to a level comparable to those for CH3 rotations. This entropic reduction in energy barriers may accelerate molecular processes requiring hydrogen bond breakage and play a kinetically important role in protein function.

Quantifying conformational dynamics using solid-state R₁ρ experiments

We demonstrate the determination of quantitative rates of molecular reorientation in the solid state with rotating frame (R(1ρ)) relaxation measurements. Reorientation of the carbon chemical shift anisotropy (CSA) tensor was used to probe site-specific conformational exchange in a model system, d(6)-dimethyl sulfone (d(6)-DMS). The CSA as a probe of exchange has the advantage that it can still be utilized when there is no dipolar mechanism (i.e. no protons attached to the site of interest). Other works have presented R(1ρ) measurements as a general indicator of dynamics, but this study extracts quantitative rates of molecular reorientation from the R(1ρ) values. Some challenges of this technique include precise knowledge of sample temperature and determining the R(2)(0) contribution to the observed relaxation rate from interactions other than molecular reorientation, such as residual dipolar couplings or fast timescale dynamics; determination of this term is necessary in order to quantify the exchange rate due to covariance between the 2 terms. Low-temperature experiments measured an R(2)(0) value of 1.8±0.2s(-1) Allowing for an additional relaxation term (R(2)(0)), which was modeled as both temperature-dependent and temperature-independent, rates of molecular reorientation were extracted from field strength-dependent R(1ρ) measurements at four different temperatures and the activation energy was determined from these exchange rates. The activation energies determined were 74.7±4.3kJ/mol and 71.7±2.9kJ/mol for the temperature-independent and temperature-dependent R(2)(0) models respectively, in excellent agreement with literature values. The results of this study suggest important methodological considerations for the application of the method to more complicated systems such as proteins, such as the importance of deuterating samples and the need to make assumptions regarding the R(2)(0) contribution to relaxation.

Sensitivity of 2H NMR spectroscopy to motional models: proteins and highly viscous liquids as examples

In order to study to what extent mechanisms of molecular motion can be unambiguously revealed by (2)H NMR spectroscopy, (2)H spectra for proteins (chicken villin protein headpiece HP36, selectively methyl-deuterated at leucine-69, C(δ) D(3)) and binary systems of high viscosity (benzene-d(6) in tricresyl phosphate) have been carefully analyzed as illustrative examples (the spectra are taken from the literature). In the first case, a model of restricted diffusion mediated by jumps between rotameric orientations has been tested against jump- and free diffusion models which describe rotational motion combined with jump dynamics. It has been found that the set of (2)H spectra of methyl-deuterated at leucine-69 chicken villin protein headpiece HP36 can be consistently explained by different motional models as well as by a gaussian distribution of correlation times assuming isotropic rotation (simple brownian diffusion model). The last finding shows that when the possible distribution of correlation times is not very broad one might not be able to distinguish between heterogeneous and homogenous (but more complex) dynamics by analyzing (2)H lineshapes. For benzene-d(6) in tricresyl phosphate, the dynamics is heterogeneous and it has been demonstrated that a gaussian distribution of correlation times reproduces well the experimental lineshapes, while for a Cole-Davidson distribution the agreement is somewhat worse. For inquires into the sensitivity of quadrupolar NMR spectral analysis (by "quadrupolar NMR spectroscopy we understand NMR spectroscopy of nuclei possessing quadrupole moment), the recently presented theoretical approach [Kruk et al., J. Chem. Phys. 135, 224511 (2011)] has been used as it allows simulating quadrupolar spectra for arbitrary motional conditions by employing the stochastic Liouville equation.

Deuterated peptides and proteins: structure and dynamics studies by MAS solid-state NMR

Perdeuteration and back substitution of exchangeable protons in microcrystalline proteins, in combination with recrystallization from D(2)O-containing buffers, significantly reduce (1)H, (1)H dipolar interactions. This way, amide proton line widths on the order of 20 Hz are obtained. Aliphatic protons are accessible either via specifically protonated precursors or by using low amounts of H(2)O in the bacterial growth medium. The labeling scheme enables characterization of structure and dynamics in the solid-state without dipolar truncation artifacts.

Solid-state 2H NMR relaxation illuminates functional dynamics of retinal cofactor in membrane activation of rhodopsin

Rhodopsin is a canonical member of the family of G protein-coupled receptors, which transmit signals across cellular membranes and are linked to many drug interventions in humans. Here we show that solid-state (2)H NMR relaxation allows investigation of light-induced changes in local ps-ns time scale motions of retinal bound to rhodopsin. Site-specific (2)H labels were introduced into methyl groups of the retinal ligand that are essential to the activation process. We conducted solid-state (2)H NMR relaxation (spin-lattice, T(1Z), and quadrupolar-order, T(1Q)) experiments in the dark, Meta I, and Meta II states of the photoreceptor. Surprisingly, we find the retinylidene methyl groups exhibit site-specific differences in dynamics that change upon light excitation--even more striking, the C9-methyl group is a dynamical hotspot that corresponds to a crucial functional hotspot of rhodopsin. Following 11-cis to trans isomerization, the (2)H NMR data suggest the β-ionone ring remains in its hydrophobic binding pocket in all three states of the protein. We propose a multiscale activation mechanism with a complex energy landscape, whereby the photonic energy is directed against the E2 loop by the C13-methyl group, and toward helices H3 and H5 by the C5-methyl of the β-ionone ring. Changes in retinal structure and dynamics initiate activating fluctuations of transmembrane helices H5 and H6 in the Meta I-Meta II equilibrium of rhodopsin. Our proposals challenge the Standard Model whereby a single light-activated receptor conformation yields the visual response--rather an ensemble of substates is present, due to the entropy gain produced by photolysis of the inhibitory retinal lock.

Retinal dynamics underlie its switch from inverse agonist to agonist during rhodopsin activation

X-ray and magnetic resonance approaches, though central to studies of G protein-coupled receptor (GPCR)-mediated signaling, cannot address GPCR protein dynamics or plasticity. Here we show that solid-state (2)H NMR relaxation elucidates picosecond-to-nanosecond-timescale motions of the retinal ligand that influence larger-scale functional dynamics of rhodopsin in membranes. We propose a multiscale activation mechanism whereby retinal initiates collective helix fluctuations in the meta I-meta II equilibrium on the microsecond-to-millisecond timescale.

Comparative dynamics of leucine methyl groups in FMOC-leucine and in a protein hydrophobic core probed by solid-state deuteron nuclear magnetic resonance over 7-324 K temperature range

Quantitative dynamics of methyl groups in 9-fluorenylmethyloxycarbonyl-leucine (FMOC-leu) have been analyzed and compared with earlier studies of methyl dynamics in chicken villin headpiece subdomain protein (HP36) labeled at L69, a key hydrophobic core position. A combination of deuteron solid-state nuclear magnetic resonance experiments over the temperature range of 7-324 K and computational modeling indicated that while the two compounds show the same modes of motions, there are marked differences in the best-fit parameters of these motions. One of the main results is that the crossover observed in the dynamics of the methyl groups in the HP36 sample at 170 K is absent in FMOC-leu. A second crossover at around 95-88 K is present in both samples. The differences in the behavior of the two compounds suggest that some of the features of methyl dynamics reflect the complexity of the protein hydrophobic core and are not determined solely by local interactions.

Dynamics of reassembled thioredoxin studied by magic angle spinning NMR: snapshots from different time scales

Solid-state NMR spectroscopy can be used to probe internal protein dynamics in the absence of the overall molecular tumbling. In this study, we report (15)N backbone dynamics in differentially enriched 1-73(U-(13)C,(15)N)/74-108(U-(15)N) reassembled thioredoxin on multiple time scales using a series of 2D and 3D MAS NMR experiments probing the backbone amide (15)N longitudinal relaxation, (1)H-(15)N dipolar order parameters, (15)N chemical shift anisotropy (CSA), and signal intensities in the temperature-dependent and (1)H T(2)'-filtered NCA experiments. The spin-lattice relaxation rates R(1) (R(1) = 1/T(1)) were observed in the range from 0.012 to 0.64 s(-1), indicating large site-to-site variations in dynamics on pico- to nanosecond time scales. The (1)H-(15)N dipolar order parameters, <S>, and (15)N CSA anisotropies, delta(sigma), reveal the backbone mobilities in reassembled thioredoxin, as reflected in the average <S> = 0.89 +/- 0.06 and delta(sigma) = 92.3 +/- 5.2 ppm, respectively. From the aggregate of experimental data from different dynamics methods, some degree of correlation between the motions on the different time scales has been suggested. Analysis of the dynamics parameters derived from these solid-state NMR experiments indicates higher mobilities for the residues constituting irregular secondary structure elements than for those located in the alpha-helices and beta-sheets, with no apparent systematic differences in dynamics between the alpha-helical and beta-sheet residues. Remarkably, the dipolar order parameters derived from the solid-state NMR measurements and the corresponding solution NMR generalized order parameters display similar qualitative trends as a function of the residue number. The comparison of the solid-state dynamics parameters to the crystallographic B-factors has identified the contribution of static disorder to the B-factors. The combination of longitudinal relaxation, dipolar order parameter, and CSA line shape analyses employed in this study provides snapshots of dynamics and a new insight on the correlation of these motions on multiple time scales.

Probing the dynamics of a protein hydrophobic core by deuteron solid-state nuclear magnetic resonance spectroscopy

With the goal of investigating dynamical features of hydrophobic cores of proteins over a wide range of temperatures, the chicken villin headpiece subdomain protein (HP36) was labeled at a "single" site corresponding to any one of the two C(delta)D(3) groups of leucine-69, which is located in a key position of the core. The main techniques employed are deuteron NMR quadrupolar echo line shape analysis, and T(1Z) (Zeeman) and T(1Q) (quadrupolar order) relaxation experiments performed at 11.7 and 17.6 T over the temperature range of 112 to 298 K. The experimental data are compared with computer simulations. The deuteron line shapes give an excellent fit to a three-mode motional model that consists of (a) fast three-site rotational jumps about the pseudo C(3) methyl spinning axis, (b) slower reorientation of the spinning axis, described by diffusion along a restricted arc, and (c) large angle jumps between traces of rotameric conformers. Relaxation behavior is described by a phenomenological distribution of activation energies for three-site hops at high temperatures that collapses to a single, distinctly smaller value for lower temperatures.

Retinal conformation and dynamics in activation of rhodopsin illuminated by solid-state H NMR spectroscopy

Solid-state NMR spectroscopy gives a powerful avenue for investigating G protein-coupled receptors and other integral membrane proteins in a native-like environment. This article reviews the use of solid-state (2)H NMR to study the retinal cofactor of rhodopsin in the dark state as well as the meta I and meta II photointermediates. Site-specific (2)H NMR labels have been introduced into three regions (methyl groups) of retinal that are crucially important for the photochemical function of rhodopsin. Despite its phenomenal stability (2)H NMR spectroscopy indicates retinal undergoes rapid fluctuations within the protein binding cavity. The spectral lineshapes reveal the methyl groups spin rapidly about their three-fold (C(3)) axes with an order parameter for the off-axial motion of SC(3) approximately 0.9. For the dark state, the (2)H NMR structure of 11-cis-retinal manifests torsional twisting of both the polyene chain and the beta-ionone ring due to steric interactions of the ligand and the protein. Retinal is accommodated within the rhodopsin binding pocket with a negative pretwist about the C11=C12 double bond. Conformational distortion explains its rapid photochemistry and reveals the trajectory of the 11-cis to trans isomerization. In addition, (2)H NMR has been applied to study the retinylidene dynamics in the dark and light-activated states. Upon isomerization there are drastic changes in the mobility of all three methyl groups. The relaxation data support an activation mechanism whereby the beta-ionone ring of retinal stays in nearly the same environment, without a large displacement of the ligand. Interactions of the beta-ionone ring and the retinylidene Schiff base with the protein transmit the force of the retinal isomerization. Solid-state (2)H NMR thus provides information about the flow of energy that triggers changes in hydrogen-bonding networks and helix movements in the activation mechanism of the photoreceptor.

Solid-state 2H NMR spectroscopy of retinal proteins in aligned membranes

Solid-state 2H NMR spectroscopy gives a powerful avenue to investigating the structures of ligands and cofactors bound to integral membrane proteins. For bacteriorhodopsin (bR) and rhodopsin, retinal was site-specifically labeled by deuteration of the methyl groups followed by regeneration of the apoprotein. 2H NMR studies of aligned membrane samples were conducted under conditions where rotational and translational diffusion of the protein were absent on the NMR time scale. The theoretical lineshape treatment involved a static axial distribution of rotating C-C2H3 groups about the local membrane frame, together with the static axial distribution of the local normal relative to the average normal. Simulation of solid-state 2H NMR lineshapes gave both the methyl group orientations and the alignment disorder (mosaic spread) of the membrane stack. The methyl bond orientations provided the angular restraints for structural analysis. In the case of bR the retinal chromophore is nearly planar in the dark- and all-trans light-adapted states, as well upon isomerization to 13-cis in the M state. The C13-methyl group at the "business end" of the chromophore changes its orientation to the membrane upon photon absorption, moving towards W182 and thus driving the proton pump in energy conservation. Moreover, rhodopsin was studied as a prototype for G protein-coupled receptors (GPCRs) implicated in many biological responses in humans. In contrast to bR, the retinal chromophore of rhodopsin has an 11-cis conformation and is highly twisted in the dark state. Three sites of interaction affect the torsional deformation of retinal, viz. the protonated Schiff base with its carboxylate counterion; the C9-methyl group of the polyene; and the beta-ionone ring within its hydrophobic pocket. For rhodopsin, the strain energy and dynamics of retinal as established by 2H NMR are implicated in substituent control of activation. Retinal is locked in a conformation that is twisted in the direction of the photoisomerization, which explains the dark stability of rhodopsin and allows for ultra-fast isomerization upon absorption of a photon. Torsional strain is relaxed in the meta I state that precedes subsequent receptor activation. Comparison of the two retinal proteins using solid-state 2H NMR is thus illuminating in terms of their different biological functions.

Methyl rotation barriers in proteins from 2H relaxation data. Implications for protein structure

Side-chain 2H and backbone 15N relaxation data have been collected at multiple temperatures in the samples of the SH3 domain from alpha-spectrin. Combined analyses of the data allowed for determination of the temperature-dependent correlation times tauf characterizing fast methyl motion. Molecular dynamics simulations confirmed that tauf are dominated by methyl rotation; the corresponding activation energies approximate methyl rotation barriers. For 33 methyl groups in the alpha-spectrin SH3 domain the average barrier height was thus determined to be 2.8 +/- 0.9 kcal/mol. This value is deemed representative of the "fluid" hydrophobic protein core where some barriers are increased and others are lowered because of the contacts with surrounding atoms, but there is no local order that could produce systematically higher (lower) barriers. For comparison, the MD simulation predicts the average barrier of 3.1 kcal/mol (calculated via the potential of mean force) or 3.4-3.5 kcal/mol (rigid barriers after appropriate averaging over multiple MD snapshots). The latter result prompted us to investigate rigid methyl rotation barriers in a series of NMR structures from the Protein Databank. In most cases the barriers proved to be higher than expected, 4-6 kcal/mol. To a certain degree, this is caused by tight packing of the side chains in the NMR structures and stems from the structure calculation procedure where the coordinates are first annealed toward the temperature of 0 K and then subjected to energy minimization. In several cases the barriers >10 kcal/mol are indicative of van der Waals violations. The notable exceptions are (i) the structures solved using the GROMOS force field where tight methyl packing is avoided (3.0-3.6 kcal/mol) and (ii) the structure solved by means of the dynamic ensemble refinement method (Lindorff-Larsen, K.; Best, R. B.; DePristo, M. A.; Dobson, C. M.; Vendruscolo, M. Nature 2005, 433, 128) (3.5 kcal/mol). These results demonstrate that methyl rotation barriers, derived from the experiments that are traditionally associated with studies of protein dynamics, can be also used in the context of structural work. This is particularly interesting in view of the recent efforts to incorporate dynamics data in the process of protein structure elucidation.

Two-dimensional solid-state NMR applied to a chimeric potassium channel

Solid-state NMR (ssNMR) represents a spectroscopic method to study membrane protein structure and dynamics in lipid bilayers. We present two-dimensional correlation experiments conducted on a fully [13C,15N] labeled version of a chimeric potassium (KcsA-Kv1.3) channel. Data obtained by using two different ion concentrations suggest a structural conservation of the selectivity filter region. SsNMR experiments conducted at two different temperatures point to differential molecular dynamics of the channel.

Conformational Flexibility of a Microcrystalline Globular Protein: Order Parameters by Solid-State NMR Spectroscopy

The majority of protein structures are determined in the crystalline state, yet few methods exist for the characterization of dynamics for crystalline biomolecules. Solid-state NMR can be used to probe detailed dynamic information in crystalline biomolecules. Recent advances in high-resolution solid-state NMR have enabled the site-specific assignment of (13)C and (15)N nuclei in proteins. With the use of multidimensional separated-local-field experiments, we report the backbone and side chain conformational dynamics of ubiquitin, a globular microcrystalline protein. The measurements of molecular conformational order parameters are based on heteronuclear dipolar couplings, and they are correlated to assigned chemical shifts, to obtain a global perspective on the sub-microsecond dynamics in microcrystalline ubiquitin. A total of 38 Calpha, 35 Cbeta and multiple side chain unique order parameters are collected, and they reveal the high mobility of ubiquitin in the microcrystalline state. In general the side chains show elevated motion in comparison with the backbone sites. The data are compared to solution NMR order parameter measurements on ubiquitin. The SSNMR measurements are sensitive to motions on a broader time scale (low microsecond and faster) than solution NMR measurements (low nanosecond and faster), and the SSNMR order parameters are generally lower than the corresponding solution values. Unlike solution NMR relaxation-based order parameters, order parameters for (13)C(1)H(2) spin systems are readily measured from the powder line shape data. These results illustrate the potential for detailed, extensive, and site-specific dynamic studies of biopolymers by solid-state NMR.

Order parameters based on (13)C(1)H, (13)C(1)H(2) and (13)C(1)H(3) heteronuclear dipolar powder patterns: a comparison of MAS-based solid-state NMR sequences

Order parameters describing conformational exchange processes on the nanosecond to microsecond timescale can be obtained from powder patterns in solid-state NMR (SSNMR) experiments. Extensions of these experiments to magic-angle spinning (MAS) based high-resolution experiments have been demonstrated, which show a great promise for site-specific probes of biopolymers. In this study, we present a detailed comparison of two pulse sequences, transverse Manfield-Rhim-Elleman-Vaughn (T-MREV) and Lee-Goldburg cross-polarization (LGCP), using experimental and simulation tools to explore their utility in the study of order parameters. We discuss systematic errors due to passively coupled (13)C or (1)H nuclei, as well as due to B(1) inhomogeneity. Both pulse sequences can provide quantitative measurements of the order parameter, but the LGCP experiment is capable of greater accuracy provided that the B(1) field is highly homogeneous. The T-MREV experiment is far better compensated for B(1) inhomogeneity, and it also performs better in situations with limited signal.

Characterization of Dynamics of Perdeuterated Proteins by MAS Solid-State NMR

We show in this communication that dynamic information for uniformly 2H,13C,15N isotopically enriched, crystalline proteins can be obtained by MAS solid-state NMR spectroscopy. The experiments make use of the deuterium quadrupolar tensor, which is the dominant interaction mechanism. Dynamic properties are accessed by measurement of the size of the quadrupolar coupling constant, Cq, and the value of the asymmetry parameter, eta, via evolution of the deuterium chemical shift, as well as by measurement of deuterium T1 relaxation times. Three-dimensional experiments are performed in order to obtain site-specific resolution. Due to proton dilution, no proton decoupling is required in the carbon evolution periods at MAS rotation frequencies of 10 kHz.

2H NMR and X-ray diffraction studies of methyl rotation in crystals of ortho-methyldibenzocycloalkanones

We have used (2)H NMR lineshape analyses and single crystal X-ray diffraction (XRD) to investigate the effects of molecular structure and crystalline environment on the rotational dynamics of methyl groups in four aromatic cycloalkanones. These include two methyl-substituted anthrones, one anthraquinone and one dibenzosuberone, which are known to undergo excited state H-atom tunneling from the ortho-methyl group to the carbonyl oxygen. With experiments conducted between 100 and 300K, samples 1,4-dimethylanthrone (DMAT) and 1,4-dimethylanthraquinone (DMAQ) were shown to enter the intermediate exchange regime (k(rot) approximately <10(7)s(-1)) at ca. 120K while samples of 1,4,10,10-tetramethylanthrone (TMAT) and 1,4-dimethyldibenzosuberone (DMDBS) remained in the fast exchange regime even at ca. 100K. Single crystal XRD analyses suggest that high intramolecular hindrance is avoided by molecular distortions, and that intermolecular contacts play an important role.

Magnetic resonance studies of the bacteriorhodopsin pump cycle

Active transport requires the alternation of substrate uptake and release with a switch in the access of the substrate binding site to the two sides of the membrane. Both the transfer and switch aspects of the photocycle have been subjects of magnetic resonance studies in bacteriorhodopsin. The results for ion transfer indicate that the Schiff base of the chromophore is hydrogen bonded before, during, and after its deprotonation. This suggests that the initial complex counterion of the Schiff base decomposes in such a way that the Schiff base carries its immediate hydrogen-bonding partner with it as it rotates during the first half of the photocycle. If so, bacteriorhodopsin acts as an inward-directed hydroxide pump rather than as an outward-directed proton pump. The studies of the access switch explore both protein-based and chromophore-based mechanisms. Combined with evidence from functional studies of mutants and other forms of spectroscopy, the results suggest that maintaining access to the extracellular side of the protein after photoisomerization involves twisting of the chromophore and that the decisive switch in access to the cytoplasmic side results from relaxation of the chromophore when the constraints on the Schiff base are released by decomposition of the complex counterion.

Dynamics and orientation of N+(CD3)3-bromoacetylcholine bound to its binding site on the nicotinic acetylcholine receptor

Dynamic and structural information has been obtained for an analogue of acetylcholine while bound to the agonist binding site on the nicotinic acetylcholine receptor (nAcChoR), using wide-line deuterium solid-state NMR. Analysis of the deuterium lineshape obtained at various temperatures from unoriented nAcChoR membranes labeled with deuterated bromoacetylcholine (BAC) showed that the quaternary ammonium group of the ligand is well constrained within the agonist binding site when compared with the dynamics observed in the crystalline solids. This motional restriction would suggest that a high degree of complementarity exists between the quaternary ammonium group of the ligand and the protein within the agonist binding site. nAcChoR membranes were uniaxially oriented by isopotential centrifugation as determined by phosphorous NMR of the membrane phospholipids. Analysis of the deuterium NMR lineshape of these oriented membranes enriched with the nAcChoR labeled with N(+)(CD(3))(3)-BAC has enabled us to determine that the angle formed between the quaternary ammonium group of the BAC and the membrane normal is 42 degrees in the desensitized form of the receptor. This measurement allows us to orient in part the bound ligand within the proposed receptor binding site.

Chromophore reorientation during the photocycle of bacteriorhodopsin: experimental methods and functional significance

Light-induced isomerization leads to orientational changes of the retinylidene chromophore of bacteriorhodopsin in its binding pocket. The chromophore reorientation has been characterized by the following methods: polarized absorption spectroscopy in the visible, UV and IR; polarized resonance Raman scattering; solid-state deuterium nuclear magnetic resonance; neutron and X-ray diffraction. Most of these experiments were performed at low temperatures with bacteriorhodopsin trapped in one or a mixture of intermediates. Time-resolved measurements at room temperature with bacteriorhodopsin in aqueous suspension can currently only be carried out with transient polarized absorption spectroscopy in the visible. The results obtained to date for the initial state and the K, L and M intermediates are presented and discussed. The most extensive data are available for the M intermediate, which plays an essential role in the function of bacteriorhodopsin. For this intermediate the various methods lead to a consistent picture: the curved all-trans polyene chain in the initial state straightens out in the M intermediate (13-cis) and the chain segment between C(5) and C(13) tilts upwards in the direction of the cytoplasmic surface. The kink at C(13) allows the positions of beta-ionone ring and Schiff base nitrogen to remain approximately fixed.

Bioenergetics of the Archaea

In the late 1970s, on the basis of rRNA phylogeny, Archaea (archaebacteria) was identified as a distinct domain of life besides Bacteria (eubacteria) and Eucarya. Though forming a separate domain, Archaea display an enormous diversity of lifestyles and metabolic capabilities. Many archaeal species are adapted to extreme environments with respect to salinity, temperatures around the boiling point of water, and/or extremely alkaline or acidic pH. This has posed the challenge of studying the molecular and mechanistic bases on which these organisms can cope with such adverse conditions. This review considers our cumulative knowledge on archaeal mechanisms of primary energy conservation, in relationship to those of bacteria and eucarya. Although the universal principle of chemiosmotic energy conservation also holds for Archaea, distinct features have been discovered with respect to novel ion-transducing, membrane-residing protein complexes and the use of novel cofactors in bioenergetics of methanogenesis. From aerobically respiring Archaea, unusual electron-transporting supercomplexes could be isolated and functionally resolved, and a proposal on the organization of archaeal electron transport chains has been presented. The unique functions of archaeal rhodopsins as sensory systems and as proton or chloride pumps have been elucidated on the basis of recent structural information on the atomic scale. Whereas components of methanogenesis and of phototrophic energy transduction in halobacteria appear to be unique to Archaea, respiratory complexes and the ATP synthase exhibit some chimeric features with respect to their evolutionary origin. Nevertheless, archaeal ATP synthases are to be considered distinct members of this family of secondary energy transducers. A major challenge to future investigations is the development of archaeal genetic transformation systems, in order to gain access to the regulation of bioenergetic systems and to overproducers of archaeal membrane proteins as a prerequisite for their crystallization.

Hydrogen bonding effects on amine rotation rates in crystalline amino acids

Rates of rotation for amines in a variety of crystalline environments are reported, and the trends are explained in terms of the strengths of local hydrogen bonding interactions. Proton spin-lattice relaxation times (T1) and deuterium broad-line NMR spectra have been measured for D-, DL- and L- aspartic acid, two polymorphs of glycine, alanine, and leucine in the temperature range from -40 to 110 degrees C. The energy barriers for amine rotation are 27 +/- 2 kJ mol-1 for D- or L-aspartic acid and 22 +/- 2 kJ mol-1 for DL-aspartic acid; these energies are slightly lower than the previously reported value for the L from based on direct proton T1 measurements at 60 MHz. The values for the alpha and gamma forms of glycine were 24 +/- 2 and 30 +/- 2 kJ mol-1 respectively, that for L-alanine was 40 +/- 2 and that for L-leucine was 49 +/- 3 kJ mol-1. These are all in rough agreement with previously reported values (although the differences for the polymorphs of glycine and for L- vs. DL-aspartic acid were not reported). Crystal structures of these amino acids indicate differences in hydrogen bonding environments around the R-NH3+ groups that are probably responsible for the different activation barriers. A molecular mechanics calculation of the rotation energy barriers for L- and DL-aspartic acid based on the crystal structures gave satisfactory agreement with experimental results if a uniform (and arbitrarily chosen) dielectric constant of 2.5 was assumed. Differences between L- and DL-aspartic acids and between two polymorphs of glycine were well represented qualitatively. Including additional neighboring molecules not involved in the hydrogen bonding or including periodic boundary conditions to describe the crystal packing did not significantly affect these results. If vacuum dielectric constants are used, the barriers are uniformly overestimated, and if the experimental macroscopic dielectric constant values are used, the barriers are generally underestimated. Dielectric constants differ substantially from one amino acid to another and significantly affect the estimated barriers; in fact, the bulk dielectric constants appear to be the major difference between the highest and the lowest values. The difficulty of accurately including dielectric relaxation into molecular mechanics calculations resulted in the disagreement between experimental measurements and theoretical calculations.

Membrane protein structure: the contribution and potential of novel solid state NMR approaches

Alternative methods for describing molecular detail for large integral membrane proteins are required in the absence of routine crystallographic approaches. Novel solid state NMR methods, devised for the study of large molecular assemblies, are now finding applications in biological systems, including integral membrane proteins. Wild-type and genetically engineered proteins can be investigated and detailed information about side chains, prosthetic groups, ligands (e.g. drugs) and binding sites can be deduced. The molecular structure and dynamics of selected parts of the proteins are accessible by a range of different solid state NMR approaches. Inter- and intra-atomic distances can be determined rather accurately (within ångströms) and the orientation of molecular bonds (within 2 degrees) can be measured in ideal cases. Here, a brief description of the methods is given and then some specific examples described with an indication of the future potential for the approaches in studying membrane proteins. It is anticipated that this emerging NMR methodology will be more widely used in the future, not only for resolving local structure, but also for more expansive descriptions of membrane protein structure at atomic resolution.