1H-15N-NMR studies of bacteriorhodopsin Halobacterium halobium. Conformational dynamics of the four-helical bundle

Direct observation of coherent oscillations in solution due to microheterogeneous environment

We report, for the first time, direct observation of coherent oscillations in the ground-state of IR775 dye due to microheterogeneous environment. Using ultrafast near-infrared degenerate pump-probe technique centered at 800 nm, we present the dynamics of IR775 in a binary mixture of methanol and chloroform at ultra-short time resolution of 30 fs. The dynamics of the dye in binary mixtures, in a time-scale of a few fs to ~740 ps, strongly varies as a function of solvent composition (volume fraction). Multi-oscillation behavior of the coherent vibration was observed, which increased with decreasing percentage of methanol in the dye mixture. Maximum number of damped oscillations were observed in 20% methanol. The observed vibrational wavepacket motion in the ground-state is periodic in nature. We needed two cosine functions to fit the coherent oscillation data as two different solvents were used. Dynamics of the dye molecule in binary mixtures can be explained by wavepacket motion in the ground potential energy surface. More is the confinement of the dye molecule in binary mixtures, more is the number of damped oscillations. The vibrational cooling time, τ₂, increases with increase in the confinement of the system. The observed wavepacket oscillations in ground-state dynamics continued until 1.6 ps.

Specificity of helix packing in transmembrane dimer of the cell death factor BNIP3: a molecular modeling study

BNIP3 is a mitochondrial 19-kDa proapoptotic protein, a member of the Bcl-2 family. It has a single COOH-terminal transmembrane (TM) alpha-helical domain, which is required for membrane targeting, proapoptotic activity, hetero- and homo-dimerization in membrane. The role and the molecular details of association of TM helices of BNIP3 are yet to be established. Here, we present a molecular modeling study of helix interactions in its membrane domain. The approach combines Monte Carlo conformational search in an implicit hydrophobic slab followed by molecular dynamics simulations in a hydrated full-atom lipid bilayer. The former technique was used for exhaustive sampling of the peptides' conformational space and for generation of putative "native-like" structures of the dimer. The latter ones were taken as realistic starting points to assess stability and dynamic behavior of the complex in explicit lipid-water surrounding. As a result, several groups of tightly packed right-handed structures of the dimer were proposed. They have almost similar helix-helix interface, which includes the motif A(176)xxxG(180)xxxG(184) and agrees well with previous mutagenesis data and preliminary NMR analysis. Molecular dynamics simulations of these structures reveal perfect adaptation of most of them to heterogeneous membrane environment. A remarkable feature of the predicted dimeric structures is the occurrence of a cluster of H-bonded histidine 173 and serines 168 and 172 on the helix interface, near the N-terminus. Because of specific polar interactions between the monomers, this part of the dimer has no such dense packing as the C-terminal one, thus allowing penetration of water from the extramembrane side into the membrane interior. We propose that the ionization state of His(173) can mediate structural and dynamic properties of the dimer. This, in turn, may be related to pH-dependent proapoptotic activity of BNIP3, which is triggering on by acidosis appearing under hypoxic conditions.

Differential dynamics in the G protein-coupled receptor rhodopsin revealed by solution NMR

G protein-coupled receptors are cell-surface seven-helical membrane proteins that undergo conformational changes on activation. The mammalian photoreceptor, rhodopsin, is the best-studied member of this superfamily. Here, we provide the first evidence that activation in rhodopsin may involve differential dynamic properties of side-chain versus backbone atoms. High-resolution NMR studies of alpha-(15)N-labeled receptor revealed large backbone motions in the inactive dark state. In contrast, indole side-chain (15)N groups of tryptophans showed well resolved, equally intense NMR signals, suggesting restriction to a single specific conformation.

Solution NMR spectroscopy of [alpha -15N]lysine-labeled rhodopsin: The single peak observed in both conventional and TROSY-type HSQC spectra is ascribed to Lys-339 in the carboxyl-terminal peptide sequence

