Biological membranes are vital for diverse cellular functions such as maintaining cell and organelle structure, selective permeability, active transport, and signaling. The surface charge of the membrane bilayer plays a critical role in these myriad processes. For most biomembranes, the surface charge of anionic phospholipids contributes to the negative surface charge density within the interfacial region of the bilayer. To quantify surface charge, it is essential to understand the proton dissociation behavior of the titratable headgroups within such lipids. We describe a protocol that uses model membranes for electrokinetic zeta potential measurements coupled with data analysis using Gouy-Chapman-Stern formalism to determine the pKa value of the component lipids. A detailed example is provided for homogeneous bilayers composed of the monoanionic lipid phosphatidylglycerol. This approach can be adapted for the measurement of bilayers with a heterogeneous lipid combination, as well as for lipids with multiple titratable sites in the headgroup (e.g., cardiolipin).