Protonation states of membrane-embedded carboxylic acid groups in rhodopsin and metarhodopsin II: a Fourier-transform infrared spectroscopy study of site-directed mutants

Generalized QM/MM Force Matching Approach Applied to the 11-cis Protonated Schiff Base Chromophore of Rhodopsin

We extended a previously developed force matching approach to systems with covalent QM/MM boundaries and describe its user-friendly implementation in the publicly available software package CPMD. We applied this approach to the challenging case of the retinal protonated Schiff base in dark state bovine rhodopsin. We were able to develop a highly accurate force field that is able to capture subtle structural changes within the chromophore that have a pronounced influence on the optical properties. The optical absorption spectrum calculated from configurations extracted from a MD trajectory using the new force field is in excellent agreement with QM/MM and experimental references.

GPCR activation: protonation and membrane potential

GPCR proteins represent the largest family of signaling membrane proteins in eukaryotic cells. Their importance to basic cell biology, human diseases, and pharmaceutical interventions is well established. Many crystal structures of GPCR proteins have been reported in both active and inactive conformations. These data indicate that agonist binding alone is not sufficient to trigger the conformational change of GPCRs necessary for binding of downstream G-proteins, yet other essential factors remain elusive. Based on analysis of available GPCR crystal structures, we identified a potential conformational switch around the conserved Asp2.50, which consistently shows distinct conformations between inactive and active states. Combining the structural information with the current literature, we propose an energy-coupling mechanism, in which the interaction between a charge change of the GPCR protein and the membrane potential of the living cell plays a key role for GPCR activation.

Precision vs flexibility in GPCR signaling

The G protein coupled receptor (GPCR) rhodopsin activates the heterotrimeric G protein transducin (Gt) to transmit the light signal into retinal rod cells. The rhodopsin activity is virtually zero in the dark and jumps by more than one billion fold after photon capture. Such perfect switching implies both high fidelity and speed of rhodopsin/Gt coupling. We employed Fourier transform infrared (FTIR) spectroscopy and supporting all-atom molecular dynamics (MD) simulations to study the conformational diversity of rhodopsin in membrane environment and extend the static picture provided by the available crystal structures. The FTIR results show how the equilibria of inactive and active protein states of the receptor (so-called metarhodopsin states) are regulated by the highly conserved E(D)RY and Yx7K(R) motives. The MD data identify an intrinsically unstructured cytoplasmic loop region connecting transmembrane helices 5 and 6 (CL3) and show how each protein state is split into conformational substates. The C-termini of the Gtγ- and Gtα-subunits (GαCT and GγCT), prepared as synthetic peptides, are likely to bind sequentially and at different sites of the active receptor. The peptides have different effects on the receptor conformation. While GγCT stabilizes the active states but preserves CL3 flexibility, GαCT selectively stabilizes a single conformational substate with largely helical CL3, as it is found in crystal structures. Based on these results we propose a mechanism for the fast and precise signal transfer from rhodopsin to Gt, which assumes a stepwise and mutual reduction of their conformational space. The mechanism relies on conserved amino acids and may therefore underlie GPCR/G protein coupling in general.

Rhodopsin Absorption from First Principles: Bypassing Common Pitfalls

Bovine rhodopsin is the most extensively studied retinal protein and is considered the prototype of this important class of photosensitive biosystems involved in the process of vision. Many theoretical investigations have attempted to elucidate the role of the protein matrix in modulating the absorption of retinal chromophore in rhodopsin, but, while generally agreeing in predicting the correct location of the absorption maximum, they often reached contradicting conclusions on how the environment tunes the spectrum. To address this controversial issue, we combine here a thorough structural and dynamical characterization of rhodopsin with a careful validation of its excited-state properties via the use of a wide range of state-of-the-art quantum chemical approaches including various flavors of time-dependent density functional theory (TDDFT), different multireference perturbative schemes (CASPT2 and NEVPT2), and quantum Monte Carlo (QMC) methods. Through extensive quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations, we obtain a comprehensive structural description of the chromophore-protein system and sample a wide range of thermally accessible configurations. We show that, in order to obtain reliable excitation properties, it is crucial to employ a sufficient number of representative configurations of the system. In fact, the common use of a single, ad hoc structure can easily lead to an incorrect model and an agreement with experimental absorption spectra due to cancelation of errors. Finally, we show that, to properly account for polarization effects on the chromophore and to quench the large blue-shift induced by the counterion on the excitation energies, it is necessary to adopt an enhanced description of the protein environment as given by a large quantum region including as many as 250 atoms.

Protein Field Effect on the Dark State of 11-cis Retinal in Rhodopsin by Quantum Monte Carlo/Molecular Mechanics

The accurate determination of the geometrical details of the dark state of 11-cis retinal in rhodopsin represents a fundamental step for the rationalization of the protein role in the optical spectral tuning in the vision mechanism. We have calculated geometries of the full retinal protonated Schiff base chromophore in the gas phase and in the protein environment using the correlated variational Monte Carlo method. The bond length alternation of the conjugated carbon chain of the chromophore in the gas phase shows a significant reduction when moving from the β-ionone ring to the nitrogen, whereas, as expected, the protein environment reduces the electronic conjugation. The proposed dark state structure is fully compatible with solid-state NMR data reported by Carravetta et al. [J. Am. Chem. Soc. 2004, 126, 3948-3953]. TDDFT/B3LYP calculations on such geometries show a blue opsin shift of 0.28 and 0.24 eV induced by the protein for S1 and S2 states, consistently with literature spectroscopic data. The effect of the geometrical distortion alone is a red shift of 0.21 and 0.16 eV with respect to the optimized gas phase chromophore. Our results open new perspectives for the study of the properties of chromophores in their biological environment using correlated methods.

High-resolution crystal structure of human protease-activated receptor 1

