Suppression of the back proton-transfer from Asp85 to the retinal Schiff base in bacteriorhodopsin: A theoretical analysis of structural elements

Bacteriorhodopsin proton-pumping mechanism: Successes and challenges in computational approaches

Bacteriorhodopsin (bR) is perhaps the best-studied proton pump. Over about four decades, research on this fascinating photocyclic light-driven protein inspired the development of key experimental and computational methodologies that are now widely used in membrane protein studies. We review here failures and successes in computational approaches that have been applied to study the bR proton-transfer steps. Conflict between experimental results pertaining to the proton transfer mechanisms in the early photocycle intermediates was resolved by detailed quantum mechanical/molecular mechanical computation, the results of which were confirmed more than a decade later. Key to this approach was the realization that, to understand how the pump works and achieves directional transfer of protons, the individual reaction steps-proton transfer and reorganization of the internal hydrogen-bond network-needed to be considered within the context of the energy landscape of the complete reaction cycle.

Graphs of protein-water hydrogen bond networks to dissect structural movies of ion-transfer microbial rhodopsins

Microbial rhodopsins are membrane proteins that use the energy absorbed by the covalently bound retinal chromophore to initiate reaction cycles resulting in ion transport or signal transduction. Thousands of distinct microbial rhodopsins are known and, for many rhodopsins, three-dimensional structures have been solved with structural biology, including as entire sets of structures solved with serial femtosecond crystallography. This sets the stage for comprehensive studies of large datasets of static protein structures to dissect structural elements that provide functional specificity to the various microbial rhodopsins. A challenge, however, is how to analyze efficiently intra-molecular interactions based on large datasets of static protein structures. Our perspective discusses the usefulness of graph-based approaches to dissect structural movies of microbial rhodopsins solved with time-resolved crystallography.

Hydroxide Ion Mechanism for Long-Range Proton Pumping in the Third Proton Transfer of Bacteriorhodopsin

In bacteriorhodopsin, representative light-driven proton pump, five proton transfers yield vectorial active proton translocation, resulting in a proton gradient in microbes. Third proton transfer occurs from Asp96 to the Schiff base on the photocycle, which is expected to be a long-range proton transfer via the Grotthuss mechanism through internal water molecules. Here, large-scale quantum molecular dynamics simulations are performed for the third proton transfer, where all the atoms (~50000 atoms) are treated quantum-mechanically. The simulations demonstrate that two reaction paths exist along the water wire, namely, via hydronium and via hydroxide ions. The free energy analysis confirms that the path via hydroxide ions is considerably favorable and consistent with the observed lifetime of the transient water wire. Therefore, the proposed hydroxide ion mechanism, as in the first proton transfer, is responsible for the third long-range proton transfer.

Light-driven proton transfers and proton transport by microbial rhodopsins - A biophysical perspective

In the last twenty years, our understanding of the rules and mechanisms for the outward light-driven proton transport (and underlying proton transfers) by microbial rhodopsins has been changing dramatically. It transitioned from a very detailed atomic-level understanding of proton transport by bacteriorhodopsin, the prototypical proton pump, to a confounding variety of sequence motifs, mechanisms, directions, and modes of transport in its newly found homologs. In this review, we will summarize and discuss experimental data obtained on new microbial rhodopsin variants, highlighting their contribution to the refinement and generalization of the ideas crystallized in the previous century. In particular, we will focus on the proton transport (and transfers) vectoriality and their structural determinants, which, in many cases, remain unidentified.

Structural Factors Determining the Absorption Spectrum of Channelrhodopsins: A Case Study of the Chimera C1C2

