Electrostatic Protein−Chromophore Interactions Promote the all-trans → 13-cis Isomerization of the Protonated Retinal Schiff Base in Bacteriorhodopsin: An ab Initio CASSCF/MRCI Study

Tuning the primary reaction of channelrhodopsin-2 by imidazole, pH, and site-specific mutations

Femtosecond time-resolved absorption measurements were performed to investigate the influence of the pH, imidazole concentration, and point mutations on the isomerization process of Channelrhodopsin-2. Apart from the typical spectral characteristics of retinal isomerization, an additional absorption feature rises for the wild-type (wt) on a timescale from tens of ps to 1 ns within the spectral range of the photoproduct and is attributed to an equilibration between different K-intermediates. Remarkably, this absorption feature vanishes upon addition of imidazole or lowering the pH. In the latter case, the isomerization is dramatically slowed down, due to protonation of negatively charged amino acids within the retinal binding pocket, e.g., E123 and D253. Moreover, we investigated the influence of several point mutations within the retinal binding pocket E123T, E123D, C128T, and D156C. For E123T, the isomerization is retarded compared to wt and E123D, indicating that a negatively charged residue at this position functions as an effective catalyst in the isomerization process. In the case of the C128T mutant, all primary processes are slightly accelerated compared to the wt, whereas the isomerization dynamics for the D156C mutant is similar to wt after addition of imidazole.

Temperature investigations of E/Z isomers in ketimines based of p-dibenzoylobenzene with aniline and 2,6-dimethylaniline by infrared spectroscopy

The reaction of aniline with p-dibenzoylobenzene (K1) can lead to Z/Z, Z/E and E/E isomers however the only Z/Z and E/E were formed. At room temperature these isomers may be separated, thus the corresponding FTIR spectra could be recorded. The observed bands were assigned and temperature investigations were lead to monitor the structural changes during heating Z and E forms of K1 from 20 to 240 degrees C. FTIR spectroscopy showed that the bigger changes of the Z form was observed with an increase of temperature. Similar experience was lead with the ketimine synthesized from 2,6-dimethylaniline and p-dibenzoylobenzene (K2) investigated as a mixture of isomers.

Electrostatic and steric interactions determine bacteriorhodopsin single-molecule biomechanics

Bacteriorhodopsin (bR) is a haloarchaeal membrane protein that converts the energy of single photons into large structural changes to directionally pump protons across purple membrane. This is achieved by a complex combination of local dynamic interactions controlling bR biomechanics at the submolecular level, producing efficient amplification of the retinal photoisomerization. Using single molecule force spectroscopy at different salt concentrations, we show that tryptophan (Trp) residues use steric specific interactions to create a rigid scaffold in bR extracellular region and are responsible for the main unfolding barriers. This scaffold, which encloses the retinal, controls bR local mechanical properties and anchors the protein into the membrane. Furthermore, the stable Trp-based network allows ion binding to two specific sites on the extracellular loops (BC and FG), which are involved in proton release and lateral transport. In contrast, the cytoplasmic side of bR is mainly governed by relatively weak nonspecific electrostatic interactions that provide the flexibility necessary for large cytoplasmic structural rearrangements during the photocycle. The presence of an extracellular Trp-based network tightly enclosing the retinal seems common to most haloarchaeal rhodopsins, and could be relevant to their exceptional efficiency.

Excited-State Charge Transfer at a Conical Intersection: Effects of an Environment

The influence of a polar and polarizable environment on charge transfer processes at a conical intersection (CI) can be described by a diabatic free energy model yielding coupled surfaces as a function of both molecular coordinates and a solvent coordinate. We extend and apply this model for the S1-S0 CI in protonated Schiff bases, representing a model for retinal isomerization (Faraday Discuss. 2004, 127, 395, 2004). A dielectric continuum description of the solvent is combined with a minimal, two-electron-two-orbital electronic structure model according to Bonacić-Koutecký, Koutecký, and Michl (Angew. Chem. 1987, 26, 170), which characterizes the charge translocation effects at the CI. The model predicts that the nonequilibrium solvent state resulting from the S0-->S1 Franck-Condon transition can entail the disappearance of the CI, such that solvent motion is necessary to reach the CI seam. The concerted evolution of the intramolecular coordinates and the solvent coordinate is illustrated by an excited-state minimum energy path.

