Cation binding to beta-casein. A comparison of electrostatic models

The binding of cations of β-casein at pH 6.6 was considered previously. Available for three sodium concentiations, I = 0.04, 0.08, or 0.16 M are: [1] proton releases between I and [2] for each I, as calcium activity is increased, correlated sequences of monomer net charge, proton release, site bound calcium and protein Solvation- Models for ion binding are examined. Critical considerations are the intrinsic binding constants between hydrogen[H], calcium[Ca] and sodium[Na] ions and phosphate[P] and caiboxyIate[C] sites, and the effects of electrostatic interaction between sites as influenced by spatial fixed charge distribution, ionic strength and dielectric constant [D]. Anticipated intrinsic binding constants are k_H,P^o = 3 × 10⁶, k_Ca,P^o = 120, k_Na,P^o = 1, k_H,C^o = 7 × 10⁴ and k_Ca,C^o = 5.6Distributed charge models, either surface or volume, are inadequate since any reasonable monomer size yields fixed charge densities requiring k_H,P^o and k_Ca,C^o which are too low when the maximum in D is 75. Also, with increasing calcium binding, calculated proton release is only 0.4 to 0.5 of that observed.Discrete charge models accept anticipated k^o and yield calculated sequences of calcium binding and proton release which are in good agreement with those observed provided that: (1) using the known amino acid sequence of the phosphate-containing acidic peptide portion of the molecule, pep tide fixed charge is distributed at the lowest I so as to minimize electrostatic free energy; (2) in the region of fixed charge, D is approximately 5; (3) the distances between peptide fixed charges decrease with increasing ionic strength or calcium binding and (4) while protein is in solution, the acidic peptide and the remainder of the molecule are essentially electrostatically independent.     

Cation binding to β-casein. a comparison of electrostatic models

The binding of cations of β-casein at pH 6.6 was considered previously. Available for three sodium concentiations, I = 0.04, 0.08, or 0.16 M are: [1] proton releases between I and [2] for each I, as calcium activity is increased, correlated sequences of monomer net charge, proton release, site bound calcium and protein Solvation- Models for ion binding are examined. Critical considerations are the intrinsic binding constants between hydrogen[H], calcium[Ca] and sodium[Na] ions and phosphate[P] and caiboxyIate[C] sites, and the effects of electrostatic interaction between sites as influenced by spatial fixed charge distribution, ionic strength and dielectric constant [D]. Anticipated intrinsic binding constants are k_H,P^o = 3 × 10⁶, k_Ca,P^o = 120, k_Na,P^o = 1, k_H,C^o = 7 × 10⁴ and k_Ca,C^o = 5.6Distributed charge models, either surface or volume, are inadequate since any reasonable monomer size yields fixed charge densities requiring k_H,P^o and k_Ca,C^o which are too low when the maximum in D is 75. Also, with increasing calcium binding, calculated proton release is only 0.4 to 0.5 of that observed.Discrete charge models accept anticipated k^o and yield calculated sequences of calcium binding and proton release which are in good agreement with those observed provided that: (1) using the known amino acid sequence of the phosphate-containing acidic peptide portion of the molecule, pep tide fixed charge is distributed at the lowest I so as to minimize electrostatic free energy; (2) in the region of fixed charge, D is approximately 5; (3) the distances between peptide fixed charges decrease with increasing ionic strength or calcium binding and (4) while protein is in solution, the acidic peptide and the remainder of the molecule are essentially electrostatically independent.     

Calbindin D(9k): a protein optimized for calcium binding at neutral pH.

The binding of calcium ions by EF-hand proteins depends strongly on the electrostatic interactions between Ca(2+) ions and negatively charged residues of these proteins. We have investigated the pH dependence of the binding of Ca(2+) ions by calbindin D(9k). This protein offers a unique possibility for interpretation of such data since the pK(a) values of all ionizable groups are known. The binding is independent of pH between 7 and 9, where maximum calcium affinity is observed. An abrupt decrease in the binding affinity is observed at pH values below 7. This decrease is due to protonation of acidic groups, leading to modification of protein charges. The pH dependence of the product of the two macroscopic Ca(2+)-binding constants can be formally described by the involvement of two acidic groups with pK(a) = 6.6. Monte Carlo calculations show that the reduction of Ca(2+) binding is strictly determined by variable electrostatic interactions due to pH-dependent changes not only in the binding sites, but also of the overall charge of the protein.     

