Reversible proton transfer dynamics in bacteriorhodopsin

In a previous study of ab initio dynamics, the proton transfer in bacteriorhodopsin from protonated asp96 in the cytoplasmic region toward the deprotonated Schiff base was investigated. A quantum mechanics/molecular mechanics model was constructed from the X-ray structure of bacteriorhodopsin E204Q mutant. In this model, asp96, asp85, and thr89 as well as most of the retinal chromophore and the Schiff base link of lys216 were treated quantum mechanically while the rest of the atoms were treated molecular mechanically. A channel was found in the X-ray structure allowing a water chain to form between the asp96 and Schiff base. In the present study, a chain of four waters from asp96 to the Schiff base N coupled with one branching water supports proton transfer as a concerted event in about 3.5 ps. With both a neutral asp85 and a branched water, the dynamics is now found to be more complicated than observed in the initial study for the transition from the photocycle late M state to the N state. Proton transfer is also observed from the Schiff base back to asp96 demonstrating that there is no effective barrier to proton transfer larger than kT in a strong H-bonded network. The binding of the branched water to the four water chains can dynamically hinder the proton transfer.     

Computational analysis of the proton translocation from Asp96 to schiff base in bacteriorhodopsin

The potential energy change during the M --> N process in bacteriorhodopsin has been evaluated by ab initio quantum chemical and advanced quantum chemical calculations following molecular dynamics (MD) simulations. Many previous experimental studies have suggested that the proton transfer from Asp96 to the Schiff base occurs under the following two conditions: (1) the hydrogen bond between Thr46 and Asp96 breaks and Thr46 is detached from Asp96 and (2) a stable chain of four water molecules spans an area from Asp96 --> Schiff base. In this work, we successfully reproduced the proton-transfer process occurring under these two conditions by molecular dynamics and quantum chemical calculations. The quantum chemical computation revealed that the proton transfer from Asp96 to Shiff base occurs in two-step reactions via an intermediate in which an H(3)O(+) appears around Ala215. The activation energy for the proton transfer in the first reaction was calculated to be 9.7 kcal/mol, which enables fast and efficient proton pump action. Further QM/MM (quantum mechanical/molecular mechanical) and FMO (fragment molecular orbital) calculations revealed that the potential energy change during the proton transfer is tightly regulated by the composition and the geometry of the surrounding amino acid residues of bacteriorhodopsin. Here, we report in detail the Asp96 --> Schiff base proton translocation mechanism of bacteriorhodopsin. Additionally, we discuss the effectiveness of combining quantum chemical calculations with truncated cluster models followed by advanced quantum chemical calculations applied to a whole protein to elucidate its reaction mechanism.     

Key role of active-site water molecules in bacteriorhodopsin proton-transfer reactions

The functional mechanism of the light-driven proton pump protein bacteriorhodopsin depends on the location of water molecules in the active site at various stages of the photocycle and on their roles in the proton-transfer steps. Here, free energy computations indicate that electrostatic interactions favor the presence of a cytoplasmic-side water molecule hydrogen bonding to the retinal Schiff base in the state preceding proton transfer from the retinal Schiff base to Asp85. However, the nonequilibrium nature of the pumping process means that the probability of occupancy of a water molecule in a given site depends both on the free energies of insertion of the water molecule in this and other sites during the preceding photocycle steps and on the kinetic accessibility of these sites on the time scale of the reaction steps. The presence of the cytoplasmic-side water molecule has a dramatic effect on the mechanism of proton transfer: the proton is channeled on the Thr89 side of the retinal, whereas the transfer on the Asp212 side is hindered. Reaction-path simulations and molecular dynamics simulations indicate that the presence of the cytoplasmic-side water molecule permits a low-energy bacteriorhodopsin conformer in which the water molecule bridges the twisted retinal Schiff base and the proton acceptor Asp85. From this low-energy conformer, proton transfer occurs via a concerted mechanism in which the water molecule participates as an intermediate proton carrier.     

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.     

细菌紫红质中氢键相互作用和质子转移机理的分子模拟

