| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

* Department of Physics & Astronomy, Rice University, Houston, Texas 77251; and
National Synchrotron Light Source, Brookhaven National Laboratory, Upton, New York 11973
Correspondence: Address reprint requests to Dr. Huey W. Huang, Dept. of Physics & Astronomy, Rice University, Houston, Texas 77251-1892. Tel.: 713-348-4899; Fax: 713-348-4150; E-mail: hwhuang{at}rice.edu.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
In the standard MAD method, one uses an equation relating three quantities; i.e., the magnitude of the normal diffraction amplitude of the host molecules, the magnitude of the normal diffraction amplitude of the anomalous atoms, and the relative phase angle between the two amplitudes, with coefficients that are wavelength dependent. A set of diffraction measurements are then recorded at wavelengths below and above the absorption edge of the anomalous atom so as to vary these coefficients as much as possible to allow the three unknowns to be determined. We previously applied this standard method to a gramicidin-lipid bilayer system (17
). This required the sample to have a well-defined thickness; hence the sample was sandwiched between two substrates, with one of them being a thin, polished beryllium plate for x-ray transmission. However, more commonly used and much more conveniently prepared lipid samples are deposited on one substrate only. Such samples can be subject to a rapid hydration change, which is important for synchrotron radiation experiment. One drawback of these open samples is that the thickness of the jellylike lipid deposition is known only approximately due to at least two reasons: 1), although the amount of lipid can be precise, the lipid is difficult to confine to a well-defined area; and 2), the thickness may vary with hydration. As a result the length of the beam path through the sample is known only approximately. Therefore the correction for the x-ray absorption can be significantly uncertain, particularly if the absorption coefficient is large and sensitive to wavelength.
To alleviate the difficulty of absorption correction, we choose our x-ray energies below the absorption edge, where the x-ray absorption coefficient is relatively small and almost independent of wavelength. We measured diffraction at eight different wavelengths so as to solve the unknowns by a straight-line fitting instead of solving nonlinear equations. Even assuming that the absorption correction is error free, the multiple solutions of the simultaneous nonlinear equations constructed from MAD measurement can be ambiguous. This is because the coefficients in the equations inevitably contain experimental errors. We can imagine that if the coefficients of the equations were precise, the correct solutions would be reproduced by overdeterminations (by using more than three different wavelengths). One could then single out the correct solutions from the non-physical ones. But this is in general not the case due to the imprecise coefficients. We found the method of straight-line fitting much more straightforward by comparison.
We will use the lamellar phase of distearoyl phosphatidylcholine with brominated chains to illustrate this procedure. Interestingly the commonly used swelling method (4
,5
) cannot resolve the phase problem of this lamellar system. We will demonstrate that the MAD method correctly determines the phases.
| EXPERIMENTAL METHODS |
|---|
|
|
|---|
100
surface, P-doped), 300-µm thick, were purchased from Virginia Semiconductor (Fredericksburg, VA).
Preparation for oriented samples followed the method described in Ludtke et al (18
). Di18:0(9,10dibromo)PC was first dissolved in a 1:1 trifluoroethanol (TFE)-chloroform solvent and then uniformly deposited onto a clean, flat silicon substrate. The organic solvent was evaporated in vacuum or open air for
1 h. The deposit was then hydrated with saturated water vapor and incubated in an oven at 35°C overnight. The result was 0.4 mg of lipid spread over an area of 10 x 10 mm2, thus with an approximately uniform thickness of 4 µm. For diffraction experiment, the sample was kept inside a humidity-temperature chamber (2
). The substrate was attached to a temperature-controlled aluminum mount by heat-sink paste. Directly facing the sample surface was a water reservoir, where the water temperature was adjusted to vary the relative humidity (RH) inside the sample chamber. A temperature transducer (AD590, Analog Devices, Norwood, MA) and a relative humidity sensor (HC-600, Ohmic, Easton, MD) were mounted close to the sample to monitor the sample condition. The outputs from the sensing elements were fed to PID feedback control circuits, which in turn powered two sets of Peltier modules (Melcor, NJ), one for heating or cooling the sample and another for heating or cooling the water reservoir. The chamber was covered by a double-layered insulating wall with Kapton windows for the passage of x ray. Between the two layers, a resistive heating coil maintained the surface temperature of the chamber above that of the sample so as to avoid water condensation on the Kapton windows.
At room temperature the lipid formed hydrated bilayers parallel to the substrate. Diffraction from the lamellar phase was measured at 25°C and 90%RH. This unusual choice of humidity level requires an explanation. It is well known that for the majority of lipids, peak broadening and progressive weakening of the reflection orders occur when the humidity level exceeds
98%RH, due to the undulation fluctuations of the membranes in water (19
). For this lipid, the diffraction pattern began to lose high orders and exhibit peak broadening above
92%RH. This was found repeatedly in experiments with freshly prepared samples, except that once or twice the disordering started around 96%RH.
X-ray experiment was performed at the beamline X21 of the National Synchrotron Light Source, Brookhaven National Laboratory (Upton, NY). The setup was similar to the one described in Yang and Huang (2
). The x-ray beam was collimated by two sets of slits before the sample chamber, resulting in a beam size of 0.5 x 0.5 mm2 at the sample. Diffraction by the lamellar phase was recorded on a MarCCD detector (Mar USA, Evanston, IL). A niobium (Nb) attenuator was used to keep the first two orders from saturating the detector. The intensity of the incident beam was monitored by a Bicron scintillation detector (Saint-Gobain Crystals, Newbury, OH) that measured the elastic scattering from a 0.9-µm-thick polyethylene film inserted in the incident beam. The detector was positioned at 90° angle from the incident beam and perpendicular to the incident polarization.
Wavelength dependence of detectors
It is imperative to measure and correct for the wavelength dependence of the detectors and the attenuator. All the diffraction intensities were normalized to a fixed incident photon flux (photon number per area per time). We used an air-filled ion chamber (with 8-µm-thick Kapton windows) as the reference whose wavelength dependence is contained in the expression
The first factor
is the absorption ratio of the x ray passing through the ion chamber, where
is the absorption coefficient for air and D the length of the ion chamber. The second factor
comes from the energy of photon, because an ion chamber is an energy detector. An empirical expression for
is given by the Center for X-Ray Optics of the Lawrence Berkeley Laboratory (20
). For D = 8 cm at 1 atm and 25°C,
The wavelength dependence of the detectors and of the Nb attenuator is included in Table 1 for reference. Note that the insignificant wavelength dependence for the Bicron detector shown in Table 1 could be misleading; the detector has much stronger wavelength dependence outside of the range of wavelength shown. (Ion chambers were not used in the diffraction experiment due to space limitation.)
|
by using the theoretical values calculated by Cromer and Libermann (21
was calculated by the dispersion relation using the CHOOCH program by Evans and Pettifer (22
= 0.5 (in the unit of electron) for each successive energy as shown in Table 1 where the values of
and
are listed.
Multi-wavelength anomalous diffraction
Diffraction by the lamellar phase was measured by rotating the substrate 1°/s (2
) from incident angle
0° to
10°, at eight x-ray wavelengths chosen above. Well-aligned lipid lamella diffract strongly (which we routinely measure on a sealed-tube x-ray diffractometer; see below), thus a moderate beam intensity was used by narrowing the slits. The beam was blocked between scans so the sample was exposed to radiation only during data collection. We first completed the scans for eight wavelengths at the same sample position. Then we displaced the substrate to a previously unexposed sample position and completed the eight scans in the reversed order of wavelengths. We found that the results were consistent with each other indicating no deterioration effect from radiation damage. After the experiment the sample was examined by thin layer chromatography as described in Yang et al. (23
); the result did not show extra spots as compared with fresh lipid. All data were also reproduced by at least two freshly prepared samples.
We also applied the swelling method to the samples by normal diffraction collected on a sealed-tube x-ray generator and a Huber four-cycle goniometer, with a line-focused (1 x 10 mm) Cu K
source (
= 1.542 Å) operating at 40 kV and 1530 mA. The characteristics of this experimental setup have been described in detail previously (24
,25
).
Reduction of the MAD data
The intensities of the diffraction peaks were integrated directly on the detector image in two ways. 1), The peaks were fit by two-dimensional Gaussian functions plus a background, and then integrated. 2), The peaks were first integrated in the direction parallel to the substrate surface over a width slightly wider than the apparent peak widths. The result was plotted along the substrate normal. On this one-dimensional profile, the background was obtained by using the intensities between the peaks and extrapolated into the peak regions. After the background removal, each peak was fit with a Gaussian and then integrated. The results of the two methods were consistent with each other.
The integrated intensities were corrected for the polarization factor, the Lorentz factor, diffraction volume, and x-ray absorption (2
). These are the standard corrections for normal lamellar diffraction. For MAD, there is an additional wavelength-dependent correction of
3 (26
). All these corrections are straightforward, except for the x-ray absorption. In general, the jellylike lipid samples made by the deposit method do not have a well-defined thickness, therefore there is an uncertainty in the absorption correction. This uncertainty is greatly magnified if the multi-wavelength measurement includes energies above the K-edge where the absorption coefficient is large and sensitive to wavelength. Below the K-edge, the absorption coefficient is almost constant (Table 1). Thus the absorption correction does not affect the relative magnitudes of the intensities measured at different wavelengths. For this reason we limited our MAD measurement to energies below the K-edge. As will be seen below, even in the region below the absorption edge there are sufficient dispersion variations in diffraction intensity for MAD analysis.
The diffraction amplitude from a system containing atoms of an element with anomalous scattering factor
is written as
![]() | (1) |
is the normal scattering factor of atom j at position
and
and
are the real and imaginary parts of the anomalous scattering factor; their values for Br are shown in Table 1. The index j includes all atoms except for the anomalous atoms. The index k includes all the anomalous atoms. Fo is the normal diffraction amplitude of the whole system. F2 is the normal diffraction amplitude of the anomalous atoms alone. We will assume that the unit cell of the lipid structure is centrosymmetric, so that both Fo and F2 are real (rather than complex) quantities and their phases are the signs of the amplitudes. The following analysis will test whether this assumption is valid. From Eq. 1, we have
![]() | (2) |
However, the second term is
1% of the first term, because at energies below the absorption edge
0.5, which is
10% of
(Table 1). Therefore we have the approximate relation
![]() | (3) |
by
since
is negative below the absorption edge.
We plotted
against
for each peak (Fig. 1). The data for every peak appear to satisfy a linear relation (i.e., approximately fit to a straight line). Because Eq. 3 is the combined result of a), both Fo and F2 are real quantities and b), the linear approximation from Eq. 2 to Eq. 3, we conclude that both the assumption of centrosymmetry and the assumption of linear approximation are valid. In each plot, we fitted the data with a straight line (Fig. 1). (The correlation coefficient of the fit is also given in the figure.) From Eq. 3, we see that the intercept of the fitted line gives
; the magnitude of the slope gives
; and the sign of the slope gives the sign of
(Fig. 1, insets).
|
's are determined, so are those of
's.
Distribution of the bromine label and electron density profile of the bilayer
We use the amplitudes
to build the Patterson function (26
) for the Br distribution
![]() | (4) |
The result is displayed in Fig. 2. The Patterson function shows pairs of peaks at position z and Dz. The pair of
is from self-correlation. Another pair at
is from intercorrelation, indicating that there are two Br peaks in the bilayer separated by a distance
The data can be modeled by a density
![]() | (5) |
the Gaussian width for each peak. The density extends beyond the unit cell by the periodic condition. It is easy to show that the Patterson function of
is
![]() | (6) |
The Patterson function of a model distribution with parameters d = 7.5 Å and
= 4 Å closely reproduces the experimental Patterson function as shown in Fig. 2. The points to note are: 1), the width of a Patterson peak is
times that of the corresponding peak in real space, and 2), the amplitudes of the intercorrelation peaks (located at
and
on the Patterson coordinate) are one-half of the self-correlation peaks at
and
|
that closely reproduces the experimental Patterson, one calculates the model amplitudes
![]() | (7) |
The phases of
determine the phases of experimental F2 and, from the sign of Fo / F2 obtained from the MAD analysis, those of Fo.
With their phases (signs) determined, F2 and Fo are used to construct the experimental electron density profiles for the Br distribution (Fig. 3 A) and for the whole lipid (Fig. 3 B), respectively. The experimental electron densities are normalized to the real densities by three parameters a, b, and
:
and
The three parameters are determined by the three conditions: 1), Br electron density vanishes in the water region, 2), the integration of
from D/2 to D/2 equals to the total number of Br electrons in two lipids divided by the lipid cross section, 3), the integration of
from D/2 to D/2 equals to the total number of electrons in two lipids plus 20 water molecules divided by the lipid cross section. The lipid cross section was obtained by the use of the phosphate-to-phosphate distance in Fig. 3 B and the relevant lipid volume (27
). (The number of water molecules used here is somewhat arbitrary. This number could be more precisely measured (27
), but was not done, because it would not affect the discussion here.)
|
92% (peak broadening and losing high orders). Thus the swelling experiment was performed from 89% to 91%RH. As one sees from Fig. 4 A, the sign of the third order cannot be clearly determined by the swelling method. In fact judging from the overall agreement between the Shannon constructions (24
for the phases, instead of
as determined by the MAD method. Also there is no obvious reason, based on a bilayer structure, to reject either one of the two electron density profiles built with either a positive or a negative third order (Fig. 4, B and C). (For example, in Fig. 4 C, the high central region relative to the water region could be explained as due to the contribution of bromines to the central region.)
|
| DISCUSSION |
|---|
|
|
|---|
determined by the MAD method for the di18:0(9,10dibromo)PC bilayers are correct, we measured a series of lipid mixtures of di18:1PC (DOPC) and di18:0(9,10dibromo)PC. Starting with pure DOPC where the third order is positive (as determined by the swelling method), the magnitude of the third order first diminished with the increasing fraction of di18:0(9,10dibromo)PC and then increased (Fig. 5). This is consistent with the sign change of the third order as the fraction of di18:0(9,10dibromo)PC increases. That explains the phases
for pure di18:0(9,10dibromo)PC. The reason for the positive third order for most lipids is that the position of the headgroup is about D/3 from the center of the bilayer. However, the bromine peaks in the middle of the chain give a negative third order. Thus the sign of the third order depends on their relative contributions. Apparently the bromine peaks dominate the sign when there are two bromines per chain.
|
In protein crystallography, heavy atoms are attached to proteins isomorphously. The purpose of MAD is to solve the structures of the native proteins. For lipids, heavy atom labels are used in a different way, because such labels often alter the property of the original lipid. For example, the property of di18:0(9,10dibromo)PC is somewhat between that of di18:0PC and of di18:1PC. By osmotic pressure it undergoes a transition from the lamellar phase to a rhombohedral phase, similar to di18:1PC but at a different phase boundary (29
), whereas di18:0PC does not exhibit such a phase transition. The mixture of cholesterol and di18:0(9,10dibromo)PC can be induced to a rhombohedral phase with or without peptides ((30
); the authors' unpublished data). In such cases, the unit cell structure could be a fusion intermediate state, called a stalk (1
,2
) or a transmembrane pore (3
). With the application of anomalous diffraction, the bromine labels serve the dual purposes of 1), phase determination thus resolving the electron density distribution within the unit cell; and 2), highlighting the lipid chains of the PC thus making it possible to distinguish the distributions of PC and cholesterol separately.
The MAD method supplements the existing methods of phase determination including swelling (4
,5
), pattern recognition (6
,7
), and methyl trough search (8
,9
). Although in one previous case a combination of swelling and pattern recognition solved the phase problem for a rhombohedral phase (1
,2
), we found that the swelling method is often not applicable to nonlamellar phases either because the range of swelling is insufficient or because swelling changes the unit cell structure (e.g., in the inverted hexagonal phase). The method of pattern recognition relies on the presumed invariance of density moments between different phases of the same lipid (6
,7
). This does not apply to the cases where the lipid components might redistribute depending on the degree of bending (31
). The use of methyl trough search (8
,9
) is limited to the cases where the mathematical property of the monolayer surface is known. In comparison, the MAD method is founded on a firmer theoretical basis (10
), as demonstrated by its general applicability in protein crystallography (10
16
). We believe that its applications to lipid systems are yet to be explored.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Submitted on December 23, 2005; accepted for publication April 11, 2006.
| REFERENCES |
|---|
|
|
|---|
2. Yang, L., and H. W. Huang. 2003. A rhombohedral phase of lipid containing a membrane fusion intermediate structure. Biophys. J. 84:18081817.
3. Yang, L., T. M. Weiss, and H. W. Huang. 2000. Crystallization of antimicrobial pores in membranes: magainin and protegrin. Biophys. J. 79:20022009.
4. Blaurock, A. E. 1971. Structure of the nerve myelin membrane: proof of the low-resolution profile. J. Mol. Biol. 56:3552.[CrossRef][Medline]
5. Torbet, J., and M. H. Wilkins. 1976. X-ray diffraction studies of lecithin bilayers. J. Theor. Biol. 62:447458.[CrossRef][Medline]
6. Luzzati, V., A. Tardieu, and D. Taupin. 1972. A pattern-recognition approach to the phase problem: application to the X-ray diffraction study of biological membranes and model systems. J. Mol. Biol. 64:269286.[CrossRef][Medline]
7. Mariani, P., V. Luzzati, and H. Delacroix. 1988. Cubic phases of lipid-containing systems. Structure analysis and biological implications. J. Mol. Biol. 204:165189.[CrossRef][Medline]
8. Harper, P. E., and S. M. Gruner. 2000. Electron density modeling and reconstruction of infinite periodic minimal surfaces (IPMS) based phases in lipid-water systems. I. Modeling IPMS-based phases. Eur. Phys. J. E. 2:217228.[CrossRef]
9. Harper, P. E., S. M. Gruner, R. N. Lewis, and R. N. McElhaney. 2000. Electron density modeling and reconstruction of infinite periodic minimal surfaces (IPMS) based phases in lipid-water systems. II. Reconstruction of D surface based phases. Eur. Phys. J. E. 2:229245.[CrossRef]
10. Karle, J. 1989. Macromolecular structure from anomalous dispersion. Phys. Today. 42:2229.
11. Guss, J. M., E. A. Merritt, R. P. Phizackerley, B. Hedman, M. Murata, K. O. Hodgson, and H. C. Freeman. 1988. Phase determination by multiple-wavelength x-ray diffraction: crystal structure of a basic "blue" copper protein from cucumbers. Science. 241:806811.
12. Hendrickson, W. A., J. L. Smith, R. P. Phizackerley, and E. A. Merritt. 1988. Crystallographic structure analysis of lamprey hemoglobin from anomalous dispersion of synchrotron radiation. Proteins. 4:7788.[CrossRef][Medline]
13. Hendrickson, W. A., J. R. Horton, H. M. Murthy, A. Pahler, and J. L. Smith. 1989. Multiwavelength anomalous diffraction as a direct phasing vehicle in macromolecular crystallography. Basic Life Sci. 51:317324.[Medline]
14. Kahn, R., R. Fourme, R. Bosshard, M. Chiadmi, J. L. Risler, O. Dideberg, and J. P. Wery. 1985. Crystal structure study of Opsanus tau parvalbumin by multiwavelength anomalous diffraction. FEBS Lett. 179:133137.[CrossRef][Medline]
15. Murthy, H. M., W. A. Hendrickson, W. H. Orme-Johnson, E. A. Merritt, and R. P. Phizackerley. 1988. Crystal structure of clostridium acidi-urici ferredoxin at 5-Å resolution based on measurements of anomalous X-ray scattering at multiple wavelengths. J. Biol. Chem. 263:1843018436.
16. Hendrickson, W. A. 1991. Determination of macromolecular structures from anomalous diffraction of synchrotron radiation. Science. 254:5158.
17. Liu, W., T. Y. Teng, Y. Wu, and H. W. Huang. 1991. Phase determination for membrane diffraction by anomalous dispersion. Acta Crystallogr. A47:553559.[CrossRef]
18. Ludtke, S., K. He, and H. Huang. 1995. Membrane thinning caused by magainin 2. Biochemistry. 34:1676416769.[CrossRef][Medline]
19. Chen, F. Y., W. C. Hung, and H. W. Huang. 1997. Critical swelling of phospholipid bilayers. Phys. Rev. Lett. 79:40264029.[CrossRef]
20. Laboratory, Center for X-ray Optics of the Lawrence Berkeley. http://www-cxro.lbl.gov/optical_constants. [Online.].
21. Cromer, D. T., and D. Libermann. 1970. Relativistic calculation of anomalous scattering factors for X rays. J. Chem. Phys. 53:18911898.[CrossRef]
22. Evans, G., and R. F. Pettifer. 2001. CHOOCH: a program for deriving anomalous-scattering factors from X-ray fluorescence spectra. J. Appl. Crystallogr. 34:8286.[CrossRef]
23. Yang, L., T. A. Harroun, T. M. Weiss, L. Ding, and H. W. Huang. 2001. Barrel-stave model or toroidal model? A case study on melittin pores. Biophys. J. 81:14751485.
24. Olah, G. A., H. W. Huang, W. H. Liu, and Y. L. Wu. 1991. Location of ion-binding sites in the gramicidin channel by X-ray diffraction. J. Mol. Biol. 218:847858.[CrossRef][Medline]
25. Weiss, T. M., L. Yang, L. Ding, A. J. Waring, R. I. Lehrer, and H. W. Huang. 2002. Two states of cyclic antimicrobial peptide RTD-1 in lipid bilayers. Biochemistry. 41:1007010076.[CrossRef][Medline]
26. Warren, B. E. 1969. X-Ray Diffraction. Section 4.1, 9.3. Dover Publications, New York.
27. Nagle, J. F., and S. Tristram-Nagle. 2000. Structure of lipid bilayers. Biochim. Biophys. Acta. 1469:159195.[Medline]
28. McIntosh, T. J., and P. W. Holloway. 1987. Determination of the depth of bromine atoms in bilayers formed from bromolipid probes. Biochemistry. 26:17831788.[CrossRef][Medline]
29. Yang, L., L. Ding, and H. W. Huang. 2003. New phases of phospholipids and implications to the membrane fusion problem. Biochemistry. 42:66316635.[CrossRef][Medline]
30. Pan, D., W. Wang, W. Liu, L. Yang, and H. W. Huang. 2006. Chain packing in the inverted hexagonal phase of phospholipids: a study by X-ray anomalous diffraction on bromine-labeled chains. J. Am. Chem. Soc. In press.
31. Ding, L., T. M. Weiss, G. Fragneto, W. Liu, L. Yang, and H. W. Huang. 2005. Distorted hexagonal phase studied by neutron diffraction: lipid components demixed in a bent monolayer. Langmuir. 21:203210.[CrossRef][Medline]
This article has been cited by other articles:
![]() |
S. Qian, W. Wang, L. Yang, and H. W. Huang Structure of the Alamethicin Pore Reconstructed by X-Ray Diffraction Analysis Biophys. J., May 1, 2008; 94(9): 3512 - 3522. [Abstract] [Full Text] [PDF] |
||||
![]() |
W.-C. Hung, M.-T. Lee, F.-Y. Chen, and H. W. Huang The Condensing Effect of Cholesterol in Lipid Bilayers Biophys. J., June 1, 2007; 92(11): 3960 - 3967. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. Wang, L. Yang, and H. W. Huang Evidence of Cholesterol Accumulated in High Curvature Regions: Implication to the Curvature Elastic Energy for Lipid Mixtures Biophys. J., April 15, 2007; 92(8): 2819 - 2830. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |