Publications

2024
Flores-Ibarra A, Maia RNA, Olasz B, Church JR, Gotthard G, Schapiro I, Heberle J, Nogly P. Light-Oxygen-Voltage (LOV)-sensing Domains: Activation Mechanism and Optogenetic Stimulation. Photosensory Receptors-Mechanisms and Effects [Internet]. 2024;436 (5) :168356. Publisher's VersionAbstract
The light-oxygen-voltage (LOV) domains of phototropins emerged as essential constituents of light-sensitive proteins, helping initiate blue light-triggered responses. Moreover, these domains have been identified across all kingdoms of life. LOV domains utilize flavin nucleotides as co-factors and undergo structural rearrangements upon exposure to blue light, which activates an effector domain that executes the final output of the photoreaction. LOV domains are versatile photoreceptors that play critical roles in cellular signaling and environmental adaptation; additionally, they can noninvasively sense and control intracellular processes with high spatiotemporal precision, making them ideal candidates for use in optogenetics, where a light signal is linked to a cellular process through a photoreceptor. The ongoing development of LOV-based optogenetic tools, driven by advances in structural biology, spectroscopy, computational methods, and synthetic biology, has the potential to revolutionize the study of biological systems and enable the development of novel therapeutic strategies.
Mackintosh MJ, Hoischen D, Martin H-D, Schapiro I, Gärtner W. Merocyanines form bacteriorhodopsins with strongly bathochromic absorption maxima. [Internet]. 2024;23 (1) :31 - 53. Publisher's VersionAbstract

There is a need to shift the absorbance of biomolecules to the optical transparency window of tissue for applications in optogenetics and photo-pharmacology. There are a few strategies to achieve the so-called red shift of the absorption maxima. Herein, a series of 11 merocyanine dyes were synthesized and employed as chromophores in place of retinal in bacteriorhodopsin (bR) to achieve a bathochromic shift of the absorption maxima relative to bR’s \({\lambda }_{\mathrm{max}}^{a}\) of 568 nm. Assembly with the apoprotein bacterioopsin (bO) led to stable, covalently bound chromoproteins with strongly bathochromic absorbance bands, except for three compounds. Maximal red shifts were observed for molecules 9, 2, and 8 in bR where the \({\lambda }_{\mathrm{max}}^{a}\) was 766, 755, and 736 nm, respectively. While these three merocyanines have different end groups, they share a similar structural feature, namely, a methyl group which is located at the retinal equivalent position 13 of the polyene chain. The absorption and fluorescence data are also presented for the retinal derivatives in their aldehyde, Schiff base (SB), and protonated SB (PSB) forms in solution. According to their hemicyanine character, the PSBs and their analogue bRs exhibited fluorescence quantum yields (Φf) several orders of magnitude greater than native bR (Φf 0.02 to 0.18 versus 1.5 × 10–5 in bR) while also exhibiting much smaller Stokes shifts than bR (400 to 1000 cm−1 versus 4030 cm−1 in bR). The experimental results are complemented by quantum chemical calculations where excellent agreement between the experimental \({\lambda }_{\mathrm{max}}^{a}\)and the calculated \({\lambda }_{\mathrm{max}}^{a}\) was achieved with the second-order algebraic-diagrammatic construction [ADC(2)] method. In addition, quantum mechanics/molecular mechanics (QM/MM) calculations were employed to shed light on the origin of the bathochromic shift of merocyanine 2 in bR compared with native bR.

2023
Dubbert J, Valtolina M, Huber A, Scherz TD, Wölper C, Daniliuc CG, Filiba O, Sen S, Schapiro I, Rizzo F, et al. Tuning the Emission Behaviour of Halogenated Bridged Ethers in Solution, as Solids and as Aggregates by Chalcogen Substitution. ChemPhotoChem [Internet]. 2023;7 (1) :e202200169. Publisher's VersionAbstract

Abstract In this contribution, we describe a set of three chlorinated bridged ethers with varying numbers of sulfur and oxygen atoms. The substitution leads to highly emissive compounds with tunable photophysical properties in relationship to their state of aggregation, i.?e. in solution, as aggregates and in the solid state. Additionally, an in-depth X-ray diffractometric analysis supported by a Hirshfeld study of non-covalent interactions and quantum chemical simulations was carried out. As the outcome, it was found that the content of sulfur in the compounds regulates the tuning of emission in solution as well as in the aggregated states as a consequence of their variation of planarity.