Protease-Activated Receptor-1 (PAR1) is the prototypical member of a family of G protein-coupled receptors that mediate cellular responses to thrombin and related proteases. Thrombin irreversibly activates PAR1 by cleaving the N-terminal exodomain of the receptor, which exposes a tethered peptide ligand that binds the receptor’s heptahelical bundle to effect G protein-activation. Here we report a 2.2Å resolution crystal structure of human PAR1 bound to vorapaxar, a PAR1 antagonist. The structure reveals an unusual mode of drug binding that explains how a small molecule binds virtually irreversibly to inhibit receptor activation by PAR1’s tethered ligand. In contrast to deep, solvent-exposed binding pockets observed in other peptide-activated GPCRs, the vorapaxar-binding pocket is superficial but has little surface exposed to the aqueous solvent. PARs are important targets for drug development. The structure reported here will aid development of improved PAR1 antagonists and discovery of antagonists to other members of this receptor family.

Structural transitions of transmembrane helix 6 in the formation of metarhodopsin I

Absorption of light by the visual pigment rhodopsin triggers a rapid cis-trans photoisomerization of its retinal chromophore and a series of conformational changes in both the retinal and protein. The largest structural change is an outward tilt of transmembrane helix H6 that increases the separation of the intracellular ends of H6 and H3 and opens up the G-protein binding site. In the dark state of rhodopsin, Glu247 at the intracellular end of H6 forms a salt bridge with Arg135 on H3 to tether H6 in an inactive conformation. The Arg135-Glu247 interaction is broken in the active state of the receptor, and Arg135 is then stabilized by interactions with Tyr223, Met257, and Tyr306 on helices H5, H6, and H7, respectively. To address the mechanism of H6 motion, solid-state NMR measurements are undertaken of Metarhodopsin I (Meta I), the intermediate preceding the active Metarhodopsin II (Meta II) state of the receptor. (13)C NMR dipolar recoupling measurements reveal an interhelical contact of (13)Cζ-Arg135 with (13)Cε-Met257 in Meta I but not with (13)Cζ-Tyr223 or (13)Cζ-Tyr306. These observations suggest that helix H6 has rotated in the formation of Meta I but that structural changes involving helices H5 and H7 have not yet occurred. Together, our results provide insights into the sequence of events leading up to the outward motion of H6, a hallmark of G protein-coupled receptor activation.

Structure and dynamics of the M3 muscarinic acetylcholine receptor

Acetylcholine (ACh), the first neurotransmitter to be identified¹, exerts many of its physiological actions via activation of a family of G protein-coupled receptors (GPCRs) known as muscarinic ACh receptors (mAChRs). Although the five mAChR subtypes (M₁-M₅) share a high degree of sequence homology, they show pronounced differences in G protein coupling preference and the physiological responses they mediate.^2–4 Unfortunately, despite decades of effort, no therapeutic agents endowed with clear mAChR subtype selectivity have been developed to exploit these differences.^5–6 We describe here the structure of the G_q/11-coupled M₃ mAChR bound to the bronchodilator drug tiotropium and identify the binding mode for this clinically important drug. This structure, together with that of the G_i/o-coupled M₂ receptor, offers new possibilities for the design of mAChR subtype-selective ligands. Importantly, the M₃ receptor structure allows the first structural comparison between two members of a mammalian GPCR subfamily displaying different G-protein coupling selectivities. Furthermore, molecular dynamics simulations suggest that tiotropium binds transiently to an allosteric site en route to the binding pocket of both receptors. These simulations offer a structural view of an allosteric binding mode for an orthosteric GPCR ligand and raise additional opportunities for the design of ligands with different affinities or binding kinetics for different mAChR subtypes. Our findings not only offer new insights into the structure and function of one of the most important GPCR families, but may also facilitate the design of improved therapeutics targeting these critical receptors.

Membrane Protein Simulations Using AMBER Force Field and Berger Lipid Parameters

AMBER force fields are among the most commonly used in molecular dynamics (MD) simulations of proteins. Unfortunately, they lack a specific set of lipid parameters, thus limiting its use in membrane protein simulations. In order to overcome this limitation we assessed whether the widely used united-atom lipid parameters described by Berger and co-workers could be used in conjunction with AMBER force fields in simulations of membrane proteins. Thus, free energies of solvation in water and in cyclohexane, and free energies of water to cyclohexane transfer, were computed by thermodynamic integration procedures for neutral amino acid side-chains employing AMBER99, AMBER03, and OPLS-AA amino acid force fields. In addition, MD simulations of three membrane proteins in a POPC lipid bilayer, the β2 adrenergic G protein-coupled receptor, Aquaporin-1, and the outer membrane protein Omp32, were performed with the aim of comparing the AMBER99SB/Berger combination of force fields with the OPLS-AA/Berger combination. We have shown that AMBER99SB and Berger force fields are compatible, they provide reliable free energy estimations relative to experimental values, and their combination properly describes both membrane and protein structural properties. We then suggest that the AMBER99SB/Berger combination is a reliable choice for the simulation of membrane proteins, which links the easiness of ligand parametrization and the ability to reproduce secondary structure of AMBER99SB force field with the largely validated Berger lipid parameters.

Quantum mechanical/molecular mechanical structure, enantioselectivity, and spectroscopy of hydroxyretinals and insights into the evolution of color vision in small white butterflies

