Optimizing and characterizing alignment of oriented lipid bilayers containing gramicidin D

Novel chelate-induced magnetic alignment of biological membranes

A phospholipid chelate complexed with ytterbium (DMPE-DTPA:Yb3+) is shown to be readily incorporated into a model membrane system, which may then be aligned in a magnetic field such that the average bilayer normal lies along the field. This so-called positively ordered smectic phase, whose lipids consist of less than 1% DMPE-DTPA:Yb3+, is ideally suited to structural studies of membrane proteins by solid-state NMR, low-angle diffraction, and spectroscopic techniques that require oriented samples. The chelate, 1,2-dimyristoyl-sn-glycero-3-phosphoethanolamine diethylenetriaminepentaacetic acid, which strongly binds the lanthanide ions and serves to orient the membrane in a magnetic field, prevents direct lanthanide-protein interactions and significantly reduces paramagnetic shifts and line broadening. Similar low-spin lanthanide chelates may have applications in field-ordered solution NMR studies of water-soluble proteins and in the design of new magnetically aligned liquid crystalline phases.

13C solid-state NMR of gramicidin A in a lipid membrane

The natural-abundance 13C NMR spectrum of gramicidin A in a lipid membrane was acquired under magic-angle spinning conditions. With fast sample spinning (15 kHz) at approximately 65 degrees C the peaks from several of the aliphatic, beta-, alpha-, aromatic, and carbonyl carbons in the peptide could be resolved. The resolution in the 13C spectrum was superior that observed with 1H NMR under similar conditions. The 13C linewidths were in the range 30-100 Hz, except for the alpha- and beta-carbons, the widths of which were approximately 350 Hz. The beta-sheet-like local structure of gramicidin A was observed as an upfield shift of the gramicidin alpha and carbonyl resonances. Under slow sample spinning (500 Hz), the intensity of the spinning sidebands from 13C in the backbone carbonyls was used to determine the residual chemical shift tensor. As expected, the elements of the residual chemical shift tensor were consistent with the single-stranded, right-handed beta6.3 helix structure proposed for gramicidin A in lipid membranes.

Solid-state nuclear magnetic resonance characterization of gramicidin channel structure

The method of using orientational constraints derived from solid-state NMR for structural characterization of polypeptides in heterogeneous environments has now been demonstrated. A very high resolution structure has been achieved that has led to greater functional understanding of this channel. Much can be done to improve this structural technique to make it more efficient and more generally applicable. Others as well as ourselves are applying this approach to membrane proteins. Although solid-phase synthesis and specific site isotopic labeling has been essential for the development described here, one of the primary challenges is to be able to use amino acid-specific and uniform labeling of peptides and proteins by biosynthetic means for isotopic incorporation. This will allow for the study of many more proteins and significantly large proteins. Unlike solution NMR structural methods, there are no intrinsic molecular weight limitations. In fact, as the molecular weight increases the molecular motion will become less and the spectroscopic properties will improve. The major limitation will be sensitivity: as the molecular weight increases the number of moles will decrease in the samples, causing sensitivity to decrease. Advances in field strength and NMR technology help to address this problem. With larger molecules and more isotopically labeled sites resolution could also be a problem; however, the two- and three-dimensional methods demonstrated by Opella and co-workers clearly show the potential for enormous resolving power. In the 15N dimension alone it is shown that the resolution is greater than in solution NMR. Although challenges such as spectral assignments have yet to be completely solved, several approaches have been described, and the prospects are excellent for solving this and other problems facing the development of this novel approach for structural elucidation. Although there is an attempt to get away from solid-phase synthesis to solve larger molecular weight structures, peptide synthesis will continue to be important for generating single- and double-site labeled model compounds for characterizations of spin interaction tensors. Such characterizations will continue to be a very important aspect of this structural approach.

Transmembrane four-helix bundle of influenza A M2 protein channel: structural implications from helix tilt and orientation

