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 Moncelli, M. R.
Right arrow Articles by Guidelli, R.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Moncelli, M. R.
Right arrow Articles by Guidelli, R.

Biophys J, May 1998, p. 2388-2397, Vol. 74, No. 5

Surface Dipole Potential at the Interface between Water and Self-Assembled Monolayers of Phosphatidylserine and Phosphatidic Acid

Maria Rosa Moncelli, Lucia Becucci, Francesco Tadini Buoninsegni, and Rolando Guidelli

Department of Chemistry, Florence University, 50121 Florence, Italy

    ABSTRACT
Top
Abstract
Introduction
Results
Discussion
Conclusions
Appendix
References

The nature and magnitude of the surface dipole potential chi  at a membrane/water interface still remain open to discussion. By combining measurements of differential capacity C and charge density sigma  at the interface between self-assembled monolayers of phosphatidylserine and phosphatidic acid supported by mercury and aqueous electrolytes of different concentration and pH, a sigmoidal dependence of chi  upon sigma  is revealed, with the inflection at sigma  = 0. This behavior is strongly reminiscent of the surface dipole potential due to reorientation of adsorbed water molecules at electrified interfaces. The small increase in C with a decrease in the frequency of the AC signal below ~80 Hz, as observed with phospholipid monolayers with partially protonated polar groups, is explained either by a sluggish collective reorientation of some polar groups of the lipid or by a sluggish movement of protons across the polar head region.

    INTRODUCTION
Top
Abstract
Introduction
Results
Discussion
Conclusions
Appendix
References

Even though the existence of an appreciable dipole potential difference between the interior of a membrane and the adjacent aqueous solution is universally accepted, the origin of this dipole potential remains obscure. Thus it may stem from the orientation of dipoles in 1) the water molecules adjacent to the membrane, 2) the polar headgroups, and/or 3) the ester linkages to the glycerol backbone (McLaughlin, 1977). The dipole potential chi  is not a thermodynamically significant quantity, and hence cannot be measured directly. In principle, however, an insight into the origin of chi  at membrane/water interfaces can be gained by measuring its changes with a change in the charge density sigma lip on the membrane surface. Incorporation of lipophilic ions into a membrane to change sigma lip is to be avoided, because these ions may easily alter the dipole potential by their very presence. A convenient procedure for altering the charge density of a phospholipid monolayer consists of changing the extent of protonation of its ionizable groups by varying the pH.

In recent years we have measured the charge density sigma lip of self-assembled monolayers of different phospholipids as a function of pH (Moncelli et al., 1994; Moncelli and Becucci, 1995), by using a biomimetic membrane that consists of a mercury electrode coated with a phospholipid monolayer (Miller, 1981; Nelson and Benton, 1986). This half-membrane provides an inherent mechanical stability and a resistance to high electric fields that are not shared by BLMs. Over the potential region of minimum capacity, which ranges from -0.15 to -0.75 V/SCE, the monolayer is impermeable to inorganic metal ions, whereas it becomes permeable outside this region. The differential capacity C of a lipid monolayer on mercury over this region is ~1.7-1.9 µF cm-2, that is, twice the value for a BLM. The charge density sigma lip of a self-assembled monolayer of phosphatidylserine (PS) supported by mercury was found to vary from slightly negative to slightly positive values as the bulk pH of the bathing solution is varied from 7 to 4 (Moncelli et al., 1994). Analogously, the charge density sigma lip of a monolayer of phosphatidic acid (PA) passes from negative to positive values as the pH is varied from 4 to 1.5 (Moncelli and Becucci, 1995). PS and PA are therefore ideal candidates for measuring chi  changes with varying pH.

In Moncelli at al. (1994) and Moncelli and Becucci (1995), the charge density sigma lip as a function of pH was determined by measuring the small changes in the overall differential capacity C of a lipid-coated mercury electrode after a change in the concentration of the electrolyte KCl from 5 × 10-3 to 0.1 M. These changes were considered to be due exclusively to a change in the differential capacity Cd of the diffuse layer, which can be regarded as being in series with the very low capacity, Clip, of the monolayer; the reciprocal of the experimental capacity was therefore set equal to 1/C = 1/Clip + 1/Cd. The change in Cd after a given change in the electrolyte concentration c is expected to decrease rapidly with an increase in the absolute value of the overall charge sigma  experienced by the ions of the diffuse layer. Thus, if we plot values of the reciprocal, 1/CdGC, of the diffuse-layer capacity calculated on the basis of the Gouy-Chapman (GC) theory at different electrolyte concentrations c and at constant charge sigma  against the corresponding values calculated at sigma  = 0, 1/Cd,0GC, we obtain roughly straight lines whose slope decreases progressively with an increase in |sigma |, and ultimately vanishes for |sigma >=  4 µC cm-2. Because the capacity Clip of the lipid monolayer is approximately independent of the electrolyte concentration, the slope, Sexp, of an experimental plot of 1/C against 1/Cd,0GC for a given range of electrolyte concentrations was regarded as a measure of the slope of the plot of the reciprocal, 1/Cd, of the experimental diffuse-layer capacity against 1/Cd,0GC. Slopes, Scalc, of plots of the reciprocal 1/CdGC of the diffuse-layer capacity against 1/Cd,0GC were therefore calculated on the basis of the GC theory at different pH values, for different sets of values of the protonation constants of the ionizable groups of the lipid. The resulting plots of Scalc versus pH were then compared with the experimental plot of Sexp versus pH for the lipid under study. Finally, the protonation constants of the lipid were ascribed the values providing the best fit between Sexp versus pH and Scalc versus pH plots. The overall charge density sigma  experienced by the diffuse-layer ions is the sum of the charge density on the lipid, sigma lip, plus the small charge density on the mercury surface, sigma M. Hence, to calculate Scalc, the charge density sigma M as a function of c and pH had also to be estimated. In doing so, we assumed that the surface dipole potential chi  was independent of the solution composition, for simplicity.

In this work the simplifying assumption of a pH-independent surface dipole potential chi  is abandoned. To draw conclusions about the pH dependence of chi , the charge density sigma M on PS- and PA-coated mercury electrodes was measured at different pH values. Moreover, the dependence of the differential capacity C upon the frequency nu  was checked over the frequency range from 2 to 500 Hz. The results of these measurements suggest a contribution to chi  from the reorientation of adsorbed water molecules; tentative explanations for the frequency dependence of C at frequencies less than 80 Hz will be provided.

    EXPERIMENTAL

The water used was obtained from light mineral water by distilling it once, and by then distilling the water so obtained from alkaline permanganate, while constantly discarding the heads. Merck Suprapur® KCl was baked at 500°C before use to remove any organic impurities. Dioleoylphosphatidylcholine (PC) and dioleoyl PS were obtained from Lipid Products (South Nutfield, Surrey, England), and dioleoyl PA was obtained from Avanti Polar Lipids (Birmingham, AL). The desired pH values were realized with Merck Suprapur® HCl over the pH range from 2 to 5, with a 1 × 10-3 M HPO42-/H2PO4- buffer over the pH range from 6.5 to 7.5, and with a 1 × 10-3 M H3BO3/NaOH buffer over the pH range from 8.5 to 9.8.

The home-made hanging mercury drop electrode (HMDE), the cell, and the procedure for the preparation of the self-assembled phospholipid monolayers are described elsewhere (Moncelli et al., 1994). Measurements of the differential capacity C at a constant frequency of 75 Hz were carried out with a Metrohm Polarecord E506 (Herisau, Switzerland). In view of the low capacity of the lipid-coated electrode (<2 µF cm-2), C was directly measured by the quadrature component of the AC current, other than at the lowest salt concentrations; in the latter case, both quadrature and in-phase components of the AC current were measured, to correct for the cell resistance. The system was calibrated using a precision capacitor (Decade Capacitor type 1412-BC; General Radio, Concord, MA). All potentials were measured versus a saturated calomel electrode (SCE). The reproducibility of the differential capacity in passing from one mercury drop to another was better than 0.05 µF cm-2. At any rate, each set of differential capacity measurements at variable KCl concentration and constant pH was carried out on the same lipid-coated mercury drop, so as to practically eliminate the effect of any slight irreproducibilities in the drop surface area or in the lipid transfer. This permitted us to estimate the changes in differential capacity after an increase in electrolyte concentration with an accuracy better than 0.02 µF cm-2. The electrolyte concentration in the cell was progressively increased by adding a deaerated solution of the concentrated electrolyte from a microsyringe (Hamilton, Reno, NV). The plunger of the syringe was fastened tightly to the rod of a digital display micrometer screw with a 0.005-mm pitch (no. 297-101-01; Mitutoyo, Tokyo, Japan). The micrometer screw was held by a movable stand that permitted the syringe needle to be lowered into the solution during the addition and raised above the solution just after the addition. After each addition the solution was stirred mildly for ~30 s with a magnetic stirrer on the bottom of the cell. The stability of the differential capacity was tested by recording it over the whole potential region of minimum capacity two or three times consecutively, interposing a mild stirring between each measurement; whenever detectable differences between these recordings were observed, the whole series of measurements was discarded. Differential capacity measurements at different frequencies were carried out with a Stanford Research 850 lock-in amplifier. To check the stability of the lipid monolayer during measurements, the frequency of the AC signal was first varied progressively from 2 to 500 Hz, and then in the opposite direction on the same lipid-coated mercury drop. Measurements were discarded whenever the difference in the capacity values at the same frequency in the two opposite runs was found to be greater than 1%.

The surface charge density sigma M at the HMDE coated with a self-assembled phospholipid monolayer was measured by a technique described elsewhere (Becucci et al., 1996). Briefly, sigma M was obtained by analogical integration of the capacitive current that flows at constant applied potential as a consequence of a slight contraction of the mercury drop. The contraction must be carried out while keeping the neck of the lipid-coated mercury drop in contact with the lipid film spread on the surface of the electrolytic solution. This procedure ensures that the monolayer maintains its properties, including its thickness, as the drop is expanded or compressed. The capacitive charge flowing during a change Delta A in the drop area, once divided by Delta A, yields directly the charge density sigma M on the metal.

    RESULTS
Top
Abstract
Introduction
Results
Discussion
Conclusions
Appendix
References

Protonation

Fig. 1 shows plots of Sexp versus pH for PS- and PA-coated mercury electrodes. Both plots show a maximum. In the case of PS, Sexp is practically zero at pH 7.5, which implies that at this pH value the polar head is negatively charged. The Sexp value approaching unity in the proximity of pH 6 indicates that the polar head of PS is practically neutral at this pH. This implies that one of the two anionic groups, either the phosphate or the carboxyl group, is almost completely protonated; in this way, the negative charge borne by the other anionic group is practically neutralized by the positive charge of the amino group, which is completely protonated over the whole pH range investigated. As the pH is decreased further from 6 to 3, the Sexp value decreases again, attaining the zero value: this implies that the polar head is now positively charged, and hence that even the further anionic group starts to be appreciably protonated. To justify the very weak acidity of at least one of the two anionic groups, we must necessarily assume that it is buried somewhere inside the polar head region of PS (Moncelli at al., 1994). In such a position an anionic group has a much lower acidity than in bulk water, because the negative average potential difference between the position of the anion and the aqueous solution attracts protons electrostatically; moreover, the reaction of the monoanion with a proton annihilates the charges of both reactants, and hence is strongly favored by the low dielectric constant of the polar head region. Of the two anionic groups, the one that is likely to be more deeply buried in the polar head region is the phosphate group, because of its closer vicinity to the hydrocarbon tails. As concerns the PA film, the -PO42- group is monoprotonated over the whole pH range investigated. The fact that Sexp becomes different from zero at pH less than 5 denotes a further protonation of the phosphate group and a resulting tendency of the polar heads to become uncharged. However, the rapid decay of Sexp in passing from pH 2 to pH 1.5 can only be explained by a tendency of the polar heads to become positively charged at these low pH values. A possible explanation is a labile binding of a proton to the >Cdouble bond O group of one of the ester groups of the lipid (Moncelli and Becucci, 1995).


View larger version (14K):
[in this window]
[in a new window]
 
FIGURE 1   Values of Sexp versus pH for PS (black-square) and PA (triangle ) self-assembled monolayers supported by mercury. The solid curves are Scalc versus pH plots calculated as described in the text for PS, with K1 = 5 × 106 M-1, K2 = 1 × 105 M-1, and gamma /epsilon gamma  = 0.1 × 10-8 cm (a), and for PA with K1 = 1 × 108 M-1, K2 = 1 × 105 M-1, K3 = 50 M-1, and gamma /epsilon gamma  = 0.2 × 10-8 cm (b).

Charge measurements

Fig. 2 shows the charge density sigma M on PS- and PA-coated mercury in 0.1 M KCl at a constant applied potential of -0.5 V as a function of pH. In the case of the PS film, sigma M attains a maximum value in the proximity of pH 5, where the PS polar head is almost uncharged, and then decreases again with a further increase in pH. A decrease in sigma M with an increase in pH is also observed with the PA film. This behavior contradicts our original assumption of a pH-independent dipole potential chi  (Moncelli et al., 1994; Moncelli and Becucci, 1995). If this assumption were correct, then the charge density sigma M at constant applied potential would shift gradually in the positive direction with an increase in pH, to compensate for the negative shift in the surface potential psi d after the progressive deprotonation of the ionizable groups. The behavior of the sigma M versus pH plots in Fig. 2 denotes an appreciable change in chi  with a change in pH.


View larger version (13K):
[in this window]
[in a new window]
 
FIGURE 2   Values of sigma M versus pH for PS (black-square) and PA (triangle ) self-assembled monolayers supported by mercury in 0.1 M KCl at -0.5 V/SCE.

Frequency dispersion

The frequency dependence of the differential capacity C of self-assembled monolayers of PS, PA, and PC in contact with aqueous solutions of 0.1 M KCl of different pH values is shown in Figs. 3 and 4 over the frequency range from 2 to 500 Hz at a bias potential of -0.5 V. With all systems investigated, C is independent of the frequency nu  for nu  >=  80 Hz. However, as the frequency is decreased below this value, the PS and PA films start to show a small but progressive increase in C that is still observed at 2 Hz. This frequency dispersion is observed at all pH values investigated, but is more pronounced at the lower pH values. The behavior of PC differs from that of PS and PA in that no frequency dispersion is observed at pH greater than 4; only at lower pH values does the differential capacity increase slightly with a decrease in frequency below 80 Hz.


View larger version (13K):
[in this window]
[in a new window]
 
FIGURE 3   Plots of C versus nu  for a PS self-assembled monolayer at -0.5 V/SCE in 0.1 M KCl-buffered solutions of pH 4.3 (black-down-triangle ), 3.2 (square ), 5.3 (black-triangle), 7.3 (open circle ), and 6.3 (bullet ).


View larger version (13K):
[in this window]
[in a new window]
 
FIGURE 4   Plots of C versus nu  for a PC self-assembled monolayer at -0.5 V/SCE in 0.1 M KCl-buffered solutions of pH 1.8 (bullet ) and 7.5 (black-down-triangle ), and for a PA self-assembled monolayer at -0.5 V/SCE in 0.1 M KCl-buffered solutions of pH 1.5 (black-triangle) and 7.2 (square ).

    DISCUSSION
Top
Abstract
Introduction
Results
Discussion
Conclusions
Appendix
References

Frequency dependence of the differential capacity

In view of the complexity of the structure of the self-assembled lipid monolayer, the frequency dispersion in Figs. 3 and 4 cannot be ascribed unambiguously to a single phenomenon. Thus, in principle, any movement of charged species or reorientation of dipoles within the lipid film that is in a condition of following the small (10 mV peak-to-peak) AC signal will oppose the corresponding external alternating field with a resulting increase in differential capacity; only when a sufficiently high increase in frequency causes this charge movement and/or dipole reorientation to lag behind the AC signal will the differential capacity attain a constant minimum value. A sluggish change in the tilt of the hydrocarbon tails of the lipid after the AC signal, with a resulting change in the thickness of the film and in its differential capacity, seems to be excluded. In fact, it would be expected not only with PS and PA, but also with PC, because all of these lipids have the same dioleoyl hydrocarbon tails. Moreover, it cannot explain the passage from a frequency independence of C to a frequency dispersion with a decrease in pH, as observed with PC films (see Fig. 4). The frequency dispersion in Figs. 3 and 4 can be tentatively explained 1) by a sluggish collective reorientation of some polar groups of the lipids after the AC signal, or 2) by a sluggish movement of protons from partially protonated ionizable groups buried inside the polar head region to the bathing solution and vice versa.

In the first case the change in orientation of the polar heads of the lipid after the small AC signal must also be very small. Thus, e.g., if the dipole moment µ of the dipole consisting of the charged carboxyl group and of the charged ammonium group in a PS molecule is estimated at 6 D, the change in chi  involved in its passage from an orientation parallel to the monolayer to a vertical orientation amounts to 4pi NNAµ/epsilon gamma , where N approx  2 × 10-10 mol cm-2 is the density of PS in the monolayer, NA is Avogadro's number, and epsilon gamma is the dielectric constant of the polar head region. If we set epsilon gamma  = 30, the change in chi  is equal to 0.09 V and involves a charge flow of about (2 µF cm-2) × 0.09 V = 0.18 µC cm-2; if the AC signal of 10 mV peak-to-peak were to produce such a drastic reorientation of the above dipole moment, it would give rise to an increase in differential capacity as high as 18 µF cm-2. A slight change in the orientation of the polar heads of a PS monolayer supported on Hg can be tentatively justified by considering that the acidity of the phosphate group in this monolayer is much lower than that normally reported in the literature for PS vesicles (Tsui et al., 1986), dispersions (MacDonald et al., 1976), monolayers (Ohki and Kurland, 1981), and BLMs (Matinyan et al., 1985), where the electric field acting on the hydrocarbon tails and the parameters related to intermolecular spacing and state of compression may be somewhat different. In Moncelli at al. (1994), this difference in behavior was explained by assuming that PS self-assembled monolayers may assume at least two different conformations of the polar head, with similar Gibbs energies but quite different acidities of the ionizable groups. Thus a conformation of the PS polar heads with two negative and one positive charge on the same plane parallel to the lipid layer is not as electrostatically favored as the conformation assumed by zwitterionic lipids such as PC: a conformation with the phosphate group deep inside the polar head region and a C-N dipole roughly parallel to the lipid plane and in direct contact with the aqueous phase may well have a comparable Gibbs energy. This interpretation of the apparently anomalous behavior of PS self-assembled monolayers is supported by the observation that in the presence of adsorbed tetraphenylphosphonium cations, these monolayers behave as though they were actually negatively charged (Moncelli et al., 1995); this behavior was explained by a conformational change in the PS polar heads induced by the tetraphenylammonium cations, leading to a deprotonation of the phosphate groups. It is therefore possible that at frequencies less than 80 Hz, the applied AC field may start to be accompanied by a modest fluctuation in the conformation of the PS polar heads.

