O to bR transition in bacteriorhodopsin occurs through a proton hole mechanism

Abstract

Extensive classical and quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations are used to establish the structural features of the O state in bacteriorhodopsin (bR) and its conversion back to the bR ground state. The computed free energy surface is consistent with available experimental data for the kinetics and thermodynamics of the O to bR transition. The simulation results highlight the importance of the proton release group (PRG, consisting of Glu194/204) and the conserved arginine 82 in modulating the hydration level of the protein cavity. In particular, in the O state, deprotonation of the PRG and downward rotation of Arg82 lead to elevated hydration level and a continuous water network that connects the PRG to the protonated Asp85. Proton exchange through this water network is shown by ∼0.1-μs semiempirical QM/MM free energy simulations to occur through the generation and propagation of a proton hole, which is relayed by Asp212 and stabilized by Arg82. This mechanism provides an explanation for the observation that the D85S mutant of bacteriorhodopsin pumps chloride ions. The electrostatics–hydration coupling mechanism and the involvement of all titration states of water are likely applicable to many biomolecules involved in bioenergetic transduction.