The Conformational Behaviour of Fluorinated Tetrahydrothiopyran.

We present a comprehensive study on the conformational behavior of diversely substituted 4-fluorotetrahydrothiopyran derivatives. Through quantum chemical simulations including DFT as well as NBO and NPA analysis, we elucidate the pivotal role of electrostatic interactions, occasionally complemented by hyperconjugative interactions, in stabilizing axial fluorine conformers. Less polar conformers were occasionally obtained, attributed to the interplay of electrostatic and hyperconjugative interactions. Experimental validation through NMR spectroscopy aligns with the computational analysis, thus providing a coherent understanding of the structural dynamics of these compounds.

RhoMax: Computational Prediction of Rhodopsin Absorption Maxima Using Geometric Deep Learning.

Microbial rhodopsins (MRs) are a diverse and abundant family of photoactive membrane proteins that serve as model systems for biophysical techniques. Optogenetics utilizes genetic engineering to insert specialized proteins into specific neurons or brain regions, allowing for manipulation of their activity through light and enabling the mapping and control of specific brain areas in living organisms. The obstacle of optogenetics lies in the fact that light has a limited ability to penetrate biological tissues, particularly blue light in the visible spectrum. Despite this challenge, most optogenetic systems rely on blue light due to the scarcity of red-shifted opsins. Finding additional red-shifted rhodopsins would represent a major breakthrough in overcoming the challenge of limited light penetration in optogenetics. However, determining the wavelength absorption maxima for rhodopsins based on their protein sequence is a significant hurdle. Current experimental methods are time-consuming, while computational methods lack accuracy. The paper introduces a new computational approach called RhoMax that utilizes structure-based geometric deep learning to predict the absorption wavelength of rhodopsins solely based on their sequences. The method takes advantage of AlphaFold2 for accurate modeling of rhodopsin structures. Once trained on a balanced train set, RhoMax rapidly and precisely predicted the maximum absorption wavelength of more than half of the sequences in our test set with an accuracy of 0.03 eV. By leveraging computational methods for absorption maxima determination, we can drastically reduce the time needed for designing new red-shifted microbial rhodopsins, thereby facilitating advances in the field of optogenetics.

The photochemistry and photophysics of benzoyl-carbazole.

Benzoyl-carbazole and its derivatives are considered a platform for exploring processes such as room temperature phosphorescence (RTP) and thermally activated delayed fluorescence (TADF). They have also been reported to exhibit dual emission, but there is a great spectral variability in the relative intensity of the emission bands reported in different studies. To better understand the fundamental photophysical properties, we set to explore BCz and its perfluorinated derivative F5BCz using spectroscopy and quantum chemical simulations. We find that the reported dual fluorescence in solution and in films results from a photochemical process (photo-Fries rearrangement), producing carbazole among other products, explaining the variation in the reported emission spectra. In addition, BCz exhibits solvent dependent TADF, which is explained by the stabilization of the charge transfer S1 state in polar solvents. F5BCz undergoes an efficient photochemical process (Mallory reaction) from its single state to produce highly fluorescent product c-F5BCz, in 40% isolated yield. This photoreactivity also proceeds in films under ambient conditions, which have significant implications on the applications of BCz-based materials for optoelectronic applications.

Merocyanines form bacteriorhodopsins with strongly bathochromic absorption maxima.

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 [Formula: see text] 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 [Formula: see text] 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 [Formula: see text] and the calculated [Formula: see text] 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.

Light-Oxygen-Voltage (LOV)-sensing Domains: Activation Mechanism and Optogenetic Stimulation.

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.

Short-Range Effects in the Special Pair of Photosystem II Reaction Centers: The Nonconservative Nature of Circular Dichroism.

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.

Disequilibrating azobenzenes by visible-light sensitization under confinement.

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.

Embedding Nonrigid Solutes in an Averaged Environment: A Case Study on Rhodopsins.

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.

Transient 2D IR Spectroscopy and Multiscale Simulations Reveal Vibrational Couplings in the Cyanobacteriochrome Slr1393-g3.

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.

Unexpected Nucleophile Masking in Acyl Transfer to Sterically Crowded and Conformationally Restricted Galactosides.

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.

Optical Activity and Spin Polarization: The Surface Effect.

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.

