The binding of charged drugs to neutral phosphatidylcholine membranes was assessed by measuring their zeta-potential values in the presence of different drug concentrations. This methodology was applied to the study of the concentration effects of two nonsteroidal antiinflammatory drugs (NSAIDs). Results revealed an intense membrane charging that was proportional to the amount of negatively charged drug in the media. A mathematical formalism was adapted and an analytical expression derived to calculate directly surface potentials from zeta-potential data. The membrane loading state, expressed as the number of molecules per unit area, was calculated for the negative and for the neutral forms of the drugs. An approach was also developed that allows the determination of the maximum number of molecules per unit area by fitting a binding isotherm to the dependence of the number molecules per unit area with the drug concentration. The calculation of the maximum mol lipid/drug ratio can also be estimated and related to the binding stoichiometry, as well as to the maximum lipid loading capacity. Furthermore, the concentration profiles for both drugs can be established in terms of the distance to the liposome surface. The developed methodology allowed for the simultaneous determination of partition coefficients (K-p) for the NSAIDs in lipid/aqueous media because zeta-potential values can be related to the drug concentration at the lipid/ aqueous media interface. Alternative independent methodologies were used to determine K-p: spectrophotometric and centrifugation assays. A mathematical relation was developed to compare the K-p values determined from the zeta-potential data with those obtained from the other techniques used because in the former case they are calculated on the basis of the number of molecules per unit area and in the latter on the basis of the total drug concentrations in solution, and the values of the partition coefficients obtained from all the techniques were found to be equal, within the experimental error. This methodology constitutes a more straightforward method than the other techniques used because partition coefficients for all drug forms (charged and noncharged) can be assessed with a minimum number of experimental determinations and it allows for a characterization of the electrostatic properties of neutral membranes upon binding of charged drugs.