Free-energy simulations of the retinal cis --> trans isomerization in bacteriorhodopsin.

Free-energy profiles for ground-state cis --> trans isomerization of retinal in vacuum, in solution, and in the protein bacteriorhodopsin are calculated using free-energy simulations. The free-energy barriers in the protein were 9 kcal/mol for ionized Asp85 and 14 kcal/mol for neutral Asp85, significantly lower than those found in solution (18 kcal/mol) or vacuum (19 kcal/mol). Therefore, bacteriorhodopsin can be said to act as a catalyst in the isomerization. The barrier in the protein is due mainly to stabilization of the transition state through favorable nonbonded interactions with the protein part of the system, with internal strain and interactions with solvent playing minor roles. The protonated Asp85 simulation models the behavior of the system in the N --> O transition. Our calculated 14 kcal/mol barrier and 4-ms relaxation time for this process are in excellent agreement with experimentally measured values of 12 kcal/mol and 5 ms, respectively. The ionized Asp85 simulation models two hypothetical processes: the N --> O transition with a proton removed from Asp85 and the initial BR568 --> L transition on the ground-state energy surface. The cis-trans isomerization barrier in this system is 9 kcal/mol, the lowest of all the studied cases. The presence of the charged carboxylate group in the ionized Asp85 system leads to strong stabilization of the transition state by interactions with the surroundings and changes the distance between Asp85 and the Schiff base proton compared to the corresponding distance in the neutral Asp85 system. This suggests that the protonation of Asp85 plays an important role in regulating access to the Schiff base proton. For both Asp85 ionization states the calculated cis-trans free-energy difference was close to 0, indicating that the protein can accommodate both retinal isomers equally well. The computed negligible difference between the N and O free-energy levels is in accord with experimental data.