[alpha-(15)N]Lysine-labeled rhodopsin, prepared by expression of a synthetic gene in HEK293 cells, was investigated both by conventional and transverse relaxation optimized spectroscopy-type heteronuclear single quantum correlation spectroscopy. Whereas rhodopsin contains 11 lysines, 8 in cytoplasmic loops and 1 each in the C-terminal peptide sequence and the intradiscal and transmembrane domains, only a single sharp peak was observed in dodecyl maltoside micelles. This result did not change when dodecyl maltoside was replaced by octyl glucoside or octyl glucoside-phospholipid-mixed micelles. Additional signals of much lower and variable intensity appeared at temperatures above 20 degrees C and under denaturing conditions. Application of the transverse relaxation optimized spectroscopy sequence resulted in sharpening of resonances but also losses of signal intensity. The single peak observed has been assigned to the C-terminal Lys-339 from the following lines of evidence. First, the signal is observed in HNCO spectra of rhodopsin, containing the labeled [(13)C]Ser-338/[(15)N]Lys-339 dipeptide. Second, addition of a monoclonal anti-rhodopsin antibody that binds to the C-terminal 8 aa of rhodopsin caused disappearance of the peak. Third, truncated rhodopsin lacking the C-terminal sequence Asp-330-Ala-348 showed no signal, whereas the enzymatically produced peptide fragment containing the above sequence showed the single peak. The results indicate motion in the backbone amide groups of rhodopsin at time scales depending on their location in the sequence. At the C terminus, conformational averaging occurs at the nanosecond time scale but varies from microsecond to millisecond in other parts of the primary sequence. The motions reflecting conformational exchange may be general for membrane proteins containing transmembrane helical bundles.

A (13)C NMR study on [3-(13)C]-, [1-(13)C]Ala-, or [1-(13)C]Val-labeled transmembrane peptides of bacteriorhodopsin in lipid bilayers: insertion, rigid-body motions, and local conformational fluctuations at ambient temperature

We have recorded (13)C NMR spectra of selectively [3-(13)C]Ala-, [1-(13)C]Ala-, or [1-(13)C]Val-labeled synthetic transmembrane peptides of bacteriorhodopsin (bR) and enzymatically cleaved C-2 fragment in the solid and dimyristoylphosphatidylcholine bilayer. It turned out that these transmembrane peptides either in hexafluoroisopropanol or cast from it take an ordinary alpha-helix (alpha(I)-helix) irrespective of their amino acid sequences with reference to the conformation-dependent (13)C chemical shifts of (Ala)(n) taking the alpha-helix form. These transmembrane peptides are not always static in the lipid bilayer as in the solid state but undergo rigid-body motions with various frequencies as estimated from suppressed peaks either by fast isotropic or large-amplitude motions (>10(8) Hz) or intermediate frequencies (10(5) or 10(3) Hz). Further, (13)C chemical shifts of the [3-(13)C]Ala-labeled peptides in the bilayer were displaced downfield by 0.3-1.1 ppm depending upon amino acid sequence with respect to those in the solid state, which were explained in terms of local conformational fluctuation (10(2) Hz) deviated from the torsion angles (alpha(II)-helix) from those of standard alpha-helix, under anisotropic environment in lipid bilayer, in addition to the above-mentioned rigid-body motions. The carbonyl (13)C peaks, on the other hand, are not sensitively displaced by such local anisotropic fluctuations, because they are more sensitive to the manner of hydrogen-bond interactions. The amino acid sequences of these peptides inserted within the bilayer were not always the same as those of intact bR, causing disposition of the transmembrane alpha-helical segment from that of intact bR. Finally, we confirmed that the (13)C NMR peak positions of the random coil form are located at the boundary between the alpha-helix and a turned structure in loop regions.

Customizing model membranes and samples for NMR spectroscopic studies of complex membrane proteins

Both solution and solid state nuclear magnetic resonance (NMR) techniques for structural determination are advancing rapidly such that it is possible to contemplate bringing these techniques to bear upon integral membrane proteins having multiple transmembrane segments. This review outlines existing and emerging options for model membrane media for use in such studies and surveys the special considerations which must be taken into account when preparing larger membrane proteins for NMR spectroscopic studies.

13C-13C rotational resonance in a transmembrane peptide: a comparison of the fluid and gel phases