Focusing of the electrostatic potential at EF-hands of calbindin D(9k): titration of acidic residues

Biological functions for a large class of calmodulin-related proteins, such as target protein activation and Ca(2+) buffering, are based on fine-tuned binding and release of Ca(2+) ions by pairs of coupled EF-hand metal binding sites. These are abundantly filled with acidic residues of so far unknown ionization characteristics, but assumed to be essential for protein function in their ionized forms. Here we describe the measurement and modeling of pK(a) values for all aspartic and glutamic acid residues in apo calbindin D(9k), a representative of calmodulin-related proteins. We point out that while all the acidic residues are ionized predominantly at neutral pH, the onset of proton uptake by Ca(2+) ligands with high pK(a) under these conditions may have functional implications. We also show that the negative electrostatic potential is focused at the bidental Ca(2+) ligand of each site, and that the potential is significantly more negative at the N-terminal binding site.     

Site-specific phosphorylation and point mutations of telokin modulate its Ca2+-desensitizing effect in smooth muscle

Forskolin and 8-bromoguanosine 3'-5'-cyclic monophosphate (8-Br-cGMP) induce phosphorylation of Ser-13 of telokin and relaxation of smooth muscle at constant calcium. Comparison with the effect of wild type with aspartate (D; to mimic phosphorylation) and alanine (A; non-phosphorylatable) mutants of telokin showed that the S13D mutant was more effective than wild type in relaxing smooth muscle at constant calcium. The efficacy of the Ser-13A, S12A, and S12D mutants was not significantly different from that of wild-type telokin. The effect of neither S13D nor Ser-13A was affected by 8-Br-cGMP, whereas the effect of wild type, S12A, and S12D was enhanced by 8-Br-cGMP, indicating the specificity of Ser-13 charge modification. Mutation of Ser-19 (a mitogen-activated protein kinase site) showed the S19A to be more effective than, and S19D to be not different from, wild-type telokin. The effect of both mutants was slightly enhanced by 8-Br-cGMP. A truncated (residues 1-142) form lacking the acidic C terminus had the same relaxant effect as wild-type telokin, whereas the C-terminal peptide (residues 142-155) had no effect. We conclude that site-specific modification of the N terminus modulates the Ca2+ -desensitizing effect of telokin on force.     

Thermodynamics of the binding of Ca2+ to porcine pancreatic phospholipase A2

The binding of Ca2+ to porcine pancreatic phospholipase A2 was studied by batch microcalorimetry. Enthalpies of binding at 25 degrees C were determined as a function of Ca2+ concentration in buffered solutions at pH 8.0 using both the Tris-HCl and Hepes-NaOH buffer systems. The calorimetric results indicate that protons are released on calcium binding and that in addition to the binding of the active-site calcium, there appears to be weak binding of a second Ca2+. Results from potentiometric titrations indicate that this proton release on binding Ca2+ arises from a change in pK of a histidine(s) functional group. The thermodynamic functions delta G0, delta H0 and delta S0 for calcium binding to phospholipase A2 have been determined. These results are compared with literature data for Ca2+ complex formation with some small molecules and also the protein troponin-C.     

Phytic acid-metal ion interactions. II. The effect of pH on Ca(II) binding

