Conformation and orientation of the retinyl chromophore in rhodopsin: a critical evaluation of recent NMR data on the basis of theoretical calculations results in a minimum energy structure consistent with all experimental data

Cellular and Molecular Mechanisms of Pathogenesis Underlying Inherited Retinal Dystrophies

Inherited retinal dystrophies (IRDs) are congenital retinal degenerative diseases that have various inheritance patterns, including dominant, recessive, X-linked, and mitochondrial. These diseases are most often the result of defects in rod and/or cone photoreceptor and retinal pigment epithelium function, development, or both. The genes associated with these diseases, when mutated, produce altered protein products that have downstream effects in pathways critical to vision, including phototransduction, the visual cycle, photoreceptor development, cellular respiration, and retinal homeostasis. The aim of this manuscript is to provide a comprehensive review of the underlying molecular mechanisms of pathogenesis of IRDs by delving into many of the genes associated with IRD development, their protein products, and the pathways interrupted by genetic mutation.

6-s-cis Conformation and polar binding pocket of the retinal chromophore in the photoactivated state of rhodopsin

The visual pigment rhodopsin is unique among the G protein-coupled receptors in having an 11-cis retinal chromophore covalently bound to the protein through a protonated Schiff base linkage. The chromophore locks the visual receptor in an inactive conformation through specific steric and electrostatic interactions. This efficient inverse agonist is rapidly converted to an agonist, the unprotonated Schiff base of all-trans retinal, upon light activation. Here, we use magic angle spinning NMR spectroscopy to obtain the (13)C chemical shifts (C5-C20) of the all-trans retinylidene chromophore and the (15)N chemical shift of the Schiff base nitrogen in the active metarhodopsin II intermediate. The retinal chemical shifts are sensitive to the conformation of the chromophore and its molecular interactions within the protein-binding site. Comparison of the retinal chemical shifts in metarhodopsin II with those of retinal model compounds reveals that the Schiff base environment is polar. In particular, the (13)C15 and (15)Nepsilon chemical shifts indicate that the C horizontal lineN bond is highly polarized in a manner that would facilitate Schiff base hydrolysis. We show that a strong perturbation of the retinal (13)C12 chemical shift observed in rhodopsin is reduced in wild-type metarhodopsin II and in the E181Q mutant of rhodopsin. On the basis of the T(1) relaxation time of the retinal (13)C18 methyl group and the conjugated retinal (13)C5 and (13)C8 chemical shifts, we have determined that the conformation of the retinal C6-C7 single bond connecting the beta-ionone ring and the retinylidene chain is 6-s-cis in both the inactive and the active states of rhodopsin. These results are discussed within the general framework of ligand-activated G protein-coupled receptors.

Molecular mechanism of long-range synergetic color tuning between multiple amino acid residues in conger rhodopsin

The synergetic effects of multiple rhodopsin mutations on color tuning need to be completely elucidated. Systematic genetic studies and spectroscopy have demonstrated an interesting example of synergetic color tuning between two amino acid residues in conger rhodopsin's ancestral pigment (p501): -a double mutation at one nearby and one distant residue led to a significant λ(max) blue shift of 13 nm, whereas neither of the single mutations at these two sites led to meaningful shifts.To analyze the molecular mechanisms of this synergetic color tuning, we performed homology modeling, molecular simulations, and electronic state calculations. For the double mutant, N195A/A292S, in silico mutation analysis demonstrated conspicuous structural changes in the retinal chromophore, whereas that of the single mutant, A292S, was almost unchanged. Using statistical ensembles of QM/MM optimized structures, the excitation energy of retinal chromophore was evaluated for the three visual pigments. As a result, the λ(max) shift of double mutant (DM) from p501 was -8 nm, while that of single mutant (SM) from p501 was +1 nm. Molecular dynamics simulation for DM demonstrated frequent isomerization between 6-s-cis and 6-s-trans conformers. Unexpectedly, however, the two conformers exhibited almost identical excitation energy, whereas principal component analysis (PCA) identified the retinal-counterion cooperative change of BLA (bond length alternation) and retinal-counterion interaction lead to the shift.

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.

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.

