Aromatic N-heterocycles have been used in electrochemical CO2 reduction, but their precise role is not yet fully understood. We used first-principles quantum chemistry to determine how the molecular sizes and substituent groups of these molecules affect their standard redox potentials involving various proton and electron transfers. We then use that data to generate molecular Pourbaix diagrams to find the electrochemical conditions at which the aromatic N-heterocycle molecules could participate in multiproton and electron shuttling in accordance with the Sabatier principle. While one electron standard redox potentials for aromatic N-heterocycles can vary significantly with molecule size and the presence of substituent groups, the two-electron and two-proton standard redox potentials depend much less on structural modifications and substituent groups. This indicates that a wide variety of aromatic N-heterocycles can participate in proton, electron, and/or hydride shuttling under suitable electrochemical conditions.