Publications by Year: 2018

Toker Y, Langeland J, Gruber E, Kjær C, Nielsen SB, Andersen LH, Borin VA, Schapiro I. Counterion-controlled spectral tuning of the protonated Schiff-base retinal. Physical Review A [Internet]. 2018;98 (4) :043428. Publisher's VersionAbstract

Color vision is based on the ability of different opsin proteins to tune the absorption band of their chromophore, the retinal protonated Schiff base (RPSB). Two main mechanisms proposed for this tunability are geometric and electrostatic. Here we probe the latter effect experimentally and by a quantum chemical calculation of the absorption by an isolated complex containing the retinal chromophore and molecules with a strong dipole moment. Betaine complexation causes an anomalously large blue shift. The shift provides direct evidence that the strong charge-transfer character of the electronic transition is the cause of the opsin shift, and shows that the electric field of the counterion is responsible for the color tuning, which allows absorption of light in the blue region of the visible spectrum by opsin proteins.

Guo Y, Wolff FE, Schapiro I, Elstner M, Marazzi M. Different hydrogen bonding environments of the retinal protonated Schiff base control the photoisomerization in channelrhodopsin-2. Physical Chemistry Chemical Physics [Internet]. 2018;20 (43) :27501 - 27509. Publisher's VersionAbstract

The first event of the channelrhodopsin-2 (ChR2) photocycle, i.e. trans-to-cis photoisomerization, is studied by means of quantum mechanics/molecular mechanics, taking into account the flexible retinal environment in the ground state. By treating the chromophore at the ab initio multiconfigurational level of theory, we can rationalize the experimental findings based on pump–probe spectroscopy, explaining the different and more complex scenario found for ChR2 in comparison to other rhodopsins. In particular, we find that depending on the hydrogen bonding pattern, different excited states are involved, hence making it possible to suggest one pattern as the most productive. Moreover, after photoisomerization the structure of the first photocycle intermediate, P5001, is characterized by simulating the infrared spectrum and compared to available experimental data. This was obtained by extensive molecular dynamics, where the chromophore is described by a semi-empirical method based on density functional theory. The results clearly identify which counterion is responsible for accepting the proton from the retinal Schiff base: the side chain of the glutamic acid E123.

Sen S, Schapiro I. A comprehensive benchmark of the XMS-CASPT2 method for the photochemistry of a retinal chromophore model. Molecular Physics [Internet]. 2018;116 (19 - 20) :2571 - 2582. Publisher's VersionAbstract

ABSTRACT The performance of the extended multi-state (XMS)-complete active space second-order perturbation theory (CASPT2) method has been assessed for the benchmark of a truncated retinal model, the penta-2,4-dieniminium cation (PSB3). This benchmark presents a challenge for multireference electronic structure methods because the wave function character is changing considerably. The assessment comprises ground and excited state pathways of the isomerisation, including transition states and conical intersection (CI) points. It also includes circular paths centred around different CIs, and 2D potential energy scans located in the branching planes. In this work, we compare the performance of the previous formulations of CASPT2, the single-state and the multi-state, with the recently developed XMS-CASPT2. Besides, we have also tested two variants of internal contraction in XMS-CASPT2, namely, the single-state single reference (SS-SR) and multi-state multireference (MS-MR) schemes. In our study, we find that XMS-CASPT2 corrects the artefacts and discontinuities present in other CASPT2 variants. The investigation of a circular loop and 2D potential energy surfaces around the surface crossing point shows that XMS-CASPT2 exhibits a smooth topology at the CI with the correct degeneracy. It also agrees better with the reference method MRCISD+Q in regions of the potential energy surfaces further away from CIs. Another observation is the close agreement between the results from the SS-SR contraction scheme and the more demanding MS-MR scheme.

Nogly P, Weinert T, James D, Carbajo S, Ozerov D, Furrer A, Gashi D, Borin V, Skopintsev P, Jaeger K, et al. Retinal isomerization in bacteriorhodopsin captured by a femtosecond x-ray laser. Science [Internet]. 2018. Publisher's VersionAbstract
Ultrafast isomerization of retinal is the primary step in photoresponsive biological functions including vision in humans and ion-transport across bacterial membranes. We studied the sub-picosecond structural dynamics of retinal isomerization in the light-driven proton pump bacteriorhodopsin using an x-ray laser. A series of structural snapshots with near-atomic spatial and temporal resolution in the femtosecond regime show how the excited all-trans retinal samples conformational states within the protein binding pocket prior to passing through a twisted geometry and emerging in the 13-cis conformation. Our findings suggest ultrafast collective motions of aspartic acid residues and functional water molecules in the proximity of the retinal Schiff base as a key ingredient for this stereo-selective and efficient photochemical reaction.
Zill JC, He Z, Tank M, Ferlez BH, Canniffe DP, Lahav Y, Bellstedt P, Alia A, Schapiro I, Golbeck JH, et al. 15N photo-CIDNP MAS NMR analysis of reaction centers of Chloracidobacterium thermophilum. Photosynthesis Research [Internet]. 2018;137 (2) :295-305. Publisher's VersionAbstract