Dynamic structure of retinylidene ligand of rhodopsin probed by molecular simulations

Rhodopsin is currently the only available atomic-resolution template for understanding biological functions of the G protein-coupled receptor (GPCR) family. The structural basis for the phenomenal dark state stability of 11-cis-retinal bound to rhodopsin and its ultrafast photoreaction are active topics of research. In particular, the beta-ionone ring of the retinylidene inverse agonist is crucial for the activation mechanism. We analyzed a total of 23 independent, 100 ns all-atom molecular dynamics simulations of rhodopsin embedded in a lipid bilayer in the microcanonical (N,V,E) ensemble. Analysis of intramolecular fluctuations predicts hydrogen-out-of-plane (HOOP) wagging modes of retinal consistent with those found in Raman vibrational spectroscopy. We show that sampling and ergodicity of the ensemble of simulations are crucial for determining the distribution of conformers of retinal bound to rhodopsin. The polyene chain is rigidly locked into a single, twisted conformation, consistent with the function of retinal as an inverse agonist in the dark state. Most surprisingly, the beta-ionone ring is mobile within its binding pocket; interactions are non-specific and the cavity is sufficiently large to enable structural heterogeneity. We find that retinal occupies two distinct conformations in the dark state, contrary to most previous assumptions. The beta-ionone ring can rotate relative to the polyene chain, thereby populating both positively and negatively twisted 6-s-cis enantiomers. This result, while unexpected, strongly agrees with experimental solid-state (2)H NMR spectra. Correlation analysis identifies the residues most critical to controlling mobility of retinal; we find that Trp265 moves away from the ionone ring prior to any conformational transition. Our findings reinforce how molecular dynamics simulations can challenge conventional assumptions for interpreting experimental data, especially where existing models neglect conformational fluctuations.

Structural analysis and dynamics of retinal chromophore in dark and meta I states of rhodopsin from 2H NMR of aligned membranes

Rhodopsin is a prototype for G protein-coupled receptors (GPCRs) that are implicated in many biological responses in humans. A site-directed (2)H NMR approach was used for structural analysis of retinal within its binding cavity in the dark and pre-activated meta I states. Retinal was labeled with (2)H at the C5, C9, or C13 methyl groups by total synthesis, and was used to regenerate the opsin apoprotein. Solid-state (2)H NMR spectra were acquired for aligned membranes in the low-temperature lipid gel phase versus the tilt angle to the magnetic field. Data reduction assumed a static uniaxial distribution, and gave the retinylidene methyl bond orientations plus the alignment disorder (mosaic spread). The dark-state (2)H NMR structure of 11-cis-retinal shows torsional twisting of the polyene chain and the beta-ionone ring. The ligand undergoes restricted motion, as evinced by order parameters of approximately 0.9 for the spinning C-C(2)H(3) groups, with off-axial fluctuations of approximately 15 degrees . Retinal is accommodated within the rhodopsin binding pocket with a negative pre-twist about the C11=C12 double bond that explains its rapid photochemistry and the trajectory of 11-cis to trans isomerization. In the cryo-trapped meta I state, the (2)H NMR structure shows a reduction of the polyene strain, while torsional twisting of the beta-ionone ring is maintained. Distortion of the retinal conformation is interpreted through substituent control of receptor activation. Steric hindrance between trans retinal and Trp265 can trigger formation of the subsequent activated meta II state. Our results are pertinent to quantum and molecular mechanics simulations of ligands bound to GPCRs, and illustrate how (2)H NMR can be applied to study their biological mechanisms of action.

G-protein coupled receptor structure

