All-cis-retinal and 7-cis,9-cis,11-cis-retinal1

A modified 15 + 5 route led to 7-cis,9-cis,11-cis-retinonitrile in high selectivity and also a lesser amount of all-cis-retinonitrile. DIBAL reduction gave the corresponding retinal isomers. The absorption maximum 287nm of all-cis-retinal reveals the distorted nature of the chromophore.     

Quantum mechanical/molecular mechanical studies on spectral tuning mechanisms of visual pigments and other photoactive proteins

The protein environments surrounding the retinal tune electronic absorption maximum from 350 to 630 nm. Hybrid quantum mechanical/molecular mechanical (QM/MM) methods can be used in calculating excitation energies of retinal in its native protein environments and in studying the molecular basis of spectral tuning. We hereby review recent QM/MM results on the phototransduction of bovine rhodopsin, bacteriorhodopsin, sensory rhodopsin II, nonretinal photoactive yellow protein and their mutants.     

Effects of water re-location and cavity trimming on the CASPT2//CASSCF/AMBER excitation energy of Rhodopsin

We have used quantum-mechanics/molecular-mechanics computations based on ab initio multiconfigurational perturbation theory to determine and rationalize the effect of the re-location of one crystallographic water molecule on the vertical excitation energy of the visual pigment rhodopsin. It is found that the re-location of one water molecule to the opposite side of the 11-cis retinal chromophore leads to a large 0.7–0.8 Å contraction in the chromophore—counterion salt-bridge distance. In spite of this structural effect, the change in excitation energy is found to be limited (< 1.5 kcal mol^?1). Through an analysis of different rhodopsin models in terms of “components” (isolated chromophore, isolated chromophore—counterion ion-pair and models deprived of the counterion charges) we show that the limited change of the excitation energy can be related to a displacement of the retinal chromophore to a different spot of the protein cavity.     

The role of the active site tyrosine in the mechanism of lytic polysaccharide monooxygenase

Catalytic breakdown of polysaccharides can be achieved more efficiently by means of the enzymes lytic polysaccharide monooxygenases (LPMOs). However, the LPMO mechanism has remained controversial, preventing full exploitation of their potential. One of the controversies has centered around an active site tyrosine, present in most LPMO classes. Recent investigations have for the first time obtained direct (spectroscopic) evidence for the possibility of chemical modification of this tyrosine. However, the spectroscopic features obtained in the different investigations are remarkably different, with absorption maximum at 420 and 490 nm, respectively. In this paper we use density functional theory (DFT) in a QM/MM formulation to reconcile these (apparently) conflicting results. By modeling the spectroscopy as well as the underlying reaction mechanism we can show how formation of two isomers (both involving deprotonation of tyrosine) explains the difference in the observed spectroscopic features. Both isomers have a [TyrO–Cu–OH]⁺ moiety with the OH in either the cis- or trans-position to a deprotonated tyrosine. Although the cis-[TyrO–Cu–OH]⁺ moiety is well positioned for oxidation of the substrate, preliminary calculations with the substrate reveal that the reactivity is at best moderate, making a protective role of tyrosine more likely.

With QM/MM, we investigate the mechanism of tyrosine deprotonation in lytic polysaccharide monooxygenases. Our results support deprotonation and our calculated UV-vis spectra show that two isomers must be formed to match recent experiments.     

Computational photochemistry of retinal proteins

High spectral tunability and quantum yield are the striking features of rhodopsin photochemistry. They rely on a strong and complex interaction of their chromophore, the protonated Schiff base of retinal, with its protein environment. In this article, we review the progress in the computational modeling of these systems, focusing on the optical properties and the excited state dynamics. While the earlier success of atomistic theoretical models was based on the breakthrough in X-ray crystallography and combined quantum mechanical molecular mechanical (QM/MM) methodology, recent advances point out the importance of high-level QM methods and the incorporation of effects that are neglected in conventional QM/MM or ONIOM schemes, like polarization and charge transfer.     

