Density Functional Theory Study of the Structure and 13C Chemical Shifts of Retinylidene Iminium Salts

Thermal isomerization rates in retinal analogues using Ab-Initio molecular dynamics

For a detailed understanding of chemical processes in nature and industry, we need accurate models of chemical reactions in complex environments. While Eyring transition state theory is commonly used for modeling chemical reactions, it is most accurate for small molecules in the gas phase. A wide range of alternative rate theories exist that can better capture reactions involving complex molecules and environmental effects. However, they require that the chemical reaction is sampled by molecular dynamics simulations. This is a formidable challenge since the accessible simulation timescales are many orders of magnitude smaller than typical timescales of chemical reactions. To overcome these limitations, rare event methods involving enhanced molecular dynamics sampling are employed. In this work, thermal isomerization of retinal is studied using tight-binding density functional theory. Results from transition state theory are compared to those obtained from enhanced sampling. Rates obtained from dynamical reweighting using infrequent metadynamics simulations were in close agreement with those from transition state theory. Meanwhile, rates obtained from application of Kramers' rate equation to a sampled free energy profile along a torsional dihedral reaction coordinate were found to be up to three orders of magnitude higher. This discrepancy raises concerns about applying rate methods to one-dimensional reaction coordinates in chemical reactions.

Opsin Effect on the Electronic Structure of the Retinylidene Chromophore in Rhodopsin

Direct examination of experimental NMR parameters combined with electronic structure analysis was used to provide a first-principle interpretation of NMR experiments and give a precise evaluation of how the electronic perturbation of the protein environment affects the electronic properties of the retinylidene chromophere in rhodopsin. To this end, we pursued a theoretical analysis using a combination of tools including quantum mechanics/molecular mechanics (QM/MM) at the Density Functional Theory (DFT) level, in conjunction with gauge independent atomic orbital (GIAO) calculations of (13)C NMR chemical shieldings and (1)J(CC) spin-spin coupling constants obtained with the Coupled Perturbed DFT (CPDFT) method. The opsin effect on the retinylidene chromophere is interpreted as an inductive effect of Glu-113 which readjusts the weighting factors of resonance substructures of the conjugated chain of the chromophere. These changes give a rationalization to the alternating effect of the (13)C chemical shifts magnitudes when comparing the retinylidene chromophere in the presence and absence of the protein environment. Conversely, perturbation of π orbitals has little to no effect over (1)J (13)C-(13)C spin-spin coupling constants, as they are mainly dominated by the Fermi contact term, and hence the counteraion effect is restricted to the vicinity of the perturbation. Thus, the apparent contradiction between experimental findings based on chemical shifts (deep penetration) and one-bond J-couplings (localized effects of the protonated Schiff base at the chain terminus) is in fact a consequence of different properties responding differently to the same external perturbation.

Solvent and Protein Effects on the Structure and Dynamics of the Rhodopsin Chromophore

The structure and dynamics of the retinal chromophore of rhodopsin are investigated systematically in different environments (vacuum, methanol solution, and protein binding pocket) and with different computational approaches (classical, quantum, and hybrid quantum mechanics/molecular mechanics (QM/MM) descriptions). Finite temperature effects are taken into account by molecular dynamics simulations. The different components that determine the structure and dynamics of the chromophore in the protein are dissected, both in the dark state and in the early photointermediates. In vacuum and in solution the chromophore displays a very high flexibility, which is significantly reduced by the protein environment. In the 11-cis chromophore, the bond-length alternation, which is correlated with the dipole moment, is found to be similar in solution and in the protein, while it differs greatly with respect to minimum-energy vacuum structures. In the model of the earliest protein photointermediate, the highly twisted chromophore shows a very reduced bond-length alternation.

(1)H and (13)C MAS NMR evidence for pronounced ligand-protein interactions involving the ionone ring of the retinylidene chromophore in rhodopsin

Rhodopsin is a member of the superfamily of G-protein-coupled receptors. This seven alpha-helix transmembrane protein is the visual pigment of the vertebrate rod photoreceptor cells that mediate dim light vision. In the active binding site of this protein the ligand or chromophore, 11-cis-retinal, is covalently bound via a protonated Schiff base to lysine residue 296. Here we present the complete (1)H and (13)C assignments of the 11-cis-retinylidene chromophore in its ligand-binding site determined with ultra high field magic angle spinning NMR. Native bovine opsin was regenerated with 99% enriched uniformly (13)C-labeled 11-cis-retinal. From the labeled pigment, (13)C carbon chemical shifts could be obtained by using two-dimensional radio frequency-driven dipolar recoupling in a solid-state magic angle spinning homonuclear correlation experiment. The (1)H chemical shifts were assigned by two-dimensional heteronuclear ((1)H-(13)C) dipolar correlation spectroscopy with phase-modulated Lee-Goldburg homonuclear (1)H decoupling applied during the t(1) period. The data indicate nonbonding interactions between the protons of the methyl groups of the retinylidene ionone ring and the protein. These nonbonding interactions are attributed to nearby aromatic acid residues Phe-208, Phe-212, and Trp-265 that are in close contact with, respectively, H-16/H-17 and H-18. Furthermore, binding of the chromophore involves a chiral selection of the ring conformation, resulting in equatorial and axial positions for CH(3)-16 and CH(3)-17.