The retinal conformation and its environment in rhodopsin in light of a new 2.2 A crystal structure

Comprehensive understanding of the role of GPER in estrogen receptor-alpha negative breast cancer

G protein-coupled estrogen receptor (GPER) plays a prominent role in facilitating the rapid, non-genomic signaling of estrogens in breast cancer cells. Herein, a comprehensive overview of the role of GPER in ER-ɑ-negative breast cancer is provided. Activation of GPER affected proliferation, metastasis and epithelial mesenchymal transition in ER-ɑ negative breast cancer cells. Clinical studies have demonstrated that GPER positivity was strongly correlated with larger tumor size and advanced clinical stage, suggesting that GPER/ERK signaling may play a role in promoting tumor progression. Strong evidence existed that environmental contaminants like bisphenol A have a carcinogenic potential mediated by GPER activation. The complexity of the cross talk between GPER and other receptors including ER-β, ER-α36, Estrogen-related receptor α (ERRα) and androgen receptor has been discussed. The potential utility of small molecules and phytoestrogens targeting GPER, adds valuable insights into its therapeutic potential. This review holds promises in advancing our understanding of GPER role in ER-ɑ-negative breast cancer. Overall, the consequences of GPER activation are still an area of active research and the implication are not entirely clear.

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.

GPER binding site detection and description: A flavonoid-based docking and molecular dynamics simulations study

Flavonoids, a phenolic compounds class widely distributed in the plant kingdom, have attracted much interest for their implications on several health and disease processes. Usually, the consumption of this type of compounds is approximately 1 g/d, primarily obtained from cereals, chocolate, and dry legumes ensuring its beneficial role in maintaining the homeostasis of the human body. In this context, in cancer disease prominent data points to the role of flavonoids as adjuvant treatment aimed at the regression of the disease. GPER, an estrogen receptor on the cell surface, has been postulated as a probable orchestrator of the beneficial effects of several flavonoids through modulation/inhibition of various mechanisms that lead to cancer progression. Therefore, applying pocket and cavity protein detection and docking and molecular dynamics simulations (MD), we generate, from a cluster composed of 39 flavonoids, crucial insights into the potential role as GPER ligands, of Puerarin, Isoquercetin, Kaempferol 3-O-glucoside and Petunidin 3-O-glucoside, aglycones whose sugar moiety delimits a new described sugar-acceptor sub-cavity into the cavity binding site on the receptor, as well as of the probable activation mechanism of the receptor and the pivotal residues involved in it. Altogether, our results shed light on the potential use of the aforementioned flavonoids as GPER ligands and for further evaluations in in vitro and in vivo assays to elucidate their probable anti-cancer activity.

Comparative quantum-classical dynamics of natural and synthetic molecular rotors show how vibrational synchronization modulates the photoisomerization quantum efficiency

We use quantum-classical trajectories to investigate the origin of the different photoisomerization quantum efficiency observed in the dim-light visual pigment Rhodopsin and in the light-driven biomimetic molecular rotor para-methoxy N-methyl indanylidene-pyrrolinium (MeO-NAIP) in methanol. Our results reveal that effective light-energy conversion requires, in general, an auxiliary molecular vibration (called promoter) that does not correspond to the rotary motion but synchronizes with it at specific times. They also reveal that Nature has designed Rhodopsin to exploit two mechanisms working in a vibrationally coherent regime. The first uses a wag promoter to ensure that ca. 75% of the absorbed photons lead to unidirectional rotations. The second mechanism ensures that the same process is fast enough to avoid directional randomization. It is found that MeO-NAIP in methanol is incapable of exploiting the above mechanisms resulting into a 50% quantum efficiency loss. However, when the solvent is removed, MeO-NAIP rotation is predicted to synchronize with a ring-inversion promoter leading to a 30% increase in quantum efficiency and, therefore, biomimetic behavior.

Graph-based algorithms to dissect long-distance water-mediated H-bond networks for conformational couplings in GPCRs

Changes in structure and dynamics elicited by agonist ligand binding at the extracellular side of G protein coupled receptors (GPCRs) must be relayed to the cytoplasmic G protein binding side of the receptors. To decipher the role of water-mediated hydrogen-bond networks in this relay mechanism, we have developed graph-based algorithms and analysis methodologies applicable to datasets of static structures of distinct GPCRs. For a reference dataset of static structures of bovine rhodopsin solved at the same resolution, we show that graph analyses capture the internal protein-water hydrogen-bond network. The extended analyses of static structures of rhodopsins and opioid receptors suggest a relay mechanism whereby inactive receptors have in place much of the internal core hydrogen-bond network required for long-distance relay of structural change, with extensive local H-bond clusters observed in structures solved at high resolution and with internal water molecules.

Glu1022.53-Mediated Early Conformational Changes in the Process of Light-Induced Green Cone Pigment Activation

Ligand-triggered activation of G protein-coupled receptors (GPCRs) relies on the phenomenon of loose allosteric coupling, which involves conformational alterations spanning from the extracellular ligand-binding domain to the cytoplasmic region, where interactions with G proteins occur. During the GPCR activation process, several intermediate and equilibrium states orchestrate the movement of the flexible and rigid transmembrane (TM) segments of the GPCR. Monitoring early conformational changes is important in unraveling the structural intricacies of the loose allosteric coupling. Here, we focus on the lumi intermediate formed by thermal relaxation from the initial photointermediate, batho in primate green cone pigment (MG), a light-sensitive GPCR responsible for color vision. Our findings from light-induced Fourier transform infrared difference spectroscopy reveal its similarity with rhodopsin, which mediates twilight vision, specifically involving the flip motion of the β-ionone ring, the relaxation of the torsional structure of the retinal, and local perturbations in the α-helix upon lumi intermediate formation. Conversely, we observe a hydrogen bond modification specific to MG's protonated carboxylic acid, identifying its origin as Glu1022.53 situated in TM2. The weakening of the hydrogen bond strength at Glu1022.53 during the transition from the batho to the lumi intermediates corresponds to a slight outward movement of TM2. Additionally, within the X-ray crystal structure of the rhodopsin lumi intermediate, we note the relocation of the Met862.53 side chain in TM2, expanding the volume of the retinal binding pocket. Consequently, the position of 2.53 emerges as the early step in the conformational shift toward light-induced activation. Moreover, given the prevalence of IR-insensitive hydrophobic amino acids at position 2.53 in many rhodopsin-like GPCRs, including rhodopsin, the hydrogen bond alteration in the C═O stretching band at Glu1022.53 of MG can be used as a probe for tracing conformational changes during the GPCR activation process.

Convergent mechanism underlying the acquisition of vertebrate scotopic vision

High sensitivity of scotopic vision (vision in dim light conditions) is achieved by the rods' low background noise, which is attributed to a much lower thermal activation rate (kth) of rhodopsin compared with cone pigments. Frogs and nocturnal geckos uniquely possess atypical rods containing noncanonical cone pigments that exhibit low kth, mimicking rhodopsin. Here, we investigated the convergent mechanism underlying the low kth of rhodopsins and noncanonical cone pigments. Our biochemical analysis revealed that the kth of canonical cone pigments depends on their absorption maximum (λmax). However, rhodopsin and noncanonical cone pigments showed a substantially lower kth than predicted from the λmax dependency. Given that the λmax is inversely proportional to the activation energy of the pigments in the Hinshelwood distribution-based model, our findings suggest that rhodopsin and noncanonical cone pigments have convergently acquired low frequency of spontaneous-activation attempts, including thermal fluctuations of the protein moiety, in the molecular evolutionary processes from canonical cone pigments, which contributes to highly sensitive scotopic vision.