The interaction of phytic acid with Ca(II) has been studied by potentiometric titration and by measurement of free Ca(II) concentrations using an ion Selective electrode. With increasing Ca(II) concentration, the titration curve of phytic acid is displaced to regions of lower pH. In the binding of calcium ions to phytic acid, there is no evidence that significant binding occurs below approximately pH 5. Above this pH, the extent of binding is dependent upon both pH and the calcium to phytic acid ratios. Maximum binding obtains at a Ca(II):phytate ratio of 6 with 4.8 mol of Ca(II) bound per mol of phytate above pH ca. 8. Binding constants are apparently very large since binding isotherms at any Ca(II):phytate ratio are a linear function of the total calcium ion concentration. In all cases, binding occurs only when one or more phosphate groups have been converted to the oxo dianion form. The apparent pK' values (curve-fit parameters) that describe the potentiometric titration data are in good agreement with the constants evaluated from the binding of Ca(II) to phytate as a function of pH. Using CPK space-filling models, structures containing six metal ions in coordinate linkage to pairs of oxo dianions have been constructed and discussed within the framework of the axial conformation of phytic acid and the order of proton removal with an increase in pH based upon NMR studies.     

Electrostatic binding of proteins to membranes. Theoretical predictions and experimental results with charybdotoxin and phospholipid vesicles.

We previously applied the Poisson-Boltzmann equation to atomic models of phospholipid bilayers and basic peptides to calculate their electrostatic interactions from first principles (Ben-Tal, N., B. Honig, R. M. Peitzsch, G. Denisov, and S. McLaughlan. 1996. Binding of small basic peptides to membranes containing acidic lipids. Theoretical models and experimental results. Biophys. J. 71:561-575). Specifically, we calculated the molar partition coefficient, K (the reciprocal of the lipid concentration at which 1/2 the peptide is bound), of simple basic peptides (e.g., pentalysine) with phospholipid vesicles. The theoretical predictions agreed well with experimental measurements of the binding, but the agreement could have been fortuitous because the structure(s) of these flexible peptides is not known. Here we use the same theoretical approach to calculate the membrane binding of two small proteins of known structure: charybdotoxin (CTx) and iberiotoxin (IbTx); we also measure the binding of these proteins to phospholipid vesicles. The theoretical model describes accurately the dependence of K on the ionic strength and mol % acidic lipid in the membrane for both CTx (net charge +4) and IbTx (net charge +2). For example, the theory correctly predicts that the value of K for the binding of CTx to a membrane containing 33% acidic lipid should decrease by a factor of 10(5) when the salt concentration increases from 10 to 200 mM. We discuss the limitations of the theoretical approach and also consider a simple extension of the theory that incorporates nonpolar interactions.     

An alpha-helical peptide model for electrostatic interactions of proteins with DNA. The N terminus of RecA

A series of synthetic peptides have been studied as models for non-specific protein-DNA interactions. In an alpha-helical conformation, the charged amino acid residues of the N-terminal 24 residues of RecA protein are asymmetrically distributed; at neutral pH there is a +4 charge on one face of the helix and a -3 charge on the other face. Modeling suggests that the positive face of the helix can bind five DNA phosphate groups by electrostatic interactions. Circular dichroism (c.d.) spectra indicate that the analogous peptide, Rec24 (AIDENKQKALAAALGQIEKQFGKG-amide), is largely unstructured in water but becomes highly helical in the presence of DNA. Peptide titrations of fluorescent etheno-DNA confirm that the changes in the c.d. spectrum of the peptide are associated with binding, although a dependence of the c.d. signal on the degree of DNA saturation is observed, indicating that peptide can be bound in more than one conformation. At saturation the peptide binds to 5.0(+/- 0.5) DNA phosphate groups as predicted and the electrostatic nature of the binding is confirmed by a strong dependence on salt concentration. A "mutant" peptide where an acidic glutamate residue replaces an alanine on the basic face of the Rec24 helix exhibits weaker binding to single-stranded DNA, also consistent with the electrostatic nature of the proposed peptide-DNA interaction. Extending Rec24 by ten amino acid residues, where the additional residues do not participate in the helical motif, does not noticeably affect binding. Thus, we show experimentally that an asymmetric charge distribution on an alpha-helix can represent an important element for binding nucleic acids.     

Highly cooperative and hysteretic response of the skeletal muscle ryanodine receptor to changes in proton concentrations.