Factors influencing the energetics of electron and proton transfers in proteins. What can be learned from calculations

A protein structure should provide the information needed to understand its observed properties. Significant progress has been made in developing accurate calculations of acid/base and oxidation/reduction reactions in proteins. Current methods and their strengths and weaknesses are discussed. The distribution and calculated ionization states in a survey of proteins is described, showing that a significant minority of acidic and basic residues are buried in the protein and that most of these remain ionized. The electrochemistry of heme and quinones are considered. Proton transfers in bacteriorhodopsin and coupled electron and proton transfers in photosynthetic reaction centers, 5-coordinate heme binding proteins and cytochrome c oxidase are highlighted as systems where calculations have provided insight into the reaction mechanism.

Optimal control of ultrafastcis-transphotoisomerization of retinal in rhodopsin via a conical intersection

Optimal control simulation is applied to the cis-trans photoisomerization of retinal in rhodopsin within a two-dimensional, two-electronic-state model with a conical intersection [S. Hahn and G. Stock, J. Phys. Chem. B 104, 1146 (2000)]. For this case study, we investigate coherent control mechanisms, in which laser pulses work cooperatively with a conical intersection that acts as a "wave-packet cannon." Optimally designed pulses largely consist of shaping subpulses that prepare a wave packet, which is localized along a reaction coordinate and has little energy in the coupling mode, through multiple electronic transitions. This shaping process is shown to be essential for achieving a high target yield although the envelopes of the calculated pulses depend on the local topography of the potential-energy surfaces around the conical intersection and the choice of target. The control mechanisms are analyzed by considering the motion of reduced wave packets in a nuclear configuration space as well as by snapshots of probability current-density maps.

Ultrafast excited state dynamics of the protonated Schiff base of all-trans retinal in solvents

We present a comparative study of the ultrafast photophysics of all-trans retinal in the protonated Schiff base form in solvents with different polarities and viscosities. Steady-state spectra of retinal in the protonated Schiff base form show large absorption-emission Stokes shifts (6500-8100 cm(-1)) for both polar and nonpolar solvents. Using a broadband fluorescence up-conversion experiment, the relaxation kinetics of fluorescence is investigated with 120 fs time resolution. The time-zero spectra already exhibit a Stokes-shift of approximately 6000 cm(-1), indicating depopulation of the Franck-Condon region in < or =100 fs. We attribute it to relaxation along skeletal stretching. A dramatic spectral narrowing is observed on a 150 fs timescale, which we assign to relaxation from the S(2) to the S(1) state. Along with the direct excitation of S(1), this relaxation populates different quasistationary states in S(1), as suggested from the existence of three distinct fluorescence decay times with different decay associated spectra. A 0.5-0.65 ps decay component is observed, which may reflect the direct repopulation of the ground state, in line with the small isomerization yield in solvents. Two longer decay components are observed and are attributed to torsional motion leading to photo-isomerization. The various decay channels show little or no dependence with respect to the viscosity or dielectric constant of the solvents. This suggests that in the protein, the bond selectivity of isomerization is mainly governed by steric effects.

Supramolecular Assemblies and Molecular Recognition of Amphiphilic Schiff Bases with Barbituric Acid in Organized Molecular Films

A bolaform Schiff base, N,N'-bis(salicylidene)-1,10-decanediamine (BSC10), has been synthesized and its interfacial hydrogen bond formation or molecular recognition with barbituric acid was investigated in comparison with that of a single chain Schiff base, 2-hydroxybenzaldehyde-octadecylamine (HBOA). It has been found that while HBOA formed a monolayer at the air/water interface, the bolaform Schiff base formed a multilayer film with ordered layer structure on water surface. When the Schiff bases were spread on the subphase containing barbituric acid, both of the Schiff bases could form hydrogen bonds with barbituric acid in situ in the spreading films. As a result, an increase of the molecular areas in the isotherms was observed. The in situ H-bonded films could be transferred onto solid substrates, and the transferred multilayer films were characterized by various methods such as UV-vis and FT-IR spectrosopies. Spectral changes were observed for the films deposited from the barbituric acid subphase, which supported the hydrogen bond formation between the Schiff bases and barbituric acid. By measuring the MS-TOF of the deposited films dissolved in CHCl3 solution, it was concluded that a 2:1 complex of HBOA with barbituric acid and a 1:2 complex of BSC10 with barbituric acid were formed. On the other hand, when the multilayer films of both Schiff bases were immersed in an aqueous solution of barbituric acid, a similar molecular recognition through the hydrogen bond occurred. A clear conformational change of the alkyl spacer in the bolaform Schiff base was observed during the complex formation with the barbituric acid.