Because of their central role in regulation of cellular function, structure/function relationships for G-protein coupled receptors (GPCR) are of vital importance, yet only recently have sufficient data been obtained to begin mapping those relationships. GPCRs regulate a wide range of cellular processes, including the senses of taste, smell, and vision, and control a myriad of intracellular signaling systems in response to external stimuli. Many diseases are linked to GPCRs. A critical need exists for structural information to inform studies on mechanism of receptor action and regulation. X-ray crystal structures of only one GPCR, in an inactive state, have been obtained to date. However considerable structural information for a variety of GPCRs has been obtained using non-crystallographic approaches. This review begins with a review of the very earliest GPCR structural information, mostly derived from rhodopsin. Because of the difficulty in crystallizing GPCRs for X-ray crystallography, the extensive published work utilizing alternative approaches to GPCR structure is reviewed, including determination of three-dimensional structure from sparse constraints. The available X-ray crystallographic analyses on bovine rhodopsin are reviewed as the only available high-resolution structures for any GPCR. Structural information available on ligand binding to several receptors is included. The limited information on excited states of receptors is also reviewed. It is concluded that while considerable basic structural information has been obtained, more data are needed to describe the molecular mechanism of activation of a GPCR.

X-ray crystallographic studies for ligand–protein interaction changes in rhodopsin

G-protein-coupled receptors constitute the largest transmembrane receptor family in human. They are generally activated on binding their specific ligands at the extracellular side of membranes. The signal carried by an agonist is then transmitted to the intracellular side through a conformational change of the receptor, which becomes competent to catalyse GDP/GTP exchange in the α-subunit of heterotrimeric G-protein. Since most of the G-protein-coupled receptors (rhodopsin-like subfamily) share a set of conserved amino acid residues in the transmembrane domain, it is probable that the ligand-triggered activation process involves a common mechanism of rearrangement of the hepta-helical transmembrane bundle. For understanding the nature of this event that is not yet characterized sufficiently, X-ray crystallographic studies of rhodopsin with or without light stimulation can provide valuable information. In rhodopsin, the initial cis–trans photoisomerization of retinal chromophore triggers the structural changes of transmembrane helices. This activation process has been characterized with some spectroscopically distinct photoreaction intermediates (batho, lumi, Meta I and Meta II). With recent advances in the conditions for crystallographic experiments, the diffraction limit of the rhodopsin crystals has been substantially extended. As a result, it becomes possible to detect small structural changes evoked after photoactivation under cryogenic conditions.

QM/MM study of energy storage and molecular rearrangements due to the primary event in vision