A comparative study of two doubly 13C labeled amphiphilic transmembrane peptides was undertaken to determine the potential of rotational resonance for measuring internuclear distances through the direct dipolar coupling in the presence of motion. The two peptides, having the sequence acetyl-K2-G-L16-K2-A-amide, differed only in the position of 13C labels. The first peptide, [1-13C]leu(11):[alpha-13C]leu(12), had labels on adjacent residues, at the carbonyl of leu(11) and the alpha carbon of leu(12). The second, [1-13C]leu(8):[alpha-(13)/C]leu(11), was labeled on consecutive turns of the alpha-helical peptide. The internuclear distance between labeled positions of the first peptide, which for an ideal alpha helix has a value of 2.48 A, is relatively independent of internal flexibility or peptide conformational change. The dipolar coupling between these two nuclei is sensitive to motional averaging by molecular reorientation, however, making this peptide ideal for investigating these motions. The internuclear distance between labels on the second peptide has an expected static ideal alpha-helix value of 4.6 A, but this is sensitive to internal flexibility. In addition, the dipolar coupling between these two nuclei is much weaker because of their larger separation, making this peptide a much more difficult test of the rotational resonance technique. The dipolar couplings between the labeled nuclei of these two peptides were measured by rotational resonance in the dry peptide powders and in multilamellar dispersions with dimyristoylphosphatidylcholine in the gel phase, at -10 degrees C, and in the fluid phase, at 40 degrees C. The results for the peptide having adjacent labels can be readily interpreted in terms of a simple model for the peptide motion. The results for the second peptide show that, in the fluid phase, the motionally averaged dipolar coupling is too small to be measured by rotational resonance. Rotational resonance, rotational echo double resonance, and related techniques can be used to obtain reliable and valuable dipolar couplings in static solid and membrane systems. The interpretation of these couplings in terms of internuclear distances is straightforward in the absence of molecular motion. These techniques hold considerable promise for membrane protein structural studies under conditions, such as at low temperatures, where molecular motion does not modulate the dipolar couplings. However, a typical membrane at physiological temperatures exhibits complex molecular motions. In the absence of an accurate and detailed description of both internal and whole body molecular motions, it is unlikely that techniques of this type, which are based on extracting distances from direct internuclear dipolar couplings, can be used to study molecular structure under these conditions. Furthermore, the reduction in the strengths of the dipolar couplings by these motions dramatically reduces the useful range of distances which can be measured.

Escherichia coli diacylglycerol kinase: a case study in the application of solution NMR methods to an integral membrane protein

Diacylglycerol kinase (DAGK) is a 13-kDa integral membrane protein that spans the lipid bilayer three times and which is active in some micellar systems. In this work DAGK was purified using metal ion chelate chromatography, and its structural properties in micelles and organic solvent mixtures studies were examined, primarily to address the question of whether the structure of DAGK can be determined using solution NMR methods. Cross-linking studies established that DAGK is homotrimeric in decyl maltoside (DM) micelles and mixed micelles. The aggregate detergent-protein molecular mass of DAGK in both octyl glucoside and DM micelles was determined to be in the range of 100-110 kDa-much larger than the sum of the molecular weights of the DAGK trimers and the protein-free micelles. In acidic organic solvent mixtures, DAGK-DM complexes were highly soluble and yielded relatively well-resolved NMR spectra. NMR and circular dichroism studies indicated that in these mixtures the enzyme adopts a kinetically trapped monomeric structure in which it irreversibly binds several detergent molecules and is primarily alpha-helical, but in which its tertiary structure is largely disordered. Although these results provide new information regarding the native oligomeric state of DAGK and the structural properties of complex membrane proteins in micelles and organic solvent mixtures, the results discourage the notion that the structure of DAGK can be readily determined at high resolution with solution NMR methods.

1H-15N backbone resonance assignments of bacteriorhodopsin

