A modeling method is presented for protein systems in which proton transport is coupled to conformational change, as in proton pumps and in motors driven by the proton-motive force. Previously developed methods for calculating pK(a) values in proteins using a macroscopic dielectric model are extended beyond the equilibrium case to a master-equation model for the time evolution of the system through states defined by ionization microstate and a discrete set of conformers. The macroscopic dielectric model supplies free energy changes for changes of protonation microstate, while the method for obtaining the energetics of conformational change and the relaxation rates, the other ingredients needed for the master equation, are system dependent. The method is applied to the photoactivated proton pump, bacteriorhodopsin, using conformational free energy differences from experiment and treating relaxation rates through three adjustable parameters. The model is found to pump protons with an efficiency relatively insensitive to parameter choice over a wide range of parameter values, and most of the main features of the known photocycle from very early M to the return to the resting state are reproduced. The boundaries of these parameter ranges are such that short-range proton transfers are faster than longer-range ones, which in turn are faster than conformational changes. No relaxation rates depend on conformation. The results suggest that an "accessibility switch", while not ruled out, is not required and that vectorial proton transport can be achieved through the coupling of the energetics of ionization and conformational states.