The transmembrane portion of the M2 protein from the Influenza A virus has been studied in hydrated dimyristroylphosphotidylcholine lipid bilayers with solid-state NMR. Orientational constraints were obtained from isotopically labeled peptide samples mechanically aligned between thin glass plates. 15N chemical shifts from single site labeled samples constrain the molecular frame with respect to the magnetic field. When these constraints are applied to the peptide, modeled as a uniform alpha-helix, the tilt of the helix with respect to the bilayer normal was determined to be 33 degrees +/- 3 degrees. Furthermore, the orientation about the helix axis was also determined within an error of +/- 30 degrees. These results imply that the packing of this tetrameric protein is in a left-handed four-helix bundle. Only with such a large tilt angle are the hydrophilic residues aligned to the channel axis.

Structure, energetics, and dynamics of lipid-protein interactions: A molecular dynamics study of the gramicidin A channel in a DMPC bilayer

The microscopic details of lipid-protein interactions are examined using molecular dynamics simulations of the gramicidin A channel embedded in a fully hydrated dimyristoyl phosphatidylcholine (DMPC) bilayer. A novel construction protocol was used to assemble the initial configurations of the membrane protein complex for the simulations. Three hundred systems were constructed with different initial lipid placement and conformations. Seven systems were simulated with molecular dynamics. One system was simulated for a total of 600 psec, four were simulated for 300 psec, and two for 100 psec. Analysis of the resulting trajectories shows that the bulk solvent-membrane interface region is much broader than traditionally pictured in simplified continuum theories: its width is almost 15 angstroms. In addition, lipid-protein interactions are far more varied, both structurally and energetically, than is usually assumed: the total interaction energy between the gramicidin A and the individual lipids varies from 0 to -50 kcal/mol. The deuterium quadrupolar splittings of the lipid acyl chains calculated from the trajectories are in good agreement with experimental data. The lipid chains in direct contact with the GA are ordered but the effect is not uniform due to the irregular surface of the protein. Energy decompositions shows that the most energetically favorable interactions between lipid and protein involve nearly equal contributions from van der Waals and electrostatic interactions. The tryptophans, located near the bulk-membrane interface, appear to be particularly important in mediating both hydrogen bonding interactions with the lipid glycerol backbone and water and also in forming favorable van der Waals contacts with the hydrocarbon chains. In contrast, the interactions of the leucine residues with the lipids, also located near the interface, are dominated by van der Waals interactions with the hydrocarbon lipid chains.

Reconstitution of membrane proteins into lipid-rich bilayered mixed micelles for NMR studies

This paper describes a study undertaken to assess the possibility and practical consequences of reconstituting integral and peripheral membrane proteins into bilayered discoidal mixed micelles ("bicelles") composed of dimyristoylphosphatidylcholine and smaller amounts of either CHAPSO or short-chain phosphatidylcholine. The amphiphilic assemblies in these mixtures are uniquely suited for use in NMR structural studies because they can be magnetically oriented with experimentally-tunable system order. The first step of this study was to test about 15 membrane-associating polypeptides and proteins for their ability to interfere with magnetic orientation of the bicellar assemblies. A variety of results were obtained ranging from no perturbation to a complete disruption of orientation. Second, the suitability of bicelles as mimics of natural bilayers was tested by reconstituting diacylglycerol kinase, an integral membrane enzyme. The kinase was observed to be functional and completely stable for at least 24 h when incubated at 38 degrees C in bicelles. Third, the NMR spectra from a number of bicelle-reconstituted proteins were examined. In some cases, 13C NMR resonances from reconstituted proteins were extremely broad and asymmetric. In other cases, resonances from reconstituted proteins were moderately broad, but much less so than resonances from proteins reconstituted into multilayers oriented by mechanical methods. In the cases of two surface-associating proteins (cytochrome c and leucine enkephalin), oriented sample 13C NMR spectra of extremely high resolution were obtained. For these proteins it was also demonstrated that the experimentally variable order of the bicellar assemblies could be exploited to provide a means of screening for detergent-specific structural perturbations, for making spectral assignments, and for measuring chemical shift anisotropies and dipolar couplings. Taken as a whole, these results indicate that bicelles may be uniquely and effectively employed as model membranes to facilitate NMR structural studies of many, but not all, membrane proteins.

Molecular dynamics simulation of the gramicidin channel in a phospholipid bilayer