Metzger TS, Batchu H, Kumar A, Fedotov DA, Goren N, Bhowmick DK, Shioukhi I, Yochelis S, Schapiro I, Naaman R, et al. Optical Activity and Spin Polarization: The Surface Effect. Journal of the American Chemical Society [Internet]. 2023;145 (7) :3972 - 3977. Publisher's VersionAbstract

Chirality (‘handedness’) is a property that underlies a broad variety of phenomena in nature. Chiral molecules appear in two forms, and each is a mirror image of the other, the two enantiomers. The chirality of molecules is associated with their optical activity, and circular dichroism is commonly applied to identify the handedness of chiral molecules. Recently, the chiral induced spin selectivity (CISS) effect was established, according to which transfer of electrons within chiral molecules depends on the electron’s spin. Which spin is preferred depends on the handedness of the chiral molecule and the direction of motion of the electron. Several experiments in the past indicated that there may be a relation between the optical activity of the molecules and their spin selectivity. Here, we show that for a molecule containing several stereogenic axes, when adsorbed on a metal substrate, the peaks in the CD spectra have the same signs for the two enantiomers. This is not the case when the molecules are adsorbed on a nonmetallic substrate or dissolved in solution. Quantum chemical simulations are able to explain the change in the CD spectra upon adsorption of the molecules on conductive and nonconductive surfaces. Surprisingly, the CISS properties are similar for the two enantiomers when adsorbed on the metal substrate, while when the molecules are adsorbed on nonmetallic surface, the preferred spin depends on the molecule handedness. This correlation between the optical activity and the CISS effect indicates that the CISS effect relates to the global polarizability of the molecule.Chirality (‘handedness’) is a property that underlies a broad variety of phenomena in nature. Chiral molecules appear in two forms, and each is a mirror image of the other, the two enantiomers. The chirality of molecules is associated with their optical activity, and circular dichroism is commonly applied to identify the handedness of chiral molecules. Recently, the chiral induced spin selectivity (CISS) effect was established, according to which transfer of electrons within chiral molecules depends on the electron’s spin. Which spin is preferred depends on the handedness of the chiral molecule and the direction of motion of the electron. Several experiments in the past indicated that there may be a relation between the optical activity of the molecules and their spin selectivity. Here, we show that for a molecule containing several stereogenic axes, when adsorbed on a metal substrate, the peaks in the CD spectra have the same signs for the two enantiomers. This is not the case when the molecules are adsorbed on a nonmetallic substrate or dissolved in solution. Quantum chemical simulations are able to explain the change in the CD spectra upon adsorption of the molecules on conductive and nonconductive surfaces. Surprisingly, the CISS properties are similar for the two enantiomers when adsorbed on the metal substrate, while when the molecules are adsorbed on nonmetallic surface, the preferred spin depends on the molecule handedness. This correlation between the optical activity and the CISS effect indicates that the CISS effect relates to the global polarizability of the molecule.

Church JR, Olsen JMH, Schapiro I. Induction effects on the absorption maxima of photoreceptor proteins. [Internet]. 2023;20 : - . Publisher's VersionAbstract

Multiscale simulations have been established as a powerful tool to calculate and predict excitation energies in complex systems such as photoreceptor proteins. In these simulations the chromophore is typically treated using quantum mechanical (QM) methods while the protein and surrounding environment are described by a classical molecular mechanics (MM) force field. The electrostatic interactions between these regions are often treated using electrostatic embedding where the point charges in the MM region polarize the QM region. A more sophisticated treatment accounts also for the polarization of the MM region. In this work, the effect of such a polarizable embedding on excitation energies was benchmarked and compared to electrostatic embedding. This was done for two different proteins, the lipid membrane-embedded jumping spider rhodopsin and the soluble cyanobacteriochrome Slr1393g3. It was found that the polarizable embedding scheme produces absorption maxima closer to experimental values. The polarizable embedding scheme was also benchmarked against expanded QM regions and found to be in qualitative agreement. Treating individual residues as polarizable recovered between 50% and 71% of the QM improvement in the excitation energies, depending on the system. A detailed analysis of each amino acid residue in the chromophore binding pocket revealed that aromatic residues result in the largest change in excitation energy compared to the electrostatic embedding. Furthermore, the computational efficiency of polarizable embedding allowed it to go beyond the binding pocket and describe a larger portion of the environment, further improving the results.