The energy storage and the molecular rearrangements due to the primary photochemical event in rhodopsin are investigated by using quantum mechanics/molecular mechanics hybrid methods in conjunction with high-resolution structural data of bovine visual rhodopsin. The analysis of the reactant and product molecular structures reveals the energy storage mechanism as determined by the detailed molecular rearrangements of the retinyl chromophore, including rotation of the (C11-C12) dihedral angle from -11 degrees in the 11-cis isomer to -161 degrees in the all-trans product, where the preferential sense of rotation is determined by the steric interactions between Ala-117 and the polyene chain at the C13 position, torsion of the polyene chain due to steric constraints in the binding pocket, and stretching of the salt bridge between the protonated Schiff base and the Glu-113 counterion by reorientation of the polarized bonds that localize the net positive charge at the Schiff-base linkage. The energy storage, computed at the ONIOM electronic-embedding approach (B3LYP/6-31G*:AMBER) level of theory and the S0-->S1 electronic-excitation energies for the dark and product states, obtained at the ONIOM electronic-embedding approach (TD-B3LYP/6-31G*//B3LYP/6-31G*:AMBER) level of theory, are in very good agreement with experimental data. These results are particularly relevant to the development of a first-principles understanding of the structure-function relations in prototypical G-protein-coupled receptors.

Conformational Changes of G Protein‐Coupled Receptors During Their Activation by Agonist Binding

The superfamily of G protein-coupled receptors (GPCRs) is the largest and most diverse group of transmembrane proteins involved in signal transduction. Many of the over 1000 human GPCRs represent important pharmaceutical targets. However, despite high interest in this receptor family, no high-resolution structure of a human GPCR has been resolved yet. This is mainly due to difficulties in obtaining large quantities of pure and active protein. Until now, only a high-resolution x-ray structure of an inactive state of bovine rhodopsin is available. Since no structure of an active state has been solved, information of the GPCR activation process can be gained only by biophysical techniques. In this review, we first describe what is known about the ground state of GPCRs to then address questions about the nature of the conformational changes taking place during receptor activation and the mechanism controlling the transition from the resting to the active state. Finally, we will also address the question to what extent information about the three-dimensional GPCR structure can be included into pharmaceutical drug design programs.

The crystallographic model of rhodopsin and its use in studies of other G protein-coupled receptors

G protein-coupled receptors (GPCRs) are integral membrane proteins that respond to environmental signals and initiate signal transduction pathways activating cellular processes. Rhodopsin is a GPCR found in rod cells in retina where it functions as a photopigment. Its molecular structure is known from cryo-electron microscopic and X-ray crystallographic studies, and this has reshaped many structure/function questions important in vision science. In addition, this first GPCR structure has provided a structural template for studies of other GPCRs, including many known drug targets. After presenting an overview of the major structural elements of rhodopsin, recent literature covering the use of the rhodopsin structure in analyzing other GPCRs will be summarized. Use of the rhodopsin structural model to understand the structure and function of other GPCRs provides strong evidence validating the structural model.

Molecular dynamics investigation of primary photoinduced events in the activation of rhodopsin

Retinal cis-trans isomerization and early relaxation steps have been studied in a 10-ns molecular dynamics simulation of a fully hydrated model of membrane-embedded rhodopsin. The isomerization, induced by transiently switching the potential energy function governing the C(11)==C(12) dihedral angle of retinal, completes within 150 fs and yields a strongly distorted retinal. The most significant conformational changes in the binding pocket are straightening of retinal's polyene chain and separation of its beta-ionone ring from Trp-265. In the following 500 ps, transition of 6s-cis to 6s-trans retinal and dramatic changes in the hydrogen bonding network of the binding pocket involving the counterion for the protonated Schiff base, Glu-113, occur. Furthermore, the energy initially stored internally in the distorted retinal is transformed into nonbonding interactions of retinal with its environment. During the following 10 ns, increased mobilities of some parts of the protein, such as the kinked regions of the helices, mainly helix VI, and the intracellular loop I2, were observed, as well as transient structural changes involving the conserved salt bridge between Glu-134 and Arg-135. These features prepare the protein for major structural transformations achieved later in the photocycle. Retinal's motion, in particular, can be compared to an opening turnstile freeing the way for the proposed rotation of helix VI. This was demonstrated by a steered molecular dynamics simulation in which an applied torque enforced the rotation of helix VI.

Rhodopsin: insights from recent structural studies

The recent report of the crystal structure of rhodopsin provides insights concerning structure-activity relationships in visual pigments and related G protein-coupled receptors (GPCRs). The seven transmembrane helices of rhodopsin are interrupted or kinked at multiple sites. An extensive network of interhelical interactions stabilizes the ground state of the receptor. The ligand-binding pocket of rhodopsin is remarkably compact, and several chromophore-protein interactions were not predicted from mutagenesis or spectroscopic studies. The helix movement model of receptor activation, which likely applies to all GPCRs of the rhodopsin family, is supported by several structural elements that suggest how light-induced conformational changes in the ligand-binding pocket are transmitted to the cytoplasmic surface. The cytoplasmic domain of the receptor includes a helical domain extending from the seventh transmembrane segment parallel to the bilayer surface. The cytoplasmic surface appears to be approximately large enough to bind to the transducin heterotrimer in a one-to-one complex. The structural basis for several unique biophysical properties of rhodopsin, including its extremely low dark noise level and high quantum efficiency, can now be addressed using a combination of structural biology and various spectroscopic methods. Future high-resolution structural studies of rhodopsin and other GPCRs will form the basis to elucidate the detailed molecular mechanism of GPCR-mediated signal transduction.