help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Contreras, L. M.
Right arrow Articles by Prieto, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Contreras, L. M.
Right arrow Articles by Prieto, M.

Biophys J, May 2001, p. 2273-2283, Vol. 80, No. 5

Interaction of alpha -Melanocyte Stimulating Hormone with Binary Phospholipid Membranes: Structural Changes and Relevance of Phase Behavior

Lellys M. Contreras,* Rodrigo F. M. de Almeida,dagger José Villalaín,* Aleksandre Fedorov,dagger and Manuel Prietodagger

 *Centro de Biología Molecular y Celular, Universidad "Miguel Hernández," E-03206 Elche-Alicante, Spain; and  dagger Centro de Química-Física Molecular, Instituto Superior Técnico, P-1049-001 Lisboa, Portugal


    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
CONCLUDING REMARKS
REFERENCES

The interaction of alpha -melanocyte stimulating hormone (alpha -MSH) with negatively charged binary membrane systems composed of either 1,2-dimyristoyl-sn-glycero-3-phosphocholine/1,2-dimyristoyl-sn-glycero-3-[phospho-rac-(1-glycerol)], (DMPC/DMPG) or DMPC/1,2-dimyristoyl-sn-glycero-3-phosphate (DMPC/DMPA), both at a 3:1 ratio, was studied using complementary techniques (differential scanning calorimetry, infrared and ultraviolet absorption spectroscopy, and steady-state and time-resolved fluorescence). The peptide structure in buffer, at medium to high concentrations, is a mixture of aggregated beta -strands and random coil, and upon increasing the temperature the random coil configuration becomes predominant. At low concentrations (micromolar) there are essentially no aggregates. When in interaction with the lipidic systems this transition is prevented and the peptide is stabilized in a specific conformation different from the one in solution. The incorporation of alpha -MSH into phosphatidic acid-containing systems produced a significant alteration of the calorimetric data. Lateral heterogeneity can be induced by the peptide in the DMPA-containing mixture, at variance with the one of DMPG. In addition, the lipid/water partition coefficient for the peptide in the presence of DMPC/DMPA is greater in the gel phase as compared to the fluid phase. From the high values of limiting anisotropies it can be concluded that the peptide presents a very reduced rotational dynamics when in interaction with the lipids, pointing out to a strong interaction. Overall, these results show that the structure and stability of alpha -MSH in a negatively charged membrane environment are substantially different from those of the peptide in solution, being stabilized in a specific conformation that could be important to eliciting its biological activity.


    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
CONCLUDING REMARKS
REFERENCES

alpha -Melanocyte stimulating hormone (alpha -MSH) is a peptide hormone known for its role in regulating skin pigmentation in vertebrates (Eberle, 1988). The primary structure of alpha -MSH is Ac-Ser-Tyr-Ser-Met-Glu-His-Phe-Arg-Trp-Gly-Lys-Pro-Val-NH2. It appears that the peptide has no preferred structure in water (Biaggi et al., 1997), being very flexible, whereas all its synthetic analogs with superpotent biological activity present a beta -turn (or sometimes another kind of turn) stabilized within their central region comprising residues 6-9, i.e., His-Phe-Arg-Trp (Sawyer et al., 1980, 1982; Al-Obeidi et al., 1989). This Trp-containing region of the peptide is furthermore the minimum melanotropic message sequence, essential for ligand binding and biological function (Hruby et al., 1987).

alpha -MSH has a net +1 charge at physiological pH, and a mean hydrophobicity and hydrophobic moment that are indicative of a receptor or specific lipid-mediated membrane interaction, and not a surface-seeking characteristic peptide (González et al., 1996). alpha -MSH should be, according to these authors, classified as a secondary amphiphilic peptide with a low penetration into the lipid monolayer and a small effect on its thermotropic properties.

It was recently demonstrated by fluorescence spectroscopy that alpha -MSH and analogs interact with lipid vesicles composed of one anionic phospholipid, as concluded from alterations of its photophysical parameters (Ito et al., 1993; Macêdo et al., 1996). The fluorescence parameters were dependent on the lipid used, but the interaction was greater in the liquid crystalline phase as compared to the gel phase. The interaction was mainly electrostatic, because changes in the peptide (Ito et al., 1993), bilayer (Biaggi et al., 1996), or peptide/monolayer (González et al., 1996) parameters could not be detected by a wide variety of techniques when using zwitterionic phospholipids. In a recent molecular dynamics study (Pascutti et al., 1999) the peptide was reported to acquire a stable structure in a low dielectric constant medium, as compared to the structure in water, where a hydrophilic core was surrounded by a hydrophobic surface, forming a beta -turn structure.

All these results suggest that the lipid matrix of the cell membrane could be a key factor for the peptide biological activity, namely by stabilizing an appropriate conformation. In line with this, we have studied in detail the interaction of alpha -MSH with two-component lipid bilayers composed of a zwitterionic phospholipid (1,2-dimyristoyl-sn-glycero-3-phosphocholine, DMPC) and a negatively charged one (1,2-dimyristoyl-sn-glycero-3-[phospho-rac-(1-glycerol)], DMPG, or 1,2-dimyristoyl-sn-glycero-3-phosphate, DMPA) by using high-sensitivity differential scanning calorimetry (DSC), Fourier-transform infrared spectroscopy (FTIR), and steady-state and time-resolved fluorescence spectroscopy. The main issues addressed were the extent of lipid/peptide interaction, the effect of the peptide on the structure and thermotropism of phospholipids, and the structural features of the peptide, both in the aqueous and lipidic phases.


    MATERIALS AND METHODS
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
CONCLUDING REMARKS
REFERENCES

Chemicals

alpha -MSH and D2O were obtained from Sigma Chemical Co. (St. Louis, MO). The lipids DMPC, 1,2-dimyristoyl-d54-sn-glycero-3-phosphocholine (DMPCd54), DMPA, and DMPG were obtained from Avanti Polar Lipids (Birmingham, AL). The buffer used was 10 mM HEPES (4-(2-hydroxyethyl)-1-piperazine ethanesulfonic acid) at pH 7.4, 20 mM NaCl, and 0.1 mM EDTA. All reagents were used as received.

Differential scanning calorimetry

