Intramolecular photo-induced charge transfer in visual retinal chromophore mimics: electron density-based indices at the TD-DFT and post-HF levels

Color Tuning in Bovine Rhodopsin through Polarizable Embedding

Bovine rhodopsin is among the most studied proteins in the rhodopsin family. Its primary activation mechanism is the photoisomerization of 11-cis retinal, triggered by the absorption of a UV-visible photon. Different mutants of the same rhodopsin show different absorption wavelengths due to the influence of the specific amino acid residues forming the cavity in which the retinal chromophore is embedded, and rhodopsins activated at different wavelengths are, for example, exploited in the field of optogenetics. In this letter, we present a procedure for systematically investigating color tuning in models of bovine rhodopsin and a set of its mutants embedded in a membrane bilayer. Vertical excitation energy calculations were carried out with the polarizable embedding potential for describing the environment surrounding the chromophore. We show that polarizable embedding outperformed regular electrostatic embedding in determining both the vertical excitation energies and associated oscillator strengths of the systems studied.

Seeking for Optimal Excited States in Photoinduced Electron-Transfer Processes─The Case Study of Brooker's Merocyanine

Material design enters an era in which control of electrons in atoms, molecules, and materials is an essential property to be predicted and thoroughly understood in view of discovering new compounds with properties optimized toward specific optical/optoelectronic applications. π-electronic delocalization and charge separation/recombination enter notably into the set of features that are highly desirable to tailor. Diverse domains are particularly relying on photoinduced electron-transfer (PET), including fields of paramount importance such as energy production through light-harvesting, efficient chemoreceptive sensors, or organic field-effect transistors. In view of completing the arsenal of strategies in this area, we selected Brooker's merocyanine─a typical [D-π-A] compound─as the case study and examined from time-dependent density functional theory the opportunity offered by selected excited states to reach a suited manipulation of the charge transfer (CT) extent. In addition to the consideration of diagnostic tools able to spot the charge amount (i.e., magnitude of electron fraction) transferred upon excitation (q_CT), the spatial extent associated with such an electronic transition or CT length (D_CT), as well as the corresponding variation in dipole moment between the ground and the excited states (μ_CT), further analysis of the excitation process was undertaken. The advantage of going beyond the above-mentioned molecular indicators─which can be considered as PET global indices─was explored on the basis of a partitioning of the electron density. Relevant insight was gained on the relation these global indices have with the evolution of (local) features characterizing either chemical bond or electron delocalization upon vertical excitations.

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.

Predicting the substituent effects in the optical and electrochemical properties of N,N'-substituted isoindigos

Isoindigo, the structural isomer of the well-known dye indigo, has seen a major revival recently because of the increasing interest of its use as a potential drug core structure and for the development of organic photovoltaic materials. Highly beneficial for diverse applications are its facile synthesis, straightforward functionalisation and the broad absorption band in the visible range. Moreover, its intrinsic electron deficiency renders isoindigo a promising acceptor structure in bulk heterojunction architectures. Here we present new insights into the substituent effects of N-functionalised isoindigos, developing a reliable and fast in silico screening approach of a library of compounds. Using experimental UV-Vis and electrochemical data increased the accuracy of the TD-DFT method employed. This procedure allowed us to accurately predict the optical and electrochemical properties of N-functionalised isoindigos and the elucidation of the relationship between substituent effects and electronic properties.

Multi-scale simulation reveals that an amino acid substitution increases photosensitizing reaction inputs in Rhodopsins

Evaluating the availability of molecular oxygen (O₂ ) and energy of excited states in the retinal binding site of rhodopsin is a crucial challenging first step to understand photosensitizing reactions in wild-type (WT) and mutant rhodopsins by absorbing visible light. In the present work, energies of the ground and excited states related to 11-cis-retinal and the O₂ accessibility to the β-ionone ring are evaluated inside WT and human M207R mutant rhodopsins. Putative O₂ pathways within rhodopsins are identified by using molecular dynamics simulations, Voronoi-diagram analysis, and implicit ligand sampling while retinal energetic properties are investigated through density functional theory, and quantum mechanical/molecular mechanical methods. Here, the predictions reveal that an amino acid substitution can lead to enough energy and O₂ accessibility in the core hosting retinal of mutant rhodopsins to favor the photosensitized singlet oxygen generation, which can be useful in understanding retinal degeneration mechanisms and in designing blue-lighting-absorbing proteic photosensitizers.

CpHMD-Then-QM/MM Identification of the Amino Acids Responsible for the Anabaena Sensory Rhodopsin pH-Dependent Electronic Absorption Spectrum

Anabaena Sensory Rhodopsin (ASR), a microbial photoactive protein featuring the retinal chromophore in two different conformations, exhibits a pH-dependent electronic absorption spectrum. Using the recently developed CpHMD-then-QM/MM multiscale protocol applied to ASR embedded in a membrane model, the pH-induced changes in its maximum absorption wavelength have been reproduced and analyzed. While the acidic tiny red-shift is essentially correlated with the deprotonation of an aspartic acid located on the ASR extracellular side, the larger blue-shift experimentally reported at pH values larger than 5 involves a cluster of titrating residues sitting on the cytoplasmic side. The ASR pH-dependent spectrum is the consequence of the competitive stabilization of retinal ground and excited states by the protein electrostatic potential.

TD-M06-2X insights into the absorption and emission spectra of dichlorvos and its molecularly imprinted recognition by methacrylic acid

The absorption and emission spectra of dichlorvos and the dichlorvos-MAA complex in methanol, water, and chloroform in the molecularly imprinted recognition were investigated systematically. The M06-2X results revealed that: 1) the hydroxyl groups in polar solvents such as methanol and water may markedly influence the weak interactions, and then alter the adsorption and emission spectra; 2) the electronic excitation in absorption spectra of dichlorvos is dominated by the configuration HOMO → LUMO, but in the most stable dichlorvos-MAA it becomes the ππ* excitation of HOMO → LUMO + 1; 3) Mulliken charges reveal that dichlorvos almost dissociates to Cl^- and a cation in its S₁ excitation state; 4) the phosphorescence spectra of dichlorvos-MAA are relatively weak. Graphical Abstract The absorption and emission spectra of dichlorvos and the dichlorvos-MAA complex in the molecularly imprinted recognition of dichlorvos were investigated systematically in methanol, water, and chloroform as solvents.