用CHARMM程序以细菌紫红质1R84晶体为模型, 模拟了在等温定容条件下细菌紫红质在1 ps过程中的变化, 分析了与质子传递相关的ASP85, ASP212和水分子与视黄醛间氢键的结构变化情况. 考虑到氨基酸残基和席夫碱质子的不同距离, 考察了EC和PC两种结构的变化情况, 探讨了紫红质中质子传递的可能途径. 模拟结果表明1R84中可能的质子连续传递的机理是质子由席夫碱向水传递, 再由水向ASP85传递. 发现Asp212在模拟过程中保持EC结构, 这样可能更有利于顺序质子传递.     

Water as a cofactor in the unidirectional light-driven proton transfer steps in bacteriorhodopsin

Recent evidence for involvement of internal water molecules in the mechanism of bacteriorhodopsin is reviewed. Water O-H stretching vibration bands in the Fourier transform IR difference spectra of the L, M and N intermediates of bacteriorhodopsin were analyzed by photoreactions at cryogenic temperatures. A broad vibrational band in L was shown to be due to formation of a structure of water molecules connecting the Schiff base to the Thr46-Asp96 region. This structure disappears in the M intermediate, suggesting that it is involved in transient stabilization of the L intermediate prior to proton transfer from the Schiff base to Asp85. The interaction of the Schiff base with a water molecule is restored in the N intermediate. We propose that water is a critical mobile component of bacteriorhodopsin, forming organized structures in the transient intermediates during the photocycle and, to a large extent, determining the chemical behavior of these transient states.     

THE TWO CONSECUTIVE M SUBSTATES IN THE PHOTOCYCLE OF BACTERIORHODOPSIN ARE AFFECTED SPECIFICALLY BY THE D85N AND D96N RESIDUE REPLACEMENTS

The photocycle of the proton pump bacteriorhodopsin contains two consecutive intermediates in which the retinal Schiff base is unprotonated; the reaction between these states, termed M1 and M2, was suggested to be the switch in the proton transport which reorients the Schiff base from D85 on the extracellular side to D96 on the cytoplasmic side (Váró and Lanyi, Biochemistry 30, 5016-5022, 1991). At pH 10 the absorption maxima of both M1 and M2 could be determined in the recombinant D96N protein. We find that M1 absorbs at 411 nm as do M1 and M2 in wild-type bacteriorhodopsin, but M2 absorbs at 404 nm. Thus, in M2 but not M1 the unprotonated Schiff base is affected by the D96N residue replacement. The connectivity of the Schiff base to D96 in the detected M2 state, but not in M1, is thereby established. On the other hand, the photostationary state which develops during illumination of D85N bacteriorhodopsin contains an M state corresponding to M1 with an absorption maximum shifted to 400 nm, suggesting that this species in turn is affected by D85. These results are consistent with the suggestion that M1 and M2 are pre-switch and post-switch states, respectively.     

Electric-field effects in dry films of D85N and D85,96N mutant bacteriorhodopsin

In the D85N mutant of the protein bacteriorhodopsin (BR), the Schiff base, by which the retinal chromophore is bound to the protein, exhibits an abnormally low proton affinity (pKa approximately 8.9). Recent experiments on thin films of this protein have shown that this causes the protonation state of the Schiff base, and thus the visible absorption spectrum, to be sensitive to external electric fields. In this paper, we explore the dependence of this effect on parameters such as pH, humidity, and film thickness. The results of these experiments point to the importance of water molecules bound in the acceptor part of the proton channel as sources and donors in field-induced proton-transfer reactions. We describe additional results obtained with the D85,96N mutant, which also exhibits a low Schiff-base pK. The similar behavior of the two mutants under applied electric fields at high pH implies that the residue Asp-96 plays no role in field-induced Schiff-base protonation.     

The proton uptake channel of bacteriorhodopsin as studied by a photoelectrochemical method

A series of the mutant proteins (D96N, D96N/D85N, D115N, L93T, T46V, V49A) where the residues are located at the cytoplasmic domain of bacteriorhodopsin (bR) were studied photoelectrochemically and their photocurrent response characteristics at the electrode/electrolyte interface were compared with those of the wild-type bR. While the wild-type bR of normal proton pumping activity yields symmetrical cathodic (positive) and anodic (negative) responses, corresponding to proton release and proton uptake, respectively, these mutants, with the exception of D115N, showed diminished amplitudes in the negative response. This indicates retardation of proton translocation from the cytoplasmic surface to the retinal Schiff base. The mutation that gave the strongest influence on the negative response was D96N while moderate influence was obtained with L93T, T46V, and V49A. These results suggest that residues other than D96 also participate in the cytoplasmic proton uptake channel, either by interacting with D96 directly or by forming a hydrogen-bonded network with water molecules. The D96N/D85N double mutant yielded little response at neutral pH, but the response was partially recovered by addition of azide, while it was fully recovered in the single mutant D96N. The D115N mutant showed the response profile that closely resembles the wild-type, indicating that D115 is not crucially involved in the event of proton transfer relay at the cytoplasmic region. It was also found that every mutant in this study releases protons prior to uptake at the other membrane surface, as does the wild-type.     

