MOLECULAR-DYNAMICS STUDY OF THE EARLY INTERMEDIATES IN THE BACTERIORHODOPSIN PHOTOCYCLE

The early stages of the bacteriorhodopsin photocycle, including the J( 625), K-590, and L(550) intermediates and the role of water molecules within the protein interior, are studied by means of molecular dynamic s simulations. Our calculations examine two model's for the excited st ate potential surface governing the observed all-trans --> 13-cis phot oisomerization: one surface hindering a C-14-C-15 single-bond corotati on and the other surface allowing such corotation. The investigations use as a starting structure a model of bacteriorhodopsin based on elec tron-microscopy studies and subsequent molecular dynamics refinement. The following questions are addressed: How does the binding site guide retinal's photoisomerization? How does the photoisomerization depend on features of the excited state potential surface? Can one recognize a J(625) intermediate? How does water participate in the early part of the pump cycle? How is the initial photoreaction affected by a loweri ng of temperature? To model the quantum yield, i.e., the dependence of the dynamics on initial conditions, 50 separate isomerization trials are completed for each potential surface, at both 300 and 77 K, the tr ials being distinguished by different initial, random velocity distrib utions. From these trials emerge, besides all-transretinal, three diff erent photoproducts as candidates for the K-590 intermediate: (1) 13-c is-retinal, with the Schiff base proton oriented toward Asp-96; (2) 13 -cis-retinal, highly twisted about the C-6-C-7 bond, with the Schiff b ase proton oriented perpendicular to the membrane normal; (3) 13,14-di cis-retinal, with the Schiff base proton oriented toward the extracell ular side. Two candidates for the K-590 intermediate, case 2 and case 3 above, were subjected to simulated annealing to determine correspond ing L(550) structures. We suggest that photoproduct 2 above most likel y represents the true K-590 intermediate. Water molecules near the Sch iff base binding site are found to play a crucial role in stabilizing the K-590 state and in establishing a pathway for proton transfer to A sp-85.