Ryanodine receptors are key molecules in excitation-contraction coupling of skeletal muscle. They form the pore of the calcium release channel, which is regulated by Ca and ATP. Multiple proton titration sites are involved in controlling the different open states of the channel, as indicated by the following: i) the channel had a biphasic response to changes in proton concentrations around neutral pH; ii) the activities of the channel were inhibited by acidic pHs in a highly cooperative manner; and iii) the channel exhibited pronounced hysteresis to changes in pH. Four distinct conductance states can be identified in the single ryanodine-activated calcium release channel. The distribution of the multiple conductance states depends on the level of [Ca], ATP, and pH in the recording solution. The data are consistent with the multimeric structure of the skeletal muscle ryanodine receptor.     

Binding of small basic peptides to membranes containing acidic lipids: theoretical models and experimental results.

We measured directly the binding of Lys3, Lys5, and Lys7 to vesicles containing acidic phospholipids. When the vesicles contain 33% acidic lipids and the aqueous solution contains 100 mM monovalent salt, the standard Gibbs free energy for the binding of these peptides is 3, 5, and 7 kcal/mol, respectively. The binding energies decrease as the mol% of acidic lipids in the membrane decreases and/or as the salt concentration increases. Several lines of evidence suggest that these hydrophilic peptides do not penetrate the polar headgroup region of the membrane and that the binding is mainly due to electrostatic interactions. To calculate the binding energies from classical electrostatics, we applied the nonlinear Poisson-Boltzmann equation to atomic models of the phospholipid bilayers and the basic peptides in aqueous solution. The electrostatic free energy of interaction, which arises from both a long-range coulombic attraction between the positively charged peptide and the negatively charged lipid bilayer, and a short-range Born or image charge repulsion, is a minimum when approximately 2.5 A (i.e., one layer of water) exists between the van der Waals surfaces of the peptide and the lipid bilayer. The calculated molar association constants, K, agree well with the measured values: K is typically about 10-fold smaller than the experimental value (i.e., a difference of about 1.5 kcal/mol in the free energy of binding). The predicted dependence of K (or the binding free energies) on the ionic strength of the solution, the mol% of acidic lipids in the membrane, and the number of basic residues in the peptide agree very well with the experimental measurements. These calculations are relevant to the membrane binding of a number of important proteins that contain clusters of basic residues.     

Conformational and thermodynamic properties of peptide binding to the human S100P protein

S100P is a member of the S100 subfamily of calcium-binding proteins that are believed to be associated with various diseases, and in particular deregulation of S100P expression has been documented for prostate and breast cancer. Previously, we characterized the effects of metal binding on the conformational properties of S100P and proposed that S100P could function as a Ca2+ conformational switch. In this study we used fluorescence and CD spectroscopies and isothermal titration calorimetry to characterize the target-recognition properties of S100P using a model peptide, melittin. Based on these experimental data we show that S100P and melittin can interact in a Ca2+-dependent and -independent manner. Ca2+-independent binding occurs with low affinity (Kd approximately 0.2 mM), has a stoichiometry of four melittin molecules per S100P dimer and is presumably driven by favorable electrostatic interactions between the acidic protein and the basic peptide. In contrast, Ca2+-dependent binding of melittin to S100P occurs with high affinity (Kd approximately 5 microM) has a stoichiometry of two molecules of melittin per S100P dimer, appears to have positive cooperativity, and is driven by hydrophobic interactions. Furthermore, Ca2+-dependent S100P-melittin complex formation is accompanied by significant conformational changes: Melittin, otherwise unstructured in solution, adopts a helical conformation upon interaction with Ca2+-S100P. These results support a model for the Ca2+-dependent conformational switch in S100P for functional target recognition.     

Characterisation of the interaction of the C-terminus of the dopamine D2 receptor with neuronal calcium sensor-1