Counterion Controlled Photoisomerization of Retinal Chromophore Models: a Computational Investigation

CASPT2//CASSCF photoisomerization path computations have been used to unveil the effects of an acetate counterion on the photochemistry of two retinal protonated Schiff base (PSB) models: the 2-cis-penta-2,4-dieniminium and the all-trans-epta-2,4,6-trieniminium cations. Different positions/orientations of the counterion have been investigated and related to (i) the spectral tuning and relative stability of the S0, S1, and S2 singlet states; (ii) the selection of the photochemically relevant excited state; (iii) the control of the radiationless decay and photoisomerization rates; and, finally, (iv) the control of the photoisomerization stereospecificity. A rationale for the results is given on the basis of a simple (electrostatic) qualitative model. We show that the model readily explains the computational results providing a qualitative explanation for different aspects of the experimentally observed "environment" dependent PSB photochemistry. Electrostatic effects likely involved in controlling retinal photoisomerization stereoselectivity in the protein are also discussed under the light of these results, and clues for a stereocontrolled electrostatically driven photochemical process are presented. These computations provide a rational basis for the formulation of a mechanistic model for photoisomerization electrostatic catalysis.

Design and photochemical characterization of a biomimetic light-driven Z/E switcher

Protonated Schiff bases (PSBs) of polyenals constitute a class of light-driven switchers selected by biological evolution that provide model compounds for the development of artificial light-driven molecular devices or motors. In the present paper, our primary target is to show, through combined computational and experimental studies, that it is possible to approach the design of artificial PSBs suitable for such applications. Below, we use the methods of computational photochemistry to design and characterize the prototype biomimetic molecular switchers 4-cyclopenten-2'-enylidene-3,4-dihydro-2H-pyrrolinium and its 5,5'-dimethyl derivative both containing the penta-2,4-dieniminium chromophore. To find support for the predicted behavior, we also report the photochemical reaction path of the synthetically accessible compound 4-benzylidene-3,4-dihydro-2H-pyrrolinium. We show that the preparation and photochemical characterization of this compound (together with three different N-methyl derivatives) provide both support for the predicted photoisomerization mechanism and information on its sensitivity to the molecular environment.

Calculation of proton transfers in Bacteriorhodopsin bR and M intermediates

Residue ionization states were calculated in nine crystal structures of bacteriorhodopsin trapped in bR, early M, and late M states by multiconformation continuum electrostatics. This combines continuum electrostatics and molecular mechanics, deriving equilibrium distributions of ionization states and polar residue and water positions. The three central cluster groups [retinal Schiff base (SB), Asp 85 and Asp 212] are ionized in bR structures while a proton has transferred from SB(+) to Asp 85(-) in late M structures matching experimental results. The proton shift in M is due to weaker SB(+)-ionized acid and more favorable SB(0)-ionized acid interactions following retinal isomerization. The proton release cluster (Glu 194 and Glu 204) binds one proton in bR, which is lost to water by pH 8 in late M. In bR the half-ionized state is stabilized by charge-dipole interactions while full ionization is disallowed by charge-charge repulsion between the closely spaced acids. In M the acids move apart, permitting full ionization. Arg 82 movement connects the proton shifts in the central and proton release clusters. Changes in total charge of the two clusters are coupled by direct long-range interactions. Separate calculations consider continuum or explicit water in internal cavities. The explicit waters and nearby polar residues can reorient to stabilize different charge distributions. Proton release to the low-pH, extracellular side of the protein occurs in these calculations where residue ionization remains at equilibrium with the medium. Thus, the key changes distinguishing the intermediates are indeed trapped in the structures.