An alternative explanation for the increase in differential capacity with a decrease in frequency below ~80 Hz consists of assuming a progressive increase in the ability of the protons to move to and fro across the polar head region after the AC signal, and hence to oppose the applied AC field. This implies a slow equilibration of the protons between the polar head region of the lipid film and the bathing solution. This interpretation contrasts with kinetic analyses of time-resolved proton-phospholipid interactions in micelles and liposomes (Nachliel and Gutman, 1988; for a review see Gutman and Nachliel, 1990), according to which the rate of proton binding to the phospholipid lies in the microsecond and submicrosecond time scale. Gutman's conclusions tend to support the "delocalized chemiosmotic theory" (Kasianowicz et al., 1987; Polle and Junge, 1989), according to which the proton movement from proton pumps to proton sinks in photosynthesis and respiration takes place in the aqueous bulk phase because of a very rapid equilibration of the protons between the lipid and the adjacent bathing solution. However, in several laboratories, evidence has also been gathered in favor of a "localized theory," according to which protons move exclusively along the membrane surface; the latter evidence relies on measurements with both biomimetic membranes (Kell, 1979; Prats et al., 1985, 1986; Teissié et al., 1985; Morgan et al., 1988; Antonenko et al., 1993) and fragments of biomembranes. Thus Heberle et al. (1994) showed that a pH sensor positioned at the surface of a purple membrane, at an average distance of 240 nm from the proton ejecting bacteriorhodopsin, detects the liberated protons eight times faster than a pH probe in the bulk aqueous phase at an average distance of only 17 nm. According to these authors, the proton's lateral motion along the membrane surface is faster than in the adjacent bulk water phase, not because of a higher diffusion coefficient of protons, but rather because of a surprisingly low rate of proton transfer from the membrane surface to the water phase, lying in the millisecond time scale; protons should therefore move within an extended Coulomb cage formed by the lipid headgroups and the proteinous amino acids. An interpretation of the frequency dispersion in Figs. 3 and 4 in terms of a slow equilibration of protons between the polar head region of the lipid monolayer and the bulk aqueous phase would therefore provide a further piece of evidence in favor of the delocalized theory.

Under the assumption that the frequency dispersion is due to sluggish protonation equilibria, the experimental Sexp versus pH plots in Fig. 1 refer to a situation in which these protonation equilibria do not follow the AC signal, because they were obtained at a frequency of 75 Hz, which practically marks the upper boundary of the region of frequency dispersion, as appears from Figs. 3 and 4. The two-capacitor model adopted in the previous work carried out in this laboratory (Moncelli et al., 1994) does not account for this situation, because it locates all ionizable groups in direct contact with the aqueous phase; moreover, the capacity C includes a finite contribution due to the rate of change, dsigma lip/dsigma M, of the overall charge density sigma lip of the polar heads of the lipid with a change in sigma M. If sigma M is shifted in the positive direction, protons are repelled electrostatically from the aqueous phase in the immediate vicinity of the lipid film, increasing the local pH there; this causes an instantaneous partial deprotonation of the ionizable groups of the lipid and a decrease in sigma lip. In practice, a protonation-deprotonation step with a relaxation time much shorter than the period, nu -1, of the AC signal causes the term dsigma lip/dsigma M to be negative. On the other hand, if its relaxation time is much longer than nu -1, the term dsigma lip/dsigma M tends to vanish, causing a decrease in the differential capacity C.

A model of the mercury/phospholipid/ water interphase

In what follows we will adopt a model of three capacitors in series, schematically depicted in Fig. 5, to account for a possible slow equilibration of protons. The model will only consider the two extreme situations in which the protonation equilibria involving the ionizable groups buried well inside the polar head region either follow the AC signal perfectly or else do not follow it at all, and hence are blocked at the bias potential. This model will serve to show that the experimental frequency dispersion can be justified by assuming that certain protonation equilibria are blocked at frequencies greater than 80 Hz, and hence that their contribution to sigma lip does not follow the fluctuations of sigma M produced by the AC signal. Fig. 5 shows a model of a lipid monolayer deposited on mercury, consisting of a hydrocarbon tail region of dielectric constant epsilon beta enclosed between the electrode surface plane x = 0 and the plane x = beta , and of a polar head region of dielectric constant epsilon gamma enclosed between x = beta  and the lipid/solution boundary x = d triple-bond  (beta  + gamma ) (Moncelli et al., 1995). For simplicity, the ionizable groups of the lipid are considered to be located either at the boundary x = beta  between these two regions, or else at x = d, that is, in direct contact with the aqueous phase. The groups at x = d experience a hydrogen ion concentration satisfying the Boltzmann distribution law, i.e., cH+ exp(-Fpsi d/RT), where cH+ is the bulk hydrogen ion concentration and psi d is the "average" electric potential at x = d, i.e., the surface potential. The hydrogen ion concentration at x = beta  is also assumed to satisfy a Boltzmann distribution law, cH+ exp(-Fphi beta /RT); in this case, however, the electric potential phi beta at x = beta  is considered to have a local character, and hence to experience discreteness-of-charge effects. These effects are considered in the framework of the "cutoff disk model," according to which an adsorbed ion is surrounded by a circular charge-free region (the exclusion disk) that is imaged infinite times in the x = 0 and x = d planes (Levine et al., 1962, 1965). The expected values for the parameters of a typical lipid monolayer are 10-20 Å for beta , 4-10 Å for gamma , ~2 for epsilon beta , and 8-50 for epsilon gamma (Flewelling and Hubbell, 1986). Moreover, the cross-sectional area of a lipid molecule is close to 60 Å2; hence, if each polar head contains only one ionizable group at x = beta , the "steric hard-core radius" between two neighboring ionizable groups is on the order of 8-9 Å. The exclusion disk radius, rho , cannot be smaller than this steric hard-core radius. With such a large value for rho , Levine's expression for the local potential phi beta , as measured with respect to the bulk solution, is satisfactorily approximated by its limiting form for rho  right-arrow infinity (Levine et al., 1972; Moncelli et al., 1995):
&phgr;<SUB>&bgr;</SUB>=<FR><NU>&ggr;/&egr;<SUB>&ggr;</SUB></NU><DE>&bgr;/&egr;<SUB>&bgr;</SUB>+&ggr;/&egr;<SUB>&ggr;</SUB></DE></FR> (&psgr;<SUB>0</SUB>−&khgr;−&psgr;<SUB><UP>d</UP></SUB>)+&psgr;<SUB><UP>d</UP></SUB> (1)
In fact, when rho  is comparable to the distance 2gamma between the discrete charges at x = beta  and their nearest-neighboring images, the screening effect of these images becomes so large as to cause the limiting behavior for rho  right-arrow infinity  to be closely approached.