Samples containing 2.6 mmol of phospholipids (DMPCd54 plus either DMPA or DMPG at a 3:1 molar ratio) dissolved in chloroform/methanol (1:1 v/v) were mixed alone or with alpha -MSH, also in the same solvent, to give a final phospholipid to peptide molar ratio of 10:1. The mixtures were dried under vacuum for 3 h to remove all traces of the organic solvents. Then, multilamellar vesicles (MLV) were formed in 1.4 ml buffer heated at a temperature ~10°C above the temperature of the gel-to-liquid-crystalline phase transition (Tm) of the mixture and frozen at -80°C, this process being repeated five times. Differential scanning calorimetry was performed in a high-resolution Microcal MC-2 calorimeter (Microcal Inc., Northampton, MA). Differences in the heat capacity between the sample and the reference cell were obtained by raising the temperature at a constant rate of 45°C/h over a temperature range of 5-60°C. The excess heat capacity functions were obtained after baseline subtraction and correction for the instrument time response. Unless otherwise stated, the third scan was used for transition temperature and enthalpy calculations.

Infrared spectroscopy

Samples containing the same quantity of lipids and peptide as described above were hydrated in the same deuterated buffer at pD = 7.4. In order to wash unbound peptide, vesicles were pelleted by centrifugation and washed with the same deuterated buffer, freeze-thawed, and centrifuged again. Approximately 20 µl of the pelleted sample were placed between two CaF2 windows separated by a 50-µm-thick Teflon spacer in a liquid demountable cell (Harrick, Ossining, NY). Peptide solutions of 12.8 mM concentration obtained under bath sonication were used for obtaining the spectra in buffer. The spectra were obtained in a Nicolet 580 FTIR spectrometer using a deuterated triglycine sulfate detector. Each spectrum was obtained by collecting 200 scans with a triangular function at a resolution of 2 cm-1. For temperature studies, samples were scanned between 5 and 70°C with 2°C intervals and a 2-min delay between each consecutive scan. Data treatment and resolution enhancement methods performed according to previously published methods (Kauppinen et al., 1981) were made using Spectra Calc or Grams software (Galactic Industries, Salem, NH). The frequency of the band due to the CH2 or CD2 symmetric stretching vibration was calculated from the center of gravity taking the top 10 points of each specific band and fitted to a Gaussian band. The criterion used for buffer subtraction in the C==O and amide regions was the removal of the band near 1210 cm-1, and a flat baseline between 1800 and 2100 cm-1.

Absorption and fluorescence studies

For fluorescence studies, large unilamellar vesicles (LUV) prepared by the extrusion technique were used (Hope et al., 1985). First, MLV were made as described for DSC, except that DMPC was used instead of DMPCd54. Extrusion cycles were then performed on MLV suspensions with Nucleopore polycarbonate filters, first with 0.4 µm pore diameter until no resistance was offered, followed by eight cycles with 0.1 µm pore diameter filters. The final concentration of total phospholipid in LUV suspensions was determined by phosphorus analysis (McClare, 1971).

Direct dissolution of the peptide in a small volume of buffer (~0.1 absorbance at 290 nm) was achieved by vortexing and bath sonication (Bandelin Sonorex RK156). Then, the samples were prepared by taking aliquots of the peptide solution and adding the same volume of lipid suspension plus buffer in different proportions. The final peptide absorbance at 290 nm was ~0.05 in all samples. This corresponds to a peptide concentration of ~30 µM (epsilon  ~ 3500 M-1 cm-1 (Wetlaufer, 1962)).

Absorbance spectra were carried out in a Shimadzu UV-3101PC absorption spectrophotometer using spectral bandwidths of 2.0 nm. The correction for light scattering was carried out subtracting the baseline from an identical lipidic suspension. Steady-state fluorescence spectra were obtained in a SLM-Aminco 8100 Series 2 spectrofluorimeter. Emission spectra were corrected using standard emission spectra of L-Tyr and L-Trp (Chen, 1967). 5 mm × 5 mm quartz cuvettes were used and temperature was controlled up to ±0.5°C in a thermostatted cuvette holder. Spectral bandwidth was 2 nm for excitation and 4 nm for emission in most measurements.

The time-resolved instrumentation (correlated single-photon timing technique) was previously described (Loura et al., 1996). Excitation was at 295 nm and emission (at 350 nm) was detected at the magic angle (54.7°) relative to the vertically polarized beam. The number of counts on the peak channel was 20,000 and the number of channels per curve used for analysis was 800, with 8.4 or 11 ps/channel at 37°C and 14 or 15.3 ps/channel at 20°C. Data analysis was carried out using a nonlinear, least-square iterative convolution method based on the algorithm of Marquardt (1963). The goodness-of-fit was judged from the reduced chi 2, weighted residuals, and autocorrelation plots.

The average lifetime <A><AC>&tgr;</AC><AC>&cjs1171;</AC></A> of a fluorophore with a complex fluorescence decay is defined as (e.g., Lakowicz, 1999),
<A><AC>&tgr;</AC><AC>&cjs1171;</AC></A>=<LIM><OP>∑</OP><LL><UP>i</UP></LL></LIM> a<SUB><UP>i</UP></SUB>&tgr;<SUP><UP>2</UP></SUP><SUB><UP>i</UP></SUB>&cjs1134;<LIM><OP>∑</OP></LIM> a<SUB><UP>i</UP></SUB>&tgr;<SUB><UP>i</UP></SUB> (1)
where ai are the normalized pre-exponentials (amplitudes) and tau i the lifetime components. For the purpose of partition coefficients determination, the lifetime-weighted quantum yield, < tau > , leads to simpler formalisms, and it will be used throughout this work
⟨&tgr;⟩=<LIM><OP>∑</OP><LL><UP>i</UP></LL></LIM> a<SUB><UP>i</UP></SUB>&tgr;<SUB><UP>i</UP></SUB> (2)
The time-resolved fluorescence anisotropies were determined using Glan-Thompson polarizers, and calculated (e.g., Lakowicz, 1999) according to:
r(t)=<FR><NU>I<SUB><UP>VV</UP></SUB>(t)−GI<SUB><UP>VH</UP></SUB>(t)</NU><DE>I<SUB><UP>VV</UP></SUB>(t)+2GI<SUB><UP>VH</UP></SUB>(t)</DE></FR> (3)
where the different intensities are the vertical (IVV) and horizontal (IVH) components of the fluorescence emission with excitation vertical to the emission axis. The intensity decays of polarized light IVV(t) and IVH(t) were obtained with the same accumulation time. The instrumental factor G = IHV/IHH for the time-resolved instrumentation is G = 1, once a scrambler is used before the detector. The limiting values of anisotropy at infinite time (rinfinity ) (Best et al., 1987), were determined directly from the plot of r(t) vs. t. For this purpose there is no need of deconvolution with the laser pulse profile.


    RESULTS
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
CONCLUDING REMARKS
REFERENCES