Multiple Ecological Axes Drive Molecular Evolution of Cone Opsins in Beloniform Fishes

Ecological and evolutionary transitions offer an excellent opportunity to examine the molecular basis of adaptation. Fishes of the order Beloniformes include needlefishes, flyingfishes, halfbeaks, and allies, and comprise over 200 species occupying a wide array of habitats-from the marine epipelagic zone to tropical rainforest rivers. These fishes also exhibit a diversity of diets, including piscivory, herbivory, and zooplanktivory. We investigated how diet and habitat affected the molecular evolution of cone opsins, which play a key role in bright light and colour vision and are tightly linked to ecology and life history. We analyzed a targeted-capture dataset to reconstruct the evolutionary history of beloniforms and assemble cone opsin sequences. We implemented codon-based clade models of evolution to examine how molecular evolution was affected by habitat and diet. We found high levels of positive selection in medium- and long-wavelength beloniform opsins, with piscivores showing increased positive selection in medium-wavelength opsins and zooplanktivores showing increased positive selection in long-wavelength opsins. In contrast, short-wavelength opsins showed purifying selection. While marine/freshwater habitat transitions have an effect on opsin molecular evolution, we found that diet plays a more important role. Our study suggests that evolutionary transitions along ecological axes produce complex adaptive interactions that affect patterns of selection on genes that underlie vision.

Picosecond quantum-classical dynamics reveals that the coexistence of light-induced microbial and animal chromophore rotary motion modulates the isomerization quantum yield of heliorhodopsin

Rhodopsins are light-responsive proteins forming two vast and evolutionary distinct superfamilies whose functions are invariably triggered by the photoisomerization of a single retinal chromophore. In 2018 a third widespread superfamily of rhodopsins called heliorhodopsins was discovered using functional metagenomics. Heliorhodopsins, with their markedly different structural features with respect to the animal and microbial superfamilies, offer an opportunity to study how evolution has manipulated the chromophore photoisomerization to achieve adaptation. One question is related to the mechanism of such a reaction and how it differs from that of animal and microbial rhodopsins. To address this question, we use hundreds of quantum-classical trajectories to simulate the spectroscopically documented picosecond light-induced dynamics of a heliorhodopsin from the archaea thermoplasmatales archaeon (TaHeR). We show that, consistently with the observations, the trajectories reveal two excited state decay channels. However, inconsistently with previous hypotheses, only one channel is associated with the -C13C14- rotation of microbial rhodopsins while the second channel is characterized by the -C11C12- rotation typical of animal rhodopsins. The fact that such -C11C12- rotation is aborted upon decay and ground state relaxation, explains why illumination of TaHeR only produces the 13-cis isomer with a low quantum efficiency. We argue that the documented lack of regioselectivity in double-bond excited state twisting motion is the result of an "adaptation" that could be completely lost via specific residue substitutions modulating the steric hindrance experienced along the isomerization motion.

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.

Adaptive Evolution of Nearctic Deepwater Fish Vision: Implications for Assessing Functional Variation for Conservation

Intraspecific functional variation is critical for adaptation to rapidly changing environments. For visual opsins, functional variation can be characterized in vitro and often reflects a species' ecological niche but is rarely considered in the context of intraspecific variation or the impact of recent environmental changes on species of cultural or commercial significance. Investigation of adaptation in postglacial lakes can provide key insight into how rapid environmental changes impact functional evolution. Here, we report evidence for molecular adaptation in vision in 2 lineages of Nearctic fishes that are deep lake specialists: ciscoes and deepwater sculpin. We found depth-related variation in the dim-light visual pigment rhodopsin that evolved convergently in these 2 lineages. In vitro characterization of spectral sensitivity of the convergent deepwater rhodopsin alleles revealed blue-shifts compared with other more widely distributed alleles. These blue-shifted rhodopsin alleles were only observed in deep clear postglacial lakes with underwater visual environments enriched in blue light. This provides evidence of remarkably rapid and convergent visual adaptation and intraspecific functional variation in rhodopsin. Intraspecific functional variation has important implications for conservation, and these fishes are of conservation concern and great cultural, commercial, and nutritional importance to Indigenous communities. We collaborated with the Saugeen Ojibway Nation to develop and test a metabarcoding approach that we show is efficient and accurate in recovering the ecological distribution of functionally relevant variation in rhodopsin. Our approach bridges experimental analyses of protein function and genetics-based tools used in large-scale surveys to better understand the ecological extent of adaptive functional variation.

Osmotic stress studies of G-protein-coupled receptor rhodopsin activation

We summarize and critically review osmotic stress studies of the G-protein-coupled receptor rhodopsin. Although small amounts of structural water are present in these receptors, the effect of bulk water on their function remains uncertain. Studies of the influences of osmotic stress on the GPCR archetype rhodopsin have given insights into the functional role of water in receptor activation. Experimental work has discovered that osmolytes shift the metarhodopsin equilibrium after photoactivation, either to the active or inactive conformations according to their molar mass. At least 80 water molecules are found to enter rhodopsin in the transition to the photoreceptor active state. We infer that this movement of water is both necessary and sufficient for receptor activation. If the water influx is prevented, e.g., by large polymer osmolytes or by dehydration, then the receptor functional transition is back shifted. These findings imply a new paradigm in which rhodopsin becomes solvent swollen in the activation mechanism. Water thus acts as an allosteric modulator of function for rhodopsin-like receptors in lipid membranes.

Xport-A functions as a chaperone by stabilizing the first five transmembrane domains of rhodopsin-1

Rhodopsin-1 (Rh1), the main photosensitive protein of Drosophila, is a seven-transmembrane domain protein, which is inserted co-translationally in the endoplasmic reticulum (ER) membrane. Biogenesis of Rh1 occurs in the ER, where various chaperones interact with Rh1 to aid in its folding and subsequent transport from the ER to the rhabdomere, the light-sensing organelle of the photoreceptors. Xport-A has been proposed as a chaperone/transport factor for Rh1, but the exact molecular mechanism for Xport-A activity upon Rh1 is unknown. Here, we propose a model where Xport-A functions as a chaperone during the biogenesis of Rh1 in the ER by stabilizing the first five transmembrane domains (TMDs) of Rh1.

Hydrophobic interaction between the TM1 and H8 is essential for rhodopsin trafficking to vertebrate photoreceptor outer segments

A major unsolved question in vertebrate photoreceptor biology is the mechanism of rhodopsin transport to the outer segment. In rhodopsin-like class A G protein-coupled receptors, hydrophobic interactions between C-terminal α-helix 8 (H8), and transmembrane α-helix-1 (TM1) have been shown to be important for transport to the plasma membrane, however whether this interaction is important for rhodopsin transport to ciliary rod outer segments is not known. We examined the crystal structures of vertebrate rhodopsins and class A G protein-coupled receptors and found a conserved network of predicted hydrophobic interactions. In Xenopus rhodopsin (xRho), this interaction corresponds to F313, L317, and L321 in H8 and M57, V61, and L68 in TM1. To evaluate the role of H8-TM1 hydrophobic interactions in rhodopsin transport, we expressed xRho-EGFP where hydrophobic residues were mutated in Xenopus rods and evaluated the efficiency of outer segment enrichment. We found that substituting L317 and M57 with hydrophilic residues had the strongest impact on xRho mislocalization. Substituting hydrophilic amino acids at positions L68, F313, and L321 also had a significant impact. Replacing L317 with M resulted in significant mislocalization, indicating that the hydrophobic interaction between residues 317 and 57 is exquisitely sensitive. The corresponding experiment in bovine rhodopsin expressed in HEK293 cells had a similar effect, showing that the H8-TM1 hydrophobic network is essential for rhodopsin transport in mammalian species. Thus, for the first time, we show that a hydrophobic interaction between H8 and TM1 is critical for efficient rhodopsin transport to the vertebrate photoreceptor ciliary outer segment.

