Abstract
The high-resolution structure of heliorhodopsin crystallized at low pH reveals the presence of a planar triangle molecule, acetate, in the inner water cavity. Here, we investigate how the acetate molecule is stabilized at the counterion Glu107 moiety, using molecular dynamics (MD) simulations and a quantum mechanical/molecular mechanical (QM/MM) approach. QM/MM calculations indicate that the density is best described as acetate among triangle acids, including nitric acid and bicarbonate. The calculated protonation state indicates that protonated acetate donates an H-bond to deprotonated Glu107 in the low-pH crystal structure. The observed red-shift of ∼30 nm in the absorption wavelength with pK(a) ≈ 4 is likely due to the His23/His80 protonation, rather than the Glu107 protonation. MD simulations also show that acetate can exist at the Glu107 moiety only when it is protonated. When ionized, acetate is released from the Glu107 moiety via Asn101 at the channel bottleneck and Arg91 on the intracellular protein surface. These observations could explain how acetate binds at low pH and releases at high pH.