Differential scanning calorimetry

DSC was used to study the effect of alpha -MSH on phospholipid mixtures composed of either DMPCd54/DMPA or DMPCd54/DMPG, 3:1. Thermograms recorded for the phospholipid mixture DMPCd54/DMPA, either in the absence or in the presence of alpha -MSH, are shown in Fig. 1 A. In the absence of alpha -MSH, the transition peak for the mixture DMPCd54/DMPA was broad and asymmetrical, with a maximum at ~22°C, in agreement with the literature (Garidel et al., 1997), and no pretransition was observed at variance with the pure phospholipid. Tm for the mixture is close to the one observed for the pure DMPCd54, as Tm is ~20°C for DMPCd54 (Muga et al., 1991) and 50°C for DMPA (Garidel et al., 1997). The asymmetric broadening can be related to a non-ideal behavior of the phospholipid mixture. Upon incorporation of alpha -MSH a small shift in the onset temperature of the gel-to-liquid-crystalline phase transition (Tc) of ~2°C is observed, as is the appearance of sharper peaks within the broad component. These peaks (at ~23°C, 27°C, and 32°C) were only apparent in the pure phospholipid mixture. Despite the significant effect on the thermogram, alpha -MSH did not significantly affect the enthalpy change of the gel-to-liquid-crystalline phase transition of the phospholipid mixture.



View larger version (14K):
[in this window]
[in a new window]
 
FIGURE 1   DCS heating thermograms for aqueous dispersions of mixtures at a 3:1 molar ratio of DMPCd54/DMPA (A) and of DMPCd54/DMPG (B), with (dotted lines) or without (solid lines) alpha -MSH at a phospholipid/peptide molar ratio of 10:1.

Thermograms recorded for the phospholipid mixture DMPCd54/DMPG, either in the absence or in the presence of alpha -MSH, are shown in Fig. 1 B. The mixture DMPCd54/DMPG has ideal behavior, and as it is found in the pure phospholipids, presents a pretransition at ~12°C and a highly cooperative main transition at ~22°C, in agreement with the literature (Muga et al., 1991). When alpha -MSH was incorporated in this mixture, the pretransition was abolished, the Tm value decreased slightly to ~21°C, and at the same time the peak broadened slightly, indicating a decrease in the cooperativity. The incorporation of alpha -MSH did not significantly affect the enthalpy change of the gel-to-liquid-crystalline phase transition of the phospholipid mixture.

Determination of the lipid/water partition coefficients of the peptide

Upon incorporation of Trp-containing peptides into membranes, an increase in Trp fluorescence lifetime usually takes place, and there is a concomitant increase in its fluorescence quantum yield. Any of these parameters can be used to quantify the extent of interaction of the peptide with different membrane model systems, as previously shown (Santos et al., 1998), but time-resolved measurements have the advantage of avoiding the artifacts introduced by light-scattering of the lipid suspension. The quantification should be based on the determination of the lipid/water partition coefficient (White et al., 1998), which can be described by
K<SUB><UP>p</UP></SUB>=<FR><NU>n<SUB><UP>L</UP></SUB>/V<SUB><UP>L</UP></SUB></NU><DE>n<SUB><UP>W</UP></SUB>/V<SUB><UP>W</UP></SUB></DE></FR> (4)
where ni stands for moles of peptide in phase i and Vi for volume of phase i. The phase is either aqueous (i = W) or lipidic (i = L). The relationship and assumptions used to determine Kp were previously described (Santos et al., 1998),
⟨&tgr;⟩=<FR><NU>⟨&tgr;⟩<SUB><UP>W</UP></SUB>+⟨&tgr;⟩<SUB><UP>L</UP></SUB>K<SUB><UP>p</UP></SUB><A><AC>V</AC><AC>&cjs1171;</AC></A><SUB><UP>L</UP></SUB>[L]</NU><DE>1+K<SUB><UP>p</UP></SUB><A><AC>V</AC><AC>&cjs1171;</AC></A><SUB><UP>L</UP></SUB>[L]</DE></FR> (5)
where the lifetime-weighted quantum yield of the peptide < tau > is calculated through Eq. 2 for each lipid concentration [L]. A nonlinear fit of Eq. 5 to the data yields the values of < tau > in phase i (i = W, L) and Kp. The values of VL used in this work (0.637 M-1 for vesicles at 20°C and 0.698 M-1 for vesicles at 37°C) are extrapolated from available data in the literature (Marsh, 1990). In Fig. 2 are shown the fluorescence lifetime-weighted quantum yields < tau > of alpha -MSH versus lipid concentration (either DMPC/DMPG or DMPC/DMPA mixtures both at a 3:1 molar ratio) at 20°C and 37°C. < tau > W and the best-fit parameters < tau > L and Kp are shown in Table 1. Although the four Kp values have the same order of magnitude, there are clear differences depending on the lipidic phase; the values for the fluid phase (37°C) are similar for both lipidic mixtures, but for the gel phase (20°C) the partition coefficient for the DMPC/DMPA system is higher as compared to the one for DMPC/DMPG.



View larger version (22K):
[in this window]
[in a new window]
 
FIGURE 2   Lifetime-weighted quantum yield < tau > of alpha -MSH (lambda exc = 295 nm) versus lipid concentration [L] of DMPC/DMPG (3:1) (A and B) and DMPC/DMPA (3:1) (C and D) in the gel phase at 20°C (A and C) and in the fluid phase at 37°C (B and D). Peptide concentration (~30 µM) was kept constant during the experiments. The solid line is the fit of Eq. 5 to the data and the best-fit parameters are shown in Table 1.


                              
View this table:
[in this window]
[in a new window]
 
