Femtosecond-to-millisecond structural changes in a light-driven sodium pump.

Light-driven sodium pumps actively transport small cations across cellular membranes1. These pumps are used by microorganisms to convert light into membrane potential and have become useful optogenetic tools with applications in neuroscience. Although the resting state structures of the prototypical sodium pump Krokinobacter eikastus rhodopsin 2 (KR2) have been solved2,3, it is unclear how structural alterations over time allow sodium to be translocated against a concentration gradient. Here, using the Swiss X-ray Free Electron Laser4, we have collected serial crystallographic data at ten pump-probe delays from femtoseconds to milliseconds. High-resolution structural snapshots throughout the KR2 photocycle show how retinal isomerization is completed on the femtosecond timescale and changes the local structure of the binding pocket in the early nanoseconds. Subsequent rearrangements and deprotonation of the retinal Schiff base open an electrostatic gate in microseconds. Structural and spectroscopic data, in combination with quantum chemical calculations, indicate that a sodium ion binds transiently close to the retinal within one millisecond. In the last structural intermediate, at 20 milliseconds after activation, we identified a potential second sodium-binding site close to the extracellular exit. These results provide direct molecular insight into the dynamics of active cation transport across biological membranes.

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.

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.

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.

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.

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.

Structural elements regulating the photochromicity in a cyanobacteriochrome.

The three-dimensional (3D) crystal structures of the GAF3 domain of cyanobacteriochrome Slr1393 (Synechocystis PCC6803) carrying a phycocyanobilin chromophore could be solved in both 15-Z dark-adapted state, Pr, λmax = 649 nm, and 15-E photoproduct, Pg, λmax = 536 nm (resolution, 1.6 and 1.86 Å, respectively). The structural data allowed identifying the large spectral shift of the Pr-to-Pg conversion as resulting from an out-of-plane rotation of the chromophore's peripheral rings and an outward movement of a short helix formed from a formerly unstructured loop. In addition, a third structure (2.1-Å resolution) starting from the photoproduct crystals allowed identification of elements that regulate the absorption maxima. In this peculiar form, generated during X-ray exposition, protein and chromophore conformation still resemble the photoproduct state, except for the D-ring already in 15-Z configuration and tilted out of plane akin the dark state. Due to its formation from the photoproduct, it might be considered an early conformational change initiating the parental state-recovering photocycle. The high quality and the distinct features of the three forms allowed for applying quantum-chemical calculations in the framework of multiscale modeling to rationalize the absorption maxima changes. A systematic analysis of the PCB chromophore in the presence and absence of the protein environment showed that the direct electrostatic effect is negligible on the spectral tuning. However, the protein forces the outer pyrrole rings of the chromophore to deviate from coplanarity, which is identified as the dominating factor for the color regulation.

The interplay between chromophore and protein determines the extended excited state dynamics in a single-domain phytochrome.

Phytochromes are a diverse family of bilin-binding photoreceptors that regulate a wide range of physiological processes. Their photochemical properties make them attractive for applications in optogenetics and superresolution microscopy. Phytochromes undergo reversible photoconversion triggered by the Z ⇄ E photoisomerization about the double bond in the bilin chromophore. However, it is not fully understood at the molecular level how the protein framework facilitates the complex photoisomerization dynamics. We have studied a single-domain bilin-binding photoreceptor All2699g1 (Nostoc sp. PCC 7120) that exhibits photoconversion between the red light-absorbing (Pr) and far red-absorbing (Pfr) states just like canonical phytochromes. We present the crystal structure and examine the photoisomerization mechanism of the Pr form as well as the formation of the primary photoproduct Lumi-R using time-resolved spectroscopy and hybrid quantum mechanics/molecular mechanics simulations. We show that the unusually long excited state lifetime (broad lifetime distribution centered at ∼300 picoseconds) is due to the interactions between the isomerizing pyrrole ring D and an adjacent conserved Tyr142. The decay kinetics shows a strongly distributed character which is imposed by the nonexponential protein dynamics. Our findings offer a mechanistic insight into how the quantum efficiency of the bilin photoisomerization is tuned by the protein environment, thereby providing a structural framework for engineering bilin-based optical agents for imaging and optogenetics applications.

Deciphering the Spectral Tuning Mechanism in Proteorhodopsin: The Dominant Role of Electrostatics Instead of Chromophore Geometry.

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.

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.

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.

The impact of twisting on the intersystem crossing in acenes: an experimental and computational study.

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.

Two-photon conversion of a bacterial phytochrome.

In nature, sensory photoreceptors underlie diverse spatiotemporally precise and generally reversible biological responses to light. Photoreceptors also serve as genetically encoded agents in optogenetics to control by light organismal state and behavior. Phytochromes represent a superfamily of photoreceptors that transition between states absorbing red light (Pr) and far-red light (Pfr), thus expanding the spectral range of optogenetics to the near-infrared range. Although light of these colors exhibits superior penetration of soft tissue, the transmission through bone and skull is poor. To overcome this fundamental challenge, we explore the activation of a bacterial phytochrome by a femtosecond laser emitting in the 1 µm wavelength range. Quantum chemical calculations predict that bacterial phytochromes possess substantial two-photon absorption cross sections. In line with this notion, we demonstrate that the photoreversible Pr ↔ Pfr conversion is driven by two-photon absorption at wavelengths between 1170 and 1450 nm. The Pfr yield was highest for wavelengths between 1170 and 1280 nm and rapidly plummeted beyond 1300 nm. By combining two-photon activation with bacterial phytochromes, we lay the foundation for enhanced spatial resolution in optogenetics and unprecedented penetration through bone, skull, and soft tissue.