NCS-1 is a member of the neuronal calcium sensor (NCS) family of EF-hand Ca(2+) binding proteins which has been implicated in several physiological functions including regulation of neurotransmitter release, membrane traffic, voltage gated Ca(2+) channels, neuronal development, synaptic plasticity, and learning. NCS-1 binds to the dopamine D2 receptor, potentially affecting its internalisation and controlling dopamine D2 receptor surface expression. The D2 receptor binds NCS-1 via a short 16-residue cytoplasmic C-terminal tail. We have used NMR and fluorescence spectroscopy to characterise the interactions between the NCS-1/Ca(2+) and D2 peptide. The data show that NCS-1 binds D2 peptide with a K(d) of ～14.3 μM and stoichiometry of peptide binding to NCS-1 of 2:1. NMR chemical shift mapping confirms that D2 peptide binds to the large, solvent-exposed hydrophobic groove, on one face of the NCS-1 molecule, with residues affected by the presence of the peptide spanning both the N and C-terminal portions of the protein. The NMR and mutagenesis data further show that movement of the C-terminal helix 11 of NCS-1 to fully expose the hydrophobic groove is important for D2 peptide binding. Molecular docking using restraints derived from the NMR chemical shift data, together with the experimentally-derived stoichiometry, produced a model of the complex between NCS-1 and the dopamine receptor, in which two molecules of the receptor are able to simultaneously bind to the NCS-1 monomer.     

Binding of calcium to lysozyme and its derivatives

Bindings of calcium to lysozyme and its derivatives were studied by UV difference spectroscopy at various pH's. The binding constant was ca. 40 m-1 at around neutral pH. The binding caused proton release from lysozyme and did not inhibit the binding of tri-N-acetylglucosamine to lysozyme. In the presence of 0.2 M Ca2+, lysozyme showed 26% of the activity of the free enzyme toward hexa-N-acetylglucosamine but the cleavage pattern was similar to that of the free enzyme. Thus, calcium was predicted to bind near the catalytic carboxyls to cause inhibition of lysozyme activity. It was found from the results of protease digestion that calcium binding shifted the native-denatured transition in lysozyme toward the native state.     

Cholecystokinin receptor occupation and cholecystokinin-induced calcium mobilization in the early phase in rat pancreatic acini.

We examined receptor occupation, calcium mobilization and amylase release for cholecystokinin octapeptide (CCK-8) within a 3-min incubation period at 37 degrees C using dispersed acini from rat pancreas. Analysis of competitive binding inhibition data obtained after a 3-min incubation revealed the presence of only a single class of CCK receptors, while two classes of CCK receptor, i.e., high-affinity and low-affinity CCK receptors, were detected when binding reached a steady-state after a 60-min incubation. The IC50 of CCK receptors calculated from the 3-min binding data was 19.0 +/- 0.5 nM (mean +/- S.D.), close to the Kd of the low-affinity CCK receptors determined by equilibrium binding studies. Exposure of fura-2-loaded acini to 10-1000 pM CCK-8 caused an immediate and dose-dependent increase in [Ca2+]i followed by a gradual decrease in [Ca2+]i. The CCK-stimulated amylase release after 3 min of incubation was biphasic; amylase release increased over the dose range of 3-300 pM CCK-8, peaked at 300 pM CCK-8 and decreased with supramaximal concentrations of CCK-8. Our data suggest that occupation of the low-affinity, but not the high-affinity, CCK receptors is more directly associated with calcium mobilization and subsequent stimulation of amylase release in rat pancreatic acini.     

Mg2+-induced proton release from Escherichia coli ribosome and ribosomal RNA

Escherichia coli ribosome released protons upon addition of Mg2+. The Mg2+-induced proton release was studied by means of the pH-stat technique. The number of protons released from a 70 S ribosome in the Mg2+ concentration range 1-20 mM was about 30 at pH 7 and 7.6, and increased to about 40 at pH 6.5. The rRNA mixture extracted from 70 S ribosome showed proton release of amount and of pH dependence similar to those of the 70 S ribosome but the ribosomal protein mixture released few. This indicates that rRNA is the main source of the protons released from ribosome. The pH titration of rRNA showed that the pKa values of nucleotide bases were downward shifted upon Mg2+ binding. This pKa shift can account for the proton release. The Scatchard plots of proton release from rRNA and ribosome were concave upward, showing that the Mg2+-binding sites leading to proton release were either heterogeneous or had a negative cooperativity. A model assuming heterogeneous Mg2+-binding sites is shown to be unable to explain the proton release. Electrostatic field effect models are proposed in which Mg2+ modulates the electrostatic field of phosphate groups and the potential change induces a shift of the pKa values of bases that leads to the proton release. These models can explain the main features of the proton release.     