TABLE 1   Lipid/water partition coefficients (Kp), and 9Trp photophysical parameters of alpha -MSH in different systems

Steady-state fluorescence spectra and anisotropy

The fluorescence spectra of the peptide under different experimental conditions are shown in Fig. 3. At the excitation wavelength of lambda exc = 290 nm, the emission is due to the single Trp residue, and its interaction with the lipidic systems can be appreciated from the spectral shifts observed. The maximum emission wavelength in buffer (lambda em = 349 nm) undergoes a blue-shift (lambda em = 346 nm) in the presence of the lipidic mixture DMPC/DMPG (3:1) in the gel phase, this effect being greater (lambda em = 343 nm) in the case of DMPC/DMPA (3:1), also in the gel phase. For the fluid phase, an identical spectral shift (Delta lambda ) is observed for both lipidic systems (Table 1).



View larger version (18K):
[in this window]
[in a new window]
 
FIGURE 3   Corrected fluorescence emission spectra (lambda exc = 290 nm) of alpha -MSH at 20°C in aqueous solution (solid line), incorporated in gel phase LUV of DMPC/DMPG (3:1) (dashed line), and DMPC/DMPA (3:1) (dotted line). The spectra were corrected using Eq. 7 as described in the text.

The incorporation of the peptide in the membrane is not quantitative. If the partition coefficients Kp are known, the fraction of peptide in the aqueous phase and in the membrane surface can be determined. Because the samples with the highest concentration of lipid still have an appreciable amount of peptide in water (e.g., in the case of DMPC/DMPG at 20°C, ~52% of the peptide), the spectra obtained do not correspond only to alpha -MSH incorporated into the membrane. The spectra shown in Fig. 3 are corrected for this effect and were calculated from the experimental spectra in water and at the highest concentration of lipid. The molar fraction of the peptide in the aqueous phase is given by
x<SUB><UP>W</UP></SUB>=<FR><NU>1</NU><DE>K<SUB><UP>p</UP></SUB><A><AC>V</AC><AC>&cjs1171;</AC></A><SUB><UP>L</UP></SUB>[L]+1</DE></FR> (6)
and in this way the spectrum in the membrane F<UP><SUB>&lgr;</SUB><SUP>L</SUP></UP> is determined from
F<SUP><UP>L</UP></SUP><SUB>&lgr;</SUB>=C · <FENCE>F<SUP><UP>L+W</UP></SUP><SUB>&lgr;</SUB>−x<SUB><UP>W</UP></SUB><FR><NU>1</NU><DE>1+⟨&tgr;⟩<SUB><UP>L</UP></SUB>/⟨&tgr;⟩<SUB><UP>W</UP></SUB></DE></FR> F<SUP><UP>W</UP></SUP><SUB>&lgr;</SUB></FENCE> (7)
where F<UP><SUB>&lgr;</SUB><SUP>i</SUP></UP> is the spectral distribution (unitary area) function of alpha -MSH experimentally determined in water (i = W), and at the highest concentration of lipid (i = L + W). C is a normalization constant. The spectral shifts presented in Table 1 are obtained after this correction.

Steady-state fluorescence anisotropy values, < r> , were determined according to Eq. 3 using the steady-state intensities of the components. An anisotropy value of < r>  = 0.0144 ± 0.0003 (mean standard error of four different experiments) was obtained at 20°C (lambda exc = 290 nm), this being identical to the one obtained in the fluid phase (37°C), < r>  = 0.0146 ± 0.0007.

alpha -MSH structure in water and in the presence of lipids by IR

The infrared spectrum of the amide I' region of alpha -MSH in D2O buffer, at pH 7.4 and at different temperatures, is shown in Fig. 4 A. The amide I' of alpha -MSH in the 1700-1600 cm-1 region is formed by different underlying components that give place to a broad and asymmetric band. Whereas at 15°C the most characteristic feature is a band with maximum intensity at 1618 cm-1, at higher temperatures there is a dominant band with a maximum at 1649 cm-1 (Fig. 4 A). Other smaller components are also discernible in the original spectra at ~1685 and 1586 cm-1. Bands appearing at ~1618 cm-1 accompanied by high-wavenumber counterparts at ~1680-1685 cm-1 have been previously reported as originated from self-associated peptides forming an intermolecular network of hydrogen-bonded beta -strands (Arrondo and Goñi, 1999; Krimm and Bandekar, 1986).



View larger version (36K):
[in this window]
[in a new window]
 
FIGURE 4   Infrared spectra in the amide I' region of alpha -MSH in solution at different temperatures as indicated (A), plot of the 1618 cm-1/1650 cm-1 band intensity ratio versus temperature (B), and two-basis spectra fitting of spectra obtained at 70°C (C) and 15°C (D). Arrows in (A) indicate the direction change of the bands upon increase in temperature.

Bands located at ~1656-1650 cm-1 are usually observed for alpha -helix in D2O solution, and bands appearing at ~1644-1642 cm-1 can be assigned in D2O to unordered (random) structures (Arrondo and Goñi, 1999). The frequency of the 1649 cm-1 band is somewhat below the typical range for most alpha -helical structures. In this way we obtained circular dichroism (CD) spectra (results not shown), which allows unequivocally assigning this band to a random coil structure, coexisting with the previously described structure (beta -strands). In addition, the CD spectra were carried out at a lower concentration (millimolar range), showing that the aggregates detected by IR are not due to any solubilization problem at the higher concentrations.

It is apparent in Fig. 4 A that upon a temperature increase, the area corresponding to the bands at 1618 and 1685 cm-1 (intermolecularly associated beta -strands) decreases, whereas at the same time the area of the 1649 cm-1 (random coil structure) component band increases. A small but noteworthy inflection can be observed in the temperature variation pattern at ~40°C (Fig. 4 B). Decomposition analysis of the infrared spectra at different temperatures showed that these spectra can be described as linear combinations of only two distinct basis spectra, which would correspond to the peptide either in aggregated form or with random coil structure (Fig. 4, C and D). Whereas the area of the aggregated component is about the same as the random coil structure at 15°C (Fig. 4 D), it decreases to 17% at 70°C accompanied by a concomitant increase of the random coil structure (Fig. 4 C). These data suggest that the structural transition observed upon increasing the temperature is a two-state process involving the contribution from two distinct conformations.