Assessment of the Electron Correlation Treatment on the Quantum-Classical Dynamics of Retinal Protonated Schiff Base Models: XMS-CASPT2, RMS-CASPT2, and REKS Methods

We compare the performance of three different multiconfigurational wave function-based electronic structure methods and two implementations of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method. The study is characterized by three features: (i) it uses a small set of quantum-classical trajectories rather than potential energy surface mapping, (ii) it focuses, exclusively, on the photoisomerization of retinal protonated Schiff base models, and (iii) it probes the effect of both methyl substitution and the increase in length of the conjugate π-system. For each tested method, the corresponding analytical gradients are used to drive the quantum-classical (Tully's FSSH method) trajectory propagation, including the recent multistate XMS-CASPT2 and RMS-CASPT2 gradients. It is shown that while CASSCF, XMS-CASPT2, and RMS-CASPT2 yield consistent photoisomerization dynamics descriptions, REKS produces, in some of these systems, qualitatively different behavior that is attributed to a flatter and topographically different excited state potential energy surface. The origin of this behavior can be traced back to the effect of the employed density functional approximation. The above studies are further expanded by benchmarking, at the CASSCF and REKS levels, the electronic structure methods using a QM/MM model of the visual pigment rhodopsin.

Similarities and Differences in Photochemistry of Type I and Type II Rhodopsins

The diversity of the retinal-containing proteins (rhodopsins) in nature is extremely large. Fundamental similarity of the structure and photochemical properties unites them into one family. However, there is still a debate about the origin of retinal-containing proteins: divergent or convergent evolution? In this review, based on the results of our own and literature data, a comparative analysis of the similarities and differences in the photoconversion of the rhodopsin of types I and II is carried out. The results of experimental studies of the forward and reverse photoreactions of the bacteriorhodopsin (type I) and visual rhodopsin (type II) rhodopsins in the femto- and picosecond time scale, photo-reversible reaction of the octopus rhodopsin (type II), photovoltaic reactions, as well as quantum chemical calculations of the forward photoreactions of bacteriorhodopsin and visual rhodopsin are presented. The issue of probable convergent evolution of type I and type II rhodopsins is discussed.

Structure-function analysis suggests that the photoreceptor LITE-1 is a light-activated ion channel

Sensation of light is essential for all organisms. The eye-less nematode Caenorhabditis elegans detects UV and blue light to evoke escape behavior. The photosensor LITE-1 absorbs UV photons with an unusually high extinction coefficient, involving essential tryptophans. Here, we modeled the structure and dynamics of LITE-1 using AlphaFold2-multimer and molecular dynamics (MD) simulations and performed mutational and behavioral assays in C. elegans to characterize its function. LITE-1 resembles olfactory and gustatory receptors from insects, recently shown to be tetrameric ion channels. We identified residues required for channel gating, light absorption, and mechanisms of photo-oxidation, involving a likely binding site for the peroxiredoxin PRDX-2. Furthermore, we identified the binding pocket for a putative chromophore. Several residues lining this pocket have previously been established as essential for LITE-1 function. A newly identified critical cysteine pointing into the pocket represents a likely chromophore attachment site. We derived a model for how photon absorption, via a network of tryptophans and other aromatic amino acids, induces an excited state that is transferred to the chromophore. This evokes conformational changes in the protein, possibly leading to a state receptive to oxidation of cysteines and, jointly, to channel gating. Electrophysiological data support the idea that LITE-1 is a photon and H2O2-coincidence detector. Other proteins with similarity to LITE-1, specifically C. elegans GUR-3, likely use a similar mechanism for photon detection. Thus, a common protein fold and assembly, used for chemoreception in insects, possibly by binding of a particular compound, may have evolved into a light-activated ion channel.

Quantum Mechanics/Molecular mechanics calculations predict A1, not A2, is present in melanopsin (Opn4m) of red-eared slider turtles (Trachemys scripta elegans)

Melanopsin is a photopigment that plays a role in non-visual, light-driven, cellular processes such as modulation of circadian rhythms, retinal vascular development, and the pupillary light reflex (PLR). In this study, computational methods were used to understand which chromophore is harbored by melanopsin in red-eared slider turtles (Trachemys scripta elegans). In mammals, the vitamin A derivative 11-cis-retinal (A1) is the chromophore, which provides functionality for melanopsin. However, in red-eared slider turtles, a member of the reptilian class, the identity of the chromophore remains unclear. Red-eared slider turtles, similar to other freshwater vertebrates, possess visual pigments that harbor a different vitamin A derivative, 11-cis-3,4-didehydroretinal (A2), making their pigments more sensitive to red-light than blue-light, therefore, suggesting the chromophore to be the A2 derivative instead of the A1. To help resolve the chromophore identity, in this work, computational homology models of melanopsin in red-eared slider turtles were first constructed. Next, quantum mechanics/molecular mechanics (QM/MM) calculations were carried out to compare how A1 and A2 derivatives bind to melanopsin. Time dependent density functional theory (TDDFT) calculations were then used to determine the excitation energy of the pigments. Lastly, calculated excitation energies were compared to experimental spectral sensitivity data from responses by the irises of red-eared sliders. Contrary to what was expected, our results suggest that melanopsin in red-eared slider turtles is more likely to harbor the A1 chromophore than the A2. Furthermore, a glutamine (Q622.56) and tyrosine (Y853.28) residue in the chromophore binding pocket are shown to play a role in the spectral tuning of the chromophore.

A blue-shifted anion channelrhodopsin from the Colpodellida alga Vitrella brassicaformis

Microbial rhodopsins, a family of photoreceptive membrane proteins containing the chromophore retinal, show a variety of light-dependent molecular functions. Channelrhodopsins work as light-gated ion channels and are widely utilized for optogenetics, which is a method for controlling neural activities by light. Since two cation channelrhodopsins were identified from the chlorophyte alga Chlamydomonas reinhardtii, recent advances in genomic research have revealed a wide variety of channelrhodopsins including anion channelrhodopsins (ACRs), describing their highly diversified molecular properties (e.g., spectral sensitivity, kinetics and ion selectivity). Here, we report two channelrhodopsin-like rhodopsins from the Colpodellida alga Vitrella brassicaformis, which are phylogenetically distinct from the known channelrhodopsins. Spectroscopic and electrophysiological analyses indicated that these rhodopsins are green- and blue-sensitive pigments (λ_max =  ~ 550 and ~ 440 nm) that exhibit light-dependent ion channeling activities. Detailed electrophysiological analysis revealed that one of them works as a monovalent anion (Cl^-, Br^- and NO₃^-) channel and we named it V. brassicaformis anion channelrhodopsin-2, VbACR2. Importantly, the absorption maximum of VbACR2 (~ 440 nm) is blue-shifted among the known ACRs. Thus, we identified the new blue-shifted ACR, which leads to the expansion of the molecular diversity of ACRs.

Improving the Modeling of Extracellular Ligand Binding Pockets in RosettaGPCR for Conformational Selection

G protein-coupled receptors (GPCRs) are the largest class of drug targets and undergo substantial conformational changes in response to ligand binding. Despite recent progress in GPCR structure determination, static snapshots fail to reflect the conformational space of putative binding pocket geometries to which small molecule ligands can bind. In comparative modeling of GPCRs in the absence of a ligand, often a shrinking of the orthosteric binding pocket is observed. However, the exact prediction of the flexible orthosteric binding site is crucial for adequate structure-based drug discovery. In order to improve ligand docking and guide virtual screening experiments in computer-aided drug discovery, we developed RosettaGPCRPocketSize. The algorithm creates a conformational ensemble of biophysically realistic conformations of the GPCR binding pocket between the TM bundle, which is consistent with a knowledge base of expected pocket geometries. Specifically, tetrahedral volume restraints are defined based on information about critical residues in the orthosteric binding site and their experimentally observed range of Cα-Cα-distances. The output of RosettaGPCRPocketSize is an ensemble of binding pocket geometries that are filtered by energy to ensure biophysically probable arrangements, which can be used for docking simulations. In a benchmark set, pocket shrinkage observed in the default RosettaGPCR was reduced by up to 80% and the binding pocket volume range and geometric diversity were increased. Compared to models from four different GPCR homology model databases (RosettaGPCR, GPCR-Tasser, GPCR-SSFE, and GPCRdb), the here-created models showed more accurate volumes of the orthosteric pocket when evaluated with respect to the crystallographic reference structure. Furthermore, RosettaGPCRPocketSize was able to generate an improved realistic pocket distribution. However, while being superior to other homology models, the accuracy of generated model pockets was comparable to AlphaFold2 models. Furthermore, in a docking benchmark using small-molecule ligands with a higher molecular weight between 400 and 700 Da, a higher success rate in creating native-like binding poses was observed. In summary, RosettaGPCRPocketSize can generate GPCR models with realistic orthosteric pocket volumes, which are useful for structure-based drug discovery applications.

Counterion at an atypical position: Investigating the mechanism of photoisomerization in jellyfish rhodopsin

The position of the counterion in animal rhodopsins plays a crucial role in maintaining visible light sensitivity and facilitating the photoisomerization of their retinal chromophore. The counterion displacement is thought to be closely related to the evolution of rhodopsins, with different positions found in invertebrates and vertebrates. Interestingly, box jellyfish rhodopsin (JelRh) acquired the counterion in transmembrane 2 (TM2) independently. This is a unique feature, as in most animal rhodopsins, the counterion is found in a different location. In this study, we used Fourier Transform Infrared spectroscopy to examine the structural changes that occur in the early photointermediate state of JelRh. We aimed to determine whether the photochemistry of JelRh is similar to that of other animal rhodopsins by comparing its spectra to those of vertebrate bovine rhodopsin (BovRh) and invertebrate squid rhodopsin (SquRh). We observed that the N-D stretching band of the retinal Schiff base was similar to that of BovRh, indicating the interaction between the Schiff base and the counterion is similar in both rhodopsins, despite their different counterion positions. Furthermore, we found that the chemical structure of the retinal in JelRh is similar to that in BovRh, including the changes in the hydrogen-out-of-plane band that indicates a retinal distortion. Overall, the protein conformational changes induced by the photoisomerization of JelRh yielded spectra that resemble an intermediate between BovRh and SquRh, suggesting a unique spectral property of JelRh, and making it the only animal rhodopsin with a counterion in TM2 and an ability to activate G_s protein.

Function and dynamics of the intrinsically disordered carboxyl terminus of β2 adrenergic receptor

Advances in structural biology have provided important mechanistic insights into signaling by the transmembrane core of G-protein coupled receptors (GPCRs); however, much less is known about intrinsically disordered regions such as the carboxyl terminus (CT), which is highly flexible and not visible in GPCR structures. The β2 adrenergic receptor's (β2AR) 71 amino acid CT is a substrate for GPCR kinases and binds β-arrestins to regulate signaling. Here we show that the β2AR CT directly inhibits basal and agonist-stimulated signaling in cell lines lacking β-arrestins. Combining single-molecule fluorescence resonance energy transfer (FRET), NMR spectroscopy, and molecular dynamics simulations, we reveal that the negatively charged β2AR-CT serves as an autoinhibitory factor via interacting with the positively charged cytoplasmic surface of the receptor to limit access to G-proteins. The stability of this interaction is influenced by agonists and allosteric modulators, emphasizing that the CT plays important role in allosterically regulating GPCR activation.

Difference FTIR Spectroscopy of Jumping Spider Rhodopsin-1 at 77 K

Animal visual rhodopsins can be classified into monostable and bistable rhodopsins, which are typically found in vertebrates and invertebrates, respectively. The former example is bovine rhodopsin (BovRh), whose structures and functions have been extensively studied. On the other hand, those of bistable rhodopsins are less known, despite their importance in optogenetics. Here, low-temperature Fourier-transform infrared (FTIR) spectroscopy was applied to jumping spider rhodopsin-1 (SpiRh1) at 77 K, and the obtained light-induced spectral changes were compared with those of squid rhodopsin (SquRh) and BovRh. Although chromophore distortion of the resting state monitored by HOOP vibrations is not distinctive between invertebrate and vertebrate rhodopsins, distortion of the all-trans chromophore after photoisomerization is unique for BovRh, and the distortion was localized at the center of the chromophore in SpiRh1 and SquRh. Highly conserved aspartate (D83 in BovRh) does not change the hydrogen-bonding environment in invertebrate rhodopsins. Thus, present FTIR analysis provides specific structural changes, leading to activation of invertebrate and vertebrate rhodopsins. On the other hand, the analysis of O-D stretching vibrations in D₂O revealed unique features of protein-bound water molecules. Numbers of water bands in SpiRh1 and SquRh were less and more than those in BovRh. The X-ray crystal structure of SpiRh1 observed a bridged water molecule between the protonated Schiff base and its counterion (E194), but strongly hydrogen-bonded water molecules were never detected in SpiRh1, as well as SquRh and BovRh. Thus, absence of strongly hydrogen-bonded water molecules is substantial for animal rhodopsins, which is distinctive from microbial rhodopsins.

Role of Monomer/Tetramer Equilibrium of Rod Visual Arrestin in the Interaction with Phosphorylated Rhodopsin

The phototransduction cascade in vertebrate rod visual cells is initiated by the photoactivation of rhodopsin, which enables the activation of the visual G protein transducin. It is terminated by the phosphorylation of rhodopsin, followed by the binding of arrestin. Here we measured the solution X-ray scattering of nanodiscs containing rhodopsin in the presence of rod arrestin to directly observe the formation of the rhodopsin/arrestin complex. Although arrestin self-associates to form a tetramer at physiological concentrations, it was found that arrestin binds to phosphorylated and photoactivated rhodopsin at 1:1 stoichiometry. In contrast, no complex formation was observed for unphosphorylated rhodopsin upon photoactivation, even at physiological arrestin concentrations, suggesting that the constitutive activity of rod arrestin is sufficiently low. UV-visible spectroscopy demonstrated that the rate of the formation of the rhodopsin/arrestin complex well correlates with the concentration of arrestin monomer rather than the tetramer. These findings indicate that arrestin monomer, whose concentration is almost constant due to the equilibrium with the tetramer, binds to phosphorylated rhodopsin. The arrestin tetramer would act as a reservoir of monomer to compensate for the large changes in arrestin concentration in rod cells caused by intense light or adaptation.

Correction of rhodopsin serial crystallography diffraction intensities for a lattice-translocation defect

Rhodopsin is a G-protein-coupled receptor that detects light and initiates the intracellular signalling cascades that underpin vertebrate vision. Light sensitivity is achieved by covalent linkage to 11-cis retinal, which isomerizes upon photo-absorption. Serial femtosecond crystallography data collected from rhodopsin microcrystals grown in the lipidic cubic phase were used to solve the room-temperature structure of the receptor. Although the diffraction data showed high completeness and good consistency to 1.8 Å resolution, prominent electron-density features remained unaccounted for throughout the unit cell after model building and refinement. A deeper analysis of the diffraction intensities uncovered the presence of a lattice-translocation defect (LTD) within the crystals. The procedure followed to correct the diffraction intensities for this pathology enabled the building of an improved resting-state model. The correction was essential to both confidently model the structure of the unilluminated state and interpret the light-activated data collected after photo-excitation of the crystals. It is expected that similar cases of LTD will be observed in other serial crystallography experiments and that correction will be required in a variety of systems.

Ultrafast structural changes direct the first molecular events of vision

Vision is initiated by the rhodopsin family of light-sensitive G protein-coupled receptors (GPCRs)1. A photon is absorbed by the 11-cis retinal chromophore of rhodopsin, which isomerizes within 200 femtoseconds to the all-trans conformation2, thereby initiating the cellular signal transduction processes that ultimately lead to vision. However, the intramolecular mechanism by which the photoactivated retinal induces the activation events inside rhodopsin remains experimentally unclear. Here we use ultrafast time-resolved crystallography at room temperature3 to determine how an isomerized twisted all-trans retinal stores the photon energy that is required to initiate the protein conformational changes associated with the formation of the G protein-binding signalling state. The distorted retinal at a 1-ps time delay after photoactivation has pulled away from half of its numerous interactions with its binding pocket, and the excess of the photon energy is released through an anisotropic protein breathing motion in the direction of the extracellular space. Notably, the very early structural motions in the protein side chains of rhodopsin appear in regions that are involved in later stages of the conserved class A GPCR activation mechanism. Our study sheds light on the earliest stages of vision in vertebrates and points to fundamental aspects of the molecular mechanisms of agonist-mediated GPCR activation.

Rhodopsin, light-sensor of vision

The light sensor of vertebrate scotopic (low-light) vision, rhodopsin, is a G-protein-coupled receptor comprising a polypeptide chain with bound chromophore, 11-cis-retinal, that exhibits remarkable physicochemical properties. This photopigment is extremely stable in the dark, yet its chromophore isomerises upon photon absorption with 70% efficiency, enabling the activation of its G-protein, transducin, with high efficiency. Rhodopsin's photochemical and biochemical activities occur over very different time-scales: the energy of retinaldehyde's excited state is stored in <1 ps in retinal-protein interactions, but it takes milliseconds for the catalytically active state to form, and many tens of minutes for the resting state to be restored. In this review, we describe the properties of rhodopsin and its role in rod phototransduction. We first introduce rhodopsin's gross structural features, its evolution, and the basic mechanisms of its activation. We then discuss light absorption and spectral sensitivity, photoreceptor electrical responses that result from the activity of individual rhodopsin molecules, and recovery of rhodopsin and the visual system from intense bleaching exposures. We then provide a detailed examination of rhodopsin's molecular structure and function, first in its dark state, and then in the active Meta states that govern its interactions with transducin, rhodopsin kinase and arrestin. While it is clear that rhodopsin's molecular properties are exquisitely honed for phototransduction, from starlight to dawn/dusk intensity levels, our understanding of how its molecular interactions determine the properties of scotopic vision remains incomplete. We describe potential future directions of research, and outline several major problems that remain to be solved.

Molecular Mechanism of Spectral Tuning by Chloride Binding in Monkey Green Sensitive Visual Pigment

The visual pigments of the cones perceive red, green, and blue colors. The monkey green (MG) pigment possesses a unique Cl^- binding site; however, its relationship to the spectral tuning in green pigments remains elusive. Recently, FTIR spectroscopy revealed the characteristic structural modifications of the retinal binding site by Cl^- binding. Herein, we report the computational structural modeling of MG pigments and quantum-chemical simulation to investigate its spectral redshift and physicochemical relevance when Cl^- is present. Our protein structures reflect the previously suggested structural changes. AlphaFold2 failed to predict these structural changes. Excited-state calculations successfully reproduced the experimental red-shifted absorption energies, corroborating our protein structures. Electrostatic energy decomposition revealed that the redshift results from the His197 protonation state and conformations of Glu129, Ser202, and Ala308; however, Cl^- itself contributes to the blueshift. Site-directed mutagenesis supported our analysis. These modeled structures may provide a valuable foundation for studying cone pigments.

Computational investigation of functional water molecules in GPCRs bound to G protein or arrestin

G protein-coupled receptors (GPCRs) are membrane proteins constituting the largest family of drug targets. The activated GPCR binds either the heterotrimeric G proteins or arrestin through its activation cycle. Water molecules have been reported to play a role in GPCR activation. Nevertheless, reported studies are focused on the hydrophobic helical bundle region. How water molecules function in GPCR bound either G protein or arrestin is rarely studied. To address this issue, we carried out computational studies on water molecules in both GPCR/G protein complexes and GPCR/arrestin complexes. Using inhomogeneous fluid theory (IFT), we locate all possible hydration sites in GPCRs binding either to G protein or arrestin. We observe that the number of water molecules on the interaction surface between GPCRs and signal proteins are correlated with the insertion depths of the α5-helix from G-protein or "finger loop" from arrestin in GPCRs. In three out of the four simulation pairs, the interfaces of Rhodopsin, M₂R and NTSR1 in the G protein-associated systems show more water-mediated hydrogen-bond networks when compared to these in arrestin-associated systems. This reflects that more functionally relevant water molecules may probably be attracted in G protein-associated structures than that in arrestin-associated structures. Moreover, we find the water-mediated interaction networks throughout the NPxxY region and the orthosteric pocket, which may be a key for GPCR activation. Reported studies show that non-biased agonist, which can trigger both GPCR-G protein and GPCR-arrestin activation signal, can result in pharmacologically toxicities. Our comprehensive studies of the hydration sites in GPCR/G protein complexes and GPCR/arrestin complexes may provide important insights in the design of G-protein biased agonists.

New simulation insights on the structural transition mechanism of bovine rhodopsin activation

Inactive rhodopsin can absorb photons, which induces different structural transitions that finally activate rhodopsin. We have examined the change in spatial configurations and physicochemical factors that result during the transition mechanism from the inactive to the active rhodopsin state via intermediates. During the activation process, many existing atomic contacts are disrupted, and new ones are formed. This is related to the movement of Helix 5, which tilts away from Helix 3 in the intermediate state in lumirhodopsin and moves closer to Helix 3 again in the active state. Similar patterns of changing atomic contacts are observed between Helices 3 and 5 of the adenosine and neurotensin receptors. In addition, residues 220-238 of rhodopsin, which are disordered in the inactive state, fold in the active state before binding to the Gα, where it catalyzes GDP/GTP exchange on the Gα subunit. Finally, molecular dynamics simulations in the membrane environment revealed that the arrestin binding region adopts a more flexible extended conformation upon phosphorylation, likely promoting arrestin binding and inactivation. In summary, our results provide additional structural understanding of specific rhodopsin activation which might be relevant to other Class A G protein-coupled receptor proteins.

Identification of Small Molecular Chaperones Binding P23H Mutant Opsin through an In Silico Structure-Based Approach

N-terminal P23H opsin mutation accounts for most of retinitis pigmentosa (RP) cases. P23H functions and folding can be rescued by small chaperone ligands, which contributes to validate mutant opsin as a suitable target for pharmacological treatment of RP. However, the lack of structural details on P23H mutant opsin strongly impairs drug design, and new chemotypes of effective chaperones of P23H opsin are in high demand. Here, a computational-boosted workflow combining homology modeling with molecular dynamics (MD) simulations and virtual screening was used to select putative P23H opsin chaperones among different libraries through a structure-based approach. In vitro studies corroborated the reliability of the structural model generated in this work and identified a number of novel chemotypes of safe and effective chaperones able to promote P23H opsin trafficking to the outer cell membrane.

Energy Transport and Its Function in Heptahelical Transmembrane Proteins

Photoproteins such as bacteriorhodopsin (bR) and rhodopsin (Rho) need to effectively dissipate photoinduced excess energy to prevent themselves from damage. Another well-studied seven transmembrane (TM) helices protein is the β₂ adrenergic receptor (β₂AR), a G protein-coupled receptor for which energy dissipation paths have been linked with allosteric communication. To study the vibrational energy transport in the active and inactive states of these proteins, a master equation approach [J. Chem. Phys.2020, 152, 045103] is employed, which uses scaling rules that allow us to calculate energy transport rates solely based on the protein structure. Despite their overall structural similarity, the three 7TM proteins reveal quite different strategies to redistribute excess energy. While bR quickly removes the energy using the TM7 helix as a "lightning rod", Rho exhibits a rather poor energy dissipation, which might eventually require the hydrolysis of the Schiff base between the protein and the retinal chromophore to prevent overheating. Heating the ligand adrenaline of β₂AR, the resulting energy transport network of the protein is found to change significantly upon switching from the active state to the inactive state. While the energy flow may highlight aspects of the inter-residue couplings of β₂AR, it seems not particularly suited to explain allosteric phenomena.

Best Practices in Constant pH MD Simulations: Accuracy and Sampling

Various approaches have been proposed to include the effect of pH in molecular dynamics (MD) simulations. Among these, the λ-dynamics approach proposed by Brooks and co-workers [Kong, X.; Brooks III, C. L. J. Chem. Phys.1996, 105, 2414−2423] can be performed with little computational overhead and hfor each typeence be used to routinely perform MD simulations at microsecond time scales, as shown in the accompanying paper [Aho, N. et al. J. Chem. Theory Comput.2022, DOI: 10.1021/acs.jctc.2c00516]. At such time scales, however, the accuracy of the molecular mechanics force field and the parametrization becomes critical. Here, we address these issues and provide the community with guidelines on how to set up and perform long time scale constant pH MD simulations. We found that barriers associated with the torsions of side chains in the CHARMM36m force field are too high for reaching convergence in constant pH MD simulations on microsecond time scales. To avoid the high computational cost of extending the sampling, we propose small modifications to the force field to selectively reduce the torsional barriers. We demonstrate that with such modifications we obtain converged distributions of both protonation and torsional degrees of freedom and hence consistent pK_a estimates, while the sampling of the overall configurational space accessible to proteins is unaffected as compared to normal MD simulations. We also show that the results of constant pH MD depend on the accuracy of the correction potentials. While these potentials are typically obtained by fitting a low-order polynomial to calculated free energy profiles, we find that higher order fits are essential to provide accurate and consistent results. By resolving problems in accuracy and sampling, the work described in this and the accompanying paper paves the way to the widespread application of constant pH MD beyond pK_a prediction.

Bidirectional Photochemistry of Antarctic Microbial Rhodopsin: Emerging Trend of Ballistic Photoisomerization from the 13-cis Resting State

The decades-long ultrafast examination of nearly a dozen microbial retinal proteins, ion pumps, and sensory photoreceptors has not identified structure-function indicators which predict photoisomerization dynamics, whether it will be sub-picosecond and ballistic or drawn out with complex curve-crossing kinetics. Herein, we report the emergence of such an indicator. Using pH control over retinal isomer ratios, photoinduced transient absorption is recorded in an inward proton pumping Antarctic microbial rhodopsin (AntR) for 13-cis and all-trans retinal resting states. The all-trans fluorescent state decays with 1 ps exponential kinetics. In contrast, in 13-cis it decays within ∼300 fs accompanied by continuous spectral evolution, indicating ballistic internal conversion. The coherent wave packet nature of 13-cis isomerization in AntR matches published results for bacteriorhodopsin (BR) and Anabaena sensory rhodopsin (ASR), which also accommodate both all-trans and 13-cis retinal resting states, marking the emergence of a first structure-photodynamics indicator which holds for all three tested pigments.

Convergent evolutionary counterion displacement of bilaterian opsins in ciliary cells

Opsins are universal photoreceptive proteins in animals. Vertebrate rhodopsin in ciliary photoreceptor cells photo-converts to a metastable active state to regulate cyclic nucleotide signaling. This active state cannot photo-convert back to the dark state, and thus vertebrate rhodopsin is categorized as a mono-stable opsin. By contrast, mollusk and arthropod rhodopsins in rhabdomeric photoreceptor cells photo-convert to a stable active state to stimulate IP3/calcium signaling. This active state can photo-convert back to the dark state, and thus these rhodopsins are categorized as bistable opsins. Moreover, the negatively charged counterion position crucial for the visible light sensitivity is different between vertebrate rhodopsin (Glu113) and mollusk and arthropod rhodopsins (Glu181). This can be explained by an evolutionary scenario where vertebrate rhodopsin newly acquired Glu113 as a counterion, which is thought to have led to higher signaling efficiency of vertebrate rhodopsin. However, the detailed evolutionary steps which led to the higher efficiency in vertebrate rhodopsin still remain unknown. Here, we analyzed the xenopsin group, which is phylogenetically distinct from vertebrate rhodopsin and functions in protostome ciliary cells. Xenopsins are blue-sensitive bistable opsins that regulate cAMP signaling. We found that a bistable xenopsin of Leptochiton asellus had Glu113 as a counterion but did not exhibit elevated signaling efficiency. Therefore, our results show that vertebrate rhodopsin and L. asellus xenopsin regulate cyclic nucleotide signaling in ciliary cells and displaced the counterion position from Glu181 to Glu113 via convergent evolution, whereas subsequently only vertebrate rhodopsin elevated its signaling efficiency by acquiring the mono-stable property.

Conformational selection guides β-arrestin recruitment at a biased G protein-coupled receptor

G protein-coupled receptors (GPCRs) recruit β-arrestins to coordinate diverse cellular processes, but the structural dynamics driving this process are poorly understood. Atypical chemokine receptors (ACKRs) are intrinsically biased GPCRs that engage β-arrestins but not G proteins, making them a model system for investigating the structural basis of β-arrestin recruitment. Here, we performed nuclear magnetic resonance (NMR) experiments on ¹³CH₃-ε-methionine-labeled ACKR3, revealing that β-arrestin recruitment is associated with conformational exchange at key regions of the extracellular ligand-binding pocket and intracellular β-arrestin-coupling region. NMR studies of ACKR3 mutants defective in β-arrestin recruitment identified an allosteric hub in the receptor core that coordinates transitions among heterogeneously populated and selected conformational states. Our data suggest that conformational selection guides β-arrestin recruitment by tuning receptor dynamics at intracellular and extracellular regions.

Rhodopsins: An Excitingly Versatile Protein Species for Research, Development and Creative Engineering

The first member and eponym of the rhodopsin family was identified in the 1930s as the visual pigment of the rod photoreceptor cell in the animal retina. It was found to be a membrane protein, owing its photosensitivity to the presence of a covalently bound chromophoric group. This group, derived from vitamin A, was appropriately dubbed retinal. In the 1970s a microbial counterpart of this species was discovered in an archaeon, being a membrane protein also harbouring retinal as a chromophore, and named bacteriorhodopsin. Since their discovery a photogenic panorama unfolded, where up to date new members and subspecies with a variety of light-driven functionality have been added to this family. The animal branch, meanwhile categorized as type-2 rhodopsins, turned out to form a large subclass in the superfamily of G protein-coupled receptors and are essential to multiple elements of light-dependent animal sensory physiology. The microbial branch, the type-1 rhodopsins, largely function as light-driven ion pumps or channels, but also contain sensory-active and enzyme-sustaining subspecies. In this review we will follow the development of this exciting membrane protein panorama in a representative number of highlights and will present a prospect of their extraordinary future potential.

Conserved hydrogen-bond motifs of membrane transporters and receptors

Membrane transporters and receptors often rely on conserved hydrogen bonds to assemble transient paths for ion transfer or long-distance conformational couplings. For transporters and receptors that use proton binding and proton transfer for function, inter-helical hydrogen bonds of titratable protein sidechains that could change protonation are of central interest to formulate hypotheses about reaction mechanisms. Knowledge of hydrogen bonds common at sites of potential interest for proton binding could thus inform and guide studies on functional mechanisms of protonation-coupled membrane proteins. Here we apply graph-theory approaches to identify hydrogen-bond motifs of carboxylate and histidine sidechains in a large data set of static membrane protein structures. We find that carboxylate-hydroxyl hydrogen bonds are present in numerous structures of the dataset, and can be part of more extended H-bond clusters that could be relevant to conformational coupling. Carboxylate-carboxyamide and imidazole-imidazole hydrogen bonds are represented in comparably fewer protein structures of the dataset. Atomistic simulations on two membrane transporters in lipid membranes suggest that many of the hydrogen bond motifs present in static protein structures tend to be robust, and can be part of larger hydrogen-bond clusters that recruit additional hydrogen bonds.

Hydration-mediated G-protein-coupled receptor activation

Although G-protein–coupled receptors (GPCRs) control vast physiological pathways, their activation remains chemically and physically enigmatic. Our osmotic stress studies of the visual receptor rhodopsin have redefined the standard model of GPCR signaling by revealing the essential role of bulk water. We show results consistent with a large number of water molecules flooding the rhodopsin interior during activation to stabilize the effector binding conformation. These results suggest a model of GPCR activation in which the receptor becomes solvent-swollen upon formation of the active state. We thus demonstrate the mechanism whereby water acts as a powerful allosteric modulator of a pharmacologically important membrane protein family.

The Rhodopsin family of G-protein–coupled receptors (GPCRs) comprises the targets of nearly a third of all pharmaceuticals. Despite structural water present in GPCR X-ray structures, the physiological relevance of these solvent molecules to rhodopsin signaling remains unknown. Here, we show experimental results consistent with the idea that rhodopsin activation in lipid membranes is coupled to bulk water movements into the protein. To quantify hydration changes, we measured reversible shifting of the metarhodopsin equilibrium due to osmotic stress using an extensive series of polyethylene glycol (PEG) osmolytes. We discovered clear evidence that light activation entails a large influx of bulk water (∼80–100 molecules) into the protein, giving insight into GPCR activation mechanisms. Various size polymer osmolytes directly control rhodopsin activation, in which large solutes are excluded from rhodopsin and dehydrate the protein, favoring the inactive state. In contrast, small osmolytes initially forward shift the activation equilibrium until a quantifiable saturation point is reached, similar to gain-of-function protein mutations. For the limit of increasing osmolyte size, a universal response of rhodopsin to osmotic stress is observed, suggesting it adopts a dynamic, hydrated sponge-like state upon photoactivation. Our results demand a rethinking of the role of water dynamics in modulating various intermediates in the GPCR energy landscape. We propose that besides bound water, an influx of bulk water plays a necessary role in establishing the active GPCR conformation that mediates signaling.

A model for how Gβγ couples Gα to GPCR

Representing ∼5% of the human genome, G-protein-coupled receptors (GPCRs) are a primary target for drug discovery; however, the molecular details of how they couple to heterotrimeric G protein subunits are incompletely understood. Here, I propose a hypothetical initial docking model for the encounter between GPCR and Gβγ that is defined by transient interactions between the cytosolic surface of the GPCR and the prenyl moiety and the tripeptide motif, asparagine-proline-phenylalanine (NPF), in the C-terminus of the Gγ subunit. Analysis of class A GPCRs reveals a conserved NPF binding site formed by the interaction of the TM1 and H8. Functional studies using differentially prenylated proteins and peptides further suggest that the intracellular hydrophobic core of the GPCR is a prenyl binding site. Upon binding TM1 and H8 of GPCRs, the propensity of the C-terminal region of Gγ to convert into an α helix allows it to extend into the hydrophobic core of the GPCR, facilitating the GPCR active state. Conservation of the NPF motif in Gγ isoforms and interacting residues in TM1 and H8 suggest that this is a general mechanism of GPCR-G protein signaling. Analysis of the rhodopsin dimer also suggests that Gγ-rhodopsin interactions may facilitate GPCR dimer transactivation.

Structural basis of GPCR coupling to distinct signal transducers: implications for biased signaling

Three classes of G-protein-coupled receptor (GPCR) partners - G proteins, GPCR kinases, and arrestins - preferentially bind active GPCRs. Our analysis suggests that the structures of GPCRs bound to these interaction partners available today do not reveal a clear conformational basis for signaling bias, which would have enabled the rational design of biased GRCR ligands. In view of this, three possibilities are conceivable: (i) there are no generalizable GPCR conformations conducive to binding a particular type of partner; (ii) subtle differences in the orientation of individual residues and/or their interactions not easily detectable in the receptor-transducer structures determine partner preference; or (iii) the dynamics of GPCR binding to different types of partners rather than the structures of the final complexes might underlie transducer bias.

Evolution of the Automatic Rhodopsin Modeling (ARM) Protocol

In recent years, photoactive proteins such as rhodopsins have become a common target for cutting-edge research in the field of optogenetics. Alongside wet-lab research, computational methods are also developing rapidly to provide the necessary tools to analyze and rationalize experimental results and, most of all, drive the design of novel systems. The Automatic Rhodopsin Modeling (ARM) protocol is focused on providing exactly the necessary computational tools to study rhodopsins, those being either natural or resulting from mutations. The code has evolved along the years to finally provide results that are reproducible by any user, accurate and reliable so as to replicate experimental trends. Furthermore, the code is efficient in terms of necessary computing resources and time, and scalable in terms of both number of concurrent calculations as well as features. In this review, we will show how the code underlying ARM achieved each of these properties.

Quantum-classical simulations of rhodopsin reveal excited-state population splitting and its effects on quantum efficiency

The activation of rhodopsin, the light-sensitive G-protein coupled receptor responsible for dim-light vision in vertebrates, is driven by an ultrafast excited-state double-bond isomerization with a quantum efficiency of almost 70%. The origin of such light sensitivity is not understood and a key question is whether in-phase nuclear motion controls the quantum efficiency value. Here, we use hundreds of quantum-classical trajectories to show that, 15 femtoseconds after light absorption, a degeneracy between the reactive excited state and a neighboring state causes the splitting of the rhodopsin population into subpopulations. These subpopulations propagate with different velocities and lead to distinct contributions to the quantum efficiency. We also show that such splitting is modulated by protein electrostatics, thus linking amino-acid sequence variations to quantum efficiency modulation. Finally, we discuss how such a linkage that in principle could be exploited to achieve higher quantum efficiencies, would simultaneously increase the receptor thermal noise leading to a trade-off that may have played a role in rhodopsin evolution.

Creation of photocyclic vertebrate rhodopsin by single amino acid substitution

Opsins are universal photoreceptive proteins in animals and can be classified into three types based on their photoreaction properties. Upon light irradiation, vertebrate rhodopsin forms a metastable active state, which cannot revert back to the original dark state via either photoreaction or thermal reaction. By contrast, after photoreception, most opsins form a stable active state which can photoconvert back to the dark state. Moreover, we recently found a novel type of opsins whose activity is regulated by photocycling. However, the molecular mechanism underlying this diversification of opsins remains unknown. In this study, we showed that vertebrate rhodopsin acquired the photocyclic and photoreversible properties upon introduction of a single mutation at position 188. This revealed that the residue at position 188 contributes to the diversification of photoreaction properties of opsins by its regulation of the recovery from the active state to the original dark state.

Lipids and Phosphorylation Conjointly Modulate Complex Formation of β2-Adrenergic Receptor and β-arrestin2

G protein-coupled receptors (GPCRs) are the largest class of human membrane proteins that bind extracellular ligands at their orthosteric binding pocket to transmit signals to the cell interior. Ligand binding evokes conformational changes in GPCRs that trigger the binding of intracellular interaction partners (G proteins, G protein kinases, and arrestins), which initiate diverse cellular responses. It has become increasingly evident that the preference of a GPCR for a certain intracellular interaction partner is modulated by a diverse range of factors, e.g., ligands or lipids embedding the transmembrane receptor. Here, by means of molecular dynamics simulations of the β₂-adrenergic receptor and β-arrestin2, we study how membrane lipids and receptor phosphorylation regulate GPCR-arrestin complex conformation and dynamics. We find that phosphorylation drives the receptor’s intracellular loop 3 (ICL3) away from a native negatively charged membrane surface to interact with arrestin. If the receptor is embedded in a neutral membrane, the phosphorylated ICL3 attaches to the membrane surface, which widely opens the receptor core. This opening, which is similar to the opening in the G protein-bound state, weakens the binding of arrestin. The loss of binding specificity is manifested by shallower arrestin insertion into the receptor core and higher dynamics of the receptor-arrestin complex. Our results show that receptor phosphorylation and the local membrane composition cooperatively fine-tune GPCR-mediated signal transduction. Moreover, the results suggest that deeper understanding of complex GPCR regulation mechanisms is necessary to discover novel pathways of pharmacological intervention.

New insights into the molecular mechanism of rhodopsin retinitis pigmentosa from the biochemical and functional characterization of G90V, Y102H and I307N mutations

Mutations in the photoreceptor protein rhodopsin are known as one of the leading causes of retinal degeneration in humans. Two rhodopsin mutations, Y102H and I307N, obtained in chemically mutagenized mice, are currently the subject of increased interest as relevant models for studying the process of retinal degeneration in humans. Here, we report on the biochemical and functional characterization of the structural and functional alterations of these two rhodopsin mutants and we compare them with the G90V mutant previously analyzed, as a basis for a better understanding of in vivo studies. This mechanistic knowledge is fundamental to use it for developing novel therapeutic approaches for the treatment of inherited retinal degeneration in retinitis pigmentosa. We find that Y102H and I307N mutations affect the inactive–active equilibrium of the receptor. In this regard, the mutations reduce the stability of the inactive conformation but increase the stability of the active conformation. Furthermore, the initial rate of the functional activation of transducin, by the I307N mutant is reduced, but its kinetic profile shows an unusual increase with time suggesting a profound effect on the signal transduction process. This latter effect can be associated with a change in the flexibility of helix 7 and an indirect effect of the mutation on helix 8 and the C-terminal tail of rhodopsin, whose potential role in the functional activation of the receptor has been usually underestimated. In the case of the Y102H mutant, the observed changes can be associated with conformational alterations affecting the folding of the rhodopsin intradiscal domain, and its presumed involvement in the retinal binding process by the receptor.

Recent Advances in Structure, Function, and Pharmacology of Class A Lipid GPCRs: Opportunities and Challenges for Drug Discovery

Great progress has been made over the past decade in understanding the structural, functional, and pharmacological diversity of lipid GPCRs. From the first determination of the crystal structure of bovine rhodopsin in 2000, much progress has been made in the field of GPCR structural biology. The extraordinary progress in structural biology and pharmacology of GPCRs, coupled with rapid advances in computational approaches to study receptor dynamics and receptor-ligand interactions, has broadened our comprehension of the structural and functional facets of the receptor family members and has helped usher in a modern age of structure-based drug design and development. First, we provide a primer on lipid mediators and lipid GPCRs and their role in physiology and diseases as well as their value as drug targets. Second, we summarize the current advancements in the understanding of structural features of lipid GPCRs, such as the structural variation of their extracellular domains, diversity of their orthosteric and allosteric ligand binding sites, and molecular mechanisms of ligand binding. Third, we close by collating the emerging paradigms and opportunities in targeting lipid GPCRs, including a brief discussion on current strategies, challenges, and the future outlook.

Phospholipid Scrambling by G Protein-Coupled Receptors

Rapid flip-flop of phospholipids across the two leaflets of biological membranes is crucial for many aspects of cellular life. The transport proteins that facilitate this process are classified as pump-like flippases and floppases and channel-like scramblases. Unexpectedly, Class A G protein-coupled receptors (GPCRs), a large class of signaling proteins exemplified by the visual receptor rhodopsin and its apoprotein opsin, are constitutively active as scramblases in vitro. In liposomes, opsin scrambles lipids at a unitary rate of >100,000 per second. Atomistic molecular dynamics simulations of opsin in a lipid membrane reveal conformational transitions that expose a polar groove between transmembrane helices 6 and 7. This groove enables transbilayer lipid movement, conceptualized as the swiping of a credit card (lipid) through a card reader (GPCR). Conformational changes that facilitate scrambling are distinct from those associated with GPCR signaling. In this review, we discuss the physiological significance of GPCR scramblase activity and the modes of its regulation in cells. Expected final online publication date for the Annual Review of Biophysics, Volume 51 is May 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Structure network-based landscape of rhodopsin misfolding by mutations and algorithmic prediction of small chaperone action

Failure of a protein to achieve its functional structural state and normal cellular location contributes to the etiology and pathology of heritable human conformational diseases. The autosomal dominant form of retinitis pigmentosa (adRP) is an incurable blindness largely linked to mutations of the membrane protein rod opsin. While the mechanisms underlying the noxious effects of the mutated protein are not completely understood, a common feature is the functional protein conformational loss. Here, the wild type and 39 adRP rod opsin mutants were subjected to mechanical unfolding simulations coupled to the graph theory-based protein structure network analysis.

A robust computational model was inferred and in vitro validated in its ability to predict endoplasmic reticulum retention of adRP mutants, a feature linked to the mutation-caused misfolding. The structure-based approach could also infer the structural determinants of small chaperone action on misfolded protein mutants with therapeutic implications.

The approach is exportable to conformational diseases linked to missense mutations in any membrane protein.