Effect of polarization on the opsin shift in rhodopsins. 2. Empirical polarization models for proteins

The explicit treatment of polarization as a many-body interaction in condensed-phase systems represents a current problem in empirical force-field development. Although a variety of efficient models for molecular polarization have been suggested, polarizable force fields are still far from common use nowadays. In this work, we consider interactive polarization models employing Thole's short-range damping scheme and assess them for application on polypeptides. Despite the simplicity of the model, we find mean polarizabilities and anisotropies of amino acid side chains in excellent agreement with MP2/cc-pVQZ benchmark calculations. Combined with restrained electrostatic potential (RESP) derived atomic charges, the models are applied in a quantum-mechanical/molecular-mechanical (QM/MM) approach. An iterative scheme is used to establish a self-consistent mutual polarization between the QM and MM moieties. This ansatz is employed to study the influence of the protein polarizability on calculated optical properties of the protonated Schiff base of retinal in rhodopsin (Rh), bacterio-rhodopsin (bR), and pharaonis sensory rhodopsin II (psRII). The shifts of the excitation energy due to the instantaneous polarization response of the protein to the charge transfer on the retinal chromophore are quantified using the high level ab initio multireference spectroscopy-oriented configuration interaction (SORCI) method. The results are compared with those of previously published QM1/QM2/MM models for bR and psRII.     

Exploring the molecular mechanism for color distinction in humans

We examine here the role of the red, green, and blue human opsin structures in modulating the absorption properties of 11-cis-retinal bonded to the protein via a protonated Schiff base (PSB). We built the three-dimensional structures of the human red, green, and blue opsins using homology modeling techniques with the crystal structure of bovine rhodopsin as the template. We then used quantum mechanics (QM) combined with molecular mechanics (MM) (denoted as QM/MM) techniques in conjunction with molecular dynamics to determine how the room temperature molecular structures of the three human color opsin proteins modulate the absorption frequency of the same bound 11-cis-retinal chromophore to account for the differences in the observed absorption spectra. We find that the conformational twisting of the 11-cis-retinal PSB plays an important role in the green to blue opsin shift, whereas the dipolar side chains in the binding pocket play a surprising role of red-shifting the blue opsin with respect to the green opsin, as a fine adjustment to the opsin shift. The dipolar side chains play a large role in the opsin shift from red to green.     

Photoaffinity labeling of bovine rhodopsin

Photoaffinity labeled (3-diazoacetoxy)-9-cis-retinal (1) and (9-methylenediazoacetoxy)-9-cis-retinal (20) were synthesized and bound to absorption maxima at 465 and 460 nm respectively. Binding studies established that synthetic retinals 1 and 2 bind to the natural binding site and that the integrity of the diazoacetoxy photoaffinity label is preserved in the process. Incorporation of 3-(O¹⁴COCHN₂)-labeled 9-cis retinal could be conveniently carried out in high yield using apomembrane solubilized in CHAPS as detergent to afford the pigment analog in a pure form. Photolysis of the diazoacetoxy group within the binding site led to 15–20%, crosslinking of rhodopsin as estimated by using radiocarbon containing labeled retinal 1 thus showing that this synthetic retinal is suitable for photoaffinity labeling of the active site in rhodopsin. Subsequent experiments to establish the site(s) of crosslinking by sequencing studies will then contribute to our knowledge of the structure of rhodopsin.     

Insight into the common mechanism of the chromophore formation in the red fluorescent proteins: the elusive blue intermediate revealed