The infrared spectra of the amide I' band of alpha -MSH in the presence of either DMPCd54/DMPA or DMPCd54/DMPG at 45°C are shown in Fig. 5, A and B, respectively, and as described in the Materials and Methods section, the spectra correspond essentially to the peptide in interaction with the phospholipid membrane. The amide I' band of the peptide when bound to DMPCd54/DMPA is similar to the amide I' envelope observed for the peptide at low temperature in aqueous solution (see Fig. 4 A), but a difference in conformation can be inferred from the spectral shifts (1618 cm-1 to 1621 cm-1 and 1685 cm-1 to 1694 cm-1). As observed in the derivative and deconvolved spectra shown in Fig. 5 A, there are additional components within the amide I' envelope, which could be assigned to different substructure components. Interestingly, the spectrum does not present a change upon temperature increase (results not shown), indicating a particularly stable entity when bound to the phospholipid. The amide I' band of alpha -MSH when bound to DMPCd54/DMPG is different from that found for the peptide in the presence of DMPCd54/DMPA, being similar to one of the components found previously in the solution structure, that corresponding to the aggregated form of the peptide (see Fig. 4 D). As noted for the other lipidic system, no temperature variation was observed.



View larger version (16K):
[in this window]
[in a new window]
 
FIGURE 5   Infrared spectra in the C==O and amide I' region of alpha -MSH in the presence of DMPCd54/DMPA (A) and DMPCd54/DMPG (B) above the gel-to-liquid-crystalline phase transition of the phospholipid mixture. The derivative and deconvolved spectra are shown above the original spectra (see text for details).

Time-resolved fluorescence spectroscopy

The alpha -MSH tryptophanyl fluorescence decay in buffer is complex (described by three exponentials), and from the best-fits (reduced chi 2 <=  1.2), a short component (tau 1 = 0.448 ns (20°C) or 0.307 ns (37°C); a1 = 0.18), an intermediate component (tau 2 = 2.00 ns (20°C) or 1.36 ns (37°C); a2 = 0.48), and a long component (tau 3 = 3.45 ns (20°C) or 2.35 ns (37°C); a3 = 0.34) are obtained. In the presence of both lipidic systems at the two temperatures studied a variation of both the lifetimes and amplitudes of the components is observed. The data for the DMPC/DMPA (3:1) mixture at 20°C is depicted in Fig. 6, and the trend of variation is similar to the one obtained for the DMPC/DMPG mixture (results not shown). The short lifetime is essentially invariant, while the intermediate one increases slightly and reaches a plateau, and a marked increase is observed for the long one, which varies from 3.45 ns in buffer up two 6.72 ns for the highest lipid concentration. Regarding the amplitudes, that of the short component shows a very slight increase; the amplitude of the intermediate component increases significantly and for the amplitude of the long one a correspondent decrease is obtained, both reaching a plateau at higher lipid concentrations. The amplitudes and the lifetimes of the components are similar to those obtained by Ito et al. (1993). From the partition coefficient study previously described (Eq. 5), < tau > L, the lifetime-weighted quantum yield of the peptide interacting with the membrane, was determined and the values for the two membrane model systems are presented in Table 1. The < tau > L values are longer in the DMPC/DMPG vesicles as compared to the DMPC/DMPA ones at the same temperature.



View larger version (16K):
[in this window]
[in a new window]
 
FIGURE 6   Fluorescence lifetime components tau i and respective normalized pre-exponentials ai for alpha -MSH (lambda exc = 295 nm, lambda em = 350 nm) versus lipid concentration [L] (DMPC/DMPA (3:1) LUV) at 20°C. The solid lines are merely guides to the eye and have no physical meaning.

The time-resolved anisotropies of alpha -MSH in water and at the highest lipid concentration were obtained and the decay for the DMPC/DMPA (3:1) mixture (37°C) is presented in Fig. 7. When in water, the Trp residue fluorescence anisotropy decays to zero, as expected, and values for the rotational correlation time Phi  = 0.520 ± 0.016 ns (20°C) and Phi  = 0.317± 0.009 ns (T = 37°C, result not shown) are obtained, considering for the fundamental anisotropy r0 = 0.24 (lambda exc = 295 nm, Valeur and Weber, 1977). When in interaction with the membrane, the anisotropy of alpha -MSH decays to a value larger than zero (rinfinity ). Considering that even at the highest lipid concentration a considerable fraction of peptide still remains in water, the limiting anisotropy should be corrected for this effect as previously described for the fluorescence spectra. This is carried out according to the equation:
r<SUP><UP>L</UP></SUP><SUB>∞</SUB>=<FENCE>1+<FR><NU>x<SUB><UP>W</UP></SUB>⟨&tgr;⟩<SUB><UP>W</UP></SUB></NU><DE>x<SUB><UP>L</UP></SUB>⟨&tgr;⟩<SUB><UP>L</UP></SUB></DE></FR></FENCE> · r<SUB>∞</SUB> (8)
where r<UP><SUB>∞</SUB><SUP>L</SUP></UP> is the limiting anisotropy of alpha -MSH when in the membrane and rinfinity is the experimentally determined value for the highest lipid concentration. The other parameters are described in Table 1 or obtained from Eq. 6. For both lipidic systems r<UP><SUB>∞</SUB><SUP>L</SUP></UP> is larger at the lower temperature, but no significant variation was observed between the two lipidic systems.



View larger version (26K):
[in this window]
[in a new window]
 
FIGURE 7   Anisotropy decays of alpha -MSH (lambda exc = 295 nm) in buffer (A) and in the presence of 5.2 mM of DMPC/DMPA (3:1) LUV at 37°C (B).

Thermotropic behavior of the lipids from the CH2(CD2) stretching vibration

In order to obtain further information on the transitions observed in the thermograms (Fig. 1), we have studied the effect of alpha -MSH on the phase transition of the phospholipid mixtures by IR (Fig. 8). Because in this study one of the phospholipids was acyl chain-perdeuterated (DMPCd54), and the other was not, it is possible to independently detect changes that might occur with each one.