Channelrhodopsins are photosensitive proteins that trigger flagella motion in single-cell algae and have been successfully utilized in optogenetic applications. In optogenetics, light is used to activate neural cells in living organisms, which can be achieved by exploiting the ion channel signaling of channelrhodopsins. Tailoring channelrhodopsins for such applications includes the tuning of the absorption maximum. In order to establish rational design and to obtain a desired spectral shift, a basic understanding of the absorption spectrum is required. We have studied the chimera C1C2 as a representative of this protein family and the first member with an available crystal structure. For this purpose, we sampled the conformations of C1C2 using quantum mechanical/molecular mechanical molecular dynamics and subjected the resulting snapshots of the trajectory to excitation energy calculations using ADC(2) and simplified time-dependent density functional theory. In contrast to previous reports, we found that different hydrogen-bonding networks-involving the retinal protonated Schiff base, the putative counterions E162 and D292, and water molecules-had only a small impact on the absorption spectrum. However, in the case of deprotonated E162, increasing the distance to the Schiff base hydrogen-bonding partner led to a systematic blue shift. The β-ionone ring rotation was identified as another important contributor. Yet the most important factors were found to be the bond length alternation and bond order alternation that were linearly correlated to the absorption maximum by up to 62 and 82%, respectively. We ascribe this novel insight into the structural basis of the absorption spectrum to our enhanced protein setup that includes membrane embedding as well as long and extensive sampling.

Protein Motifs for Proton Transfers That Build the Transmembrane Proton Gradient

Biological membranes are barriers to polar molecules, so membrane embedded proteins control the transfers between cellular compartments. Protein controlled transport moves substrates and activates cellular signaling cascades. In addition, the electrochemical gradient across mitochondrial, bacterial and chloroplast membranes, is a key source of stored cellular energy. This is generated by electron, proton and ion transfers through proteins. The gradient is used to fuel ATP synthesis and to drive active transport. Here the mechanisms by which protons move into the buried active sites of Photosystem II (PSII), bacterial RCs (bRCs) and through the proton pumps, Bacteriorhodopsin (bR), Complex I and Cytochrome c oxidase (CcO), are reviewed. These proteins all use water filled proton transfer paths. The proton pumps, that move protons uphill from low to high concentration compartments, also utilize Proton Loading Sites (PLS), that transiently load and unload protons and gates, which block backflow of protons. PLS and gates should be synchronized so PLS proton affinity is high when the gate opens to the side with few protons and low when the path is open to the high concentration side. Proton transfer paths in the proteins we describe have different design features. Linear paths are seen with a unique entry and exit and a relatively straight path between them. Alternatively, paths can be complex with a tangle of possible routes. Likewise, PLS can be a single residue that changes protonation state or a cluster of residues with multiple charge and tautomer states.

Dynamic water bridging and proton transfer at a surface carboxylate cluster of photosystem II

A hydrogen-bond cluster at a negatively-charged protein interface with a bound protein and long-lived waters might be a proton storage site.

Mechanism by which water and protein electrostatic interactions control proton transfer at the active site of channelrhodopsin

Channelrhodopsins are light-sensitive ion channels whose reaction cycles involve conformation-coupled transfer of protons. Understanding how channelrhodopsins work is important for applications in optogenetics, where light activation of these proteins triggers changes in the transmembrane potential across excitable membranes. A fundamental open question is how the protein environment ensures that unproductive proton transfer from the retinal Schiff base to the nearby carboxylate counterion is avoided in the resting state of the channel. To address this question, we performed combined quantum mechanical/molecular mechanical proton transfer calculations with explicit treatment of the surrounding lipid membrane. The free energy profiles computed for proton transfer to the counterion, either via a direct jump or mediated by a water molecule, demonstrate that, when retinal is all-trans, water and protein electrostatic interactions largely favour the protonated retinal Schiff base state. We identified a conserved lysine group as an essential structural element for the proton transfer energetics in channelrhodopsins.

Catalysis of Ground State cis[Formula: see text] trans Isomerization of Bacteriorhodopsin's Retinal Chromophore by a Hydrogen-Bond Network