Induction effects on the absorption maxima of photoreceptor proteins.

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.

Photoisomerization of phytochrome chromophore models: an XMS-CASPT2 study.

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.

Retinal chromophore charge delocalization and confinement explain the extreme photophysics of Neorhodopsin.

The understanding of how the rhodopsin sequence can be modified to exactly modulate the spectroscopic properties of its retinal chromophore, is a prerequisite for the rational design of more effective optogenetic tools. One key problem is that of establishing the rules to be satisfied for achieving highly fluorescent rhodopsins with a near infrared absorption. In the present paper we use multi-configurational quantum chemistry to construct a computer model of a recently discovered natural rhodopsin, Neorhodopsin, displaying exactly such properties. We show that the model, that successfully replicates the relevant experimental observables, unveils a geometrical and electronic structure of the chromophore featuring a highly diffuse charge distribution along its conjugated chain. The same model reveals that a charge confinement process occurring along the chromophore excited state isomerization coordinate, is the primary cause of the observed fluorescence enhancement.

The involvement of triplet states in the isomerization of retinaloids.

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.

Rhodopsin-bestrophin fusion proteins from unicellular algae form gigantic pentameric ion channels.

Many organisms sense light using rhodopsins, photoreceptive proteins containing a retinal chromophore. Here we report the discovery, structure and biophysical characterization of bestrhodopsins, a microbial rhodopsin subfamily from marine unicellular algae, in which one rhodopsin domain of eight transmembrane helices or, more often, two such domains in tandem, are C-terminally fused to a bestrophin channel. Cryo-EM analysis of a rhodopsin-rhodopsin-bestrophin fusion revealed that it forms a pentameric megacomplex (~700 kDa) with five rhodopsin pseudodimers surrounding the channel in the center. Bestrhodopsins are metastable and undergo photoconversion between red- and green-absorbing or green- and UVA-absorbing forms in the different variants. The retinal chromophore, in a unique binding pocket, photoisomerizes from all-trans to 11-cis form. Heterologously expressed bestrhodopsin behaves as a light-modulated anion channel.

Ultrafast proton-coupled isomerization in the phototransformation of phytochrome.

The biological function of phytochromes is triggered by an ultrafast photoisomerization of the tetrapyrrole chromophore biliverdin between two rings denoted C and D. The mechanism by which this process induces extended structural changes of the protein is unclear. Here we report ultrafast proton-coupled photoisomerization upon excitation of the parent state (Pfr) of bacteriophytochrome Agp2. Transient deprotonation of the chromophore's pyrrole ring D or ring C into a hydrogen-bonded water cluster, revealed by a broad continuum infrared band, is triggered by electronic excitation, coherent oscillations and the sudden electric-field change in the excited state. Subsequently, a dominant fraction of the excited population relaxes back to the Pfr state, while ~35% follows the forward reaction to the photoproduct. A combination of quantum mechanics/molecular mechanics calculations and ultrafast visible and infrared spectroscopies demonstrates how proton-coupled dynamics in the excited state of Pfr leads to a restructured hydrogen-bond environment of early Lumi-F, which is interpreted as a trigger for downstream protein structural changes.

Signal transduction in light-oxygen-voltage receptors lacking the active-site glutamine.

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.

Serial crystallography captures dynamic control of sequential electron and proton transfer events in a flavoenzyme.

Flavin coenzymes are universally found in biological redox reactions. DNA photolyases, with their flavin chromophore (FAD), utilize blue light for DNA repair and photoreduction. The latter process involves two single-electron transfers to FAD with an intermittent protonation step to prime the enzyme active for DNA repair. Here we use time-resolved serial femtosecond X-ray crystallography to describe how light-driven electron transfers trigger subsequent nanosecond-to-microsecond entanglement between FAD and its Asn/Arg-Asp redox sensor triad. We found that this key feature within the photolyase-cryptochrome family regulates FAD re-hybridization and protonation. After first electron transfer, the FAD•- isoalloxazine ring twists strongly when the arginine closes in to stabilize the negative charge. Subsequent breakage of the arginine-aspartate salt bridge allows proton transfer from arginine to FAD•-. Our molecular videos demonstrate how the protein environment of redox cofactors organizes multiple electron/proton transfer events in an ordered fashion, which could be applicable to other redox systems such as photosynthesis.