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

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

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.     

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.     

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.     

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.     

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.     

Dynamics of proton transfer in bacteriorhodopsin

Proton transfer in bacteriorhodopsin from the cytoplasm to the extracellular side is initiated from protonated asp96 in the cytoplasmic region toward the deprotonated Schiff base. This occurs in the transition from the photocycle late M state to the N state. To investigate this proton-transfer process, a quantum mechanics/molecular mechanics (QM/MM) model is constructed from the bacteriorhodopsin E204Q mutant crystal structure. Three residues, asp96, asp85, and thr89, as well as most of the retinal chromophore and the Schiff base link of lys216 are treated quantum mechanically and connected to the remaining classical protein through linker atom hydrogens. Structural transformation in the M state results in the formation of a water channel between the Schiff base and asp96. Since a part of this channel is lined with hydrophobic residues, there has been a question on the mechanism of proton transfer in a hydrophobic channel. Ab initio dynamics using the CHARMM/GAMESS methodology is used to simulate the transfer of the proton through a partially hydrophobic channel. Once sufficient water molecules are added to the channel to allow the formation of a single chain of waters from asp96 to the Schiff base, the transfer occurs as a fast (less than a picosecond) concerted event irrespective of the protonation state of asp85. Dynamic transfer of the proton from asp96 to the nearest water initiates the organization of a strongly bonded water chain conducive to the transfer of the proton to the Schiff base nitrogen.     

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.     

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.     

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.     

Parallel proton transfer pathways in aqueous acid-base reactions

We study the mechanism of proton transfer (PT) between the photoacid 8-hydroxy-1,3, 6-pyrenetrisulfonic acid (HPTS) and the base chloroacetate in aqueous solution. We investigate both proton and deuteron transfer reactions in solutions with base concentrations ranging from 0.25 M to 4 M. Using femtosecond midinfrared spectroscopy, we probe the vibrational responses of HPTS, its conjugate photobase, the hydrated proton/deuteron, and chloroacetate. The measurement of these four resonances allows us to follow the sequence of proton departure from the acid, its uptake by the water solvent, and its arrival at the base. In recent studies it was shown that proton transfer to carboxylate bases proceeds via Grotthuss conduction through a water wire connecting the acid and the base [Mohammed et al., Science 310, 83 (2005);Agnew. Chem. Int. Ed. 46, 1458 (2007);Siwick and Bakker, J. Am. Chem. Soc. 129, 13412 (2007); J. Phys. Chem. B 112, 378 (2008)]. Here we show that, for the weaker base chloroacetate, an alternative channel for proton transfer arises. In this channel the proton is first transferred to the water solvent and only later taken up from the water by the base. We study the base concentration dependence of the two competing channels.     

The lifetimes of Pharaonis phoborhodopsin signaling states depend on the rates of proton transfers--effects of hydrostatic pressure and stopped flow experiments

Pharaonis phoborhodopsin (ppR), a negative phototaxis receptor of Natronomonas pharaonis, undergoes photocycle similar to the light-driven proton pump bacteriorhodopsin (BR), but the turnover rate is much slower due to much longer lifetimes of the M and O intermediates. The M decay was shown to become as fast as it is in BR in the L40T/F86D mutant. We examined the effects of hydrostatic pressure on the decay of these intermediates. For BR, pressure decelerated M decay but slightly affected O decay. In contrast, with ppR and with its L40T/F86D mutant, pressure slightly affected M decay but accelerated O decay. Clearly, the pressure-dependent factors for M and O decay are different in BR and ppR. In order to examine the deprotonation of Asp75 in unphotolyzed ppR we performed stopped flow experiments. The pH jump-induced deprotonation of Asp75 occurred with 60 ms, which is at least 20 times slower than deprotonation of the equivalent Asp85 in BR and about 10-fold faster than the O decay of ppR. These data suggest that proton transfer is slowed not only in the cytoplasmic channel but also in the extracellular channel of ppR and that the light-induced structural changes in the O intermediate of ppR additionally decrease this rate.     

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.     

Proton transfer reaction of a new orthohydroxy Schiff base at room temperature and 77 K

One new orthohydroxy Schiff base, 2-(N-benzyl-alpha-iminoethyl)naphthol (BEIN) has been synthesized. The proton transfer reaction of BEIN has been investigated by means of absorption, steady state and time resolved fluorescence spectroscopy in different solvents at room temperature and 77K. The behavior of BEIN in ethanol and water, has been studied in neutral, acidic and basic conditions. Excited state intramolecular proton transfer (ESIPT) is evidenced by a large Stokes shifted ( approximately 11,000 cm-1) fluorescence in solid crystalline media at room temperature. We present the observation of phosphorescence both in non-polar and protic solvents at 77K. The observed decay dynamics of the phosphorescence and delayed fluorescence indicates that the triplet state can be attributed to the cis-keto form. The molecular structures are determined by B3LYP/6-31G** calculation. From theoretical study it is suggested that the strengthening of hydrogen bond result from the steric repulsion of the phenyl ring. The presence of benzene ring increases the proton transfer barrier in case of BEIN compared to previously studied 7-ethylsalicylidenebenzylamine (ESBA).     