Since Vogt's discovery of A(3)-retinal or 3-hydroxyretinal in insects in 1983 and Matsui's discovery of A(4)-retinal or 4-hydroxyretinal in firefly squid in 1988, hydroxyretinal-protein interactions mediating vision have remained largely unexplored. In the present study, A(3)- and A(4)-retinals are theoretically incorporated into squid and bovine visual pigments by use of the hybrid quantum mechanics/molecular mechanics [SORCI+Q//B3LYP/6-31G(d):Amber96] method, and insights into structure, enantioselectivity, and spectroscopy are gathered and presented for the first time. Contrary to general perception, our findings rule out the formation of a hydrogen bond between the hydroxyl-bearing β-ionone ring portion of retinal and opsin. Compared to A(1)-pigments, A(3)- and A(4)-pigments exhibit slightly blue-shifted absorption maxima due to increase in bond-length alternation of the hydroxyretinal. We suggest that (i) the binding site of firefly squid (Watasenia scintillans) opsin is very similar to that of the Japanese common squid (Todarodes pacificus) opsin; (ii) the molecular mechanism of spectral tuning in small white butterflies involve sites S116 and T185 and breaking of a hydrogen bond between sites E180 and T185; and finally (iii) A(3)-retinal may have occurred during the conversion of A(1)- to A(2)-retinal and insects may have acquired them, in order to absorb light in the blue-green wavelength region and to speed up the G-protein signaling cascade.

Photoexcitation of the blue light using FAD photoreceptor AppA results in ultrafast changes to the protein matrix

Photoexcitation of the flavin chromophore in the BLUF photosensor AppA results in a conformational change that leads to photosensor activation. This conformational change is mediated by a hydrogen-bonding network that surrounds the flavin, and photoexcitation is known to result in changes in the network that include a strengthening of hydrogen bonding to the flavin C4═O carbonyl group. Q63 is a key residue in the hydrogen-bonding network, and replacement of this residue with a glutamate results in a photoinactive mutant. While the ultrafast time-resolved infrared (TRIR) spectrum of Q63E AppA(BLUF) is characterized by flavin carbonyl modes at 1680 and 1650 cm(-1), which are similar in frequency to the analogous modes from the light activated state of the wild-type protein, a band is also observed in the TRIR spectrum at 1724 cm(-1) that is unambiguously assigned to the Q63E carboxylic acid based on U-(13)C labeling of the protein. Light absorption instantaneously (<100 fs) bleaches the 1724 cm(-1) band leading to a transient absorption at 1707 cm(-1). Because Q63E is not part of the isoalloxazine electronic transition, the shift in frequency must arise from a sub picosecond perturbation to the flavin binding pocket. The light-induced change in the frequency of the Q63E side chain is assigned to an increase in hydrogen-bond strength of 3 kcal mol(-1) caused by electronic reorganization of the isoalloxazine ring in the excited state, providing direct evidence that the protein matrix of AppA responds instantaneously to changes in the electronic structure of the chromophore and supporting a model for photoactivation of the wild-type protein that involves initial tautomerization of the Q63 side chain.

The roles of transmembrane domain helix-III during rhodopsin photoactivation

BACKGROUND

Rhodopsin, the prototypic member of G protein-coupled receptors (GPCRs), undergoes isomerization of 11-cis-retinal to all-trans-retinal upon photoactivation. Although the basic mechanism by which rhodopsin is activated is well understood, the roles of whole transmembrane (TM) helix-III during rhodopsin photoactivation in detail are not completely clear.

PRINCIPAL FINDINGS

We herein use single-cysteine mutagenesis technique to investigate conformational changes in TM helices of rhodopsin upon photoactivation. Specifically, we study changes in accessibility and reactivity of cysteine residues introduced into the TM helix-III of rhodopsin. Twenty-eight single-cysteine mutants of rhodopsin (P107C-R135C) were prepared after substitution of all natural cysteine residues (C140/C167/C185/C222/C264/C316) by alanine. The cysteine mutants were expressed in COS-1 cells and rhodopsin was purified after regeneration with 11-cis-retinal. Cysteine accessibility in these mutants was monitored by reaction with 4, 4'-dithiodipyridine (4-PDS) in the dark and after illumination. Most of the mutants except for T108C, G109C, E113C, I133C, and R135C showed no reaction in the dark. Wide variation in reactivity was observed among cysteines at different positions in the sequence 108-135 after photoactivation. In particular, cysteines at position 115, 119, 121, 129, 131, 132, and 135, facing 11-cis-retinal, reacted with 4-PDS faster than neighboring amino acids. The different reaction rates of mutants with 4-PDS after photoactivation suggest that the amino acids in different positions in helix-III are exposed to aqueous environment to varying degrees.

SIGNIFICANCE

Accessibility data indicate that an aqueous/hydrophobic boundary in helix-III is near G109 and I133. The lack of reactivity in the dark and the accessibility of cysteine after photoactivation indicate an increase of water/4-PDS accessibility for certain cysteine-mutants at Helix-III during formation of Meta II. We conclude that photoactivation resulted in water-accessible at the chromophore-facing residues of Helix-III.

Glutamic acid 181 is negatively charged in the bathorhodopsin photointermediate of visual rhodopsin

Assignment of the protonation state of the residue Glu-181 is important to our understanding of the primary event, activation processes and wavelength selection in rhodopsin. Despite extensive study, there is no general agreement on the protonation state of this residue in the literature. Electronic assignment is complicated by the location of Glu-181 near the nodal point in the electrostatic charge shift that accompanies excitation of the chromophore into the low-lying, strongly allowed ππ* state. Thus, the charge on this residue is effectively hidden from electronic spectroscopy. This situation is resolved in bathorhodopsin, because photoisomerization of the chromophore places Glu-181 well within the region of negative charge shift following excitation. We demonstrate that Glu-181 is negatively charged in bathorhodopsin on the basis of the shift in the batho absorption maxima at 10 K [λ(max) band (native) = 544 ± 2 nm, λ(max) band (E181Q) = 556 ± 3 nm] and the decrease in the λ(max) band oscillator strength (0.069 ± 0.004) of E181Q relative to that of the native protein. Because the primary event in rhodopsin does not include a proton translocation or disruption of the hydrogen-bonding network within the binding pocket, we may conclude that the Glu-181 residue in rhodopsin is also charged.

Structure and function of an irreversible agonist-β(2) adrenoceptor complex

G-protein-coupled receptors (GPCRs) are eukaryotic integral membrane proteins that modulate biological function by initiating cellular signalling in response to chemically diverse agonists. Despite recent progress in the structural biology of GPCRs, the molecular basis for agonist binding and allosteric modulation of these proteins is poorly understood. Structural knowledge of agonist-bound states is essential for deciphering the mechanism of receptor activation, and for structure-guided design and optimization of ligands. However, the crystallization of agonist-bound GPCRs has been hampered by modest affinities and rapid off-rates of available agonists. Using the inactive structure of the human β(2) adrenergic receptor (β(2)AR) as a guide, we designed a β(2)AR agonist that can be covalently tethered to a specific site on the receptor through a disulphide bond. The covalent β(2)AR-agonist complex forms efficiently, and is capable of activating a heterotrimeric G protein. We crystallized a covalent agonist-bound β(2)AR-T4L fusion protein in lipid bilayers through the use of the lipidic mesophase method, and determined its structure at 3.5 Å resolution. A comparison to the inactive structure and an antibody-stabilized active structure (companion paper) shows how binding events at both the extracellular and intracellular surfaces are required to stabilize an active conformation of the receptor. The structures are in agreement with long-timescale (up to 30 μs) molecular dynamics simulations showing that an agonist-bound active conformation spontaneously relaxes to an inactive-like conformation in the absence of a G protein or stabilizing antibody.

Structure and function of an irreversible agonist-β(2) adrenoceptor complex

G protein-coupled receptors (GPCRs) are eukaryotic integral membrane proteins that modulate biological function by initiating cellular signaling in response to chemically diverse agonists. Despite recent progress in the structural biology of GPCRs1, the molecular basis for agonist binding and allosteric modulation of these proteins is poorly understood. Structural knowledge of agonist-bound states is essential for deciphering the mechanism of receptor activation, and for structure-guided design and optimization of ligands. However, the crystallization of agonist-bound GPCRs has been hampered by modest affinities and rapid off-rates of available agonists. Using the inactive structure of the human β₂ adrenergic receptor (β₂AR) as a guide, we designed a β₂AR agonist that can be covalently tethered to a specific site on the receptor through a disulfide bond. The covalent β₂AR-agonist complex forms efficiently, and is capable of activating a heterotrimeric G protein. We crystallized a covalent agonist-bound β₂AR-T4L fusion protein in lipid bilayers through the use of the lipidic mesophase method2, and determined its structure at 3.5 Å resolution. A comparison to the inactive structure and an antibody-stabilized active structure (companion paper3) shows how binding events at both the extracellular and intracellular surfaces are required to stabilize an active conformation of the receptor. The structures are in agreement with long-timescale (up to 30 μs) molecular dynamics simulations showing that an agonist-bound active conformation spontaneously relaxes to an inactive-like conformation in the absence of a G protein or stabilizing antibody.

SEIRA spectroscopy on a membrane receptor monolayer using lipoprotein particles as carriers

Surface-enhanced infrared absorption (SEIRA) difference spectroscopy can probe reactions in a protein monolayer tethered to a nanostructured gold surface. SEIRA studies of membrane proteins, however, remain challenging due to sample stability, effects of the metal surface on function, and the need for a membrane-mimicking environment. Here we demonstrate and characterize a model system for membrane receptor investigations using SEIRA spectroscopy. The system employs nanoscale apolipoprotein bound bilayer (NABB) particles, similar to discoidal high-density lipoprotein particles, as soluble carriers for the G-protein-coupled receptor rhodopsin. The His-tag of the engineered apolipoprotein allows for selective binding of the NABBs to a Ni-NTA modified surface, while the lipid environment of the particle ensures stability and protection of the embedded receptor. Using SEIRA spectroscopy, we followed specific binding of rhodopsin-loaded NABB particles to the surface and formation of a membrane protein monolayer. Functionality of the photoreceptor in the immobilized NABBs was probed by SEIRA difference spectroscopy confirming protein conformational changes associated with photoactivation. Orientation of the immobilized NABB particles was assessed by comparing SEIRA data with polarized attenuated total reflection-Fourier-transform infrared spectroscopy. Thus, SEIRA difference spectroscopy supported by the NABB technology provides a promising approach for further functional studies of transmembrane receptors.

The protonation state of Glu181 in rhodopsin revisited: interpretation of experimental data on the basis of QM/MM calculations

The structure and spectroscopy of rhodopsin have been intensely studied in the past decade both experimentally and theoretically; however, important issues still remain unresolved. Of central interest is the protonation state of Glu181, where controversial and contradictory experimental evidence has appeared. While FTIR measurements indicate this residue to be unprotonated, preresonance Raman and UV-vis spectra have been interpreted in favor of a protonated Glu181. Previous computational approaches were not able to resolve this issue, providing contradicting data as well. Here, we perform hybrid QM/MM calculations using DFT methods for the electronic ground state, MRCI methods for the electronically excited states, and a polarization model for the MM part in order to investigate this issue systematically. We constructed various active-site models for protonated as well as unprotonated Glu181, which were evaluated by computing NMR, IR, Raman, and UV-vis spectroscopic data. The resulting differences in the UV-vis and Raman spectra between protonated and unprotonated models are very subtle, which has two major consequences. First, the common interpretation of prior Raman and UV-vis experiments in favor of a neutral Glu181 appears questionable, as it is based on the assumption that a charge at the Glu181 location would have a sizable impact. Second, also theoretical results should be interpreted with care. Spectroscopic differences between the structural models must be related to modeling uncertainties and intrinsic methodological errors. Despite a detailed comparison of various rhodopsins and mutants and consistently favorite results with charged Glu181 models, we find merely weak evidence from UV-vis and Raman calculations. On the contrary, difference FTIR and NMR chemical shift measurements on Rh mutants are indicative of the protonation state of Glu181. Supported by our results, they provide strong and independent evidence for a charged Glu181.

FTIR study of the photoreaction of bovine rhodopsin in the presence of hydroxylamine

In bovine rhodopsin, 11-cis-retinal forms a Schiff base linkage with Lys296. The Schiff base is not reactive to hydroxylamine in the dark, which is consistent with the well-protected retinal binding site. In contrast, under illumination it easily forms all-trans retinal oxime, resulting in the loss of color. This suggests that activation of rhodopsin creates a specific reaction channel for hydroxylamine or loosens the chromophore binding pocket. In the present study, to extract structural information on the Schiff base vicinity and to understand the changes upon activation of rhodopsin, we compared light-induced FTIR difference spectra of bovine rhodopsin in the presence and absence of hydroxylamine under physiological pH (approximately 7). Although the previous FTIR study did not observe the complex formation between rhodopsin and G-protein transducin in hydrated films, the present study clearly shows that hydrated films can be used for studies of the interaction between rhodopsin and hydroxylamine. Hydroxylamine does not react with the Schiff base of Meta-I intermediate trapped at 240 K, possibly because of decreased conformational motions under the frozen environment, while FTIR spectroscopy showed that hydroxylamine affects the hydrogen bonds of the Schiff base and water molecules in Meta-I. In contrast, formation of the retinal oxime was clearly observed at 280 K, the characteristic temperature of Meta-II accumulation in the absence of hydroxylamine, and time-dependent formation of retinal oxime was observed from Meta-II at 265 K as well. The obtained difference FTIR spectra of retinal oxime and opsin are different from that of Meta-II. It is likely that the antiparallel beta-sheet constituting a part of the retinal binding pocket at the extracellular surface is structurally disrupted in the presence of hydroxylamine, which allows the hydrolysis of the Schiff base into retinal oxime.

Sequential rearrangement of interhelical networks upon rhodopsin activation in membranes: the Meta II(a) conformational substate

Photon absorption by rhodopsin is proposed to lead to an activation pathway that is described by the extended reaction scheme Meta I <==>Meta II(a) <==> Meta II(b) <==> Meta II(b)H(+), where Meta II(b)H(+) is thought to be the conformational substate that activates the G protein transducin. Here we test this extended scheme for rhodopsin in a membrane bilayer environment by investigating lipid perturbation of the activation mechanism. We found that symmetric membrane lipids having two unsaturated acyl chains, such as 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), selectively stabilize the Meta II(a) substate in the above mechanism. By combining FTIR and UV-visible difference spectroscopy, we characterized the structural and functional changes involved in the transition to the Meta II(a) intermediate, which links the inactive Meta I intermediate with the Meta II(b) states formed by helix rearrangement. Besides the opening of the Schiff base ionic lock, the Meta II(a) substate is characterized by an activation switch in a conserved water-mediated hydrogen-bonded network involving transmembrane helices H1/H2/H7, which is sensed by its key residue Asp83. On the other hand, movement of retinal toward H5 and its interaction with another interhelical H3/H5 network mediated by His211 and Glu122 is absent in Meta II(a). The latter rearrangement takes place only in the subsequent transition to Meta II(b), which has been previously associated with movement of H6. Our results imply that activating structural changes in the H1/H2/H7 network are triggered by disruption of the Schiff base salt bridge and occur prior to other chromophore-induced changes in the H3/H5 network and the outward tilt of H6 in the activation process.