Des-(232-248)-bacteriorhodopsin was solubilized in a membrane mimicking environment of methanol-chloroform (1:1) containing 0.1 M 2HCO2N2H4 and 1H-15N backbone resonance assignment was obtained using 2D HMQC, 3D NOESY-HMQC, 3D TOCSY-HMQC and 3D HMQC-NOESY-HMQC NMR experiments. 87 cross-peaks out of 117 present in the HMQC spectrum were assigned to particular residues in 1-73 and 195-231 parts of the protein. For these residues also signals of C alpha H and C beta H protons were assigned.

Structure and dynamics of bacteriophage IKe major coat protein in MPG micelles by solution NMR

The structure and dynamics of the 53-residue filamentous bacteriophage IKe major coat protein in fully protonated myristoyllysophosphatidylglycerol (MPG) micelles were characterized using multinuclear solution NMR spectroscopy. Detergent-solubilized coat protein [sequence: see text] mimics the membrane-bound "assembly intermediate" form of the coat protein which occurs during part of the phage life cycle. NMR studies of the IKe coat protein show that the coat protein is largely alpha-helical, exhibiting a long amphipathic surface. helix (Asn 4 to Ser 26) and a shorter "micelle-spanning" C-terminal helix which begins at TRP 29 and continues at least to Phe 48. Pro 30 likely occurs in the first turn of the C-terminal helix, where it is ideally situated given the hydrogen bonding and steric restrictions imposed by this residue. The similarity of 15N relaxation values (T1, T2, and NOE and 500 MHz and T2 at 600 MHz) among much of the N-terminal helix and all of the TM helix indicates that the N-terminal helix is as closely associated with the micelle as the TM helix. The description of the protein in the micelle is supported by the observation of NOEs between lysolipid protons and protein amide protons between asn 8 and Ser 50. The N-terminal and TM helices exhibit substantial mobility on the microsecond to second time scale, which likely reflects changes in the orientation between the two helices. The overall findings serve to clarify the role of individual residues in the context of a TM alpha-helix and provide an understanding of the secondary structure, dynamics, and aqueous and micellar environments of the coat protein.

High resolution 1H nuclear magnetic resonance of a transmembrane peptide

Although the strong 1H-1H dipolar interaction is known to result in severe homogeneous broadening of the 1H nuclear magnetic resonance (NMR) spectra of ordered systems, in the fluid phase of biological and model membranes the rapid, axially symmetric reorientation of the molecules about the local bilayer normal projects the dipolar interaction onto the motional symmetry axis. Because the linewidth then scales as (3 cos2 theta-1)/2, where theta is the angle between the local bilayer normal and the magnetic field, the dipolar broadening has been reduced to an "inhomogeneous" broadening by the rapid axial reorientation. It is then possible to obtain high resolution 1H-NMR spectra of membrane components by using magic angle spinning (MAS). Although the rapid axial reorientation effectively eliminates the homogeneous dipolar broadening, including that due to n = 0 rotational resonances, the linewidths observed in both lipids and peptides are dominated by low frequency motions. For small peptides the most likely slow motions are either a "wobble" or reorientation of the molecular diffusion axis relative to the local bilayer normal, or the reorientation of the local bilayer normal itself through surface undulations or lateral diffusion over the curved surface. These motions render the peptide 1H-NMR lines too broad to be observed at low spinning speeds. However, the linewidths due to these slow motions are very sensitive to spinning rate, so that at higher speeds the lines become readily visible. The synthetic amphiphilic peptide K2GL20K2A-amide (peptide-20) has been incorporated into bilayers of 1,2-di-d 27-myristoyl-sn-glycero-3-phosphocholine (DMPC-d54) and studied by high speed 1H-MAS-NMR. The linewidths observed for this transbilayer peptide, although too broad to be observable at spinning rates below -5 kHz, are reduced to 68 Hz at a spinning speed of 14 kHz (at 500C). Further improvements in spinning speed and modifications in sample composition designed to reduce the effectiveness of the slow motions responsible for the linewidth should result in significant further reduction in peptide linewidths. With this technique, there is now the potential for the use of 1H-MAS-NMR for the study of conformation, folding, and dynamics of small membrane peptides and protein fragments.

Backbone dynamics of (1-71)- and (1-36)bacterioopsin studied by two-dimensional (1)H- (15)N NMR spectroscopy