Water molecules in the schiff base region of bacteriorhodopsin

The Schiff base region of bacteriorhodopsin (BR), a light-driven proton pump, contains a pentagonal cluster, being composed of three water molecules and one oxygen each of Asp85 and Asp212. Asp85 and Asp212 are located at similar distances from the retinal Schiff base, whereas the Schiff base proton is transferred only to Asp85 during the pump function. The present FTIR study experimentally established the stretching vibration of water402 hydrating with Asp85 by use of various BR mutants, whose frequency (2171 cm-1 as the O-D stretch) indicates very strong hydrogen bond.     

The mechanism of bacteriorhodopsin functioning. II. Ejection and injection of proton

Transmission of vibrational excitation energy conserved in cis-conformation of retinal to the outlet proton channel is considered from the perspective of quantum theory. A distribution of vibrational excitations in the channel is found; it allowed to calculate the magnitude of the directed drift proton current. The differences between velocities of proton movement in active and passive channels are considered. A transition of retinal from cis- to all -trans- conformation with the subsequent capture of proton by Schiff base out of the inlet channel is described. The lack of proton in this channel, i.e., in the H-bonded chain, is eliminated at the expense of the capture of a proton out of the cytoplasmic water enviroment. The correspondence between theoretically established states and spectroscopically identified forms of bacteriorhodopsin (inter-mediates L, M, N, and O) is proposed.     

PROTON TRANSFER IN THE PRIMARY PROCESS OF VISION

Abstract— Proton transfer was theoretically examined as a possible primary process of vision. The motion of protons in the adiabatic potential of the Schiff base hydrogen bond was investigated in terms of quantum mechanics. The probability of proton transfer from the Schiff base nitrogen (i.e. the unprotonation of Schiff base) was found to increase as the retinal rotated around 11–12. double bond by 90°. The results also suggested that the proton transfer can take place before or during the transition from the excited to ground state (excited state proton transfer). We proposed that such excited state proton transfer is one of the elementary processes in primary visual photochemistry, and this process leads to the unprotonated visual pigment, hyposorhodopsin, which has been experimentally verified as one of the primary photoproducts of rhodopsin. The probability of this process could be comparable to the conventional process leading to the protonated intermediate, bathorhodopsin. The relation of these results with the recent experimental data is discussed.     

A role for internal water molecules in proton affinity changes in the Schiff base and Asp85 for one-way proton transfer in bacteriorhodopsin

Light-induced proton pumping in bacteriorhodospin is carried out through five proton transfer steps. We propose that the proton transfer to Asp85 from the Schiff base in the L-to-M transition is accompanied by the relocation of a water cluster on the cytoplasmic side of the Schiff base from a site close to the Schiff base in L to the Phe219-Thr46 region in M. The water cluster present in L, formed at 170 K, is more rigid than that at room temperature. This may be responsible for blocking the conversion of L to M at 170 K. In the photocycle at room temperature, this water cluster returns to the site close to the Schiff base in N, with a rigid structure similar to that of L at 170 K. The increase in the proton affinity of Asp85, which is a prerequisite for the one-way proton transfer in the M-to-N transition, is suggested to be facilitated by a structural change which disrupts interactions between Asp212 and the Schiff base, and between Asp212 and Arg82. We propose that this liberation of Asp212 is accompanied by a rearrangement of the structure of water molecules between Asp85 and Asp212, stabilizing the protonated Asp85 in M.     

Water and carboxyl group environments in the dehydration blueshift of bacteriorhodopsin