Structure and activation of the visual pigment rhodopsin

Rhodopsin is a specialized G protein-coupled receptor (GPCR) found in vertebrate rod cells. Absorption of light by its 11-cis retinal chromophore leads to rapid photochemical isomerization and receptor activation. Recent results from protein crystallography and NMR spectroscopy show how structural changes on the extracellular side of rhodopsin induced by retinal isomerization are coupled to the motion of membrane-spanning helices to create a G protein binding pocket on the intracellular side of the receptor. The signaling pathway provides a comprehensive explanation for the conservation of specific amino acids and structural motifs across the class A family of GPCRs, as well as for the conservation of selected residues within the visual receptor subfamily. The emerging model of activation indicates that, rather than being unique, the visual receptors provide a basis for understanding the common structural and dynamic elements in the class A GPCRs.

Mass spectrometry combined with oxidative labeling for exploring protein structure and folding

This review discusses various mass spectrometry (MS)-based approaches for exploring structural aspects of proteins in solution. Electrospray ionization (ESI)-MS, in particular, has found fascinating applications in this area. For example, when used in conjunction with solution-phase hydrogen/deuterium exchange (HDX), ESI-MS is a highly sensitive tool for probing conformational dynamics. The main focus of this article is a technique that is complementary to HDX, that is, the covalent labeling of proteins by hydroxyl radicals. The reactivity of individual amino acid side chains with *OH is strongly affected by their degree of solvent exposure. Thus, analysis of the oxidative labeling pattern by peptide mapping and tandem mass spectrometry provides detailed structural information. A convenient method for *OH production is the photolysis of H(2)O(2) by a pulsed UV laser, resulting in oxidative labeling on the microsecond time scale. Selected examples demonstrate the use of this technique for structural studies on membrane proteins, and the combination with rapid mixing devices for characterizing the properties of short-lived protein (un)folding intermediates.

