Proteins are polyelectrolytes with acidic and basic amino acids Asp, Glu, Arg, Lys, and His, making up ≈25% of the residues. The protonation state of residues, cofactors, and ligands defines a “protonation microstate”. In an ensemble of proteins some residues will be ionized and others neutral, leading to a mixture of protonation microstates rather than in a single one as is often assumed. The microstate distribution changes with pH. The protein environment also modifies residue proton affinity so microstate distributions change in different reaction intermediates or as ligands are bound. Particular protonation microstates may be required for function, while others exist simply because there are many states with similar energy. Here, the protonation microstates generated in Monte Carlo sampling in MCCE are characterized in HEW lysozyme as a function of pH and bacterial photosynthetic reaction centers (RCs) in different reaction intermediates. The lowest energy and highest probability microstates are compared. The ΔG, ΔH, and ΔS between the four protonation states of Glu35 and Asp52 in lysozyme are shown to be calculated with reasonable precision. At pH 7 the lysozyme charge ranges from 6 to 10, with 24 accepted protonation microstates, while RCs have ≈50,000. A weighted Pearson correlation analysis shows coupling between residue protonation states in RCs and how they change when the quinone in the QB site is reduced. Protonation microstates can be used to define input MD parameters and provide insight into the motion of protons coupled to reactions.

The photosynthetic bacterial reaction centers from purple non-sulfur bacteria use light energy to drive the transfer of electrons from cytochrome c to ubiquinone. Ubiquinone bound in the QA site cycles between quinone, QA, and anionic semiquinone, QA •-, being reduced once and never binding protons. In the QB site, ubiquinone is reduced twice by QA •-, binds two protons and is released into the membrane as the quinol, QH2. The network of hydrogen bonds formed in a molecular dynamics trajectory was drawn to investigate proton transfer pathways from the cytoplasm to each quinone binding site. QA is isolated with no path for protons to enter from the surface. In contrast, there is a complex and tangled network requiring residues and waters that can bring protons to QB. There are three entries from clusters of surface residues centered around HisH126, GluH224, and HisH68. The network is in good agreement with earlier studies, Mutation of key nodes in the network, such as SerL223, were previously shown to slow proton delivery. Mutational studies had also shown that double mutations of residues such as Asp M17 and AspL210 along multiple paths in the network presented here slow the reaction, while single mutations do not. Likewise, mutation of both HisH126 and HisH128, which are at the entry to two paths reduce the rate of proton uptake.