A molecular dynamics simulation of the gramicidin A channel in an explicit dimyristoyl phosphatidylcholine bilayer was generated to study the details of lipid-protein interactions at the microscopic level. Solid-state NMR properties of the channel averaged over the 500-psec trajectory are in excellent agreement with available experimental data. In contrast with the assumptions of macroscopic models, the membrane/solution interface region is found to be at least 12 A thick. The tryptophan side chains, located within the interface, are found to form hydrogen bonds with the ester carbonyl groups of the lipids and with water, suggesting their important contribution to the stability of membrane proteins. Individual lipid-protein interactions are seen to vary from near 0 to -50 kcal/mol. The most strongly interacting conformations are short-lived and have a nearly equal contribution from both van der Waals and electrostatic energies. This approach for performing molecular dynamics simulations of membrane proteins in explicit phospholipid bilayers should help in studying the structure, dynamics, and energetics of lipid-protein interactions.

Qualitative comparison of the bilayer-associated structures of diacylglycerol and a fluorinated analog based upon oriented sample NMR data

sn-1,2-Dimyristoylglycerol (DMDAG) and sn-1,2-dimyristoyl-3-fluoropropanediol (DMFPD) were synthesized in carbonyl 13C-labeled and acyl chain perdeuterated forms. These compounds were reconstituted at low levels into both randomly dispersed dimyristoylphosphatidylcholine (DMPC) bilayers and magnetically orientable DMPC media. Samples were subjected to NMR analysis, leading to a substantial body of 2H quadrupolar splitting, 13C-13C and 13C-19F dipolar coupling and 13C chemical shift anisotropy data for both compounds. A number of measurements were also made for DMPC. The data acquired from magnetically oriented samples were found to be undesirably affected by the presence of artifacts related to the experimental use of the oriented lipid media. However, it was possible to draw a number of qualitative conclusions from the data and to correct the data so that they may ultimately prove useful in quantitative structural analyses of bilayer-associated DMDAG and DMFPD. Comparison of the DMDAG data with corresponding measurements of DMPC suggests a high degree of similarity between the two compounds within bilayers composed primarily of L alpha phase phosphatidylcholine, consistent with previous work (S.O. Smith et al. (1992) Biochemistry 31, 11660-11664). Comparison of the data for DMDAG and DMFPD suggests that the hydroxyl proton of DMDAG is involved in hydrogen bonding interactions, which appear to be largely responsible for maintaining the orientational inequivalence of its two acyl chains. Bilayer-associated DMFPD also appears to exhibit a higher degree of whole-molecule disorder than DMDAG, again suggestive of an important structural role for the hydroxyl proton of DMDAG. Finally, comparison of 13C-NMR data for oriented DMPC in the presence and absence of low levels of DMDAG and DMFPD indicated that neither compound induces a significant change in the averaged conformational state of DMPC comprising the bilayer matrix of the oriented samples.

The structure of an integral membrane peptide: a deuterium NMR study of gramicidin

Solid state deuterium NMR was employed on oriented multilamellar dispersions consisting of 1,2-dilauryl-sn-glycero-3-phosphatidylcholine and deuterium (2H) exchange-labeled gramicidin D, at a lipid to protein molar ratio (L/P) of 15:1, in order to study the dynamic structure of the channel conformation of gramicidin in a liquid crystalline phase. The corresponding spectra were used to discriminate between several structural models for the channel structure of gramicidin (based on the left- and right-handed beta 6.3 LD helix) and other models based on a structure obtained from high resolution NMR. The oriented spectrum is complicated by the fact that many of the doublets, corresponding to the 20 exchangeable sites, partially overlap. Furthermore, the asymmetry parameter, eta, of the electric field gradient tensor of the amide deuterons is large (approximately 0.2) and many of the amide groups are involved in hydrogen bonding, which is known to affect the quadrupole coupling constant. In order to account for these complications in simulating the spectra in the fast motional regime, an ab initio program called Gaussian 90 was employed, which permitted us to calculate, by quantum mechanical means, the complete electric field gradient tensor for each residue in gramicidin (using two structural models). Our results indicated that the left-handed helical models were inconsistent with our observed spectra, whereas a model based on the high-resolution structure derived by Arseniev and coworkers, but relaxed by a simple energy minimization procedure, was consistent with our observed spectra. The molecular order parameter was then estimated from the motional narrowing assuming the relaxed (right-handed) Arseniev structure. Our resultant order parameter of SZZ = 0.91 translates into an rms angle of 14 degrees, formed by the helix axis and the local bilayer normal. The strong resemblance between our spectra (and also those reported for gramicidin in 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC) multilayers) and the spectra of the same peptide incorporated in a lyotropic nematic phase, suggests that the lyotropic nematic phase simulates the local environment of the lipid bilayer.