A lipid pathway for ligand binding is necessary for a cannabinoid G protein-coupled receptor

Recent isothiocyanate covalent labeling studies have suggested that a classical cannabinoid, (-)-7'-isothiocyanato-11-hydroxy-1',1'dimethylheptyl-hexahydrocannabinol (AM841), enters the cannabinoid CB2 receptor via the lipid bilayer (Pei, Y., Mercier, R. W., Anday, J. K., Thakur, G. A., Zvonok, A. M., Hurst, D., Reggio, P. H., Janero, D. R., and Makriyannis, A. (2008) Chem. Biol. 15, 1207-1219). However, the sequence of steps involved in such a lipid pathway entry has not yet been elucidated. Here, we test the hypothesis that the endogenous cannabinoid sn-2-arachidonoylglycerol (2-AG) attains access to the CB2 receptor via the lipid bilayer. To this end, we have employed microsecond time scale all-atom molecular dynamics (MD) simulations of the interaction of 2-AG with CB2 via a palmitoyl-oleoyl-phosphatidylcholine lipid bilayer. Results suggest the following: 1) 2-AG first partitions out of bulk lipid at the transmembrane alpha-helix (TMH) 6/7 interface; 2) 2-AG then enters the CB2 receptor binding pocket by passing between TMH6 and TMH7; 3) the entrance of the 2-AG headgroup into the CB2 binding pocket is sufficient to trigger breaking of the intracellular TMH3/6 ionic lock and the movement of the TMH6 intracellular end away from TMH3; and 4) subsequent to protonation at D3.49/D6.30, further 2-AG entry into the ligand binding pocket results in both a W6.48 toggle switch change and a large influx of water. To our knowledge, this is the first demonstration via unbiased molecular dynamics that a ligand can access the binding pocket of a class A G protein-coupled receptor via the lipid bilayer and the first demonstration via molecular dynamics of G protein-coupled receptor activation triggered by a ligand binding event.