The backbone dynamics of uniformly (15)N-labelled fragments (residues 1-71 and 1-36) of bacterioopsin, solubilized in two media (methanol-chloroform (1:1), 0.1 M (2)HCO(2)NH(4), or SDS micelles) have been investigated using 2D proton-detected heteronuclear (1)H-(15)N NMR spectroscopy at two spectrometer frequencies, 600 and 400 MHz. Contributions of the conformational exchange to the transverse relaxation rates of individual nitrogens were elucidated using a set of different rates of the CPMG spin-lock pulse train and were essentially suppressed by the high-frequency CPMG spin-lock. We found that most of the backbone amide groups of (1-71)bacterioopsin in SDS micelles are involved in the conformational exchange process over a rate range of 10(3) to 10(4) s(-1). This conformational exchange is supposed to be due to an interaction between two α-helixes of (1-71)bacterioopsin, since the hydrolysis of the peptide bond in the loop region results in the disappearance of exchange line broadening. (15)N relaxation rates and (1)H-(15)N NOE values were interpreted using the model-free approach of Lipari and Szabo [Lipari, G. and Szabo, A. (1982) J. Am. Chem. Soc., 104, 4546-4559]. In addition to overall rotation of the molecule, the backbone N-H vectors of the peptides are involved in two types of internal motions: fast, on a time scale <20 ps, and intermediate, on a time scale close to 1 ns. The intermediate dynamics in the α-helical stretches was mostly attributed to bending motions. A decrease in the order parameter of intermediate motions was also observed for residues next to Pro(50), indicating an anisotropy of the overall rotational diffusion of the molecule. Distinctly mobile regions are identified by a large decrease in the order parameter of intermediate motions and correspond to the N- and C-termini, and to a loop connecting the α-helixes of (1-71)bacterioopsin. The internal dynamics of the α-helixes on the millisecond and nanosecond time scales should be taken into account in the development of a model of the functioning bacteriorhodopsin.

Spectroscopic studies of bacteriorhodopsin fragments dissolved in organic solution

Fourier transform infrared and UV fourth-derivative spectroscopies were used to study the secondary structure of bacteriorhodopsin and its chymotryptic and one of the sodium borohydride fragments dissolved in chloroform-methanol (1:1, v/v), 0.1 M LiClO4. The C1 fragment (helices C, D, E, F, and G) showed an alpha-helical content of about 53%, whereas C2 (helices A and B) had about 60%, and B2 (helices F and G) about 65% alpha-helix. The infrared main band indicated differences in alpha-helical properties between these fragments. These techniques were also used to obtain information on the interactions among helices. According to the results obtained from the hydrogen/deuterium exchange kinetics, about 40% of the amide protons of C2 are particularly protected against exchange, whereas for the C1 fragment this process is unexpectedly fast. UV fourth-derivative spectra of these samples were used to obtain information about the environment of Trp side chains. The results showed that the Trp residues of C2 are more shielded from the solvent than those of C1 or B2. The results of this work indicate that the specific interactions existing between the transmembrane segments induce different types of helical conformations in native bacteriorhodopsin.

Backbone dynamics of (1-71)bacterioopsin studied by two-dimensional 1H-15N NMR spectroscopy

