Motivated by biological membrane-containing organelles in plants and photosynthetic bacteria, we study charge regulation in a model membrane stack. Considering (de)protonation as the simplest mechanism of charge equilibration between the membranes and with the bathing environment, we uncover a symmetrybroken charge state in the stack with a quasiperiodic effective charge sequence. In the case of a monovalent bathing salt solution our model predicts complex, inhomogeneous charge equilibria depending on the strength of the (de)protonation reaction, salt concentration, intermembrane separation, and their number in the stack. Our results shed light on the basic reorganization mechanism of thylakoid membrane stacks. * majee@is.mpg.de arXiv:1909.00217v2 [cond-mat.soft]