Conserved water-mediated hydrogen bond network between TM-I, -II, -VI, and -VII in 7TM receptor activation

Five highly conserved polar residues connected by a number of structural water molecules together with two rotamer micro-switches, TrpVI:13 and TyrVII:20, constitute an extended hydrogen bond network between the intracellular segments of TM-I, -II, -VI, and -VII of 7TM receptors. Molecular dynamics simulations showed that, although the fewer water molecules in rhodopsin were relatively movable, the hydrogen bond network of the beta2-adrenergic receptor was fully loaded with water molecules that were surprisingly immobilized between the two rotamer switches, both apparently being in their closed conformation. Manipulations of the rotamer state of TyrVII:20 and TrpVI:13 demonstrated that these residues served as gates for the water molecules at the intracellular and extracellular ends of the hydrogen bond network, respectively. TrpVI:13 at the bottom of the main ligand-binding pocket was shown to apparently function as a catching trap for water molecules. Mutational analysis of the beta2-adrenergic receptor demonstrated that the highly conserved polar residues of the hydrogen bond network were all important for receptor signaling but served different functions, some dampening constitutive activity (AsnI:18, AspII:10, and AsnVII:13), whereas others (AsnVII:12 and AsnVII:16) located one helical turn apart and sharing a water molecule were shown to be essential for agonist-induced signaling. It is concluded that the conserved water hydrogen bond network of 7TM receptors constitutes an extended allosteric interface between the transmembrane segments being of crucial importance for receptor signaling and that part of the function of the rotamer micro-switches, TyrVII:20 and TrpVI:13, is to gate or trap the water molecules.

A conserved protonation-induced switch can trigger "ionic-lock" formation in adrenergic receptors

The mechanism of signal transduction in G-protein-coupled receptors (GPCRs) is a crucial step in cell signaling. However, the molecular details of this process are still largely undetermined. Carrying out submicrosecond molecular dynamics simulations of beta-adrenergic receptors, we found that cooperation between a number of highly conserved residues is crucial to alter the equilibrium between the active state and the inactive state of diffusible ligand GPCRs. In particular, "ionic-lock" formation in beta-adrenergic receptors is directly correlated with the protonation state of a highly conserved aspartic acid residue [Asp(2.50)] even though the two sites are located more than 20 A away from each other. Internal polar residues, acting as local microswitches, cooperate to propagate the signal from Asp(2.50) to the G-protein interaction site at the helix III-helix VI interface. Evolutionarily conserved differences between opsin and non-opsin GPCRs in the surrounding of Asp(2.50) influence the acidity of this residue and can thus help in rationalizing the differences in constitutive activity of class A GPCRs.

Light activation of rhodopsin: insights from molecular dynamics simulations guided by solid-state NMR distance restraints

Structural restraints provided by solid-state NMR measurements of the metarhodopsin II intermediate are combined with molecular dynamics simulations to help visualize structural changes in the light activation of rhodopsin. Since the timescale for the formation of the metarhodopsin II intermediate (>1 ms) is beyond that readily accessible by molecular dynamics, we use NMR distance restraints derived from 13C dipolar recoupling measurements to guide the simulations. The simulations yield a working model for how photoisomerization of the 11-cis retinylidene chromophore bound within the interior of rhodopsin is coupled to transmembrane helix motion and receptor activation. The mechanism of activation that emerges is that multiple switches on the extracellular (or intradiscal) side of rhodopsin trigger structural changes that converge to disrupt the ionic lock between helices H3 and H6 on the intracellular side of the receptor.

Vertebrate rhodopsin adaptation to dim light via rapid meta-II intermediate formation

Rhodopsin is a photoreceptive protein present in vertebrate rod photoreceptor cells, which are responsible for scotopic vision. Recent molecular studies have shown that several aquatic vertebrate species have independently acquired rhodopsin containing Asp83Asn, Glu122Gln, and Ala292Ser substitutions, causing a blue shift in the rhodopsin absorption spectra for adaptation to the blue-green photic environment in deep water. Here, we provide new evidence for the evolutionary and functional relevance of the Asp83Asn substitution. Spectroscopic and kinetic analyses of rhodopsins in six cichlid fishes from the East African Great Lakes using charge-coupled device spectrophotometer revealed that the Asp83Asn substitution accelerated the formation of meta-II, a rhodopsin intermediate crucial for activation of the G-protein transducin. Because rapid formation of meta-II likely results in effective transduction of photic signals, it is reasonable to assume that deep-water cichlid species have acquired rhodopsin containing Asn83 to adapt to dim lighting. Remarkably, rhodopsin containing Asn83 has been identified in terrestrial vertebrates such as bats, and these rhodopsin variants also exhibit accelerated meta-II formation. Our results indicated that the Asp83Asn substitution observed in a variety of animal species was acquired independently in many different lineages during vertebrate evolution for adaptation to dimly lit environments.

Phospholipids are needed for the proper formation, stability, and function of the photoactivated rhodopsin-transducin complex

