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

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.

Hydrogen-Bonding and Hydrophobic Interaction Networks as Structural Determinants of Microbial Rhodopsin Function

Microbial pump rhodopsins are highly versatile light-driven membrane proteins that couple protein conformational dynamics with ion translocation across the cell membranes. Understanding how microbial pump rhodopsins use specific amino acid residues at key functional sites to control ion selectivity and ion pumping direction is of general interest for membrane transporters, and could guide site-directed mutagenesis for optogenetics applications. To enable direct comparisons between proteins with different sequences we implement, for the first time, a unique numbering scheme for the microbial pump rhodopsin residues, NS-mrho. We use NS-mrho to show that distinct microbial pump rhodopsins typically have hydrogen-bond networks that are less conserved than anticipated from the amino acid residue conservation, whereas their hydrophobic interaction networks are largely conserved. To illustrate the role of the hydrogen-bond networks as structural elements that determine the functionality of microbial pump rhodopsins, we performed experiments, atomic-level simulations, and hydrogen bond network analyses on GR, the outward proton pump from Gloeobacter violaceus, and KR2, the outward sodium pump from Krokinobacter eikastus. The experiments indicate that multiple mutations that recover KR2 amino acid residues in GR not only fail to convert it into a sodium pump, but completely inactivate GR by abolishing photoisomerization of the retinal chromophore. This observation could be attributed to the drastically altered hydrogen-bond interaction network identified with simulations and network analyses. Taken together, our findings suggest that functional specificity could be encoded in the collective hydrogen-bond network of microbial pump rhodopsins.

Thermal isomerization rates in retinal analogues using Ab-Initio molecular dynamics

For a detailed understanding of chemical processes in nature and industry, we need accurate models of chemical reactions in complex environments. While Eyring transition state theory is commonly used for modeling chemical reactions, it is most accurate for small molecules in the gas phase. A wide range of alternative rate theories exist that can better capture reactions involving complex molecules and environmental effects. However, they require that the chemical reaction is sampled by molecular dynamics simulations. This is a formidable challenge since the accessible simulation timescales are many orders of magnitude smaller than typical timescales of chemical reactions. To overcome these limitations, rare event methods involving enhanced molecular dynamics sampling are employed. In this work, thermal isomerization of retinal is studied using tight-binding density functional theory. Results from transition state theory are compared to those obtained from enhanced sampling. Rates obtained from dynamical reweighting using infrequent metadynamics simulations were in close agreement with those from transition state theory. Meanwhile, rates obtained from application of Kramers' rate equation to a sampled free energy profile along a torsional dihedral reaction coordinate were found to be up to three orders of magnitude higher. This discrepancy raises concerns about applying rate methods to one-dimensional reaction coordinates in chemical reactions.

Existence of two substates in the O intermediate of the bacteriorhodopsin photocycle

The proton pumping cycle of bacteriorhodopsin (bR) is initiated when the retinal chromophore with the 13-trans configuration is photo-isomerized into the 13-cis configuration. To understand the recovery processes of the initial retinal configuration that occur in the late stage of the photocycle, we have performed a comprehensive analysis of absorption kinetics data collected at various pH levels and at different salt concentrations. The result of analysis revealed the following features of the late stages of the trans photocycle. i) Two substates occur in the O intermediate. ii) The visible absorption band of the first substate (O1) appears at a much shorter wavelength than that of the late substate (O2). iii) O1 is in rapid equilibrium with the preceding state (N), but O1 becomes less stable than N when an ionizable residue (X₁) with a pKa value of 6.5 (in 2 M KCl) is deprotonated. iv) At a low pH and at a low salt concentration, the decay time constant of O2 is longer than those of the preceding states, but the relationship between these time constants is altered when the medium pH or the salt concentration is increased. On the basis of the present observations and previous studies on the structure of the chromophore in O, we suspect that the retinal chromophore in O1 takes on a distorted 13-cis configuration and the O1-to-O2 transition is accompanied by cis-to-trans isomerization about C13C14 bond.

Mechanisms of long-distance allosteric couplings in proton-binding membrane transporters

Membrane transporters that use proton binding and proton transfer for function couple local protonation change with changes in protein conformation and water dynamics. Changes of protein conformation might be required to allow transient formation of hydrogen-bond networks that bridge proton donor and acceptor pairs separated by long distances. Inter-helical hydrogen-bond networks adjust rapidly to protonation change, and ensure rapid response of the protein structure and dynamics. Membrane transporters with known three-dimensional structures and proton-binding groups inform on general principles of protonation-coupled protein conformational dynamics. Inter-helical hydrogen bond motifs between proton-binding carboxylate groups and a polar sidechain are observed in unrelated membrane transporters, suggesting common principles of coupling protonation change with protein conformational dynamics.

Three-Step Isomerization of the Retinal Chromophore during the Anion Pumping Cycle of Halorhodopsin

The anion pumping cycle of halorhodopsin from Natronomonas pharaonis ( pHR) is initiated when the all- trans/15- anti isomer of retinal is photoisomerized into the 13- cis/15- anti configuration. A recent crystallographic study suggested that a reaction state with 13- cis/15- syn retinal occurred during the anion release process, i.e., after the N state with the 13- cis/15- anti retinal and before the O state with all- trans/15- anti retinal. In this study, we investigated the retinal isomeric composition in a long-living reaction state at various bromide ion concentrations. It was found that the 13- cis isomer (^csHR'), in which the absorption spectrum was blue-shifted by ∼8 nm compared with that of the trans isomer (^taHR), accumulated significantly when a cold suspension of pHR-rich claret membranes in 4 M NaBr was illuminated with continuous light. Analysis of flash-induced absorption changes suggested that the branching of the trans photocycle into the 13- cis isomer (^csHR') occurs during the decay of an O-like state (O') with 13- cis/15- syn retinal; i.e., O' can decay to either ^csHR' or O with all- trans/15- anti retinal. The efficiency of the branching reaction was found to be dependent on the bromide ion concentration. At a very high bromide ion concentration, the anion pumping cycle is described by the scheme ^taHR -( hν) → K → L_1a ↔ L_1b ↔ N ↔ N' ↔ O' ↔ ^csHR' ↔ ^taHR. At a low bromide ion concentration, on the other hand, O' decays into ^taHR via O.