Spectral tuning of chlorophylls in proteins - electrostatics vs. ring deformation.

In photosynthetic complexes, tuning of chlorophyll light-absorption spectra by the protein environment is crucial to their efficiency and robustness. Recombinant type II water soluble chlorophyll-binding proteins from Brassicaceae (WSCPs) are useful for studying spectral tuning mechanisms due to their symmetric homotetramer structure, and the ability to rigorously modify the chlorophyll's protein surroundings. Our previous comparison of the crystal structures of two WSCP homologues suggested that protein-induced chlorophyll ring deformation is the predominant spectral tuning mechanism. Here, we implement a more rigorous analysis based on hybrid quantum mechanics and molecular mechanics calculations to quantify the relative contributions of geometrical and electrostatic factors to the absorption spectra of WSCP-chlorophyll complexes. We show that when considering conformational dynamics, geometry distortions such as chlorophyll ring deformation accounts for about one-third of the spectral shift, whereas the direct polarization of the electron density accounts for the remaining two-thirds. From a practical perspective, protein electrostatics is easier to manipulate than chlorophyll conformations, thus, it may be more readily implemented in designing artificial protein-chlorophyll complexes.

Correction: Histidine protonation controls structural heterogeneity in the cyanobacteriochrome AnPixJg2.

Correction for 'Histidine protonation controls structural heterogeneity in the cyanobacteriochrome AnPixJg2' by Aditya G. Rao et al., Phys. Chem. Chem. Phys., 2021, DOI: 10.1039/d0cp05314g.

Histidine protonation controls structural heterogeneity in the cyanobacteriochrome AnPixJg2.

Cyanobacteriochromes are compact and spectrally diverse photoreceptor proteins that bind a linear tetrapyrrole as a chromophore. They show photochromicity by having two stable states that can be interconverted by the photoisomerization of the chromophore. These photochemical properties make them an attractive target for biotechnological applications. However, their application is impeded by structural heterogeneity that reduces the yield of the photoconversion. The heterogeneity can originate either from the chromophore structure or the protein environment. Here, we study the origin of the heterogeneity in AnPixJg2, a representative member of the red/green cyanobacteriochrome family, that has a red absorbing parental state and a green absorbing photoproduct state. Using molecular dynamics simulations and umbrella sampling we have identified the protonation state of a conserved histidine residue as a trigger for structural heterogeneity. When the histidine is in a neutral form, the chromophore structure is homogenous, while in a positively charged form, the chromophore is heterogeneous with two different conformations. We have identified a correlation between the protonation of the histidine and the structural heterogeneity of the chromophore by detailed characterization of the interactions in the protein binding site. Our findings reconcile seemingly contradicting spectroscopic studies that attribute the heterogeneity to different sources. Furthermore, we predict that circular dichroism can be used as a diagnostic tool to distinguish different substates.

Bending versus Twisting Acenes - A Computational Study.

Polycyclic aromatic hydrocarbons (PAHs) are widely used in organic electronic devices. The electronic, magnetic, and optical properties of PAHs can be tuned by structural modifications to the aromatic backbone to introduce an inherent distortion from planarity, such as bending or twisting. However, it remains difficult to isolate and control the effects of such distortions. Here, we sought to understand how backbone twisting and bending affect the electronic properties of acenes, as models for larger PAHs. We found that, even when highly distorted from planarity (30° per ring), acenes maintain their aromatic character and π orbital delocalization with minor mixing of the σ and π orbitals. In addition, the energy gap between the HOMO and LUMO decreases with increasing twist, while the gap is hardly affected by bending, since the energy of both orbitals increase to a similar extent. For bent acenes in the triplet state, the spin becomes more localized with increasing bend, whereas twisting produces an evenly distributed spin delocalization. These findings can guide the synthesis of PAHs with tailored properties.

The Impact of Retinal Configuration on the Protein-Chromophore Interactions in Bistable Jumping Spider Rhodopsin-1.

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.

Modulation of Light Energy Transfer from Chromophore to Protein in the Channelrhodopsin ReaChR.

The function of photoreceptors relies on efficient transfer of absorbed light energy from the chromophore to the protein to drive conformational changes that ultimately generate an output signal. In retinal-binding proteins, mainly two mechanisms exist to store the photon energy after photoisomerization: 1) conformational distortion of the prosthetic group retinal, and 2) charge separation between the protonated retinal Schiff base (RSBH+) and its counterion complex. Accordingly, energy transfer to the protein is achieved by chromophore relaxation and/or reduction of the charge separation in the RSBH+-counterion complex. Combining FTIR and UV-Vis spectroscopy along with molecular dynamics simulations, we show here for the widely used, red-activatable Volvox carteri channelrhodopsin-1 derivate ReaChR that energy storage and transfer into the protein depends on the protonation state of glutamic acid E163 (Ci1), one of the counterions of the RSBH+. Ci1 retains a pKa of 7.6 so that both its protonated and deprotonated forms equilibrate at physiological conditions. Protonation of Ci1 leads to a rigid hydrogen-bonding network in the active-site region. This stabilizes the distorted conformation of the retinal after photoactivation and decelerates energy transfer into the protein by impairing the release of the strain energy. In contrast, with deprotonated Ci1 or removal of the Ci1 glutamate side chain, the hydrogen-bonded system is less rigid, and energy transfer by chromophore relaxation is accelerated. Based on the hydrogen out-of-plane (HOOP) band decay kinetics, we determined the activation energy for these processes in dependence of the Ci1 protonation state.