The proton channels of the bacteriorhodopsin (BR) proton pump contain bound water molecules. The channels connect the purple membrane surfaces with the protonated retinal Schiff base at the membrane center. Films of purple membrane equilibrated at low relative humidity display a shift of the 570 nm retinal absorbance maximum to 528 nm, with most of the change occurring below 15% relative humidity. Purple membrane films were dehydrated to defined humidities between about 50 and 4.5% and examined by Fourier transform infrared difference spectroscopy. In spectra of dehydrated-minus-hydrated purple membrane, troughs are observed at 3645 and 3550 cm-1, and peaks are observed at 3665 and 3500 cm-1. We attribute these changes to water dissociation from the proton uptake channel and the resulting changes in hydrogen bonding of water that remains bound. Also, in the carboxylic acid spectral region, a trough was observed at 1742 cm-1 and a peak at 1737 cm-1. The magnitude of the trough to peak difference between 1737 and 1742 cm-1 correlates linearly with the extent of the 528 nm pigment. This suggests that a carboxylic acid group or groups is undergoing a change in environment as a result of dehydration, and that this change is linked to the appearance of the 528 nm pigment. Dehydration difference spectra with BR mutants D96N and D115N show that the 1737-1742 cm-1 change is due to Asp 96 and Asp 115. A possible mechanism is suggested that links dissociation of water in the proton uptake channel to the environmental change at the Schiff base site.     

Excitation of the M Intermediates of Bacteriorhodopsin†

Protein electric response signals (PERS) of the M intermediates of wild‐type bacteriorhodopsin (bR) were recorded. Contrary to earlier findings reporting on a single‐phase response upon excitation of the M intermediates, a kinetic analysis of the signals revealed the existence of three components, the fastest and the slowest ones of negative, while the middle one of positive sign with respect to the normal direction of proton pumping. Based on proton motion indicator experiments and molecular dipole calculations, the components were assigned to proton transfer steps and conformational changes driving the bR molecule back from the M to the ground state upon blue light excitation. The fastest, negative pump component was assigned to the proton transfer from D85 to the Schiff base. The subsequent positive component was attributed to rearrangements in the protein core (in the vicinity of the retinal molecule), triggered by the primary proton transfer process. The slowest component was established to reflect charge rearrangements associated with proton uptake by the protein from the bulk.     

IDENTIFICATION OF THE PROTON ACCEPTOR OF SCHIFF BASE DEPROTONATION IN BACTERIORHODOPSIN: A FOURIER-TRANSFORM-INFRARED STUDY OF THE MUTANT ASP85 → GLU IN ITS NATURAL LIPID ENVIRONMENT

Abstract— In order to assign the proton acceptor for Schiff base deprotonation in bacteriorhodopsin to a specific Asp residue, the photoreaction of the Asp85 → Glu mutant, as expressed in Halobacterium sp . GRB, was investigated by static low-temperature and time-resolved infrared difference spec-troscopy. Measurements were also performed on the mutant protein labeled with [4-¹³C]Asp which allowed discrimination between Asp and Glu residues. 14,15-di¹³C-retinal was incorporated to distinguish amide-II absorbance changes from changes of the ethylenic mode of the chromophore. In agreement with earlier UV-VIS measurements, our data show that from both the 540 and 610 nm species present in a pH-dependent equilibrium, intermediates similar to K and L can be formed. The 14 ms time-resolved spectrum of the 540 nm species shows that a glutamic acid becomes protonated in the M-like intermediate, whereas the comparable difference spectrum of the 610 nm species demonstrates that in the initial state a glutamic acid is already protonated. In conjunction with earlier observations of protonation of an Asp residue in wild-type M, the data provide direct evidence that the proton acceptor in the deprotonation reaction of the Schiff base is Asp85.     

The mechanism of bacteriorhodopsin functioning. I. The light-induced proton throw-over by retinal

The path of proton motion in protein bacteriorhodopsin is constructed: Tyr 43 → Tyr 26 → Asp 212 → Tyr 83 → Tyr 57 → Tyr 79 → Glu 204 → Tyr 64 → Glu 74. From the position of quantum theory the light-induced isomerization of retinal, the process carrying out the proton and throw-over from the inlet H-bond chain (Tyr 43 → Tyr 26 → Lys 216) to the outlet one (Lys 216 → Asp 212 → … → Glu 74) are investigated. The kinetics of the transition of the excitation from the π-electron subsystem to the vibrational one is studied; the function of the nonequilibrial excitation distribution of the collective intramolecular retinal vibrations is found. The mechanism of deprotonation of Schiff base (transmission of proton to the outlet channel) is considered. The correspondence between the theoretically estimated states and spectroscopically identified forms of bacteriorhodosin (initial form bR, intermediates J, K) is proposed.     