Understanding the chromophore maturation process in fluorescent proteins is important for the design of proteins with improved properties. Here, we present the results of electronic structure calculations identifying the nature of a blue intermediate, a key species in the process of the red chromophore formation in DsRed, TagRFP, fluorescent timers, and PAmCherry. The chromophore of the blue intermediate has a structure in which the π-system of the imidazole ring is extended by the acylimine bond, which can be represented by the model N-[(5-hydroxy-1H-imidazole-2yl)methylidene]acetamide (HIMA) compound. Ab initio and QM/MM calculations of the isolated model and protein-bound (mTagBFP) chromophores identify the anionic form of HIMA as the only structure that has absorption that is consistent with the experiment and is stable in the protein binding pocket. The anion and zwitterion are the only protonation forms of HIMA whose absorption (421 and 414 nm, or 2.95 and 3.00 eV) matches the experimental spectrum of the blue form in DsRed (the absorption maximum is 408 nm or 3.04 eV) and mTagBFP (400 nm or 3.10 eV). The QM/MM optimization of the protein-bound anionic form results in a structure that is close to the X-ray one, whereas the zwitterionic chromophore is unstable in the protein binding pocket and undergoes prompt proton transfer. The computed excitation energy of the protein-bound anionic form of the mTagBFP-like chromophore (3.04 eV) agrees with the experimental absorption spectrum of the protein. The DsRed-like chromophore formation in red fluorescent proteins is revisited on the basis of ab initio results and verified by directed mutagenesis revealing a key role of the amino acid residue 70, which is the second after the chromophore tripeptide, in the formation process.     

Modeling of phytochrome absorption spectra

Phytochromes constitute one of the six well‐characterized families of photosensory proteins in Nature. From the viewpoint of computational modeling, however, phytochromes have been the subject of much fewer studies than most other families of photosensory proteins, which is likely a consequence of relevant high‐resolution structural data becoming available only in recent years. In this work, hybrid quantum mechanics/molecular mechanics (QM/MM) methods are used to calculate UV‐vis absorption spectra of Deinococcus radiodurans bacteriophytochrome. We investigate how the choice of QM/MM methodology affects the resulting spectra and demonstrate that QM/MM methods can reproduce the experimental absorption maxima of both the Q and Soret bands with an accuracy of about 0.15 eV. Furthermore, we assess how the protein environment influences the intrinsic absorption of the bilin chromophore, with particular focus on the Q band underlying the primary photochemistry of phytochromes. © 2013 Wiley Periodicals, Inc.     

Matrix-isolation infrared and ultraviolet spectroscopic studies of less stable conformers of 1,3,5-hexatriene

The infrared spectra of the 3-trans and 3-cis isomers of 1,3,5-hexatriene in low-temperature Ar matrices deposited from a high-temperature nozzle or deposited under irradiation of the Hg 253.7 nm light show a number of new bands. Correspondingly, the ultraviolet spectra of the 3-trans isomer observed under similar experimental conditions show a new absorption at 276 nm. Most of these new infrared bands and the new ultraviolet absorption are attributable to less stable isomers which have the planar s-cis conformation [or gauche conformation(s) close to the planar s-cis] about either one or both of the two CC bonds.     

Probing the rhodopsin cavity with reduced retinal models at the CASPT2//CASSCF/AMBER level of theory

We show that the ab initio CASPT2//CASSCF strategy previously used to investigate the ground and excited states of the chromophore of the vision receptor rhodopsin (Rh) in vacuo can be successfully implemented in a QM/MM scheme allowing for CASPT2//CASSCF/AMBER geometry optimization and excited state property evaluation in proteins. Two receptor models (Rh-1 and Rh-2) incorporating different reduced chromophores are investigated. It is shown that Rh-2 features a chromophore equilibrium structure with the correct helicity and a lambdamax that is only 52 nm blue-shifted from the observed value. This result should open the way to a qualitatively correct ab initio QM/MM modeling of the early excited state transient species involved in the vision process.     

Unraveling the similarity of the photoabsorption of deprotonated p-coumaric acid in the gas phase and within the photoactive yellow protein

Using advanced QM/MM methods, the surprisingly negligible shift of the lowest-lying bright electronic excitation of the deprotonated p-coumaric acid (pCA(-)) within the photoactive yellow protein (PYP) is shown to stem from a subtle balance between hypsochromic and bathochromic effects. More specifically, it is found that the change in the excitation energy as a consequence of the disruption of the planarity of pCA(-) inside PYP is nearly canceled out by the shift induced by the intermolecular interactions of the chromophore and the protein as a whole. These results provide important insights about the primary absorption and the tuning of the chromophore by the protein environment in PYP.     