Sukhran Y, Alshanski I, Filiba O, Mackintosh MJ, Schapiro I, Hurevich M. Unexpected Nucleophile Masking in Acyl Transfer to Sterically Crowded and Conformationally Restricted Galactosides. The Journal of Organic Chemistry [Internet]. 2023;88 (13) :9313 - 9320. Publisher's VersionAbstract

Design and synthesis of orthogonally protected monosaccharide building blocks are crucial for the preparation of well-defined oligosaccharides in a stereo- and regiocontrolled manner. Selective introduction of protecting groups to partially protected monosaccharides is nontrivial due to the often unpredictable electronic, steric, and conformational effects of the substituents. Abolished reactivity toward a commonly used Lewis base-catalyzed acylation of O-2 was observed in conformationally restricted 4,6-O-benzylidene-3-O-Nap galactoside. Investigation of analogous systems, crystallographic characterization, and quantum chemical calculations highlighted the overlooked conformational and steric considerations, the combination of which produces a unique passivity of the 2-OH nucleophile. Evaluating the role of electrophile counterion and auxiliary base in the acylation of the sterically crowded and conformationally restricted galactoside system revealed an alternative Brønsted base-driven reaction pathway via nucleophilic activation. Insights gained from this model system were utilized to access the target galactoside intermediate within the envisioned synthetic route. The acylation strategy described herein can be implemented in future syntheses of key monomeric building blocks with unique protecting group hierarchies.Design and synthesis of orthogonally protected monosaccharide building blocks are crucial for the preparation of well-defined oligosaccharides in a stereo- and regiocontrolled manner. Selective introduction of protecting groups to partially protected monosaccharides is nontrivial due to the often unpredictable electronic, steric, and conformational effects of the substituents. Abolished reactivity toward a commonly used Lewis base-catalyzed acylation of O-2 was observed in conformationally restricted 4,6-O-benzylidene-3-O-Nap galactoside. Investigation of analogous systems, crystallographic characterization, and quantum chemical calculations highlighted the overlooked conformational and steric considerations, the combination of which produces a unique passivity of the 2-OH nucleophile. Evaluating the role of electrophile counterion and auxiliary base in the acylation of the sterically crowded and conformationally restricted galactoside system revealed an alternative Brønsted base-driven reaction pathway via nucleophilic activation. Insights gained from this model system were utilized to access the target galactoside intermediate within the envisioned synthetic route. The acylation strategy described herein can be implemented in future syntheses of key monomeric building blocks with unique protecting group hierarchies.

Buhrke D, Lahav Y, Rao A, Ruf J, Schapiro I, Hamm P. Transient 2D IR Spectroscopy and Multiscale Simulations Reveal Vibrational Couplings in the Cyanobacteriochrome Slr1393-g3. Journal of the American Chemical Society [Internet]. 2023;145 (29) :15766 - 15775. Publisher's VersionAbstract