View larger version (15K):
[in this window]
[in a new window]
 
FIGURE 5   Schematic picture of the model for a lipid monolayer deposited on mercury. The dashed curve schematically represents the profile of the average potential against the distance from the mercury surface. The diffuse-layer thickness has been compressed with respect to the monolayer thickness, for ease of representation.

For the sake of generality, let us denote the set of charge densities due to the ionizable groups at x = beta  by {sigma i}, and that due to the ionizable groups at x = d by {sigma j}, where the subscripts i and j refer to the different groups. The average potential difference psi 0 across the whole interface will then be given by
&psgr;<SUB>0</SUB>=4&pgr; <FR><NU>&bgr;</NU><DE>&egr;<SUB>&bgr;</SUB></DE></FR> &sfgr;<SUB><UP>M</UP></SUB>+4&pgr; <FR><NU>&ggr;</NU><DE>&egr;<SUB>&ggr;</SUB></DE></FR> (&sfgr;<SUB><UP>M</UP></SUB>+∑<SUB><UP>i</UP></SUB>&sfgr;<SUB><UP>i</UP></SUB>)+&khgr;+&psgr;<SUB><UP>d</UP></SUB> (2)
Here the first term is the average potential difference across the hydrocarbon tails, whereas the second is that across the polar head region; according to the GC theory, the potential difference psi d across the diffuse layer is a function of c, sigma M, and the overall charge density of the lipid, sigma lip = Sigma isigma i + Sigma jsigma j. If all protonation equilibria are perfectly mobile, differentiation of Eq. 2 with respect to sigma M yields the following expression for the reciprocal of the differential capacity C:
<FR><NU>1</NU><DE>C</DE></FR>=<FR><NU><UP>d</UP>&psgr;<SUB>0</SUB></NU><DE><UP>d</UP>&sfgr;<SUB><UP>M</UP></SUB></DE></FR>=4&pgr;<FENCE><FR><NU>&bgr;</NU><DE>&egr;<SUB>&bgr;</SUB></DE></FR>+<FR><NU>&ggr;</NU><DE>&egr;<SUB>&ggr;</SUB></DE></FR></FENCE>+4&pgr; <FR><NU>&ggr;</NU><DE>&egr;<SUB>&ggr;</SUB></DE></FR> <FR><NU><UP>d</UP>(∑<SUB><UP>i</UP></SUB>&sfgr;<SUB><UP>i</UP></SUB>)</NU><DE><UP>d</UP>&sfgr;<SUB><UP>M</UP></SUB></DE></FR>+<FR><NU><UP>d</UP>&psgr;<SUB><UP>d</UP></SUB></NU><DE><UP>d</UP>(&sfgr;<SUB><UP>M</UP></SUB>+&sfgr;<SUB><UP>lip</UP></SUB>)</DE></FR><FENCE>1+<FR><NU><UP>d</UP>&sfgr;<SUB><UP>lip</UP></SUB></NU><DE><UP>d</UP>&sfgr;<SUB><UP>M</UP></SUB></DE></FR></FENCE> (3)
where dchi /dsigma M is now regarded as negligibly small. At frequencies high enough to block the movement of protons across the polar head region (beta  < x < d), the derivatives dsigma i/dsigma M vanish, and Eq. 3 becomes
<FR><NU>1</NU><DE>C</DE></FR>=4&pgr;<FENCE><FR><NU>&bgr;</NU><DE>&egr;<SUB>&bgr;</SUB></DE></FR>+<FR><NU>&ggr;</NU><DE>&egr;<SUB>&ggr;</SUB></DE></FR></FENCE>+<FR><NU><UP>d</UP>&psgr;<SUB><UP>d</UP></SUB></NU><DE><UP>d</UP>(&sfgr;<SUB><UP>M</UP></SUB>+&sfgr;<SUB><UP>lip</UP></SUB>)</DE></FR> <FENCE>1+<FR><NU><UP>d</UP>(∑<SUB><UP>j</UP></SUB>&sfgr;<SUB><UP>j</UP></SUB>)</NU><DE><UP>d</UP>&sfgr;<SUB><UP>M</UP></SUB></DE></FR></FENCE> (4)
For an uncharged lipid, the differential capacity is approximately given by [4pi (beta /epsilon beta  + gamma /epsilon gamma )]-1, once we neglect the small contribution from the diffuse layer. This quantity can be accurately estimated at 1.7 µF cm-2, which corresponds to the differential capacity of an uncharged PC monolayer. Because all of the features of the lipid monolayer, apart from the protonation constants, depend exclusively upon the beta /epsilon beta and gamma /epsilon gamma ratios, only one of these two parameters is adjustable, whereas the other is obtained from the relation [4pi (beta /epsilon beta  + gamma /epsilon gamma )]-1 approx  1.7 µF cm-2. In particular, if we ascribe to beta  and epsilon beta the reasonable values 10 Å and 2, the gamma /epsilon gamma ratio turns out to be equal to 0.21 × 10-8 cm in electrostatic CGS units. A treatment of the model is outlined in the Appendix.

The solid curve a in Fig. 1 shows the Scalc versus pH plot for PS in best agreement with the corresponding Sexp versus pH plot, as calculated on the basis of the model by assuming that the phosphate and carboxyl groups are located at x = beta , with gamma /epsilon gamma = 0.1 × 10-8 cm, and by setting the protonation constants of these two groups equal to K1 = 5 × 106 M-1 and K2 = 1 × 105 M-1. The plot was calculated by regarding the protonation equilibria of these groups as blocked at the bias potential E = -0.5 V; in other words, the protons of the phosphate and carboxyl groups were considered to be unable to follow the 75-Hz AC signal. Incidentally, over the pH range investigated, the amino group of PS is fully protonated and therefore does not contribute to the movement of protons after the AC signal. Curve a in Fig. 6 shows a plot of Delta C versus pH, where Delta C is the difference between the differential capacity values estimated for the two extreme situations in which the protonation equilibria follow the AC signal perfectly or else are blocked at the bias potential; naturally, when we assumed that all protonation equilibria are perfectly mobile through the use of Eq. 4, the differential capacity was calculated by using a different set of protonation constants, i.e., the set that provides the best agreement with experiment under these assumptions. In practice, Delta C measures the maximum frequency dispersion resulting from the lack of proton equilibration. As expected, the maximum frequency dispersion is attained in the proximity of the pH values corresponding to the log K values of the ionizable groups buried in the polar head region, namely at those pH values at which the concentrations of the protonated and deprotonated forms of these groups are comparable. The model predicts the correct order of magnitude of the experimental frequency dispersion shown in Figs. 3 and 4.