Heterotrimeric G proteins become activated after they form a catalytically active complex with activated G protein-coupled receptors (GPCRs) and GTP replaces GDP on the G protein alpha-subunit. This transient coupling can be stabilized by nucleotide depletion, resulting in an empty-nucleotide G protein-GPCR complex. Efficient and reproducible formation of conformationally homogeneous GPCR-Gt complexes is a prerequisite for structural studies. Herein, we report isolation conditions that enhance the stability and preserve the activity and proper stoichiometry of productive complexes between the purified prototypical GPCR, rhodopsin (Rho), and the rod cell-specific G protein, transducin (Gt). Binding of purified Gt to photoactivated Rho (Rho*) in n-dodecyl beta-D-maltoside (DDM) examined by gel filtration chromatography was generally modest, and purified complexes provided heterogeneous ratios of protein components, most likely because of excess detergent. Rho*-Gt complex stability and activity were greatly increased by addition of phospholipids such as DOPC, DOPE, and DOPS and asolectin to detergent-containing solutions of these proteins. In contrast, native Rho*-Gt complexes purified directly from light-exposed bovine ROS membranes by sucrose gradient centrifugation exhibited improved stability and the expected 2:1 stoichiometry between Rho* and Gt. These results strongly indicate a lipid requirement for stable complex formation in which the likely oligomeric structure of Rho provides a superior platform for coupling to Gt, and phospholipids likely form a matrix to which Gt can anchor through its myristoyl and farnesyl groups. Our findings also demonstrate that the choice of detergent and purification method is critical for obtaining highly purified, stable, and active complexes with appropriate stoichiometry between GPCRs and G proteins needed for structural studies.

Observation of "ionic lock" formation in molecular dynamics simulations of wild-type beta 1 and beta 2 adrenergic receptors

G protein coupled receptors (GPCRs) are a large family of integral membrane proteins involved in signal transduction pathways, making them appealing drug targets for a wide spectrum of diseases. The recently crystallized structures of two engineered adrenergic receptors have opened new avenues for the understanding of the molecular mechanisms of action of GPCRs. Taking the two crystal structures as a starting point, we carried out submicrosecond molecular dynamics simulations of wild-type beta(1) and beta(2) adrenergic receptors in a lipid bilayer under physiological conditions. These simulations give access to structural and dynamic properties of the receptors in pseudo in vivo conditions. For both systems the overall fold properties of the transmembrane region as well as the binding pocket remain close to the crystal structure of the engineered systems, thus suggesting that the ligand binding mode is not affected by the introduced modifications. Both simulations indicate the presence of one or two internal water molecules absent in both crystal structures and essential for the stabilization of the binding pocket at the interface between transmembrane helices III, IV, and V. The different interactions arising from the substitution of Tyr308 in beta(2)AR into Phe325 in beta(1)AR induce different conformations of the homologous Asn(6.55) inside the binding pockets of the two receptors, suggesting a possible origin of receptor specificity in agonist binding. The equilibrated structures of both receptors recover all of the previously suggested features of inactive GPCRs including formation of a salt bridge between the cytoplasmatic moieties of helices III and VI ("ionic lock") that is absent in the crystal structures.

Rhodopsin activation switches in a native membrane environment

The elucidation of structure-function relationships of membrane proteins still poses a considerable challenge due to the sometimes profound influence of the lipid bilayer on the functional properties of the protein. The visual pigment rhodopsin is a prototype of the family of G protein-coupled transmembrane receptors and a considerable part of our knowledge on its activation mechanisms has been derived from studies on detergent-solubilized proteins. This includes in particular the events associated with the conformational transitions of the receptor from the still inactive Meta I to the Meta II photoproduct states, which are involved in signaling. These events involve disruption of an internal salt bridge of the retinal protonated Schiff base, movement of helices and proton uptake from the solvent by the conserved cytoplasmic E(D)RY network around Glu134. As the equilibria associated with these events are considerably altered by the detergent environment, we set out to investigate these equilibria in the native membrane environment and to develop a coherent thermodynamic model of these activating steps using UV-visible and Fourier-transform infrared spectroscopy as complementary techniques. Particular emphasis is put on the role of protonation of Glu134 from the solvent, which is a thermodynamic prerequisite for full receptor activation in membranes, but not in detergent. In view of the conservation of this carboxylate group in family A G protein-coupled receptors, it may also play a similar role in the activation of other family members.

Conserved waters mediate structural and functional activation of family A (rhodopsin-like) G protein-coupled receptors

G protein-coupled receptors with seven transmembrane alpha-helices (GPCRs) comprise the largest receptor superfamily and are involved in detecting a wide variety of extracellular stimuli. The availability of high-resolution crystal structures of five prototypical GPCRs, bovine and squid rhodopsin, engineered A(2A)-adenosine, beta(1)- and beta(2)-adrenergic receptors, permits comparative analysis of features common to these and likely all GPCRs. We provide an analysis of the distribution of water molecules in the transmembrane region of these GPCR structures and find conserved contacts with microdomains demonstrated to be involved in receptor activation. Colocalization of water molecules associating with highly conserved and functionally important residues in several of these GPCR crystal structures supports the notion that these waters are likely to be as important to proper receptor function as the conserved residues. Moreover, in the absence of large conformational changes in rhodopsin after photoactivation, we propose that ordered waters contribute to the functional plasticity needed to transmit activation signals from the retinal-binding pocket to the cytoplasmic face of rhodopsin and that fundamental features of the mechanism of activation, involving these conserved waters, are shared by many if not all family A receptors.

Helix movement is coupled to displacement of the second extracellular loop in rhodopsin activation

The second extracellular loop (EL2) of rhodopsin forms a cap over the binding site of its photoreactive 11-cis retinylidene chromophore. A critical question has been whether EL2 forms a reversible gate that opens upon activation or acts as a rigid barrier. Distance measurements using solid-state ¹³C NMR spectroscopy between the retinal chromophore and the β4 strand of EL2 show the loop is displaced from the retinal binding site upon activation, and there is a rearrangement in the hydrogen-bonding networks connecting EL2 with the extracellular ends of transmembrane helices H4, H5 and H6. NMR measurements further reveal that structural changes in EL2 are coupled to the motion of helix H5 and breaking of the ionic lock that regulates activation. These results provide a comprehensive view of how retinal isomerization triggers helix motion and activation in this prototypical G protein-coupled receptor.

Location of the retinal chromophore in the activated state of rhodopsin*