Cyanobacteriochromes are bistable photoreceptor proteins with desirable photochemical properties for biotechnological applications, such as optogenetics or fluorescence microscopy. Here, we investigate Slr1393-g3, a cyanobacteriochrome that reversibly photoswitches between a red-absorbing (Pr) and green-absorbing (Pg) form. We applied advanced IR spectroscopic methods to track the sequence of intermediates during the photocycle over many orders of magnitude in time. In the conversion from Pg to Pr, we have revealed a new intermediate with distinct spectroscopic features in the IR, which precedes Pr formation using transient IR spectroscopy. In addition, stationary and transient 2D IR experiments measured the vibrational couplings between different groups of the chromophore and the protein in these intermediate states, as well as their structural disorder. Anharmonic QM/MM calculations predict spectra in good agreement with experimental 2D IR spectra of the initial and final states of the photocycle. They facilitate the assignment of the IR spectra that serve as a basis for the interpretation of the spectroscopic results and suggest structural changes of the intermediates along the photocycle.Cyanobacteriochromes are bistable photoreceptor proteins with desirable photochemical properties for biotechnological applications, such as optogenetics or fluorescence microscopy. Here, we investigate Slr1393-g3, a cyanobacteriochrome that reversibly photoswitches between a red-absorbing (Pr) and green-absorbing (Pg) form. We applied advanced IR spectroscopic methods to track the sequence of intermediates during the photocycle over many orders of magnitude in time. In the conversion from Pg to Pr, we have revealed a new intermediate with distinct spectroscopic features in the IR, which precedes Pr formation using transient IR spectroscopy. In addition, stationary and transient 2D IR experiments measured the vibrational couplings between different groups of the chromophore and the protein in these intermediate states, as well as their structural disorder. Anharmonic QM/MM calculations predict spectra in good agreement with experimental 2D IR spectra of the initial and final states of the photocycle. They facilitate the assignment of the IR spectra that serve as a basis for the interpretation of the spectroscopic results and suggest structural changes of the intermediates along the photocycle.

Ricardi N, González-Espinoza CE, Adam S, Church JR, Schapiro I, Wesołowski TA. Embedding Nonrigid Solutes in an Averaged Environment: A Case Study on Rhodopsins. Journal of Chemical Theory and Computation [Internet]. 2023;19 (15) :5289 - 5302. Publisher's VersionAbstract

Many simulation methods concerning solvated molecules are based on the assumption that the solvated species and the solvent can be characterized by some representative structures of the solute and some embedding potential corresponding to this structure. While the averaging of the solvent configurations to obtain an embedding potential has been studied in great detail, this hinges on a single solute structure representation. This assumption is re-examined and generalized for conformationally flexible solutes and tested on 4 nonrigid systems. In this generalized approach, the solute is characterized by a set of representative structures and the corresponding embedding potentials. The representative structures are identified by means of subdividing the statistical ensemble, which in this work is generated by a constant-temperature molecular dynamics simulation. The embedding potential defined in the Frozen-Density Embedding Theory is used to characterize the average effect of the solvent in each subensemble. The numerical examples concern the vertical excitation energies of protonated retinal Schiff bases in protein environments. It is comprehensively shown that subensemble averaging leads to huge computational savings compared with explicit averaging of the excitation energies in the whole ensemble while introducing only minor errors in the case of the systems examined.Many simulation methods concerning solvated molecules are based on the assumption that the solvated species and the solvent can be characterized by some representative structures of the solute and some embedding potential corresponding to this structure. While the averaging of the solvent configurations to obtain an embedding potential has been studied in great detail, this hinges on a single solute structure representation. This assumption is re-examined and generalized for conformationally flexible solutes and tested on 4 nonrigid systems. In this generalized approach, the solute is characterized by a set of representative structures and the corresponding embedding potentials. The representative structures are identified by means of subdividing the statistical ensemble, which in this work is generated by a constant-temperature molecular dynamics simulation. The embedding potential defined in the Frozen-Density Embedding Theory is used to characterize the average effect of the solvent in each subensemble. The numerical examples concern the vertical excitation energies of protonated retinal Schiff bases in protein environments. It is comprehensively shown that subensemble averaging leads to huge computational savings compared with explicit averaging of the excitation energies in the whole ensemble while introducing only minor errors in the case of the systems examined.

Gemen J, Church JR, Ruoko T-P, Durandin N, Białek MJ, Weißenfels M, Feller M, Kazes M, Odaybat M, Borin VA, et al. Disequilibrating azobenzenes by visible-light sensitization under confinement. Science [Internet]. 2023;381 (6664) :1357 - 1363. Publisher's VersionAbstract