The backbone dynamics of a uniformly 15N-labelled proteolytic fragment (residues 1-71) of bacteriorhodopsin, solubilized in two media [methanol/chloroform (1:1), 0.1 M 2HCO2NH4 and SDS micelles] have been investigated using two-dimensional proton-detected heteronuclear 1H-15N NMR spectroscopy. A set of longitudinal and transverse relaxation rates of 15N nuclei and 1H-15N NOE were obtained for 61 backbone amide groups. The contribution of the conformational exchange to transverse relaxation rates of individual nitrogens was elucidated using a set of different rates of the Carr-Purcell-Meiboom-Gill (CPMG) spin-lock pulse train. We found that most of the backbone amide groups are involved in the co-operative exchange process over the rate range 10(3)-10(4) s-1, with the chemical-shift dispersion near 1 ppm. Contributions of conformational exchange to the measured transverse relaxation were essentially suppressed by the 3-kHz (spin-echo period tau = 0.083 ms) CPMG spin-lock. Under these conditions, the measured longitudinal, transverse relaxation rates and NOE values were interpreted using the model-free approach of Lipari and Szabo [Lipari, G. & Szabo, A. (1982) J. Am. Chem. Soc. 104, 4546-4559]. In both media used, the protein exhibits very similar dynamic properties, and has overall rotational correlation times of 7.0 ns and 6.6 ns in organic mixture and in SDS micelles, respectively. In addition to overall rotation of the molecule, the backbone N-H vectors are involved in two types of internal motions; fast, on a time scale of < 20 ps, and intermediate, close to 1 ns. Distinctly mobile regions are identified by a large decrease in the overall order parameter and correspond to N-terminal residues (residues 1-7 both for organic solvent and micelles), C-terminal residues (residues 65-71 and 69-71 for organic solvent and micelles, respectively) and residues connecting alpha helices (residues 33-41 and 33-38, for organic solvent and micelles, respectively). A decrease in the order parameter was also observed for residues next to Pro50, indicating a higher flexibility in this region. Thus, backbone dynamic parameters of (1-71)bacterioopsin are in good correspondence with its spatial structure [Pervushin, K. V., Orekhov, V. Yu., Popov, A., Musina, L. Yu., Arseniev, A. S., (1994) Eur. J. Biochem., in the press]. The observed conformational exchange behavior of alpha helices seems to be induced by the flickering helix-helix interaction and could be important for the functioning of bacteriorhodopsin.

Three-dimensional structure of (1-71)bacterioopsin solubilized in methanol/chloroform and SDS micelles determined by 15N-1H heteronuclear NMR spectroscopy

Spatial structures of a chymotryptic fragment C2 (residues 1-71) of bacterioopsin from Halobacterium halobium, solubilized in a mixture of methanol/chloroform (1:1, by vol.) and 0.1 M 2HCO2NH4, or in perdeuterated sodium (2H)dodecyl sulfate (SDS) micelles in the presence of perdeuterated (2,2,2-2H)trifluoroethanol, were determined by two-dimensional and three-dimensional heteronuclear 15N-1H NMR techniques. The influence of (2,2,2-2H)trifluoroethanol on the conformational dynamics of C2 in micelles and the effect of the salt (organic mixture) were studied. Under the best conditions, 1H and 15N resonances of 15N-uniformly enriched protein were assigned in both milieus by homonuclear two-dimensional NOE (NOESY) and two-dimensional total-correlated (TOCSY) spectra and heteronuclear three-dimensional NOESY-multiple-quantum-correlation (HMQC) and TOCSY-HMQC spectra. 651 (organic mixture) and 520 (micelles) interproton-distance constraints, derived from volumes of cross-peaks in two-dimensional NOESY and three-dimensional NOESY-HMQC spectra, along with deuterium exchange rates of amide groups measured in both milieus and 51 HN-C alpha H coupling constants obtained in the case of the organic mixture, were used in the construction of C2 spatial structures. Obtained structures are similar in both milieus and have two right-handed alpha-helical regions stretching from Pro8 to Met32 and Phe42 to Tyr64 (organic mixture), and from Pro8 to Met32 and Ala39 to Leu62 (micelles). In micelles, the second alpha helix is terminated by C-cap Gly63, adopting a conformation characteristic of a left-handed helix. Residues Gly65 to Thr67 from the turn of a right-handed helix. In the isotropic medium of the organic mixture, the C-terminal region of residues 65-71 lacks an ordered structure. Torsion angles chi 1 were unequivocally determined for 18 alpha-helical residues in both milieus. In the isotropic organic mixture and anisotropic micellar system, C2 remains a compact structure with a characteristic size of 3.0-3.5 nm. C2 seems to be present in at least two conformational states, packed and unpacked. Using NMR data, along with the electron cryomicroscopy model of bacteriorhodopsin [Henderson, R., Baldwin, J. M., Ceska, T. A., Zemlin, F., Beckman, E. & Downing, K. H. (1990) J. Mol. Biol. 213, 899-929], we suggested a model for the conformation of C2 in this putative close-packed state. However, no NOE contact between alpha helices was found in either milieu.