A simple electrostatic model can explain the effect of pH upon the force-pCa relation of skinned frog skeletal muscle fibers.

The relative force-pCa relation of skinned frog skeletal muscle fibers is shifted along the pCa axis by changes in pH. This shift has been interpreted as arising from competition between H+ and Ca2+ for a binding site on troponin. Unfortunately, binding studies have been unable to confirm such competition. Alternatively, however, the data fit a model where H+ influences the degree of dissociation of ionizable groups on the surface of the thin filaments, thus altering the electrostatic potential surrounding the filaments. Alterations in the potential will, in turn, change the concentration of Ca2+ near the troponin binding sites in accordance with the Boltzmann relation. A simple model, based upon the Gouy-Chapman relation between surface potential and charge density, provides a quantitative explanation for the shift of the relative force-pCa curve with pH, given a reasonable estimate of the surface charge density on the thin filament. A best fit is obtained when the ionizable groups giving rise to the potential have a log proton ionization constant (pKa) of 6.1, similar to that for the imidazole group on histidine, and when the density of these groups is near that estimated from amino acid analysis of thin filament proteins and from filament geometry. In preliminary experiments, reaction of skinned frog fibers with diethylpyrocarbonate (DEP) at pH 6 shifted the force-pCa curve toward lower Ca2+. This would be expected in the model since DEP at pH 6 is reported to specifically react with histidine imidazole groups and to irreversibly decrease their pKa, which would increase the net negative charge of the filaments.     

Calculating proton uptake/release and binding free energy taking into account ionization and conformation changes induced by protein-inhibitor association: application to plasmepsin, cathepsin D and endothiapepsin-pepstatin complexes

The protein-inhibitor binding energies of enzymes are often pH dependent, and binding induces either proton uptake or proton release. The proton uptake/release and the binding energy for three complexes with available experimental data were numerically studied: pepstatin-cathepsin D, pepstatin-plasmepsin II and pepstatin-endothiapepsin. Very good agreement with the experimental data was achieved when conformational changes were taken into account. The role of the desolvation energy and the conformational changes was revealed by modeling the complex, the separated molecules in the complex conformation and the free molecules. It was shown that the conformational changes induced by the complex formation are as important for the proton transfer as the loss of solvation energy caused by the burial of interface residues. The residues responsible for the proton transfer were identified and their contribution to the proton uptake/release calculated. These residues were found to be scattered along the whole protein rather than being localized only at the active site. In the case of cathepsin D, these residues were found to be highly conserved among the cathepsin D sequences of other species. It was shown that conformation and ionization changes induced by the complex formation are critical for the correct calculation of the binding energy. Taking into account the electrostatics and the van der Waals (vdW) energies within the Boltzmann distribution of energies and allowing ionization and conformation changes to occur makes the calculated binding energy more realistic and closer to the experimental value. The interplay between electrostatic and vdW forces makes the pH dependence of the binding energy smoother, because the vdW force acts in reaction to the changes of the electrostatic energy. It was found that a small fraction of the ionizable groups remain uncharged in both the free and complexed molecules. The sequence and structural position of these groups aligns well within the three proteases, suggesting that these may have specific role.     

Ion-binding to phospholipids. Interaction of calcium with phosphatidylserine.