Isomerization mechanism of the HcRed fluorescent protein chromophore

To understand how the protein achieves fluorescence, the isomerization mechanism of the HcRed chromophore is studied both under vacuum and in the solvated red fluorescent protein. Quantum mechanical (QM) and quantum mechanical/molecular mechanical (QM/MM) methods are applied both for the ground and the first excited state. The photoinduced processes in the chromophore mainly involve torsions around the imidazolinone-bridge bond (τ) and the phenoxy-bridge bond (φ). Under vacuum, the isomerization of the cis-trans chromophore essentially proceeds by τ twisting, while the radiationless decay requires φ torsion. By contrast, the isomerization of the cis-trans chromophore in HcRed occurs via simultaneous τ and φ twisting. The protein environment significantly reduces the barrier of this hula twist motion compared with vacuum. The excited-state isomerization barrier via the φ rotation of the cis-coplanar conformer in HcRed is computed to be significantly higher than that of the trans-non-coplanar conformer. This is consistent with the experimental observation that the cis-coplanar-conformation of the chromophore is related to the fluorescent properties of HcRed, while the trans-non-planar conformation is weakly fluorescent or non-fluorescent. Our study shows how the protein modifies the isomerization mechanism, notably by interactions involving the nearby residue Ile197, which keeps the chromophore coplanar and blocks the twisting motion that leads to photoinduced radiationless decay.     

QM/MM Studies on Photoisomerization Dynamics of Azobenzene Chromophore Tethered to a DNA Duplex: Local Unpaired Nucleobase Plays a Crucial Role

The photoresponsive azobenzene‐tethered DNAs have received growing experimental attention because of their potential applications in biotechnology and nanotechnology; however, little is known about the initial photoisomerization of azobenzene in these systems. Herein we have employed quantum mechanics/molecular mechanics (QM/MM) methods to explore the photoisomerization dynamics of an azobenzene‐tethered DNA duplex. We find that in the S₁ state the trans–cis photoisomerization path is much steeper in DNA than in vacuo, which makes the photoisomerization much faster in the DNA environment. This acceleration is primarily caused by complex steric interactions between azobenzene and the nearby unpaired thymine nucleobase, which also change the photoisomerization mechanism of azobenzene in the DNA duplex.     

Spectral tuning in visual pigments: an ONIOM(QM:MM) study on bovine rhodopsin and its mutants

We have investigated geometries and excitation energies of bovine rhodopsin and some of its mutants by hybrid quantum mechanical/molecular mechanical (QM/MM) calculations in ONIOM scheme, employing B3LYP and BLYP density functionals as well as DFTB method for the QM part and AMBER force field for the MM part. QM/MM geometries of the protonated Schiff-base 11- cis-retinal with B3LYP and DFTB are very similar to each other. TD-B3LYP/MM excitation energy calculations reproduce the experimental absorption maximum of 500 nm in the presence of native rhodopsin environment and predict spectral shifts due to mutations within 10 nm, whereas TD-BLYP/MM excitation energies have red-shift error of at least 50 nm. In the wild-type rhodopsin, Glu113 shifts the first excitation energy to blue and accounts for most of the shift found. Other amino acids individually contribute to the first excitation energy but their net effect is small. The electronic polarization effect is essential for reproducing experimental bond length alternation along the polyene chain in protonated Schiff-base retinal, which correlates with the computed first excitation energy. It also corrects the excitation energies and spectral shifts in mutants, more effectively for deprotonated Schiff-base retinal than for the protonated form. The protonation state and conformation of mutated residues affect electronic spectrum significantly. The present QM/MM calculations estimate not only the experimental excitation energies but also the source of spectral shifts in mutants.     

7-Cis isomers of retinal via 7-cis- and 7,9-dicis-β-C18-tetraene ketones : New geometric isomers of vitamin A and carotenoids—I