Photoisomerization of azobenzenes from their stable E isomer to the metastable Z state is the basis of numerous applications of these molecules. However, this reaction typically requires ultraviolet light, which limits applicability. In this study, we introduce disequilibration by sensitization under confinement (DESC), a supramolecular approach to induce the E-to-Z isomerization by using light of a desired color, including red. DESC relies on a combination of a macrocyclic host and a photosensitizer, which act together to selectively bind and sensitize E-azobenzenes for isomerization. The Z isomer lacks strong affinity for and is expelled from the host, which can then convert additional E-azobenzenes to the Z state. In this way, the host?photosensitizer complex converts photon energy into chemical energy in the form of out-of-equilibrium photostationary states, including ones that cannot be accessed through direct photoexcitation. Chemists often strive to push reactions metaphorically uphill toward less energetically favorable products. The challenge is to keep those products from rolling right back down. Gemen et al. report a clever tactic for twisting azobenzene into its higher-energy Z conformation. Specifically, they lured the more stable E isomer into a supramolecular host, along with a photosensitizer. When visible light injects energy to induce the twist, the Z isomer no longer fits in the cavity, so it gets pushed out before more light can twist it back. ?Jake S. Yeston Spatial constraints in a supramolecular host selectively convert azobenzenes to their metastable state under visible light.Photoisomerization of azobenzenes from their stable E isomer to the metastable Z state is the basis of numerous applications of these molecules. However, this reaction typically requires ultraviolet light, which limits applicability. In this study, we introduce disequilibration by sensitization under confinement (DESC), a supramolecular approach to induce the E-to-Z isomerization by using light of a desired color, including red. DESC relies on a combination of a macrocyclic host and a photosensitizer, which act together to selectively bind and sensitize E-azobenzenes for isomerization. The Z isomer lacks strong affinity for and is expelled from the host, which can then convert additional E-azobenzenes to the Z state. In this way, the host?photosensitizer complex converts photon energy into chemical energy in the form of out-of-equilibrium photostationary states, including ones that cannot be accessed through direct photoexcitation. Chemists often strive to push reactions metaphorically uphill toward less energetically favorable products. The challenge is to keep those products from rolling right back down. Gemen et al. report a clever tactic for twisting azobenzene into its higher-energy Z conformation. Specifically, they lured the more stable E isomer into a supramolecular host, along with a photosensitizer. When visible light injects energy to induce the twist, the Z isomer no longer fits in the cavity, so it gets pushed out before more light can twist it back. ?Jake S. Yeston Spatial constraints in a supramolecular host selectively convert azobenzenes to their metastable state under visible light.

Li Manni G, Fdez. Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, et al. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry. Journal of Chemical Theory and Computation [Internet]. 2023;19 (20) :6933 - 6991. Publisher's VersionAbstract

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

Gemeinhardt FG, Lahav Y, Schapiro I, Noy D, Müh F, Lindorfer D, Renger T. Short-Range Effects in the Special Pair of Photosystem II Reaction Centers: The Nonconservative Nature of Circular Dichroism. The Journal of Physical Chemistry Letters [Internet]. 2023;14 (51) :11758 - 11767. Publisher's VersionAbstract

Photosystem II reaction centers extract electrons from water, providing the basis of oxygenic life on earth. Among the light-sensitive pigments of the reaction center, a central chlorophyll a dimer, known as the special pair, so far has escaped a complete theoretical characterization of its excited state properties. The close proximity of the special pair pigments gives rise to short-range effects that comprise a coupling between local and charge transfer (CT) excited states as well as other intermolecular quantum effects. Using a multiscale simulation and a diabatization technique, we show that the coupling to CT states is responsible for 45% of the excitonic coupling in the special pair. The other short-range effects cause a nonconservative nature of the circular dichroism spectrum of the reaction center by effectively rotating the electric transition dipole moments of the special pair pigments inverting and strongly enhancing their intrinsic rotational strength.Photosystem II reaction centers extract electrons from water, providing the basis of oxygenic life on earth. Among the light-sensitive pigments of the reaction center, a central chlorophyll a dimer, known as the special pair, so far has escaped a complete theoretical characterization of its excited state properties. The close proximity of the special pair pigments gives rise to short-range effects that comprise a coupling between local and charge transfer (CT) excited states as well as other intermolecular quantum effects. Using a multiscale simulation and a diabatization technique, we show that the coupling to CT states is responsible for 45% of the excitonic coupling in the special pair. The other short-range effects cause a nonconservative nature of the circular dichroism spectrum of the reaction center by effectively rotating the electric transition dipole moments of the special pair pigments inverting and strongly enhancing their intrinsic rotational strength.