TIME-RESOLVED FOURIER TRANSFORM INFRARED SPECTROSCOPY OF THE BACTERIORHODOPSIN MUTANT TYR-185→E: ASP-96 REPROTONATES DURING O FORMATION; ASP-85 AND ASP-212 DEPROTONATE DURING O DECAY

The protonation state of key aspartic acid residues in the O intermediate of bacteriorhodopsin (bR) has been investigated by time-resolved Fourier transform infrared (FTIR) difference spectroscopy and site-directed mutagenesis. In an earlier study (Bousché et al., J. Biol Chem. 266, 11063-11067, 1991) we found that Asp-96 undergoes a deprotonation during the M-->N transition, confirming its role as a proton donor in the reprotonation pathway leading from the cytoplasm to the Schiff base. In addition, both Asp-85 and Asp-212, which protonate upon formation of the M intermediate, remain protonated in the N intermediate. In this study, we have utilized the mutant Tyr-185-->Phe (Y185F), which at high pH and salt concentrations exhibits a photocycle similar to wild type bR but has a much slower decay of the O intermediate. Y185F was expressed in native Halobacterium halobium and isolated as intact purple membrane fragments. Time-resolved FTIR difference spectra and visible difference spectra of this mutant were measured from hydrated multilayer films. A normal N intermediate in the photocycle of Y185F was identified on the basis of characteristic chromophore and protein vibrational bands. As N decays, bands characteristic of the all-trans O chromophore appear in the time-resolved FTIR difference spectra in the same time range as the appearance of a red-shifted photocycle intermediate absorbing near 640 nm. Based on our previous assignment of the carboxyl stretch bands to the four membrane embedded Asp groups: Asp-85, Asp-96, Asp-115 and Asp-212, we conclude that during O formation: (i) Asp-96 undergoes reprotonation. (ii) Asp-85 may undergo a small change in environment but remains protonated. (iii) Asp-212 remains partially protonated. In addition, reisomerization of the chromophore during the N-->O transition is accompanied by a major reversal of protein conformational changes which occurred during the earlier steps in the photocycle. These results are discussed in terms of a proposed mechanism for proton transport.     

Excited state intramolecular proton transfer in a new o-hydroxy Schiff base in non polar solvents at room temperature and 77 K

A new orthohydroxy Schiff base, 7-ethylsalicylidenebenzylamine (ESBA) has been synthesised. The excited state intramolecular proton transfer (ESIPT) processes have been investigated by means of absorption, emission and nanosecond spectroscopy at room temperature and at 77 K in non polar solvents. The ESIPT is evidenced by a large Stokes shifted emission (11 000 cm⁻¹) only at 77 K. From fluorescence and excitation spectra it is suggested that at least three different species are present in the excited state at room temperature. Our theoretical calculation at AM1 level confirm the cis-isomer to be the only viable form in the ground state.     

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 uptake of rhodobacter capsulatus reaction center mutants modified in the primary quinone environment

Flash-induced absorbance spectroscopy was used to analyze the proton uptake and electron transfer properties of photosynthetic reaction centers (RC) of Rhodobacter capsulatus that have been genetically modified near the primary quinone electron acceptor (Q(A)). M246Ala and M247Ala, which are symmetry-related to the positions of two acidic groups, L212Glu and L213Asp, in the secondary quinone electron acceptor (QB) protein environment, have been mutated to Glu and Asp, respectively. The pH dependence of the stoichiometry of proton uptake upon formation of the P+Q(A)- (H+/P+Q(A)-) and PQ(A) (H+/Q(A)-) (P is the primary electron donor, a noncovalently linked bacteriochlorophyll dimer) states have been measured in the M246Ala --> Glu and the M247Ala --> Asp mutant RC, in the M246Ala-M247Ala --> Glu-Asp double mutant and in the wild type (WT). Our results show that the introduction of an acidic group (Glu or Asp) in the QA protein region induces notable additional proton uptake over a large pH region (approximately 6-9), which reflects a delocalized response of the protein to the formation of Q(A)-. This may indicate the existence of a widely spread proton reservoir in the cytoplasmic region of the protein. Interestingly, the pH titration curves of the proton release caused by the formation of P+ (H+/P+: difference between H+/P+Q(A)- and H+/PQ(A)- curves) are nearly superimposable in the WT and the M246Ala --> Glu mutant RC, but substantial additional proton release is detected between pH 7 and 9 in the M247Ala --> Asp mutant RC. This effect can be accounted for by an increased proton release by the P+ environment in the M247Ala --> Asp mutant. The M247Ala --> Asp mutation reveals the existence of an energetic and conformational coupling between donor and acceptor sides of the RC at a distance of nearly 30A.     

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.