1A combination of molecular dynamics and quantum chemistry techniques have been employed to study the electronic excitation and conformational potential surface of retinal in the binding site of bacteriorhodopsin (bR). The CASSCF(6,9)/6-31G level of ab initio calculations (within Gaussian92) has been used for the treatment of both the ground (S0) and excited (S1) states of retinal. Charges of all atoms in the protein are represented by spherical Gaussians and explicitly included in the electronic Hamiltonian of retinal. Spectral properties have been analyzed for the native bR pigment as well as for its D85N mutant.The calculated relative shift in the absorption maxima between the two pigments is in better agreement with experiment than the computed absolute parameters of the absorption line shapes. The dark adaptation processes in bacteriorhodopsin (which involves rotation around the 13-14 and the 15-N retinal double bonds) has been modelled by following the pre-defined reaction coordinate. Our simulations support the notion that the isomerization process is catalyzed by the protonation of an aspartic acid (Asp85) side group of bacteriorhodopsin.
2
IntroductionBacteriorhodopsin (bR) spans the cell membrane of Halobacterium halobium and functions as a light-driven proton pump. bR contains seven α-helices which enclose the prosthetic group, all-trans retinal, bound via a protonated Schiff base linkage to Lys-216. Figure 1a shows the chemical structure of retinal and its conventional numbering scheme. The bR structure is presented in Fig. 2a. Retinal absorbs light and undergoes a photoisomerization process; the thermal reversal of this reaction is coupled to transfer of a proton from the cytoplasmic side (top in Fig. 2a) to the extracellular side (bottom in Fig. 2a) of the protein.Recent reviews that discuss the structure and function of bR are [1,2,3,4,5,6,7]. Figure 1 here Bacteriorhodopsin accomplishes its function through a cyclic process initiated by absorption of a photon. This photon triggers an isomerization of retinal, which proceeds then through several intermediate states identified by their absorption spectra. An accepted kinetic scheme for this cycle is an unbranched series of intermediates, shown in Fig. 3. Photoisomerization occurs in the bR 568 → J 625 transition. During the L 550 → M 412 transition the Schiff base proton is transferred to Asp-85 and, subsequently, to the outside of the cell [8,9,10,11,12,13]. During the M 412 → N 520 transition, a proton is transferred to the Schiff base from Asp-96, which then takes up a proton from the cytoplasmic environment [13]. The reaction cycle is completed as the protein returns to bR 568 via the O 640 intermediate. Figure 2 here When bR is allowed to equilibrate in the dark, it converts within an hour to a 2:1 mixture containing 13-cis retinal and all-trans retinal bR [14]. The protein containing the 13-cis retinal isomer of bRis referred to as the dark adapted (DA) pigment of bR, bR 548 , which absorbs at 548 nm. bR 548 contains, actually, retinal in a 1...