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.