View larger version (21K):
[in this window]
[in a new window]
 
FIGURE 8   Temperature dependence of the frequency of the symmetrical CH2 stretching (A and C) and symmetrical CD2 stretching (B and D) of mixtures containing DMPCd54/DMPA (A and B) and DMPCd54/DMPG (C and D) in the absence (black-square) and in the presence (down-triangle) of alpha -MSH.

The temperature dependence of the CH2/CD2 symmetric stretching frequency is shown in Fig. 8 A for DMPA and in Fig. 8 B for DMPCd54 in the DMPCd54/DMPA mixture. There, a single broad transition starting at ~22°C and ending at ~35°C is observed (mean values of Tm = 29.0°C (DMPA) and Tm = 28.8°C (DMPCd54)), quite similar to the gel-to-liquid-crystalline phase transition found by DSC for the same lipidic system (see Fig. 1 A). When alpha -MSH was present, a small shift of the mean Tm to higher temperatures (Tm = 30.2°C (DMPA) and Tm = 30.4°C(DMPCd54)) was observed, again in agreement with the DSC data. It should be emphasized that the presence of alpha -MSH induced a shift of the DMPA frequency to higher values, but only at temperatures above Tm, i.e., when the phospholipid mixture is in the liquid-crystalline phase (Fig. 8 A). This effect, not so apparent in DMPCd54 (Fig. 8 B), might indicate the presence of lipid aggregates above Tm, enriched in DMPA bound to alpha -MSH.

The temperature dependence of the CH2 symmetric frequency of the phospholipid mixture DMPCd54/DMPG is shown in Fig. 8 C for DMPG and in Fig. 8 D for DMPCd54. Single cooperative transitions at ~Tm = 21.8°C (DMPCd54) and Tm = 21.5°C (DMPG) were observed, corresponding to the gel-to-liquid-crystalline phase transition of the phospholipid mixture (see the DSC data presented in Fig. 1 B). When alpha -MSH was present in the mixture, similar transitions were observed, but they showed small shifts to lower temperatures (Tm = 20.7°C (DMPCd54) and Tm = 20.3°C (DMPG)), again in agreement with the DSC data.


    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
CONCLUDING REMARKS
REFERENCES

Thermotropic behavior of the phospholipids in the presence of alpha -MSH by calorimetry and infrared spectroscopy

As shown in Fig. 1, the DMPCd54/DMPA system displayed a non-ideal behavior because at least two peaks could be observed in the DSC thermogram. Upon incorporation of alpha -MSH, a significant change was observed. This effect was related to that observed by IR (Fig. 8), as alpha -MSH induced a significant change in the CH2 stretching frequency of DMPA as compared to the CD2 stretching frequency of DMPC, suggesting a preferential interaction of the peptide with DMPA-enriched domains. The fact that no significant effect is apparent in the gel phase is probably due to the peptide incorporation in the line defects between the gel domains (Loura et al., 1996), and so there is not such a strong lipid perturbation in that case. It is interesting that a smooth transition is observed in all cases, and there is no evidence of, e.g., shoulders that would correspond to the peaks observed in DSC. This would mean that the different lipidic domains would have no significant variations in the stretching frequencies depending on their composition, or otherwise they would superimpose preventing the observation of inflections.

In contrast to the DMPCd54/DMPA system, the DMPCd54/DMPG behaves as an ideal system, where only one peak could be discerned. Upon alpha -MSH binding to this system, a broadening of the transition peak was observed, as was a slight shift of Tc to lower temperatures, and by IR we found no significant differences in the CH2 and CD2 stretching frequencies.

In conclusion, alpha -MSH establishes an interaction with the phospholipid molecules, probably via the intercalation of the peptide molecule between the phospholipid headgroups and perturbing, although slightly, the cooperative behavior of the phospholipid acyl chains. As it was noted before, alpha -MSH has a net +1 charge at physiological pH and a low hydrophobicity, implying an essentially electrostatic interaction between the peptide and the anionic phospholipids. This interaction is more effective with the DMPA lipid containing system as compared to DMPG.

Lipid/water partition coefficients of the peptide

Several studies indicate that the peptide does not interact with neutral membranes composed of only zwitterionic phospholipids (Ito et al., 1993; de Kroon et al., 1991; Biaggi et al., 1996), i.e., the hydrophobic contribution for the partition coefficient is not significant. In fact, all the partition coefficients found in this work are of the same order of magnitude because in all cases the global charge of the lipid vesicles is the same (25% anionic phospholipid molecules). However, the different values obtained for the partition coefficients in each case reveal that other characteristics of the system are important besides its global charge. In the case of DMPC/DMPG the result is not surprising, because previous studies point out to a stronger interaction of the peptide with phospholipid vesicles in the liquid-crystalline rather than in the gel phase (Ito et al., 1993), although no quantified values (i.e., Kp or binding constants) were presented for the gel phase. The opposite behavior observed with DMPC/DMPA is unexpected. Although DMPC forms a homogeneous mixture with DMPG, DMPA presents a clear evidence of domain formation (phase separation) from the DSC thermograms (Fig. 1). It is possible that the domains enriched in the anionic phospholipid induce a stronger electrostatic interaction with the peptide, and as described before the interfaces between these domains would facilitate the peptide incorporation.

Preliminary results of both steady-state and time-resolved fluorescence on the interaction of alpha -MSH with a mixture of DMPC/1,2-dimyristoyl-sn-glycero-3-phosphoserine (for which a phase diagram similar to the DMPC/DMPA mixture is reported (Silvius and Gagné, 1984; Graham et al., 1985)) also indicate a stronger interaction in the gel phase rather than in the fluid (unpublished results).

Steady-state fluorescence spectra and anisotropy

The maximum fluorescence emission for alpha -MSH in buffer is at 349 nm, a value very close to the maximum determined for Trp in buffer (350 nm). This result implies that the Trp residue in alpha -MSH is strongly solvated, as expected for a small peptide with no tertiary structure. The interaction of the peptide with the lipidic system is clearly shown from its blue-shifts (3-6 nm, depending on temperature and lipidic system, see Table 1), pointing out to a more hydrophobic environment of the fluorophore (Lakowicz, 1999).