Role of hydrogen-bond network in energy storage of bacteriorhodopsin's light-driven proton pump revealed by ab initio normal-mode analysis

Vibrational modes of the hydrogen-bond network in the binding site of bacteriorhodopsin (bR), a protein in halobacteria functioning as a light-driven proton pump, were investigated by an ab initio quantum mechanical/molecular mechanical (QM/MM) method. Normal-mode analysis calculations for O-D and N-D stretching modes of internal water molecules and the Schiff base of the retinal chromophore in the early intermediate state, K, reproduced well experimentally observed vibrational spectra. Supported by agreement with observed spectra, the QM/MM calculation suggests that weakened hydrogen bonds upon photoisomerization of the chromophore are an important means of energy storage in bR.     

RESONANCE RAMAN SPECTRA OF BACTERIORHODOPSIN MUTANTS WITH SUBSTITUTIONS AT ASP-85, ASP-96, AND ARG-82

Detergent solubilized bacteriorhodopsin (BR) proteins which contain alterations made by site-directed mutagenesis (Asp-96----Asn, D96N; Asp-85----Asn, D85N; and Arg-82----Gln, R82Q) have been studied with resonance Raman spectroscopy. Raman spectra of the light-adapted (BRLA) and M species in D96N are identical to those of native BR, indicating that this residue is not located near the chromophore. The BRLA states of D85N and especially R82Q contain more of the 13-cis, C = N syn (BR555) species under ambient illumination compared to solubilized native BR. Replacement of Asp-85 with Asn causes a 25 nm red-shift of the absorption maximum and a frequency decrease in both the ethylenic (-7 cm-1) and the Schiff base C = NH+ (-3 cm-1) stretching modes of BRLA. These changes indicate that Asp-85 is located close to the protonated retinal Schiff base. The BRLA spectrum of R82Q exhibits a slight perturbation of the C = NH+ band, but its M spectrum is unperturbed. The Raman spectra and the absorption properties of D85N and R82Q suggest that the protein counterion environment involves the residues Asp-85-, Arg-82+ and presumably Asp-212-. These data are consistent with a model where the strength of the protein-chromophore interaction and hence the absorption maximum depends on the overall charge of the Schiff base counterion environment.     

Influence of the membrane potential on the protonation of bacteriorhodopsin: insights from electrostatic calculations into the regulation of proton pumping

Proton binding and release are elementary steps for the transfer of protons within proteins, which is a process that is crucial in biochemical catalysis and biological energy transduction. Local electric fields in proteins affect the proton binding energy compared to aqueous solution. In membrane proteins, also the membrane potential affects the local electrostatics and can thus be crucial for protein function. In this paper, we introduce a procedure to calculate the protonation probability of titratable sites of a membrane protein in the presence of a membrane potential. In the framework of continuum electrostatics, we use a modified Poisson-Boltzmann equation to include the influence of the membrane potential. Our method considers that in a transmembrane protein each titratable site is accessible for protons from only one side of the membrane depending on the hydrogen bond pattern of the protein. We show that the protonation of sites receiving their protons from different sides of the membrane is differently influenced by the membrane potential. In addition, the effect of the membrane potential is combined with the effect of the pH gradient resulting from proton pumping. Our method is applied to bacteriorhodopsin, a light-activated proton pump. We find that the proton pumping of this protein might be regulated by Asp115, a conserved residue for which no function has been identified yet. According to our calculations, the interaction of Asp115 with Asp85 leads to the protonation of the latter if the pH gradient or the membrane potential becomes too large. Since Asp85 is the primary proton acceptor in the photocycle, bacteriorhodopsin molecules in which Asp85 is protonated cannot pump protons. Furthermore, we estimate how the membrane potential affects the energetics of the individual proton-transfer reactions of the photocycle. Most reactions, except the initial proton transfer from the Schiff base to Asp85, are influenced. Our calculations give new insights into the mechanism with which bacteriorhodopsin senses the membrane potential and the pH gradient and how the proton pumping is regulated by these parameters.