Photochemically induced dynamic nuclear polarization (photo-CIDNP) has been observed in the homodimeric, type-1 photochemical reaction centers (RCs) of the acidobacterium, Chloracidobacterium (Cab.) thermophilum, by 15N magic-angle spinning (MAS) solid-state NMR under continuous white-light illumination. Three light-induced emissive (negative) signals are detected. In the RCs of Cab. thermophilum, three types of (bacterio)chlorophylls have previously been identified: bacteriochlorophyll a (BChl a), chlorophyll a (Chl a), and Zn-bacteriochlorophyll a′ (Zn-BChl a′) (Tsukatani et al. in J Biol Chem 287:5720–5732, 2012). Based upon experimental and quantum chemical 15N NMR data, we assign the observed signals to a Chl a cofactor. We exclude Zn-BChl because of its measured spectroscopic properties. We conclude that Chl a is the primary electron acceptor, which implies that the primary donor is most likely Zn-BChl a′. Chl a and 81-OH Chl a have been shown to be the primary electron acceptors in green sulfur bacteria and heliobacteria, respectively, and thus a Chl a molecule serves this role in all known homodimeric type-1 RCs.

Schnedermann C, Yang X, Liebel M, Spillane KM, Lugtenburg J, Fernández I, Valentini A, Schapiro I, Olivucci M, Kukura P, et al. Evidence for a vibrational phase-dependent isotope effect on the photochemistry of vision. Nature Chemistry [Internet]. 2018;10 (4) :449 - 455. Publisher's VersionAbstract

Vibronic coupling is key to efficient energy flow in molecular systems and a critical component of most mechanisms invoking quantum effects in biological processes. Despite increasing evidence for coherent coupling of electronic states being mediated by vibrational motion, it is not clear how and to what degree properties associated with vibrational coherence such as phase and coupling of atomic motion can impact the efficiency of light-induced processes under natural, incoherent illumination. Here, we show that deuteration of the H11–C11=C12–H12 double-bond of the 11-cis retinal chromophore in the visual pigment rhodopsin significantly and unexpectedly alters the photoisomerization yield while inducing smaller changes in the ultrafast isomerization dynamics assignable to known isotope effects. Combination of these results with non-adiabatic molecular dynamics simulations reveals a vibrational phase-dependent isotope effect that we suggest is an intrinsic attribute of vibronically coherent photochemical processes.

Sutar R L, Sen S, Eivgi O, Segalovich G, Schapiro I, Reany O, Lemcoff NG. Guiding a divergent reaction by photochemical control: bichromatic selective access to levulinates and butenolides. Chemical Science [Internet]. 2018;9 (5) :2041 - 6520. Publisher's VersionAbstract

Allylic and acrylic substrates may be efficiently transformed by a sequential bichromatic photochemical process into derivatives of levulinates or butenolides with high selectivity when phenanthrene is used as a regulator. Thus, UV-A photoinduced cross-metathesis (CM) couples the acrylic and allylic counterparts and subsequent UV-C irradiation initiates E–Z isomerization of the carbon–carbon double bond, followed by one of two competing processes; namely, cyclization by transesterification or a 1,5-H shift and tautomerization. Quantum chemical calculations demonstrate that intermediates are strongly blue-shifted for the cyclization while red-shifted for the 1,5-H shift reaction. Hence, delaying the double bond migration by employing UV-C absorbing phenanthrene, results in a selective novel divergent all-photochemical pathway for the synthesis of fundamental structural motifs of ubiquitous natural products.

Borin V, Wiebeler C, Schapiro I. A QM/MM study of the initial excited state dynamics of green-absorbing Proteorhodopsin. Faraday Discussions [Internet]. 2018;207 :137 - 152. Publisher's VersionAbstract

The primary photochemical reaction of the green-absorbing Proteorhodopsin is studied by means of a hybrid quantum mechanics/molecular mechanics (QM/MM) approach. The simulations are based on a homology model derived from the blue-absorbing Proteorhodopsin crystal structure. The geometry of retinal and the surrounding sidechains in the protein binding pocket were optimized using the QM/MM method. Starting from this geometry the isomerization was studied with a relaxed scan along the C13=C14 dihedral. It revealed an "aborted bicycle pedal" mechanism of isomerization that was originally proposed by Warshel for bovine rhodopsin and bacteriorhodopsin. However, the isomerization involved the concerted rotation about C13=C14 and C15=N, with the latter being highly tiwsted but not isomerized. Further, the simulation showed an increased steric interaction between the hydrogen at the C14 of the isomerizing bond and the hydroxyl group at the neighbouring tyrosine 200. In addition, we have simulated a nonadiabatic trajectory which showed the timing of the isomerization. In the first 20 fs upon excitation the order of the conjugated double and single bonds is inverted, consecutively the C13=C14 rotation is activated for 200 fs until the S1-S0 transition is detected. However, the isomerization is reverted due to the specific interaction with the tyrosine as observed along the relaxed scan calculation. Our simulations indicate that the retinal - tyrosine 200 interaction plays an important role in the outcome of the photoisomerization.