The results for the system with DMPG are in agreement with the studies of Ito et al. (1993), in which the shift was greater when the peptide was interacting with membranes in the liquid crystalline phase. However, when the anionic lipid is DMPA this is not the case. The information we can derive from this observation is that the hydrophobicity of the environment surrounding the Trp residue of alpha -MSH is greater in DMPC/DMPA vesicles in the gel phase than in DMPC/DMPA and DMPC/DMPG vesicles in the fluid phase. A more efficient fluorophore solvation happens for DMPC/DMPG in the gel phase. Once the hydration of the lipid bilayers is greater in the fluid phase (Cevc and Marsh, 1987), for the DMPC/DMPG system we can conclude from the relative spectral shifts that alpha -MSH penetrates deeper into the membrane in the liquid-crystalline phase. Regarding the DMPC/DMPA model membrane, and as the shift is only 1 nm or less for the fluid phase, no information about its location can be obtained. However, IR data are indicative of a deeper penetration in the fluid phase, as concluded from the greater disordering of DMPA acyl chains.

Information about the eventual aggregation of peptide in water, at the low concentrations used in fluorescence experiments, can be obtained from the anisotropy data.

If a spherical geometry is assumed for the peptide, the steady-state anisotropy can be determined from the Perrin equation (Lakowicz, 1999),
⟨r⟩=r<SUB>0</SUB>/(1+&tgr;/&PHgr;) (9)
where r0 is the fundamental anisotropy, tau  the fluorescence lifetime, and Phi  = eta V/RT is the rotational correlation time, eta  and V being the solvent viscosity and the rotating unity volume, respectively.

The peptide volume can be estimated as V = 1821 Å3 (Zamyatnin, 1972), and for water (eta  = 1.00 cP (20°C); eta  = 0.69 cP (37°C)), values of phi  = 0.45 ns (20°C) and phi  = 0.29 ns (37°C) are obtained.

From the average fluorescence lifetime (Eq. 1) of the peptide (<A><AC>&tgr;</AC><AC>&cjs1171;</AC></A> = 2.71 ns (20°C) and <A><AC>&tgr;</AC><AC>&cjs1171;</AC></A> = 1.85 ns (37°C)) and the fundamental anisotropy r0 = 0.105 at lambda exc = 290 nm (Valeur and Weber, 1977), values of < r>  = 0.015 (20°C) and < r>  = 0.014 (37°C) are obtained. These are very close to the experimental ones, and in this way we can conclude that the peptide is essentially monomeric at the physiological concentrations used in the fluorescence measurements (micromolar). In fact, for a dimer (2 V), or a larger aggregate (e.g., 5 V), values of < r>  = 0.026 (20°C) and < r>  = 0.047 (20°C), respectively, would be theoretically expected. The problems arising at the concentrations used in the other techniques (DSC, IR) will be discussed later.

Peptide structure by infrared spectroscopy

We found that the amide I band of alpha -MSH in buffer is dependent on temperature: a mixture of two structures, random coil and hydrogen-bonded beta -strand at low temperature, and nearly all random coil at high temperatures. At variance, this alteration is not observed in the presence of the lipid vesicles under all conditions tested. This behavior is in contrast with other systems, such as gramicidin S in the presence of anionic phospholipids (Lewis et al., 1999). The conformation of alpha -MSH does not change with the physical state of the host lipid bilayer, i.e., when the phospholipid changes from the gel to the liquid-crystalline phase. In this way alpha -MSH would be largely excluded from the lipid bilayer and associated with the lipid polar headgroups at the bilayer surface. Given the relatively small size of alpha -MSH, such an observation could be attributed to the formation of some type of extended intermolecular aggregates at the phospholipid surface (Thiaudière et al., 1991). This result would be in consonance with the lipid/water partition coefficients and the small spectral shifts obtained by fluorescence spectroscopy (Table 1).

It should also be noted that there are differences in the amide I band of alpha -MSH depending on the lipidic system because, as shown in Fig. 5, the proportion of random coil is greater in the presence of DMPC/DMPA than in the presence of DMPC/DMPG. These differences in spectroscopic behavior are also correlated with differences in calorimetric behavior and fluorescence lifetime-weighted quantum yields (see below), and all point to specific interactions with the lipid headgroup. In case that a strong interaction with the lipid acyl chains would exist, significant differences would be apparent for the peptide interaction with gel and fluid phases.

Due to its intrinsic lower sensitivity, it is well known that the concentrations used in the DSC and IR experiments are in general much higher than the physiological ones (submicromolar or lower). In the present work, we showed from the anisotropy data that in all studies involving fluorescence (e.g., partition constants), the monomeric species is the one under study. At variance, both from CD and IR, evidence for aggregate formation was obtained at lower temperatures, and in this way the physiological relevance of these data deserves detailed attention.

It should be stressed that at the highest temperature (70°C), the peptide structure is essentially a random coil in water, i.e., similar to the physiological monomeric species, so the different structures observed in the presence of the two lipids can be safely concluded at this temperature. In addition it was observed that the peptide structure in the presence of lipids is temperature-invariant, at variance with the studies in water. This means that at lower temperatures the interactions with the lipid and not peptide-peptide interactions are the controlling factor. Interestingly, this effect is greater for DMPC/DMPA as compared to DMPC/DMPG, in agreement with the information derived from the fluorescence data.

Time-resolved fluorescence spectroscopy

The lifetime-weighted quantum yield, < tau > , of alpha -MSH in water at 20°C is very close to the value for Trp in aqueous solution at the same temperature (< tau >  = 2.25 ns) (Lakowicz, 1999) and just slightly higher than the value for alpha -MSH reported in the literature (< tau >  = 2.09 ns) (Ito et al., 1993). If we extrapolate the values at different temperatures reported in that work, the fluorescence lifetime at 37°C (< tau >  = 1.51 ns) is in very good agreement with our data. The complex decay obtained for peptides in aqueous solution is presently assigned to different rotational conformers of the indole ring (Willis and Szabo, 1992). Upon interaction with the lipidic system the peptide undergoes alterations of its secondary structure, and these can be appreciated from the trend of variation of both lifetime components and pre-exponentials (Fig. 6).