View larger version (12K):
[in this window]
[in a new window]
 
FIGURE 6   Plots of Delta C versus pH for PS (black-square) and PA (triangle ) self-assembled monolayers supported by mercury, calculated as described in the text.

The solid curve b in Fig. 1 shows the Scalc versus pH plot for the PA film in best agreement with the corresponding Sexp versus pH plot, as calculated by assuming that the proton loosely bound to the >Cdouble bond O group of one of the ester groups of the lipid is located at x = beta  and does not follow the 75-Hz AC signal. Conversely, the phosphate group is located at x = d and follows the AC signal. The protonation constant K3 of the group buried inside the polar head region that provides the best agreement with experiment equals 50 M-1, whereas those for the two consecutive protonation equilibria of the phosphate group are equal to K1 = 1 × 108 M-1 and K2 = 1 × 105 M-1. Curve b in Fig. 6 shows the plot of Delta C versus pH for the PA film, where Delta C is the increment in the calculated value of C if the protonation equilibrium of the group located at x = beta  were perfectly mobile.

The surface dipole potential

The plot of chi -psi 0 = (chi  + const.) versus sigma  resulting from the use of the model in Fig. 5 is shown in Fig. 7; it exhibits a sigmoidal shape, with the maximum slope lying in the proximity of sigma  = 0. This plot is reminiscent of the surface dipole potential chi w due to the water molecules adsorbed at a metal/water interface as a function of the charge density sigma M on the metal. For comparison, the dashed curve in Fig. 7 is a plot of chi w versus sigma M, as calculated by Damaskin and Frumkin (1974) on the basis of a simple model of the metal/water interface. Moreover, the magnitude of the maximum change of chi  in the (chi  + const.) versus sigma  plots of Fig. 7 is comparable with that estimated by Trasatti (1975) (~300 mV) for the interface between the hydrophilic liquid gallium and water on the basis of several pieces of experimental evidence combined with a minimum of modelistic assumptions. This strongly suggests that the change in chi  with varying charge density of the PS and PA monolayers is mainly to be ascribed to the reorientation of the water molecules in contact with the polar heads; this conclusion is supported by the consideration that the only dipoles that experience the whole charge sigma lip on the lipid must lie outside the lipid film. Similar conclusions as to the molecular origin of the surface dipole potential in lipid films were drawn by Gawrisch et al. (1992) and by Zheng and Vanderkooi (1992). The surface dipole potential associated with the ester linkages to the glycerol backbone, which has been regarded as responsible for the higher permeability in lipid bilayers of lipophilic anions with respect to cations (McLaughlin, 1977; Honig et al., 1986), is apparently unaffected by a change in sigma lip. Naturally, some caution must be used in transferring these conclusions to biological membranes, which incorporate integral proteins protruding for 10 Å or so outside the lipid leaflet. Nonetheless, over the patches of the lipid leaflet free from proteins, the contribution of water reorientation to the surface dipole potential is expected to be appreciable.


View larger version (13K):
[in this window]
[in a new window]
 
FIGURE 7   Plots of (chi -psi 0 + const) versus sigma  for PS (black-square) and PA (triangle ) self-assembled monolayers at -0.5 V/SCE in buffered solutions of 0.1 M KCl, calculated as described in the text. The dashed curve is a chi w versus sigma M plot calculated by Damaskin and Frumkin (1974) for a metal/water interphase.

It should be noted that a plot of (chi  + const.) versus sigma  with a sigmoidal shape and the maximum slope lying in the proximity of sigma  = 0 are also obtained by using the crude two-capacitor model adopted in Moncelli et al. (1994), in which the ionizable groups are assumed to be in direct contact with the aqueous phase and their protonation equilibria are perfectly mobile; naturally, with this model the protonation constants providing the best fit between the Scalc versus pH plots and the corresponding Sexp versus pH plots of Fig. 1 assume different values. Hence the sigmoidal dependence of the surface dipole potential upon the charge is not subordinated to the assumption of a lack of proton equilibration.

    CONCLUSIONS
Top
Abstract
Introduction
Results
Discussion
Conclusions
Appendix
References

A completely unambiguous explanation cannot be found for the slight increase in differential capacity as the frequency of the AC signal is decreased below 80 Hz (see Figs. 3 and 4), although it may be justified either by a sluggish collective reorientation of some polar groups of the lipids after the AC signal, or else by a sluggish movement of protons from partially protonated ionizable groups buried inside the polar head region to the bathing solution, and vice versa. The lack of frequency dispersion shown by PC monolayers over a broad pH range from 4 to 9 (see Fig. 4) can be justified equally well on the basis of any of the above two tentative arguments. Thus, over this pH range, the PC film is uncharged and does not contain partially protonated ionizable groups such as to justify a movement of protons. On the other hand, over this pH range the conformation of the PC polar heads with the P-N dipoles aligned head to tail in the directions parallel to the monolayer is the most energetically favored arrangement from an electrostatic viewpoint, such as to resist changes after the AC signal. As the pH is decreased below 3, the incipient protonation of the phosphate groups begins to convert a number of P-N zwitterions into -N(CH3)3+ cations, thus undermining the network of parallel P-N dipoles. As a result, the residual P-N dipoles will tend to assume a tilted orientation that, by creating a favorable potential difference across the polar head region, will cause the protons to be attracted toward the innermost portion of this region and to protonate the phosphate groups there (Moncelli et al., 1994). Hence, at pH less than 3, the AC signal may cause either a very small fluctuation in the tilt of the P-N dipoles or a movement of protons from the partially protonated phosphate groups to the aqueous phase, and vice versa: either of these two movements may be sluggish enough to lag behind the AC signal at frequencies greater than 80 Hz, justifying the slight frequency dispersion shown by PC at pH 1.8 (see Fig. 4).

    APPENDIX