Simulation of NMR data from oriented membrane proteins: practical information for experimental design

Several hundred solid state NMR dipolar couplings and chemical shift anisotropies were simulated for the polytopic membrane protein, bacteriorhodopsin, and for an idealized transmembrane peptide conforming to several different secondary structures (alpha- and 3(10)-helices and parallel and antiparallel beta-sheets), each at several tilt angles with respect to the bilayer normal. The use of macroscopically oriented samples was assumed. The results of these simulations suggest: (i) Because of the r-3 dependence of dipolar coupling, it is likely to prove difficult to successfully execute uniform isotopic enrichment strategies to generate large numbers of quantitatively interpretable structural measurements in oriented sample NMR studies of membrane proteins. (ii) There are a number of readily implementable specific isotopic labeling schemes which can yield data patterns sufficient to identify local secondary structure for transmembrane segments of idealized proteins which are tilted by < 10 degrees with respect to the bilayer normal. (iii) The measurement of dipolar coupling constants between 13C-, 19F-, and/or 3H-labeled side chains of proximal residues may prove effective as routes to long range tertiary structural data constraints.

High-resolution structure and dynamic implications for a double-helical gramicidin A conformer

The high-resolution structure of a dimeric conformer of gramicidin A, a 15-residue polypeptide, has been determined in the mixed-solvent system of benzene and ethanol by 2D NMR techniques. NOEs, coupling constants and hydrogen-bond information were used to generate 744 experimental constraints for the dimer. Stereoassignment of most beta-methylene groups was achieved by analysis of 3J alpha beta, d alpha beta(i,i), dN beta(i,i) and dN beta(i + 1,i) distances, and consideration of the initial backbone structure determinations. Stereoassignment of several leucine methyl groups was accomplished via a distance geometry/simulated annealing routine, used for structure determination and refinement. The relatively static backbone structure was determined first and held rigid while side-chain conformations were calculated. This procedure is evaluated versus standard NMR structure determination protocols. The backbone is an antiparallel intertwined double helix, with 5.6-5.7 residues per turn, a total dimer length of 36-37 A, and a pore width of 2.5-3.0 A (van der Waals to van der Waals). The structure and dynamics of the side chains are discussed in depth, with careful attention for both the convergence of structures and the residual constraint violations per residue. Side-chain positions impart substantial amphipathic character to the helix, which could influence the conformational change that takes place upon membrane insertion of this channel-forming polypeptide.

2H NMR lineshapes of immobilized uniaxially oriented membrane proteins

As a method for the structure determination of integral membrane proteins or other large macromolecular complexes, a solid state 2H NMR approach is presented, capable of measuring the orientations of individual chemical bond vectors. In an immobilized uniaxially oriented sample, the bond angle of a deuterium-labelled methyl group relative to the axis of ordering can be calculated from the quadrupole splitting in the "zero-tilt" spectrum where the sample normal is aligned parallel to the spectrometer field direction. However, since positive and negative values of this splitting cannot be distinguished, there may appear to be two solutions, of which only one describes the correct molecular geometry. We show that it is possible to determine the bond angle uniquely between 0 degree and 90 degrees, by analysing the lineshapes of a tilt series of spectra acquired over different sample inclinations. The lineshape equation describing such oriented 2H NMR spectra will be derived (for asymmetry parameter eta = 0) and discussed, with an illustration of the various linebroadening effects from which the orientational distribution function in the macroscopically ordered system can be determined. This strategy is then applied to specifically deuterium-labelled retinal in dark-adapted bacteriorhodopsin, prepared in a uniaxially oriented sample from purple membrane fragments. From the quadrupole splitting in the zero-tilt spectrum and by lineshape simulations, the deuteromethyl group at C20 on retinal is found to make an angle of 32 degrees +/- 1 degree with the membrane normal, and the sample mosaic spread to be around +/- 8 degrees. The resulting orientation of retinal is in excellent agreement with its known structure in bacteriorhodopsin, and together with the results on other methyl groups it will be possible to construct a detailed picture of the chromophore in the protein binding pocket.