The binding of Ca2+ to monolayers and bilayers of phosphatidylserine has been investigated as a function of pH, ionic strength (NaCl concentration) and Ca2+ concentration using surface and colloid chemical techniques. The molar ratio of lipid to bound calcium decreases to 2 as the Ca2+ concentration is increased to about 0.1 mM. At [Ca2+] greater than 0.1 mM a 1:1 complex is formed. The apparent binding constant Ka ranges from about approximately 10(6) - 10(4) l/mol depending on the Ca2+ concentration. After allowing for electrostatic effects and neighbour group interactions, the intrinsic binding constant Ki of the phosphorylserine polar group at pH 7 (I = 0.01 M), where it carries a net negative charge of one, is approximately 10(4) l/mol; consistent values for Ki were obtained using several independent approaches. Ka for Ca2+ binding decreases with increasing NaCl concentration because the monovalent cations compete with Ca2+ for the same binding site. Na+ and K+ are equally effective in displacing 45Ca2+ adsorbed to monolayers of phosphatidylserine, both with respect to the kinetics and the equilibrium of the displacement. Ka for the reaction between phosphatidylserine and monovalent cations is about 10(3)-fold smaller than that of Ca2+. An investigation of the binding of Mn2+ to phosphatidylserine by both surface chemical and nuclear magnetic resonance methods shows that this cation has a similar binding constant to that of Ca2+. The Ca2+-binding capabilities of monolayers containing only carboxyl groups (i.e. arachidic acid) and phosphodiester groups (i.e. dicetyl phosphate) have also been determined; the apparent pK for the - COOH group in monolayers is larger than or equal to 9 and that for the phosphodiester group is less than 4. Since these groups do not retain the same pK values when they are in close proximity in the phosphorylserine group, the relative contributions of the two groups to the binding of Ca2+ to phosphatidylserine is not obvious.     

Calcium binding to mixed phosphatidylglycerol-phosphatidylcholine bilayers as studied by deuterium nuclear magnetic resonance

The binding of calcium to bilayer membranes composed of mixtures, in various proportions, of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) plus 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG) was investigated by using atomic absorption spectroscopy and deuterium nuclear magnetic resonance. The number of bound calcium ions, X2, was determined in the low calcium concentration range (up to 100 mM) via atomic absorption spectroscopy. Simultaneous measurements of the deuterium magnetic resonance spectra of POPC, specifically deuteriated at the alpha-methylene segment of the choline head group, revealed a linear relationship between the quadrupole splitting, delta vQ, and X2 for each particular proportion of POPC-POPG. The amount of bound calcium was then determined at much greater calcium concentrations, where the atomic absorption spectroscopy measurements were unreliable, using deuterium magnetic resonance. At low Ca2+ concentrations, the amount of bound Ca2+ increased linearly with increasing proportion of POPG, demonstrating an electrostatic contribution to Ca2+ binding. At high Ca2+ concentrations, the calcium binding isotherms exhibited saturation behavior with a maximum binding capacity of 0.5 Ca2+ and 1.0 Ca2+ per phospholipid for pure POPC and mixtures of POPC-POPG, respectively. Simultaneous deuteriation of POPG and POPC showed that both lipids remained in a fluidlike lipid bilayer at all Ca2+ concentrations tested. Any phase separation of quasi-crystalline Ca2+-POPG clusters could be excluded. The residence time of Ca2+ at an individual head group binding site was shorter than 10(-6)-10(-5) s. Thus, Ca2+ ions accumulate near the negatively charged POPG-POPC membrane surface but move freely in a "trough" of the electrical potential. The effective surface charge density, sigma, could be determined from the measured amount of bound Ca2+. Subsequently, the surface potential, psi 0, and the concentration of free Ca2+ ions at the plane of ion binding could be calculated by employing the Gouy-Chapman theory. The availability of these parameters allowed a rigorous evaluation of various models for the chemical contribution to Ca2+ binding. For mixed POPC-POPG bilayers, a simple Langmuir adsorption model yielded the best fit to the experimental data, and the binding constants were 19.5 and 18.8 M-1 for POPG contents of 20 and 50 mol %, respectively. Sodium binding was comparatively weak with a binding constant of 0.6-0.85 M-1.(ABSTRACT TRUNCATED AT 400 WORDS)