Top
Abstract
Introduction
Results
Discussion
Conclusions
Appendix
References

The sum of the charge density sigma M on the metal surface and that of the lipid, sigma lip = Sigma isigma i + Sigma jsigma j, is equal and opposite that of the diffuse-layer ions. According to the GC theory, this statement is expressed by the following implicit function:
&sfgr;<SUB><UP>M</UP></SUB>+<LIM><OP>∑</OP><LL><UP>i</UP></LL></LIM> &sfgr;<SUB><UP>i</UP></SUB>+<LIM><OP>∑</OP><LL><UP>j</UP></LL></LIM> &sfgr;<SUB><UP>j</UP></SUB>−A(1/y−y) ≡ g=0 (A1)
with
A ≡ <RAD><RCD><FR><NU>RT&egr;(c+c<SUB><UP>H<SUP>+</SUP></UP></SUB>)</NU><DE>2&pgr;</DE></FR></RCD></RAD>;  y ≡ <UP>exp</UP><FENCE><UP>−</UP><FR><NU>F&psgr;<SUB><UP>d</UP></SUB></NU><DE>2RT</DE></FR></FENCE> (A2)
Here c is the bulk concentration of the 1,1-valent electrolyte used to vary the ionic strength, cH+ is the bulk concentration of hydrogen ions, and epsilon  = 78 is the dielectric constant of the solvent.

Let us assume that the protons bind to the ionizable groups of the lipid according to a Langmuir isotherm. Moreover, let us denote by sigma max,i (sigma max,j) the maximum charge density attainable by the ith (jth) ionizable group. The charge density sigma i of the ith ionizable group located at x = beta  will then be given by
&sfgr;<SUB><UP>i</UP></SUB>=<FR><NU>&sfgr;<SUB><UP>max,i</UP></SUB></NU><DE>1+K<SUB><UP>i</UP></SUB>[<UP>H<SUP>+</SUP></UP>]<SUB>&bgr;</SUB></DE></FR>  <UP>or</UP>  &sfgr;<SUB><UP>i</UP></SUB>=<FR><NU>&sfgr;<SUB><UP>max,i</UP></SUB>K<SUB><UP>i</UP></SUB>[<UP>H<SUP>+</SUP></UP>]<SUB>&bgr;</SUB></NU><DE>1+K<SUB><UP>i</UP></SUB>[<UP>H<SUP>+</SUP></UP>]<SUB>&bgr;</SUB></DE></FR> (A3)
depending on whether sigma max, i is negative or positive. Here Ki is the protonation constant of the ith group, and [H+]beta is the proton concentration at x = beta  as affected by the local potential phi beta . In view of the expression of Eq. 1 for phi beta , we have
[<UP>H<SUP>+</SUP></UP>]<SUB>&bgr;</SUB>=c<SUB><UP>H<SUP>+</SUP></UP></SUB> <UP>exp</UP><FENCE><UP>−</UP><FR><NU>F&phgr;<SUB>&bgr;</SUB></NU><DE>RT</DE></FR></FENCE>=c<SUB><UP>H<SUP>+</SUP></UP></SUB> <UP>exp</UP><FENCE><UP>−</UP><FR><NU>F</NU><DE>RT</DE></FR> <FR><NU>&ggr;′</NU><DE>&ggr;′+&bgr;′</DE></FR> &psgr;′<SUB>0</SUB></FENCE>y<SUP>2&bgr;′/(&ggr;′<UP>+</UP>&bgr;′)</SUP> (A4)
where we have set
&psgr;′<SUB>0</SUB> ≡ &psgr;<SUB>0</SUB>−&khgr;;  &bgr;′ ≡ &bgr;/&egr;<SUB>&bgr;</SUB>;  &ggr;′ ≡ &ggr;/&egr;<SUB>&ggr;</SUB> (A5)
to simplify notations. By analogy with Eq. A3, the charge density sigma j of the jth ionizable group at the boundary x = d of the lipid film with the aqueous phase is given by
&sfgr;<SUB><UP>j</UP></SUB>=<FR><NU>&sfgr;<SUB><UP>max,j</UP></SUB></NU><DE>1+K<SUB><UP>j</UP></SUB>c<SUB><UP>H<SUP>+</SUP></UP></SUB>y<SUP>2</SUP></DE></FR>  <UP>or</UP>  &sfgr;<SUB><UP>j</UP></SUB>=<FR><NU>&sfgr;<SUB><UP>max,j</UP></SUB>K<SUB><UP>j</UP></SUB>c<SUB><UP>H<SUP>+</SUP></UP></SUB>y<SUP>2</SUP></NU><DE>1+K<SUB><UP>j</UP></SUB>c<SUB><UP>H<SUP>+</SUP></UP></SUB>y<SUP>2</SUP></DE></FR> (A6)
depending on whether sigma max,j is negative or positive; Kj is the protonation constant of the jth group, and cH+ exp(-Fpsi d/RT) cH+y2 is the hydrogen ion concentration at x = d. Strictly speaking, Eqs. A3 and A6 do not apply to the consecutive protonations of a multiply charged ionizable group, such as the -PO42- group of PA. However, it can be readily shown that the above equations are still valid, provided that the first protonation constant is both much greater than the second and much greater than 1/cH+ over the pH range investigated; these two requirements are actually satisfied by the -PO42- group of PA over the pH range covered by our measurements.

A further relationship is provided by Eq. 2, which can be written in the implicit form:
&psgr;′<SUB>0</SUB>−4&pgr;&bgr;′&sfgr;<SUB><UP>M</UP></SUB>−4&pgr;&ggr;′(&sfgr;<SUB><UP>M</UP></SUB>+∑<SUB><UP>i</UP></SUB>&sfgr;<SUB><UP>i</UP></SUB>)+2 <FR><NU>RT</NU><DE>F</DE></FR> <UP>ln</UP> y ≡ h=0 (A7)
The expressions g = 0 and h = 0 of Eqs. A1 and A7 are functions of sigma M, {sigma i}, {sigma j}, and y, which are not independent variables but are functions of only two independent variables. The real roots of the two equations A1 and A7 were obtained by the Newton-Raphson iterative procedure. To ensure a rapid convergence, it was found convenient to choose as independent variables y and the average potential difference across the polar-head region:
&psgr;<SUB>1</SUB> ≡ 4&pgr;&ggr;′<FENCE>&sfgr;<SUB><UP>M</UP></SUB>+<LIM><OP>∑</OP><LL><UP>i</UP></LL></LIM> &sfgr;<SUB><UP>i</UP></SUB></FENCE>=&psgr;′<SUB>0</SUB>−4&pgr;&bgr;′&sfgr;<SUB><UP>M</UP></SUB>+2 <FR><NU>RT</NU><DE>F</DE></FR> <UP>ln</UP> y (A8)
Moreover, it was generally necessary to mix successive approximations to y and psi 1 before they were used as input at the next level of iteration.