For the photocycle of the membrane protein bacteriorhodopsin to proceed efficiently, the thermal 13-cis to all-trans back-isomerization of the retinal chromophore must return the protein to its resting state on a time-scale of milliseconds. Here, we report on quantum mechanical/molecular mechanical energy calculations examining the structural and energetic determinants of the retinal cis-trans isomerization in the protein environment. The results suggest that a hydrogen-bonded network consisting of the retinal Schiff base, active site amino acid residues, and water molecules can stabilize the twisted retinal, thus reducing the intrinsic energy cost of the cis-trans thermal isomerization barrier.

Density functional tight binding: values of semi-empirical methods in an ab initio era

Semi-empirical (SE) methods are derived from Hartree-Fock (HF) or Density Functional Theory (DFT) by neglect and approximation of electronic integrals. Thereby, parameters are introduced which have to be determined from reference calculations and/or by fitting to available experimental data. This leads to computational methods that are about 2-3 orders of magnitude faster than the standard HF/DFT methods using medium sized basis sets while being about 3 orders of magnitude slower than empirical force field methods (Molecular Mechanics: MM). Therefore, SE methods are most appropriate for a specific range of applications. These include the study of systems that contain a large number of atoms and therefore being too large for ab initio or DFT methods and also problems where dynamic or entropic effects are particularly important. In the latter case, the errors made by considering a very limited number of molecular structures or neglecting entropic contributions can be much larger than the accuracy lost due to the use of SE methods. Another area where SE methods are attractive concerns the analysis of systems for which reliable MM models are not readily available. Therefore, even in an era when rapid progress is being made in ab initio methods, there is considerable interest in further developing SE methods. We illustrate this point by focusing on the discussion of recent development and application of the Density Functional Tight Binding method.

Schiff base switch II precedes the retinal thermal isomerization in the photocycle of bacteriorhodopsin

In bacteriorhodopsin, the order of molecular events that control the cytoplasmic or extracellular accessibility of the Schiff bases (SB) are not well understood. We use molecular dynamics simulations to study a process involved in the second accessibility switch of SB that occurs after its reprotonation in the N intermediate of the photocycle. We find that once protonated, the SB C15 = NZ bond switches from a cytoplasmic facing (13-cis, 15-anti) configuration to an extracellular facing (13-cis, 15-syn) configuration on the pico to nanosecond timescale. Significantly, rotation about the retinal’s C13 = C14 double bond is not observed. The dynamics of the isomeric state transitions of the protonated SB are strongly influenced by the surrounding charges and dielectric effects of other buried ions, particularly D96 and D212. Our simulations indicate that the thermal isomerization of retinal from 13-cis back to all-trans likely occurs independently from and after the SB C15 = NZ rotation in the N-to-O transition.

Structural, proton-transfer, thermodynamic and nonlinear optical studies of (E)-2-((2-hydroxyphenyl)iminiomethyl)phenolate

Recently, the study of imine-bridged organics is interested in proton-transfer and photo-responsive material fields. Herein, we make a investigation on the structural, thermodynamic and nonlinear optical properties of (E)-2-((2-hydroxyphenyl)iminiomethyl)phenolate (HPIMP). The structural varieties of the studied compound are characterized by the X-ray single crystal diffraction and vibrational spectral techniques, as well as the vibrational spectral bands are precisely ascribed to the studied structure with the aid of DFT theoretical calculations. The experimental results of the FT-IR and X-ray measurements supply good proofs to reveal the proton-transfer procedures of HPIMP, and exhibit that the studied compound is a good proton-transfer model. In addition, the thermodynamic properties are obtained from the theoretical vibrations of the optimized HPIMP. The linear polarizability (α(0)) and first-order hyperpolarizabilities (β(0)) respectively present the values of 26.28 Å(3) and 7.41×10(-30) cm(5)/esu predicated theoretically by the DFT-B3lYP method at 6-31G(d) level, which indicates that the studied compound is a promising nonlinear optical material candidate.

Role of Arg82 in the early steps of the bacteriorhodopsin proton-pumping cycle