Based catalyzed condensation of 7-cis- or 7,9-dicis-β-ionylideneacetaldehyde with acetone gave either 7-cis or 7,9-dicis-β-C₁₈-tetraene ketone. Reaction of the tetraene ketone mixture with diethyl cyanomethylphosphonate gave a mixture of 4 isomers of retinonitrile, all believed to have the 7-cis geometry. Partial redution of the retinonitrile mixture with di-isobutylaluminum hydride gave a mixture of retinal isomers. Repeated column chromatography resulted in partial separation of the retinal isomers. Based on NMR data the geometry of the isomers prepared are believed to be 7-cis, 7,13-dicis, and 7,9,13-tricis.     

Hybrid density functional theory/molecular mechanics calculations of two-photon absorption of dimethylamino nitro stilbene in solution

The dimethylamino nitro stilbene (DANS) molecule is studied for exploring solvent effects on two-photon absorption using the quantum mechanical/molecular mechanical (QM/MM) response theory approach, where the quantum part is represented by density functional theory. We have explored the role of geometrical change of the chromophore in solution, the importance of taking a dynamical average over the sampled structures and the role of a granular representation of the polarization and electrostatic interactions with the classically described medium. The line shape function was simulated by the QM/MM technique thereby allowing for non-empirical prediction of the absolute two-photon cross section. We report a maximum in the TPA cross section for a medium of intermediate solvent polarity (i.e. in chloroform) and provide the grounds for an explanation of this effect which recently has been experimentally observed for a series of charge transfer species in solvents of different polarity. The calculations of absorption energies reproduce well the positive solvatochromic behavior of DANS and are in good agreement with experimental spectra available for the chloroform and DMSO solvents. In line with recent development of the QM/MM response technique for color modeling, we find this methodology to offer a versatile tool to predict and analyze two-photon absorption phenomena taking place within a medium.     

Polyethers: Structural analysis of the 1,4,5,8-tetraoxadecalins and 2,2'-bis(1,3-dioxolane) by photoelectron spectroscopy,molecular mechanics and molecular orbital calculations

The photoelectron (PE) spectrum of cis-1,4,5,8-tetraoxadecalin exhibits in contrast to most polyoxa compounds an extremely well resolved low energy pan due to large ‘through-bond’ interactions between the oxygen lone pairs and the ideally oriented C-C bonds. The ‘through-bond’ interactions of the cis- as well as the unknown trans-1,4,5,8-tetraoxadecalin are discussed based on simple perturbation molecular orbital theory and PRDDO molecular orbital calculations. Additionally the PE spectrum of 2,2'-bis (1,3-dioxolane) is reported. Molecular mechanics (MM1 and MM2) and molecular orbital (PRDDO) calculations of the different conformations of the 1,4,5,8-telraoxadecalin and 2,2'-bis(1,3-dioxolane) systems allow the absence of the trans-1, 4,5,8-tetraoxadecalin, and the cis-trans energy difference in the series decalin, 1,8-dioxadecalin and 1,4,5,8-tetraoxadecalin systems to be explained.     

Preparation of sterically hindered geometric isomers of 7-cis-β-ionyl and β-ionylidene derivatives in the vitamin A series

The preparation of the relatively unknown 7-cis isomers of β-ionyl and β-ionylidene derivatives via one way sensitized geometric isomerization is described. These highly sterically crowded isomers, once prepared do not readily thermally revert back to the trans. But the trienes, at temperatures ≥ 100° undergo irreversible cyclization to cyclohexadienes. 7-cis Isomers of β-ionylideneacetone and acetaldehyde as well as higher tetraenes and pentaenes in this series could not be prepared by sensitized isomerization. Partial reduction of 7-cis-β-ionylideneacetonitrile isomers led to the preparation of 7-cis- and 7,9-dicis-β-ionylideneacetaldehyde; and methyl Grignard reaction with the same nitrile gave the triene-methyl ketones.