When the electrolyte concentration c is varied while pH is kept constant, the surface dipole potential chi  can be regarded as satisfactorily constant; hence the same is true for the quantity psi '0 = (psi 0 - chi ), because the applied potential E is also kept constant. However, psi '0 is unknown, whereas the charge density sigma M at c = 0.1 M is known, having been measured at each pH value. The first run of the Newton-Raphson iterative procedure at constant pH and variable c was therefore carried out for c = 0.1 M, using the known value of sigma M and regarding g and h as functions of y and psi '0. The partial derivatives of g and h with respect to y and psi '0 to be used in the iterative procedure were therefore obtained from g in the form of Eq. A1 and from h in the form of Eq. A7, with sigma M constant and sigma i and sigma j expressed by Eqs. A3, A4, and A6. The psi '0 value obtained from the first run was then employed in the subsequent runs for the same pH and the remaining concentrations c. In these further runs, g and h were regarded as functions of psi 1 and y, by writing them in the form
g=<FR><NU>&psgr;′<SUB>0</SUB>−&psgr;<SUB>1</SUB>+2(RT/F) <UP>ln</UP> y</NU><DE>4&pgr;&bgr;′</DE></FR>
+∑<SUB><UP>i</UP></SUB>&sfgr;<SUB><UP>i</UP></SUB>+∑<SUB><UP>j</UP></SUB>&sfgr;<SUB><UP>j</UP></SUB>−A(1/y−y)=0 (9)
and
h=<FR><NU>&bgr;′+&ggr;′</NU><DE>&bgr;′</DE></FR> &psgr;<SUB>1</SUB>−<FR><NU>&ggr;′</NU><DE>&bgr;′</DE></FR><FENCE>&psgr;′<SUB>0</SUB>+2 <FR><NU>RT</NU><DE>F</DE></FR> <UP>ln</UP> y</FENCE>−4&pgr;&ggr;′∑<SUB><UP>i</UP></SUB>&sfgr;<SUB><UP>i</UP></SUB>=0 (A10)
and differentiating them with respect to psi 1 and y at constant psi '0. This procedure was repeated for all pH values.

To determine the set of parameters {Ki}, {Kj} and gamma ' (or beta ') that provide the best fit between Scalc versus pH plots and Sexp versus pH plots, it is necessary to calculate the differential capacity C at different c and pH values. An expression for 1/C is obtained from Eqs. A7 and A8:
<FR><NU>1</NU><DE>C</DE></FR>=<FR><NU><UP>d</UP>&psgr;′<SUB>0</SUB></NU><DE><UP>d</UP>&sfgr;<SUB><UP>M</UP></SUB></DE></FR>=4&pgr;&bgr;′+<FR><NU><UP>d</UP>&psgr;<SUB>1</SUB></NU><DE><UP>d</UP>&sfgr;<SUB><UP>M</UP></SUB></DE></FR>−2 <FR><NU>RT</NU><DE>Fy</DE></FR> <FR><NU><UP>d</UP>y</NU><DE><UP>d</UP>&sfgr;<SUB><UP>M</UP></SUB></DE></FR> (A11)
C was determined by considering that, in differential capacity measurements, psi '0 varies sinusoidally about the constant value obtained at each pH by the Newton-Raphson procedure. Hence g and h must now be regarded as functions of sigma M, psi 1, and y. If we denote the partial derivatives of the functions g and h with respect to sigma M, psi 1, and y by the subscripts M, 1, and y, from the rules for the derivation of implicit functions we get
<FR><NU>1</NU><DE>C</DE></FR>=4&pgr;&bgr;′+<FR><NU>g<SUB><UP>y</UP></SUB>h<SUB><UP>M</UP></SUB>−h<SUB><UP>y</UP></SUB>g<SUB><UP>M</UP></SUB></NU><DE>g<SUB>1</SUB>h<SUB><UP>y</UP></SUB>−h<SUB>1</SUB>g<SUB><UP>y</UP></SUB></DE></FR>−2 <FR><NU>RT</NU><DE>Fy</DE></FR> <FR><NU>g<SUB><UP>M</UP></SUB>h<SUB>1</SUB>−h<SUB><UP>M</UP></SUB>g<SUB>1</SUB></NU><DE>g<SUB>1</SUB>h<SUB><UP>y</UP></SUB>−h<SUB>1</SUB>g<SUB><UP>y</UP></SUB></DE></FR> (A12)
The expression of g as a function of sigma M, psi 1, and y is provided by Eq. A1, in which sigma i and sigma j are given by Eqs. A3 and A6, and [H+]beta takes the form
[<UP>H<SUP>+</SUP></UP>]<SUB>&bgr;</SUB>=c<SUB><UP>H<SUP>+</SUP></UP></SUB>y<SUP>2</SUP> <UP>exp</UP><FENCE><UP>−</UP><FR><NU>F</NU><DE>RT</DE></FR> <FR><NU>&ggr;′</NU><DE>&bgr;′+&ggr;′</DE></FR>(4&pgr;&bgr;′&sfgr;<SUB><UP>M</UP></SUB>+&psgr;<SUB>1</SUB>)</FENCE> (A13)
Analogously, the expression of h as a function of sigma M, psi 1, and y is
h=&psgr;<SUB>1</SUB>−4&pgr;&ggr;′(&sfgr;<SUB><UP>M</UP></SUB>+∑<SUB><UP>i</UP></SUB>&sfgr;<SUB><UP>i</UP></SUB>) (A14)
where sigma i is given by Eqs. A3 and A13.

    ACKNOWLEDGMENTS

Thanks are due to ENEA, Italy, for a Ph.D. fellowship to FTB, during the tenure of which the present results were obtained. The financial support of the Ministero dell'Università e della Ricerca Scientifica e Tecnologica and of the Consiglio Nazionale delle Ricerche is gratefully acknowledged.

    FOOTNOTES