Proton-transfer reactions in the bacteriorhodopsin light-driven proton pump are coupled with structural rearrangements of protein amino acids and internal water molecules. It is generally thought that the first proton-transfer step from retinal Schiff base to the nearby Asp85 is coupled with movement of the Arg82 side chain away from Asp85 and toward the extracellular proton release group. This movement of Arg82 likely triggers the release of the proton from the proton release group to the extracellular bulk. The exact timing of the movement of Arg82 and how this movement is coupled with proton transfer are still not understood in molecular detail. Here, we address these questions by computing the free energy for the movement of the Arg82 side chain. The calculations indicate that protonation of Asp85 leads to a fast reorientation of the Arg82 side chain toward the extracellular proton release group.

Crystal structures of different substrates of bacteriorhodopsin's M intermediate at various pH levels

The hexagonal P622 crystal of bacteriorhodopsin, which is made up of stacked membranes, is stable provided that the precipitant concentration in the soaking solution is higher than a critical value (i.e., 1.5 M ammonium sulfate). Diffraction data showed that the crystal lattice shrank linearly with increasing precipitant concentration, due primarily to narrowing of intermembrane spaces. Although the crystal shrinkage did not affect the rate of formation of the photoreaction M intermediate, its lifetime increased exponentially with the precipitant concentration. It was suggested that the energetic barrier of the M-to-N transition becomes higher when the motional freedom of the EF loop is reduced by crystal lattice force. As a result of this property, the M state accumulated predominantly when the crystal that was soaked at a high precipitant concentration was illuminated at room temperature. Structural data obtained at various pH levels showed that the overall structure of M is not strongly dependent on pH, except that Glu194 and Glu204 in the proton release complex are more separated at pH 7 than at pH 4.4. This result suggests that light-induced disruption of the paired structure of Glu194 and Glu204 is incomplete when external pH is lower than the pK(a) value of the proton release group in the M state.

Long-distance proton transfer with a break in the bacteriorhodopsin active site

Bacteriorhodopsin is a proton-pumping membrane protein found in the plasma membrane of the archaeon Halobacterium salinarium. Light-induced isomerization of the retinal chromophore from all-trans to 13-cis leads to a sequence of five conformation-coupled proton transfer steps and the net transport of one proton from the cytoplasmic to the extracellular side of the membrane. The mechanism of the long-distance proton transfer from the primary acceptor Asp85 to the extracellular proton release group during the O --> bR is poorly understood. Experiments suggest that this long-distance transfer could involve a transient state [O] in which the proton resides on the intermediate carrier Asp212. To assess whether the transient protonation of Asp212 participates in the deprotonation of Asp85, we performed hybrid Quantum Mechanics/Molecular Mechanics proton transfer calculations using different protein structures and with different retinal geometries and active site water molecules. The structural models were assessed by computing UV-vis excitation energies and C=O vibrational frequencies. The results indicate that a transient [O] conformer with protonated Asp212 could indeed be sampled during the long-distance proton transfer to the proton release group. Our calculations suggest that, in the starting proton transfer state O, the retinal is strongly twisted and at least three water molecules are present in the active site.

Calcium binding to the purple membrane: A molecular dynamics study

The purple membrane (PM) is a specialized membrane patch found in halophilic archaea, containing the photoreceptor bacteriorhodopsin (bR). It is long known that calcium ions bind to the PM, but their position and role remain elusive to date. Molecular dynamics simulations in conjunction with a highly detailed model of the PM have been used to investigate the stability of calcium ions placed at three proposed cation binding sites within bR, one near the Schiff base, one in the region of the proton release group, and one near Glu9. The simulations suggest that, of the sites investigated, the binding of calcium ions was most likely at the proton release group. Binding in the region of the Schiff base, while possible, was associated with significant changes in local geometry. Calcium ions placed near Glu9 in the interior of bR (simultaneously to a Ca(2+) near the Schiff base and another one near the Glu194-Glu204 site) were not stable. The results obtained are discussed in relation to recent experimental observations and theoretical considerations.

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.