Orientational constraints as three-dimensional structural constraints from chemical shift anisotropy: the polypeptide backbone of gramicidin A in a lipid bilayer

Chemical shifts observed from samples that are uniformly aligned with respect to the magnetic field can be used as very high-resolution structural constraints. This constraint takes the form of an orientational constraint rather than the more familiar distance constraint. The accuracy of these constraints is dependent upon the quality of the tensor characterization. Both tensor element magnitudes and tensor orientations with respect to the molecular frame need to be considered. Here these constraints have been used to evaluate models for the channel conformation of gramicidin A. Of the three models used, the one experimentally derived model of gramicidin in sodium dodecyl sulfate micelles fits the data least well.

On the supramolecular organization of gramicidin channels. The elementary conducting unit is a dimer

The question, whether the conducting channels formed by the linear gramicidins are dimers (as is generally believed) or tetramers (as has been recently proposed [Stark G., M. Strässle, and Z. Takacz. 1986. J. Membr. Biol. 89:23-37; Strässle, M., G. Stark, M. Wilhelm, P. Daumas, F. Heitz, and R. Lazaro. 1989. Biochim. Biophys. Acta. 980:305-314]) has been addressed in single-channel experiments. The experimental approach was based on the ability of electrophysiological (single-channel) experiments to resolve the number of hybrid channel types that could form between gramicidin A or C and O-pyromellityl-gramicidin A or C (in which a pyromellitic acid residue has been esterified to the ethanolamine-OH group [Apell, H.-J., E. Bamberg, H. Alpes, and P. Läuger. 1977. J. Membr. Biol. 31:171-188]). The presence of the bulky, negatively charged pyromellityl group at the channel entrances endows the hybrid channels with characteristically different features and thus facilitates the resolution of the different hybrid channel types. Only two hybrid channel types were detected, indicating that the conducting channels are membrane-spanning dimers. There was likewise no evidence for lateral association between conducting channels and nonconducting monomers. These results can be reconciled with those of Stark et al. (op. cit.) if gramicidin channel formation involves a (slow) folding into beta 6.3-helical monomers followed by the dimerization step.

Experimental determination of torsion angles in the polypeptide backbone of the gramicidin A channel by solid state nuclear magnetic resonance

An analytical method for the determination of torsion angles from solid state 15N nuclear magnetic resonance (n.m.r.) spectroscopic data is demonstrated. Advantage is taken of the 15N-1H and 15N-13C dipolar interactions as well as the 15N chemical shift interaction in oriented samples. The membrane-bound channel conformation of gramicidin A has eluded an atomic resolution structure determination by more traditional approaches. Here, the torsion angles for the Ala3 site are determined by obtaining the n.m.r. data for both the Gly2-Ala3 and Ala3-Leu4 peptide linkages. Complete utilization of the orientational constraints derived from these orientation-dependent nuclear spin interactions in restricting the conformational space is most effectively achieved by utilizing spherical trigonometry. Two possible sets of torsion angles for the Ala3 site are obtained (phi, psi = -129 degrees, 153 degrees and -129 degrees, 122 degrees), both of which are consistent with a right-handed beta-helix. Other functional and computational evidence strongly supports the set for which the carbonyl oxygen atom of the Ala3-Leu4 linkage is rotated into the channel lumen.

Solid-state nuclear magnetic resonance derived model for dynamics in the polypeptide backbone of the gramicidin A channel