2022
Church JR, Olsen JMH, Schapiro I. The impact of retinal configuration on the protein–chromophore interactions in bistable jumping spider rhodopsin-1. Molecules [Internet]. 2022;27 (1). Publisher's VersionAbstract
Bistable rhodopsins have two stable forms that can be interconverted by light. Due to their ability to act as photoswitches, these proteins are considered as ideal candidates for applications such as optogenetics. In this work, we analyze a recently crystalized bistable rhodopsin, namely the jumping spider rhodopsin-1 (JSR1). This rhodopsin exhibits identical absorption maxima for the parent and the photoproduct form, which impedes its broad application. We performed hybrid QM/MM simulations to study three isomers of the retinal chromophore: the 9-cis, 11-cis and all-trans configurations. The main aim was to gain insight into the specific interactions of each isomer and their impact on the absorption maximum in JSR1. The absorption spectra were computed using sampled snapshots from QM/MM molecular dynamics trajectories and compared to their experimental counterparts. The chromophore–protein interactions were analyzed by visualizing the electrostatic potential of the protein and projecting it onto the chromophore. It was found that the distance between a nearby tyrosine (Y126) residue plays a larger role in the predicted absorption maximum than the primary counterion (E194). Geometric differences between the isomers were also noted, including a structural change in the polyene chain of the chromophore, as well as changes in the nearby hydrogen bonding network. © 2021 by the authors. Licensee MDPI, Basel, Switzerland.
Malakar P, Borin V, Bedi A, Schapiro I, Gidron O, Ruhman S. The impact of twisting on the intersystem crossing in acenes: An experimental and computational study. Physical Chemistry Chemical Physics [Internet]. 2022;24 (4) :2357 - 2362. Publisher's VersionAbstract
Due to their unique excited state dynamics, acenes play a dominant role in optoelectronic and light-harvesting applications. Their optical and electronic properties are typically tailored by side-group engineering, which often result in distortion of the acene core from planarity. However, the effect of such distortion on their excited state dynamics is not clear. In this work, we investigate the effect of twisting on the photophysics of acenes, which are helically locked to a defined twist angle by tethers of different lengths. Ultrafast transient absorption and time resolved fluorescence show a clear dependence of the rate of intersystem crossing with twisting. This trend is explained using quantum chemical calculations, showing an increase of spin-orbit coupling (SOC). At much earlier times, structural reorganization in S1, including coherent vibrational wave packet motions, is reflected in transient spectral changes. As predicted by theory, decreasing the length of diagonal tether induces enhanced activity and frequency blue-shifting of a normal vibration consisting of anthracene twisting against restraint of the tethering chain. Overall, these results serve as design principles for tuning photophysical properties of acenes via controlled twisting of their aromatic core. © 2022 the Owner Societies.
Dietler J, Gelfert R, Kaiser J, Borin V, Renzl C, Pilsl S, Ranzani AT, García de Fuentes A, Gleichmann T, Diensthuber RP, et al. Signal transduction in light-oxygen-voltage receptors lacking the active-site glutamine. Nature Communications [Internet]. 2022;13 (1). Publisher's VersionAbstract
In nature as in biotechnology, light-oxygen-voltage photoreceptors perceive blue light to elicit spatiotemporally defined cellular responses. Photon absorption drives thioadduct formation between a conserved cysteine and the flavin chromophore. An equally conserved, proximal glutamine processes the resultant flavin protonation into downstream hydrogen-bond rearrangements. Here, we report that this glutamine, long deemed essential, is generally dispensable. In its absence, several light-oxygen-voltage receptors invariably retained productive, if often attenuated, signaling responses. Structures of a light-oxygen-voltage paradigm at around 1 Å resolution revealed highly similar light-induced conformational changes, irrespective of whether the glutamine is present. Naturally occurring, glutamine-deficient light-oxygen-voltage receptors likely serve as bona fide photoreceptors, as we showcase for a diguanylate cyclase. We propose that without the glutamine, water molecules transiently approach the chromophore and thus propagate flavin protonation downstream. Signaling without glutamine appears intrinsic to light-oxygen-voltage receptors, which pertains to biotechnological applications and suggests evolutionary descendance from redox-active flavoproteins. © 2022, The Author(s).
Hou Y-N, Höppner A, Rao AG, Lahav Y, Kumar Das P, Ding W-L, Jiang X-X, Hu J-L, Schapiro I, Noy D, et al. Control of a far-red/near-infrared spectral switch in an artificial fluorescent biliprotein derived from allophycocyanin. Protein Science [Internet]. 2022;31 (9). Publisher's VersionAbstract
The molecular structure of mBDFP, a far-red fluorescent protein (FPs) derived from an allophycocyanin homolog was resolved to 2.52 Å. Its biliverdin chromophore was found to be attached to the protein in an unusual way that was never observed in natural phycobiliproteins, and only once in a sub-population of artificial bacteriophytochrome-derived FPs. One of the biliverdin's vinyl groups had two cysteine residues covalently bound to its two carbon atoms. This reduces the conjugation length of the biliverdin π-electron system, which shifts the absorption and emission spectra by about 40 nm, from the near-infrared to the far-red region of the spectrum. By spectrally characterizing a set of mBDFP mutants, we show that such spectral shifts can be induced by modifying a single residue in either one of two critical positions in the vicinity of the binding cysteines. This changes the reactivity of biliverdin and the cysteine's thiols towards forming one, or two thioether bonds to the vinyl group. The ability to control the spectral properties of BDFP by specific point mutations opens many possibilities for rational design of far-red and near-infrared FPs that are of great interest to the development of fluorescence markers for bioimaging since most biological tissues are transparent in this spectral window. © 2022 The Protein Society.
Church JR, Amoyal GS, Borin VA, Adam S, Olsen JMH, Schapiro I. Deciphering the Spectral Tuning Mechanism in Proteorhodopsin: The Dominant Role of Electrostatics Instead of Chromophore Geometry. Chemistry - A European Journal [Internet]. 2022;28 (28). Publisher's VersionAbstract
Proteorhodopsin (PR) is a photoactive proton pump found in marine bacteria. There are two phenotypes of PR exhibiting an environmental adaptation to the ocean's depth which tunes their maximum absorption: blue-absorbing proteorhodopsin (BPR) and green-absorbing proteorhodopsin (GPR). This blue/green color-shift is controlled by a glutamine to leucine substitution at position 105 which accounts for a 20 nm shift. Typically, spectral tuning in rhodopsins is rationalized by the external point charge model but the Q105L mutation is charge neutral. To study this tuning mechanism, we employed the hybrid QM/MM method with sampling from molecular dynamics. Our results reveal that the positive partial charge of glutamine near the C14−C15 bond of retinal shortens the effective conjugation length of the chromophore compared to the leucine residue. The derived mechanism can be applied to explain the color regulation in other retinal proteins and can serve as a guideline for rational design of spectral shifts. © 2022 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.
Mous S, Gotthard G, Ehrenberg D, Sen S, Weinert T, Johnson PJM, James D, Nass K, Furrer A, Kekilli D, et al. Dynamics and mechanism of a light-driven chloride pump. Science [Internet]. 2022;375 (6583) :845 - 851. Publisher's VersionAbstract
Chloride transport by microbial rhodopsins is an essential process for which molecular details such as the mechanisms that convert light energy to drive ion pumping and ensure the unidirectionality of the transport have remained elusive. We combined time-resolved serial crystallography with time-resolved spectroscopy and multiscale simulations to elucidate the molecular mechanism of a chloride-pumping rhodopsin and the structural dynamics throughout the transport cycle. We traced transient anion-binding sites, obtained evidence for how light energy is used in the pumping mechanism, and identified steric and electrostatic molecular gates ensuring unidirectional transport. An interaction with the p-electron system of the retinal supports transient chloride ion binding across a major bottleneck in the transport pathway. These results allow us to propose key mechanistic features enabling finely controlled chloride transport across the cell membrane in this light-powered chloride ion pump. © 2022 American Association for the Advancement of Science. All rights reserved.
Sen S, Kar RK, Borin VA, Schapiro I. Insight into the isomerization mechanism of retinal proteins from hybrid quantum mechanics/molecular mechanics simulations. Wiley Interdisciplinary Reviews: Computational Molecular Science [Internet]. 2022;12 (1). Publisher's VersionAbstract
The photoisomerization of retinal is a unifying primary event in the rhodopsin protein family. In vertebrate rhodopsins it is the first step in the vision process, while in the microbial rhodopsins it activates the transport of ions across the cell-membrane. This reaction is highly optimized in the protein, which is ultrafast, selective, and efficient. A great effort was directed to elucidate the mechanism due to the overall complexity of the process inside the protein. The classical one-bond-flip is too demanding in space for the confined protein cavity. Therefore, various space saving mechanisms based on the rotation of multiple double bonds have been proposed. The hybrid quantum mechanics/molecular mechanics (QM/MM) method played an important role in the elucidation of the mechanism inside the tight protein environment. It allows to take the entire protein into account while describing the ground and excited states of retinal. The predicted mechanisms include full isomerization of two or three double bonds, a simultaneous isomerization of a single and a double bond as well as the partial rotation of bonds adjacent to the central isomerization. This review summarizes mechanistic studies in the literature and compares them. This article is categorized under: Structure and Mechanism > Computational Biochemistry and Biophysics Electronic Structure Theory > Combined QM/MM Methods Software > Molecular Modeling. © 2021 Wiley Periodicals LLC.
Filiba O, Borin VA, Schapiro I. The involvement of triplet states in the isomerization of retinaloids. Physical Chemistry Chemical Physics [Internet]. 2022;114 (1). Publisher's VersionAbstract
Rhodopsins form a family of photoreceptor proteins which utilize the retinal chromophore for light energy conversion. Upon light absorption the retinal chromophore undergoes a photoisomerization. This reaction involves a non-radiative relaxation through a conical intersection between the singlet excited state and the ground state. In this work we studied the possible involvement of triplet states in the photoisomerization of retinaloids using the extended multistate (XMS) version of CASPT2. To this end, truncated models of three retinaloids were considered: protonated Schiff base, deprotonated Schiff base and the aldehyde form. The optimized geometries of the reactant, the product and the conical intersection were connected by a linear interpolation of internal coordinates to describe the isomerization. The energetic position of the low-lying singlet and triplet states as well as their spin-orbit coupling matrix elements (SOCME) were calculated along the isomerization profile. The SOCME values peaked in vicinity of the conical intersection for all the retinaloids. Furthermore, the magnitude of SOCME is invariant to the number of double bonds in the model. The SOCME for the protonated Schiff base is negligible (1.5 cm−1) which renders the involvement of the triplet state as improbable. However, the largest SOCME value of 30 cm−1 was found for the aldehyde form, followed by 15 cm−1 for the deprotonated Schiff base. © 2022 The Royal Society of Chemistry.
Rao AG, Schapiro I. Photoisomerization of phytochrome chromophore models: an XMS-CASPT2 study. Physical Chemistry Chemical Physics [Internet]. 2022;24 (48) :29393 - 29405. Publisher's VersionAbstract
Phytochromes are a superfamily of photoreceptors that harbor linear tetrapyrroles as chromophores. Upon light illumination, the linear tetrapyrrole chromophore undergoes a double bond isomerization which starts a photocycle. In this work, we studied the photoisomerization of chromophore models designed based on the C- and D-rings of the phycocyanobilin (PCB) chromophore. In total, five different models with varying substitutions were investigated. Firstly, the vertical excitation energies were benchmarked using different computational methods to establish the relative order of the excited states. Based on these calculations, we computed the photoisomerization profiles using the extended multi-state (XMS) version of the CASPT2 method. The profiles were obtained for both the clockwise and counterclockwise rotations of the C15 = C16 bond in the Z and E isomers using a linear interpolation of internal coordinates between the Franck-Condon and MECI geometries. In the minimal chromophore model that lacks the substitutions at the pyrrole rings, the isomerization involves both C14-C15 and C15 C16 bonds of the methine bridge between the C- and D-rings, resembling the hula-twist motion. The MECIs are characterized by a partial charge transfer between the two pyrrole rings pointing towards a twisted intramolecular charge transfer. Systematic introduction of substituents leads to an increase in the steric repulsion between the two pyrrole rings causing a pretwist of the dihedral around the C15 C16 bond, which creates a preference for the counterclockwise isomerization. An introduction of the carbonyl group at the D-ring increases the extent of charge transfer which changes the isomerization mechanism from hula-twist to one-bond flip. © 2022 The Royal Society of Chemistry.

Pages