Rhodopsin is a highly specialized G protein-coupled receptor (GPCR) that is activated by the rapid photochemical isomerization of its covalently bound 11-cis-retinal chromophore. Using two-dimensional solid-state NMR spectroscopy, we defined the position of the retinal in the active metarhodopsin II intermediate. Distance constraints were obtained between amino acids in the retinal binding site and specific (13)C-labeled sites located on the beta-ionone ring, polyene chain, and Schiff base end of the retinal. We show that the retinal C20 methyl group rotates toward the second extracellular loop (EL2), which forms a cap on the retinal binding site in the inactive receptor. Despite the trajectory of the methyl group, we observed an increase in the C20-Gly(188) (EL2) distance consistent with an increase in separation between the retinal and EL2 upon activation. NMR distance constraints showed that the beta-ionone ring moves to a position between Met(207) and Phe(208) on transmembrane helix H5. Movement of the ring toward H5 was also reflected in increased separation between the Cepsilon carbons of Lys(296) (H7) and Met(44) (H1) and between Gly(121) (H3) and the retinal C18 methyl group. Helix-helix interactions involving the H3-H5 and H4-H5 interfaces were also found to change in the formation of metarhodopsin II reflecting increased retinal-protein interactions in the region of Glu(122) (H3) and His(211) (H5). We discuss the location of the retinal in metarhodopsin II and its interaction with sequence motifs, which are highly conserved across the pharmaceutically important class A GPCR family, with respect to the mechanism of receptor activation.

Two protonation switches control rhodopsin activation in membranes

Activation of the G protein-coupled receptor (GPCR) rhodopsin is initiated by light-induced isomerization of the retinal ligand, which triggers 2 protonation switches in the conformational transition to the active receptor state Meta II. The first switch involves disruption of an interhelical salt bridge by internal proton transfer from the retinal protonated Schiff base (PSB) to its counterion, Glu-113, in the transmembrane domain. The second switch consists of uptake of a proton from the solvent by Glu-134 of the conserved E(D)RY motif at the cytoplasmic terminus of helix 3, leading to pH-dependent receptor activation. By using a combination of UV-visible and FTIR spectroscopy, we study the activation mechanism of rhodopsin in different membrane environments and show that these 2 protonation switches become partially uncoupled at physiological temperature. This partial uncoupling leads to approximately 50% population of an entropy-stabilized Meta II state in which the interhelical PSB salt bridge is broken and activating helix movements have taken place but in which Glu-134 remains unprotonated. This partial activation is converted to full activation only by coupling to the pH-dependent protonation of Glu-134 from the solvent, which stabilizes the active receptor conformation by lowering its enthalpy. In a membrane environment, protonation of Glu-134 is therefore a thermodynamic rather than a structural prerequisite for activating helix movements. In light of the conservation of the E(D)RY motif in rhodopsin-like GPCRs, protonation of this carboxylate also may serve a similar function in signal transduction of other members of this receptor family.

Structural basis for ligand binding and specificity in adrenergic receptors: implications for GPCR-targeted drug discovery

Crystal structures of engineered human beta 2-adrenergic receptors (ARs) in complex with an inverse agonist ligand, carazolol, provide three-dimensional snapshots of the disposition of seven transmembrane helices and the ligand-binding site of an important G protein-coupled receptor (GPCR). As expected, beta 2-AR shares substantial structural similarities with rhodopsin, the dim-light photoreceptor of the rod cell. However, although carazolol and the 11- cis-retinylidene moiety of rhodopsin are situated in the same general binding pocket, the second extracellular (E2) loop structures are quite distinct. E2 in rhodopsin shows beta-sheet structure and forms part of the chromophore-binding site. In the beta 2-AR, E2 is alpha-helical and seems to be distinct from the receptor's active site, allowing a potential entry pathway for diffusible ligands. The structures, together with extensive structure-activity relationship (SAR) data from earlier studies, provide insight about possible structural determinants of ligand specificity and how the binding of agonist ligands might alter receptor conformation. We review key features of the new beta 2-AR structures in the context of recent complementary work on the conformational dynamics of GPCRs. We also report 600 ns molecular dynamics simulations that quantified beta 2-AR receptor mobility in a membrane bilayer environment and show how the binding of an agonist ligand, adrenaline (epinephrine), causes conformational changes to the ligand-binding pocket and neighboring helices.

Functional role of the "ionic lock"--an interhelical hydrogen-bond network in family A heptahelical receptors

Activation of family A G-protein-coupled receptors involves a rearrangement of a conserved interhelical cytoplasmic hydrogen bond network between the E(D)RY motif on transmembrane helix 3 (H3) and residues on H6, which is commonly termed the cytoplasmic "ionic lock." Glu134(3.49) of the E(D)RY motif also forms an intrahelical salt bridge with neighboring Arg135(3.50) in the dark-state crystal structure of rhodopsin. We examined the roles of Glu134(3.49) and Arg135(3.50) on H3 and Glu247(6.30) and Glu249(6.32) on H6 on the activation of rhodopsin using Fourier transform infrared spectroscopy of wild-type and mutant pigments reconstituted into lipid membranes. Activation of rhodopsin is pH-dependent with proton uptake during the transition from the inactive Meta I to the active Meta II state. Glu134(3.49) of the ERY motif is identified as the proton-accepting group, using the Fourier transform infrared protonation signature and the absence of a pH dependence of activation in the E134Q mutant. Neutralization of Arg135(3.50) similarly leads to pH-independent receptor activation, but with structural alterations in the Meta II state. Neutralization of Glu247(6.30) and Glu249(6.32) on H6, which are involved in interhelical interactions with H3 and H7, respectively, led to a shift toward Meta II in the E247Q and E249Q mutants while retaining the pH sensitivity of the equilibrium. Disruption of the interhelical interaction of Glu247(6.30) and Glu249(6.32) on H6 with H3 and H7 plays its role during receptor activation, but neutralization of the intrahelical salt bridge between Glu134(3.49) and Arg135(3.50) is considerably more critical for shifting the photoproduct equilibrium to the active conformation. These conclusions are discussed in the context of recent structural data of the beta(2)-adrenergic receptor.