The dynamics of the backbone of the gramicidin A transmembrane cation channel in dimyristoylphosphatidylcholine bilayers have been investigated using solid state 15N nuclear magnetic resonance (n.m.r.) spectroscopy. With the temperature-dependent fluidity of the bilayer, the rates of motions in the helical gramicidin channel can be modulated. It is shown that in the gel phase, all substantial motions of the channel are slow on the timescale of the n.m.r. experiment (3.5 kHz). The use of oriented samples in which the axis of global channel rotation is aligned parallel to the magnetic field enables separation of global and local dynamics. Spectra obtained from oriented bilayer samples containing single-site 15N-labeled gramicidin at 8 degrees C are analyzed to yield a spatial model for local backbone motion. This model includes the axis of motion, the mean orientation, and the maximum amplitude of displacement for individual peptide planes. Specific sites in the first turn of the amino terminus were investigated, with emphasis on the Ala3 and Leu4 linkages, for which the orientation of the 15N chemical shift tensor with respect to the molecular frame has been determined. The effect of two well-characterized bilayer defect structures, parabolic focal conics and oily streaks, is included in the spectral simulations. It is found that only relatively small amplitude motions are possible at the two sites, with amplitudes of not more than +/- 8 degrees and +/- 15 degrees for the Ala3 and Leu4 sites, respectively. Detailed characterization of the bilayer surface geometry in the oriented samples is presently the major limiting factor in the use of this technique for probing the spatial extent of local motions in integral membrane proteins.

Trypsin inhibitor paradoxically stabilizes trypsin activity in sodium dodecyl sulfate, facilitating proteolytic fingerprinting

Normally trypsin has negligible activity after being dissolved in sodium dodecyl sulfate (SDS), and so it has had little utility for proteolytic fingerprinting during gel electrophoresis. Here it is demonstrated that trypsin retained activity in SDS if it was first complexed to either of two soybean-derived protease inhibitors: trypsin inhibitor (Kunitz) or trypsin-chymotrypsin inhibitor (Bowman-Birk). The inhibitors alone did not cause proteolysis. Heating or acidification in SDS inactivated the inhibitor-dependent tryptic activity, as did prior treatment with tosyl lysine chloromethyl ketone, a covalent affinity reagent for trypsin. Quenching of samples with acid at intervals prior to gel electrophoresis revealed that proteolysis did not occur in sample buffer (pH 6.8), but only at higher pH and during gel electrophoresis. Exposure of trypsin to SDS prior to addition of trypsin inhibitor resulted in an irreversible loss of activity with a half-life of about 10 s. It is proposed that the trypsin inhibitors stabilize trypsin by retarding its denaturation in SDS. The substrate for these experiments was the alpha subunit of the Na,K-ATPase. The same pattern of Na,K-ATPase fragments was obtained with bovine and porcine trypsin and with rat and porcine Na,K-ATPases. Different fragments resulted when chymotrypsin or elastase were substituted for trypsin; these proteases were active in the absence of an inhibitor, and were not markedly stabilized by interaction with soybean trypsin-chymotrypsin inhibitor (Bowman-Birk).

Magnetically orientable phospholipid bilayers containing small amounts of a bile salt analogue, CHAPSO

Buffered mixtures of the detergent 3-(cholamidopropyl)dimethylammonio-2-hydroxy-1-propanesulfonate (CHAPSO) and dimyristoylphosphatidylcholine (DMPC) orient in the presence of a strong magnetic field over a wide range of water contents (at least 65-85%) and CHAPSO:DMPC molar ratios (typically 1:10-1:3). 31P NMR studies show that the phospholipid in such mixtures is oriented with its director axis perpendicular to the magnetic field. 31P and 2H NMR results also suggest that the structure and dynamics of the DMPC molecules are similar to that of pure phospholipids existing in the liquid crystalline (L alpha) bilayer phase. The ability of 1:5 CHAPSO:DMPC samples to orient is highly tolerant of large changes in temperature, pH, and ionic strength, as well as to the addition of substantial amounts of charged amphiphiles or soluble protein. However, 2H NMR studies of deuterated beta-dodecyl melibiose (DD-MB) solubilized in the system indicate the head group conformation and/or dynamics of this glycolipid analogue is dependent upon the CHAPSO concentration. Despite the latter results, the orientational versatility of the system, together with the nondenaturing properties of CHAPSO, makes this system useful in spectroscopic studies of membrane-associated phenomena.