The structure of bacteriorhodopsin based on electron microscopy (EM) studies, as provided in Henderson et al. (1990), is refined using molecular dynamics simulations. The work is based on a previously refined and simulated structure which had added the interhelical loops to the EM model of bR. The present study applies an all-atom description to this structure and constraints to the original Henderson model, albeit with helix D shifted. Sixteen waters are then added to the protein, six in the retinal Schiff base region, four in the retinal-Asp-96 interstitial space, and six near the extracellular side. The root mean square deviation between the resulting structure and the Henderson et al. (1990) model measures only 1.8 A. Further simulations of retinal analogues for substitutions at the 2- and 4-positions of retinal and an analogue without a beta-ionone ring agree well with observed spectra. The resulting structure is characterized in view of bacteriorhodopsin's function; key features are (1) a retinal Schiff base-counterion complex which is formed by a hydrogen bridge network involving six water molecules, Asp-85, Asp-212, Tyr-185, Tyr-57, Arg-82, and Thr-89, and which exhibits Schiff base nitrogen-Asp-85 and -Asp-212 distances of 6 and 4.6 A; (2) retinal assumes a corkscrew twist as one views retinal along its backbone; and (3) a deviation from the usual alpha-helical structure of the cytoplasmic side of helix G.
The primary all-trans --> 13-cis photoisomerization of retinal in bacteriorhodopsin has been investigated by means of quantum chemical and combined classical/quantum mechanical simulations employing the density matrix evolution method. Ab initio calculations on an analog of a protonated Schiff base of retinal in vacuo reveal two excited states S1 and S2, the potential surfaces of which intersect along the reaction coordinate through an avoided crossing, and then exhibit a second, weakly avoided, crossing or a conical intersection with the ground state surface. The dynamics governed by the three potential surfaces, scaled to match the in situ level spacings and represented through analytical functions, are described by a combined classical/quantum mechanical simulation. For a choice of nonadiabatic coupling constants close to the quantum chemistry calculation results, the simulations reproduce the observed photoisomerization quantum yield and predict the time needed to pass the avoided crossing region between S1 and S2 states at tau1 = 330 fs and the S1 --> ground state crossing at tau2 = 460 fs after light absorption. The first crossing follows after a 30 degrees torsion on a flat S1 surface, and the second crossing follows after a rapid torsion by a further 60 degrees. tau1 matches the observed fluorescence lifetime of S1. Adjusting the three energy levels to the spectral shift of D85N and D212N mutants of bacteriorhodospin changes the crossing region of S1 and S2 and leads to an increase in tau1 by factors 17 and 10, respectively, in qualitative agreement with the observed increase in fluorescent lifetimes.
The availability of the structure of bacteriorhodopsin from electron microscopy studies has opened up the possibility of exploring the proton pump mechanism of this protein by means of molecular dynamics simulations. In this review we summarize earlier theoretical investigations of the photocycle of bacteriorhodopsin including relevant quantum chemistry studies of retinal, structure refinement, molecular dynamics simulations, and evaluation of pK a values. We then review a series of recent modeling efforts which refined the structure of bacteriorhodopsin adding internal water, and which studied the nature of the J intermediate and the likely geometry of the K S90 and L sso intermediates (strongly distorted 13-cis) as well as the sequence of retinal geometry and protein conformational transitions which are conventionally summarized as the M 4 12 intermediate. We also review simulations of the photocycle of light-adapted bacteriorhodopsin at T=77 K and of the photocycle of dark-adapted bacteriorhodopsin, both cycles differing from the conventional photocycle through a nonfunctional (pure 13-cis) retinal geometry of the corresponding K S90 and L sso states. The simulations demonstrate a potentially critical role of water and of minute reorientations of retinal's Schiff base nitrogen in controlling proton pumping in bR s6s ; the simulations also indicate the existence of heterogeneous photocycles. The results exemplify the important role of molecular dynamics simulations in extending investigations on bacteriorhodopsin to a level of detail which is presently beyond experimental resolution, but which needs to be known to resolve the pump mechanism of bacteriorhodopsin. Finally, we outline the major existing challenges in the field of bacteriorhodopsin modeling. 'Permanent address:
Molecular dynamics simulations and quantum chemistry calculations have been combined to describe the dark adaptation in bacteriorhodopsin (bR). The process involves the reversible thermally activated transformation of retinal from an all-trans to a 13-cis,15-syn configuration. The potential surface governing the thermal isomerization of retinal around two (13−14, 15−N) double bonds has been determined for representative protein configurations taken from molecular dynamics trajectories. CASSCF(8,8)/6-31G level ab initio calculations (within Gaussian94) were carried out for this purpose. The charge distributions of all atoms in the protein are represented by partial point charges and explicitly included in the electronic Hamiltonian. Placement of retinal into bR is found to reduce the calculated isomerization barrier. Thermal fluctuations of the protein lead to a further effective reduction of this barrier. The isomerization process is shown to be catalyzed by the protonation of an aspartic acid (Asp85) side group of bacteriorhodopsin.
The structure and the photocycle of bacteriorhodopsin (bR) containing 13-cis,15-syn retinal, so-called bR548, has been studied by means of molecular dynamics simulations performed on the complete protein. The simulated structure of bR548 was obtained through isomerization of in situ retinal around both its C13-C14 and its C15-N bond starting from the simulated structure of bR568 described previously, containing all-trans,15-anti retinal. After a 50-ps equilibration, the resulting structure of bR548 was examined by replacing retinal by analogues with modified beta-ionone rings and comparing with respective observations. The photocycle of bR548 was simulated by inducing a rapid 13-cis,15-anti-->all-trans,15-syn isomerization through a 1-ps application of a potential that destabilizes the 13-cis isomer. The simulation resulted in structures consistent with the J, K, and L intermediates observed in the photocycle of bR548. The results offer an explanation of why an unprotonated retinal Schiff base intermediate, i.e., an M state, is not formed in the bR548 photocycle. The Schiff base nitrogen after photoisomerization of bR548 points to the intracellular rather than to the extracellular site. The simulations suggest also that leakage from the bR548 to the bR568 cycle arises due to an initial 13-cis,15-anti-->all-trans,15-anti photoisomerization.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.