It is difficult to infer conclusions about the secondary structure of the peptide by the fluorescence decays in the different environments mainly because there are still very few studies made on peptides in which the Trp residue is in a position of known structure and, e.g., Dahms and Szabo (1995) point to subtleties that require caution when predicting structures. However, from the trend of variation of the time-resolved data, it can be concluded that upon membrane incorporation, alpha -MSH undergoes strong structural changes, which is in agreement with the IR data. However, it should be stressed that this pattern of fluorescence time-resolved data only reports the peptide structure in the vicinity of the Trp residue, whereas global information about all the peptide structure is obtained from IR. Thus, there is no contradiction when the very same pattern of variation was obtained for both lipids, at variance with IR data, which allowed concluding that intermolecular aggregates are the dominant structural feature in DMPG-containing vesicles, while in DMPA ones some random structure is present.

The lifetime-weighted quantum yields of the peptide interacting with the different membranes, < tau > L, are longer for DMPC/DMPG as compared to the DMPC/DMPA mixture (Table 1), probably due to the structural differences found by IR spectroscopy. The fact that < tau > L is longer in DMPG-containing vesicles but the spectral shift is smaller does not represent any contradiction. In fact, there is no clear correlation between quantum yield and the wavelength of maximum emission in proteins with a single Trp residue (Lakowicz, 1999).

The values obtained for the rotational correlation times Phi , assuming a spherical rotor, are in reasonable agreement with the ones expected from the Perrin equation (Eq. 9), once that an average lifetime is used to describe the complex decay.

Our study of time-resolved fluorescence anisotropy in membranes was restricted to the limiting anisotropies rinfinity because the correction for the fraction of peptide in the aqueous phase would be too critical to allow the determination of the dynamics of the system contained at the earlier times of the decay. The limiting anisotropies obtained are in all cases very high, considering that the fluorophore is not a typical hydrophobic molecule incorporated into the membrane core with a serious restriction to its motion. This means that in all cases the peptide is strongly adsorbed at the membrane/water interface and the wobbling motion of the Trp is very restricted, even in the fluid phase. Considering that most of the studies (e.g., Kp, Delta lambda , DSC) carried out for the DMPC/DMPA mixture in the gel phase gave peculiar results due to the phase separation observed by DSC, we could also anticipate a greater value for rinfinity , which is not observed. It can then be concluded that for DMPC/DMPG less peptide incorporation is obtained (lower KP), but the peptide in the membrane has a hindered rotation similar to the one observed for the DMPC/DMPA mixture. Irrespective of the structure adopted by the peptide in the vesicles, not only it is thermally stabilized, it is also very rigid. The limiting anisotropies rinfinity are related to the order parameter S, and the fundamental anisotropy r0 through the following relationship (e.g., Lakowicz (1999)):
S<SUP>2</SUP>=<FR><NU>r<SUB>∞</SUB></NU><DE>r<SUB>0</SUB></DE></FR>=<FENCE><FR><NU>3 <UP>cos</UP><SUP>2</SUP>&thgr;−1</NU><DE>2</DE></FR></FENCE><SUP>2</SUP> (10)
where theta  is the angular rotation of the emission transition moment and the angle brackets indicate an average over all the fluorophore population. If we consider a fundamental anisotropy for Trp excitation at 295 nm of r0 = 0.3 (Valeur and Weber, 1977), an angle theta  = 29° is obtained for the gel phase, or ~30-33° for the liquid-crystalline phase. The Trp containing central 6-9 region of alpha -MSH is the minimal melanotropic message sequence (Hruby et al., 1987). In a recent molecular dynamics study (Pascutti et al., 1999) the peptide acquired, in a low dielectric constant medium, a packed conformation that remained stabilized for >7.0 ns of the simulation. A common feature of various superpotent synthetic analogs of the native hormone is the stabilization of a turn involving the Trp residue (Sawyer et al., 1980, 1982). The behavior of the fluorescence decay components indicates a similar structure in the region of the Trp residue region, and the limiting anisotropies reveal an effective immobilization of the indole ring in both lipidic systems, above and below the Tm.


    CONCLUDING REMARKS
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
CONCLUDING REMARKS
REFERENCES

The interaction of alpha -MSH with (3:1) lipid mixtures containing the same net charge (DMPC/DMPG or DMPC/DMPA) was studied by complementary methods.

It was concluded that alpha -MSH in buffer at higher concentrations undergoes a transition from an aggregated plus random coil structures to essentially random coil upon temperature increase, this transition being prevented when in interaction with negatively charged vesicles. These can stabilize specific structures different from the ones in solution, as concluded from IR and fluorescence data. A strong interaction of the peptide with both lipidic systems was also concluded from the high limiting anisotropies observed both above and below Tm.

The peptide can induce lateral heterogeneity in the DMPC/DMPA lipidic systems, as seen from the strong alterations observed in the calorimetric thermograms, and in addition a greater partition coefficient is obtained for this last system in the gel phase as compared to the fluid phase. Analogs of alpha -MSH with superpotent biological activity possess a defined and stable structure. It is reasonable to predict that the interaction of alpha -MSH with negatively charged membranes induces and stabilizes a specific conformation, probably involving the Trp-containing message region, which could be similar to the one found for the superpotent analogs, and therefore necessary for its biological activity.

    ACKNOWLEDGMENTS

This work was supported by grants PM98-0100 from DGESIC (Spain), and Portuguese-Spanish Joint Research Program, Project PRAXIS/C/SAU/14025/98 from FCT (Portugal). A.F. acknowledges an INVOTAN grant (Portugal), and L.M.C. and R.F.M. de A. pre-doctoral fellowships from CONICIT (Venezuela) and FCT (Portugal).

A referee is acknowledged for drawing our attention to the problems associated with the high concentrations used in DSC and IR techniques.

J.V. and M.P. contributed equally to this work.

    FOOTNOTES

Received for publication 17 April 2000 and in final form 12 February 2001.

Address reprint requests to Dr. Manuel Prieto, Centro de Química-Física Molecular, Instituto Superior Técnico, P-1049-001 Lisboa, Portugal. Tel.: +35-121-841-9219; Fax: +35-121-846-4455; E-mail: prieto{at}alfa.ist.utl.pt.


    REFERENCES
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
CONCLUDING REMARKS
REFERENCES