Key role of electrostatic interactions in bacteriorhodopsin proton transfer

The first proton transport step following photon absorption in bacteriorhodopsin is from the 13-cis retinal Schiff base to Asp85. Configurational and energetic determinants of this step are investigated here by performing quantum mechanical/molecular mechanical minimum-energy reaction-path calculations. The results suggest that retinal can pump protons when in the 13-cis, 15-anti conformation but not when 13-cis, 15-syn. Decomposition of the proton transfer energy profiles for various possible pathways reveals a conflict between the effect of the intrinsic proton affinities of the Schiff base and Asp85, which favors the neutral, product state (i.e., with Asp85 protonated), with the mainly electrostatic interaction between the protein environment with the reacting partners, which favors the ion pair reactant state (i.e., with retinal protonated). The rate-limiting proton-transfer barrier depends both on the relative orientations of the proton donor and acceptor groups and on the pathway followed by the proton; depending on these factors, the barrier may arise from breaking and forming of hydrogen bonds involving the Schiff base, Asp85, Asp212, and water w402, and from nonbonded interactions involving protein groups that respond to the charge rearrangements in the Schiff base region.