Structural insights into Rab21 GTPase activation mechanism by molecular dynamics simulations

Rab proteins belong to the family of monomeric GTPases which are involved in the cellular membrane trafficking. Rab21 protein exists in interchangeable GTP- and GDP-bound states. Rabs switch between two active and inactive conformations like other GTPases. The inactive form of Rab is bound to GDP while its active form is bounded with the GTP. Interexchange between active and inactive form is mediated by the GDP/GTP exchange factor (GEF) which catalyses the conversion from GDP-bound to GTP-bound form, thereby activating the Rab. While the GTP hydrolysis of Rabs is regulated by a GTPase-activating protein (GAP) which causes Rab inactivation. Here, we report the structural flexibility of the Rab21-GTP and Rab21-GDP complexes by docking and molecular dynamics (MD) simulations. Structural analysis of exchange mechanisms of the co-factors complexed with Rab21 reveals that Cys29, Thr33, His48, Gln78 and Lys133 are essentially important in the activation of proteins. Furthermore, a significant change in the orientation of the interacting co-factors, with slight variation in the free energy of binding was observed. Complexation of GEF with Rab21-GTP and Rab21-GDP reveal a flipping of the switchable residues. Finally, 50 ns MD simulations confirm that the GTP-bound Rab21 complex is thermodynamically more favoured than the corresponding GDP-bound complex. This study provides a detailed understanding of the structural elements involved in the conformational changes of Rab21.     

The role of GTP in the activation of adenylate cyclase.

An attempt is made to integrate the knowledge on the role of hormones and guanyl nucleotides in regulating adenylate cyclase into a single molecular model. It is suggested that the hormone catalyzes the activation of the enzyme adenylate cyclase by facilitating the conversion of the enzyme from its inactive state to its active form. The hormone is also responsible for the termination of the signal namely the deactivation of the enzyme by inducing the hydrolysis of GTP at its regulatory site. The relative rates of these two processes determine the steady state concentration of the active form of the enzyme. The model also explains the difference in behaviour between GTP and its non-hydrolyzable analogs GppNHp and GTPγS.     

Relation between the conformational heterogeneity and reaction cycle of Ras: molecular simulation of Ras

Ras functions as a molecular switch by cycling between the active GTP-bound state and the inactive GDP-bound state. It is known experimentally that there is another GTP-bound state called state 1. We investigate the conformational changes and fluctuations arising from the difference in the coordinations between the switch regions and ligands in the GTP- and GDP-bound states using a total of 830 ns of molecular-dynamics simulations. Our results suggest that the large fluctuations among multiple conformations of switch I in state 1 owing to the absence of coordination between Thr-35 and Mg²⁺ inhibit the binding of Ras to effectors. Furthermore, we elucidate the conformational heterogeneity in Ras by using principal component analysis, and propose a two-step reaction path from the GDP-bound state to the active GTP-bound state via state 1. This study suggests that state 1 plays an important role in signal transduction as an intermediate state of the nucleotide exchange process, although state 1 itself is an inactive state for signal transduction.     

How does activation loop phosphorylation modulate catalytic activity in the cAMP-dependent protein kinase: a theoretical study

Phosphorylation mediates the function of many proteins and enzymes. In the catalytic subunit of cAMP-dependent protein kinase, phosphorylation of Thr 197 in the activation loop strongly influences its catalytic activity. In order to provide theoretical understanding about this important regulatory process, classical molecular dynamics simulations and ab initio QM/MM calculations have been carried out on the wild-type PKA-Mg(2) ATP-substrate complex and its dephosphorylated mutant, T197A. It was found that pThr 197 not only facilitates the phosphoryl transfer reaction by stabilizing the transition state through electrostatic interactions but also strongly affects its essential protein dynamics as well as the active site conformation.     

Coarse-grained dynamics of the receiver domain of NtrC: fluctuations, correlations and implications for allosteric cooperativity

Receiver domains are key molecular switches in bacterial signaling. Structural studies have shown that the receiver domain of the nitrogen regulatory protein C (NtrC) exists in a conformational equilibrium encompassing both inactive and active states, with phosphorylation of Asp54 allosterically shifting the equilibrium towards the active state. To analyze dynamical fluctuations and correlations in NtrC as it undergoes activation, we have applied a coarse-grained dynamics algorithm using elastic network models. Normal mode analysis reveals possible dynamical pathways for the transition of NtrC from the inactive state to the active state. The diagonalized correlation between the inactive and the active (phosphorylated) state shows that most correlated motions occur around the active site of Asp54 and in the region Thr82 to Tyr101. This indicates a coupled correlation of dynamics in the "Thr82-Tyr101" motion. With phosphorylation inducing significant flexibility changes around the active site and alpha3 and alpha4 helices, we find that this activation makes the active-site region and the loops of alpha3/beta4 and alpha4/beta5 more stable. This means that phosphorylation entropically favors the receiver domain in its active state, and the induced conformational changes occur in an allosteric manner. Analyses of the local flexibility and long-range correlated motion also suggest a dynamics criterion for determining the allosteric cooperativity of NtrC, and may be applicable to other proteins.     

Ligand based conformational space studies of the μ-opioid receptor

BackgroundG protein-coupled receptors (GPCRs) comprise a family of membrane proteins that can be activated by a variety of external factors. The μ-opioid receptor (MOR), a class A GPCR, is the main target of morphine. Recently, enhanced sampling molecular dynamics simulations of a constitutively active mutant of MOR in its apo form allowed us to capture the novel intermediate states of activation, as well as the active state. This prompted us to apply the same techniques to wild type MOR in complex with ligands, in order to explore their contributions to the receptor conformational changes in the activation process.MethodsMOR was modeled in complex with agonists (morphine, BU72), a partial agonist (naloxone benzoylhydrazone) and an antagonist (naloxone). Replica exchange with solute tempering (REST2) molecular dynamics simulations were carried out for all systems. Trajectory frames were clustered, and the activation state of each cluster was assessed by two different methods.ResultsCluster sizes and activation indices show that while agonists stabilized structures in a higher activation state, the antagonist behaved oppositely. Morphine tends to drive the receptor towards increasing R165-T279 distances, while naloxone tends to increase the NPxxYA motif conformational change.ConclusionsDespite not observing a full transition between inactive and active states, an important conformational change of transmembrane helix 5 was observed and associated with a ligand-driven step of the process.General significanceThe activation process of GPCRs is widely studied but still not fully understood. Here we carried out a step forward in the direction of gaining more details of this process.     

Molecular insight into conformational transformation of human glucokinase: conventional and targeted molecular dynamics

Abstract

Glucokinase (GK) plays a key role in the regulation of hepatic glucose metabolism. An unusual mechanism of positive cooperativity of monomeric GK containing only a single binding site for glucose is very interesting and still unclear. The activation process of GK is associated with a large-scale conformational change from the inactive to the active state. Here, conventional and targeted molecular dynamics simulations were used to study the conformational dynamics of GK in the stable configurations and in the transition from active to inactive state. Three phases of the structural reorganization of GK were detected. The first step is a transformation of GK from the active state to the intermediate structure, where the cleft between the domains is open, but alpha helix 13 is still inside the small domain. From this point, there are two alternative paths. One path leads to the inactive state through the release of helix 13 from the inside of small domain to the outside. Other path goes back to the active state. Simulation results reveal the critical role of helix 13 in the transformation of GK from the open state to inactive one and the influence of the loop 2 on the protein transformation between the open and the closed active states. Principal component analysis and covariance matrix analysis were carried out to analyze the dynamics of protein. Importance of hydrogen bonds in the stability of the closed conformation is shown. Overall, our simulations provide new information about the dynamics of GK and its structural transformation.

Communicated by Ramaswamy H. Sarma     

Mechanism of activation of a G protein-coupled receptor, the human cholecystokinin-2 receptor

G protein-coupled receptors (GPCRs) represent a major focus in functional genomics programs and drug development research, but their important potential as drug targets contrasts with the still limited data available concerning their activation mechanism. Here, we investigated the activation mechanism of the cholecystokinin-2 receptor (CCK2R). The three-dimensional structure of inactive CCK2R was homology-modeled on the basis of crystal coordinates of inactive rhodopsin. Starting from the inactive CCK2R modeled structure, active CCK2R (namely cholecystokinin-occupied CCK2R) was modeled by means of steered molecular dynamics in a lipid bilayer and by using available data from other GPCRs, including rhodopsin. By comparing the modeled structures of the inactive and active CCK2R, we identified changes in the relative position of helices and networks of interacting residues, which were expected to stabilize either the active or inactive states of CCK2R. Using targeted molecular dynamics simulations capable of converting CCK2R from the inactive to the active state, we delineated structural changes at the atomic level. The activation mechanism involved significant movements of helices VI and V, a slight movement of helices IV and VII, and changes in the position of critical residues within or near the binding site. The mutation of key amino acids yielded inactive or constitutively active CCK2R mutants, supporting this proposed mechanism. Such progress in the refinement of the CCK2R binding site structure and in knowledge of CCK2R activation mechanisms will enable target-based optimization of nonpeptide ligands.     

Vav regulates activation of Rac but not Cdc42 during FcgammaR-mediated phagocytosis

Phagocytosis is the process whereby cells direct the spatially localized, receptor-driven engulfment of particulate materials. It proceeds via remodeling of the actin cytoskeleton and shares many of the core cytoskeletal components involved in adhesion and migration. Small GTPases of the Rho family have been widely implicated in coordinating actin dynamics in response to extracellular signals and during diverse cellular processes, including phagocytosis, yet the mechanisms controlling their recruitment and activation are not known. We show herein that in response to ligation of Fc receptors for IgG (FcgammaR), the guanine nucleotide exchange factor Vav translocates to nascent phagosomes and catalyzes GTP loading on Rac, but not Cdc42. The Vav-induced Rac activation proceeds independently of Cdc42 function, suggesting distinct roles for each GTPase during engulfment. Moreover, inhibition of Vav exchange activity or of Cdc42 activity does not prevent Rac recruitment to sites of particle attachment. We conclude that Rac is recruited to Fcgamma membrane receptors in its inactive, GDP-bound state and that Vav regulates phagocytosis through subsequent catalysis of GDP/GTP exchange on Rac.     

Consequences of lysine 72 mutation on the phosphorylation and activation state of cAMP-dependent kinase

General strategies to obtain inactive kinases have utilized mutation of key conserved residues in the kinase core, and the equivalent Lys72 in cAMP-dependent kinase has often been used to generate a "dead" kinase. Here, we have analyzed the consequences of this mutation on kinase structure and function. Mutation of Lys72 to histidine (K72H) generated an inactive enzyme, which was unphosphorylated. Treatment with an exogenous kinase (PDK-1) resulted in a mutant that was phosphorylated only at Thr197 and remained inactive but nevertheless capable of binding ATP. Ser338 in K72H cannot be autophosphorylated, nor can it be phosphorylated in an intermolecular process by active wild type C-subunit. The Lys72 mutant, once phosphorylated on Thr197, can bind with high affinity to the RIalpha subunits. Thus a dead kinase can still act as a scaffold for binding substrates and inhibitors; it is only phosphoryl transfer that is defective. Using a potent inhibitor of C-subunit activity, H-89, Escherichia coli-expressed C-subunit was also obtained in its unphosphorylated state. This protein is able to mature into its active form in the presence of PDK-1 and is able to undergo secondary autophosphorylation on Ser338. Unlike the H-89-treated wild type protein, the mutant protein (K72H) cannot undergo the subsequent cis autophosphorylation following phosphorylation at Thr197. Using these two substrates and mammalian-expressed PDK-1, we can elucidate a possible two-step process for the activation of the C-subunit: initial phosphorylation on the activation loop at Thr197 by PDK-1, or a PDK-1-like enzyme, followed by second cis autophosphorylation step at Ser338.     

An Autoinhibited Noncanonical Mechanism of GTP Hydrolysis by Rheb Maintains mTORC1 Homeostasis

Rheb, an activator of mammalian target of rapamycin (mTOR), displays low intrinsic GTPase activity favoring the biologically activated, GTP-bound state. We identified a Rheb mutation (Y35A) that increases its intrinsic nucleotide hydrolysis activity ～10-fold, and solved structures of both its active and inactive forms, revealing an unexpected mechanism of GTP hydrolysis involving Asp65 in switch II and Thr38 in switch I. In the wild-type protein this noncanonical mechanism is markedly inhibited by Tyr35, which constrains the active site conformation, restricting the access of the catalytic Asp65 to the nucleotide-binding pocket. Rheb Y35A mimics the enthalpic and entropic changes associated with GTP hydrolysis elicited by the GTPase-activating protein (GAP) TSC2, and is insensitive to further TSC2 stimulation. Overexpression of Rheb Y35A impaired the regulation of mTORC1 signaling by growth factor availability. We demonstrate that the opposing functions of Tyr35 in the intrinsic and GAP-stimulated GTP catalysis are critical for optimal mTORC1 regulation.     

Pepsinogen-like activation intermediate of plasmepsin II revealed by molecular dynamics analysis

Plasmepsins are pharmaceutically relevant aspartic proteases involved in haemoglobin degradation by the malaria causing parasites Plasmodium spp. They are translated as inactive proenzymes, with an elongated prosegment. On prosegment cleavage, plasmepsins undergo a series of hitherto unresolved conformational changes before becoming active. Here, the flexibility of plasmepsin and proplasmepsin and the activation process are investigated by multiple explicit water molecular dynamics simulations. The large N-terminal displacement and the interdomain shift from the proenzyme structure to active plasmepsin are promoted by essential dynamics sampling. An intermediate, stabilized by electrostatic interactions between the catalytic dyad and the N-terminus of mature plasmepsin, is observed along all activation trajectories. Notably, the stabilizing interactions in the activation intermediate of plasmepsin are similar to those in the X-ray structure of pepsinogen. In particular, the catalytic aspartates act as hydrogen bond acceptors for the N-terminal amino group and the Ser2 hydroxyl in plasmepsin, and the side chains of Lys36pro and Tyr9 in pepsinogen. The simulation results are used to suggest in vitro experiments to test the conformational transitions involved in the maturation of plasmepsin, and design small-molecule inhibitors.     

The activity of the extracellular signal-regulated kinase 2 is regulated by differential phosphorylation in the activation loop

The mitogen-activated protein kinases (MAP kinases) play a central role in signaling pathways initiated by extracellular stimuli such as growth factors, cytokines, and various forms of environmental stress. Full activation of the MAP kinases requires dual phosphorylation of the Thr and Tyr residues in the TXY motif of the activation loop by MAP kinase kinases. Interestingly, down-regulation of MAP kinase activity can be initiated by multiple Ser/Thr phosphatases, Tyr-specific phosphatases, and dual-specificity phosphatases. This would inevitable lead to the formation of monophosphorylated MAP kinases. However, in much of the literature investigating MAP kinase signaling, there has been the implicit assumption that the monophosphorylated forms are inactive. Thus, the significance for the need of multiple phosphatases in regulating MAP kinase activity is not clear, and the biological functions of these monophosphorylated MAP kinases are currently unknown. We have prepared extracellular signal-regulated protein kinase 2 (ERK2) in all phosphorylated forms and kinetically characterized them using two proteins (the myelin basic protein and Elk-1) and ATP as substrates. Our results revealed that a single phosphorylation in the activation loop of ERK2 produces an intermediate activity state. Thus, the catalytic efficiencies of the monophosphorylated ERK2/pY and ERK2/pT (ERK2 phosphorylated on Tyr-185 and Thr-183, respectively) are approximately 2-3 orders of magnitude higher than that of the unphosphorylated ERK2 and are only 1-2 orders of magnitude lower than that of the fully active bisphosphorylated ERK2/pTpY. This raises the possibility that the monophosphorylated ERK2s may have distinct biological roles in vivo. Different phosphorylation states in the activation loop could be linked to graded effects on a single ERK2 function. Alternatively, they could be linked to distinct ERK2 functions. Although less active than the bisphosphorylated species, the monophosphorylated ERK2s may differentially phosphorylate pathway components.     

Constitutively active mutants of the alpha 1B-adrenergic receptor: role of highly conserved polar amino acids in receptor activation.

Site-directed mutagenesis and molecular dynamics simulations of the alpha 1B-adrenergic receptor (AR) were combined to explore the potential molecular changes correlated with the transition from R (inactive state) to R (active state). Using molecular dynamics analysis we compared the structural/dynamic features of constitutively active mutants with those of the wild type and of an inactive alpha 1B-AR to build a theoretical model which defines the essential features of R and R. The results of site-directed mutagenesis were in striking agreement with the predictions of the model supporting the following hypothesis. (i) The equilibrium between R and R depends on the equilibrium between the deprotonated and protonated forms, respectively, of D142 of the DRY motif. In fact, replacement of D142 with alanine confers high constitutive activity to the alpha 1B-AR. (ii) The shift of R143 of the DRY sequence out of a conserved 'polar pocket' formed by N63, D91, N344 and Y348 is a feature common to all the active structures, suggesting that the role of R143 is fundamental for mediating receptor activation. Disruption of these intramolecular interactions by replacing N63 with alanine constitutively activates the alpha 1B-AR. Our findings might provide interesting generalities about the activation process of G protein-coupled receptors.     

Crystal structures of bacterial peptidoglycan amidase AmpD and an unprecedented activation mechanism

AmpD is a cytoplasmic peptidoglycan (PG) amidase involved in bacterial cell-wall recycling and in induction of β-lactamase, a key enzyme of β-lactam antibiotic resistance. AmpD belongs to the amidase_2 family that includes zinc-dependent amidases and the peptidoglycan-recognition proteins (PGRPs), highly conserved pattern-recognition molecules of the immune system. Crystal structures of Citrobacter freundii AmpD were solved in this study for the apoenzyme, for the holoenzyme at two different pH values, and for the complex with the reaction products, providing insights into the PG recognition and the catalytic process. These structures are significantly different compared with the previously reported NMR structure for the same protein. The NMR structure does not possess an accessible active site and shows the protein in what is proposed herein as an inactive “closed” conformation. The transition of the protein from this inactive conformation to the active “open” conformation, as seen in the x-ray structures, was studied by targeted molecular dynamics simulations, which revealed large conformational rearrangements (as much as 17 Å) in four specific regions representing one-third of the entire protein. It is proposed that the large conformational change that would take the inactive NMR structure to the active x-ray structure represents an unprecedented mechanism for activation of AmpD. Analysis is presented to argue that this activation mechanism might be representative of a regulatory process for other intracellular members of the bacterial amidase_2 family of enzymes.     

Identification of novel phosphorylation sites required for activation of MAPKAP kinase-2.

MAP kinase-activated protein (MAPKAP) kinase-2 is activated in vivo by reactivating kinase (RK), a MAP kinase (MAPK) homologue stimulated by cytokines and cellular stresses. Here we show that in vitro RK phosphorylates human GST-MAPKAP kinase-2 at Thr25 in the proline-rich N-terminal region Thr222 and Ser272 in the catalytic domain and Thr334 in the C-terminal domain. Using novel methodology we demonstrate that activation of MAPKAP kinase-2 requires the phosphorylation of any two of the three residues Thr222, Ser272 and Thr334. Ser9, Thr25, Thr222, Ser272, Thr334 and Thr338 became 32P-labelled in stressed KB cells and labelling was prevented by the specific RK inhibitor SB 203580, establishing that RK phosphorylates Thr25, Thr222, Ser272 and Thr334 in vivo. The 32P-labelling of Thr338 is likely to result from autophosphorylation. GST-MAPKAP kinase-2 lacking the N-terminal domain was inactive, but activated fully when phosphorylated at Thr222, Ser272 and Thr334 by p42 MAPK or RK. In contrast, full-length GST-MAPKAP kinase-2 was phosphorylated at Thr25 (and not activated) by p42 MAPK, suggesting a role for the N-terminal domain in specifying activation by RK in vivo. The mutant Asp222/Asp334 was 20% as active as phosphorylated MAPKAP kinase-2, and this constitutively active form may be useful for studying its physiological roles.     

Long‐range molecular dynamics show that inactive forms of Protein Kinase A are more dynamic than active forms

Many protein kinases are characterized by at least two structural forms corresponding to the highest level of activity (active) and low or no activity, (inactive). Further, protein dynamics is an important consideration in understanding the molecular and mechanistic basis of enzyme function. In this work, we use protein kinase A (PKA) as the model system and perform microsecond range molecular dynamics (MD) simulations on six variants which differ from one another in terms of active and inactive form, with or without bound ligands, C‐terminal tail and phosphorylation at the activation loop. We find that the root mean square fluctuations in the MD simulations are generally higher for the inactive forms than the active forms. This difference is statistically significant. The higher dynamics of inactive states has significant contributions from ATP binding loop, catalytic loop, and αG helix. Simulations with and without C‐terminal tail show this differential dynamics as well, with lower dynamics both in the active and inactive forms if C‐terminal tail is present. Similarly, the dynamics associated with the inactive form is higher irrespective of the phosphorylation status of Thr 197. A relatively stable stature of active kinases may be better suited for binding of substrates and detachment of the product. Also, phosphoryl group transfer from ATP to the phosphosite on the substrate requires precise transient coordination of chemical entities from three different molecules, which may be facilitated by the higher stability of the active state.     

Mechanism of the chemical step for the guanosine triphosphate (GTP) hydrolysis catalyzed by elongation factor Tu

Elongation factor Tu (EF-Tu), the protein responsible for delivering aminoacyl-tRNAs (aa-tRNAs) to ribosomal A site during translation, belongs to the group of guanosine-nucleotide (GTP/GDP) binding proteins. Its active 'on'-state corresponds to the GTP-bound form, while the inactive 'off'-state corresponds to the GDP-bound form. In this work we focus on the chemical step, GTP+H(2)O-->GDP+Pi, of the hydrolysis mechanism. We apply molecular modeling tools including molecular dynamics simulations and the combined quantum mechanical-molecular mechanical calculations for estimates of reaction energy profiles for two possible arrangements of switch II regions of EF-Tu. In the first case we presumably mimic binding of the ternary complex EF-Tu.GTP.aa-tRNA to the ribosome and allow the histidine (His85) side chain of the protein to approach the reaction active site. In the second case, corresponding to the GTP hydrolysis by EF-Tu alone, the side chain of His85 stays away from the active site, and the chemical reaction GTP+H(2)O-->GDP+Pi proceeds without participation of the histidine but through water molecules. In agreement with the experimental observations which distinguish rate constants for the fast chemical reaction in EF-Tu.GTP.aa-tRNA.ribosome and the slow spontaneous GTP hydrolysis in EF-Tu, we show that the activation energy barrier for the first scenario is considerably lower compared to that of the second case.     

Stabilization of an intermediate activation state for transducin by a fluorescent GTP analogue

The GTP-binding protein (G protein), transducin, serves as a key molecular switch in vertebrate vision through the tight regulation of its GTP-binding (activation)/GTP hydrolytic (deactivation) cycle by the photoreceptor rhodopsin. To better understand the structure-function characteristics of transducin activation, we have set out to identify spectroscopic probes that bind to the guanine nucleotide-binding site of this G protein and maintain its ability to interact with its specific cellular target/effector, the cyclic GMP phosphodiesterase (PDE). In this study, we describe the characterization of a fluorescently labeled GTP analogue, BODIPY-FL GTPgammaS (BOD-GTPgammaS), that binds to the alpha subunit of transducin (alpha(T)) in a rhodopsin- and Gbetagamma-dependent manner, similar to the binding of GTP or GTPgammaS, with an apparent dissociation constant of 100 nM. The rhodopsin-dependent binding of BOD-GTPgammaS to alpha(T) is slow, relative to the rate of binding of GTPgammaS, particularly under conditions where rhodopsin must act catalytically to stimulate the exchange of BOD-GTPgammaS for GDP on multiple alpha(T) subunits. This reflects a slower rate of dissociation of rhodopsin and Gbetagamma from alpha(T)-BOD-GTPgammaS complexes, relative to their rates of dissociation from alpha(T)-GTPgammaS. The binding of BOD-GTPgammaS occurs without a change in the intrinsic tryptophan fluorescence of alpha(T), indicating that only a subtle movement of the Switch 2 domain on alpha(T) accompanies the binding of this GTPgammaS analogue. Nevertheless, the BOD-GTPgammaS-bound alpha(T) subunit is able to bind with high affinity to the recombinant, purified gamma subunit of PDE (gamma(PDE)) labeled with 5-((((2-iodoacetyl)amino)ethyl)amino)naphthalene-1-sulfonic acid (IAEDANS (K(d) approximately 13 nM)), as well as bind to and stimulate the activity of PDE, albeit less efficiently compared to alpha(T)-GTPgammaS. Taken together, these findings suggest that the binding of BOD-GTPgammaS to transducin causes it to adopt a distinct conformation that appears to be intermediate between the inactive and fully active states of alpha(T), and this fluorescent nucleotide analogue can be used as a reporter group to characterize the interactions of alpha(T) in this conformational state with its biological target/effector.     

Retinal dynamics during light activation of rhodopsin revealed by solid-state NMR spectroscopy

Rhodopsin is a canonical member of class A of the G protein-coupled receptors (GPCRs) that are implicated in many of the drug interventions in humans and are of great pharmaceutical interest. The molecular mechanism of rhodopsin activation remains unknown as atomistic structural information for the active metarhodopsin II state is currently lacking. Solid-state ²H NMR constitutes a powerful approach to study atomic-level dynamics of membrane proteins. In the present application, we describe how information is obtained about interactions of the retinal cofactor with rhodopsin that change with light activation of the photoreceptor. The retinal methyl groups play an important role in rhodopsin function by directing conformational changes upon transition into the active state. Site-specific ²H labels have been introduced into the methyl groups of retinal and solid-state ²H NMR methods applied to obtain order parameters and correlation times that quantify the mobility of the cofactor in the inactive dark state, as well as the cryotrapped metarhodopsin I and metarhodopsin II states. Analysis of the angular-dependent ²H NMR line shapes for selectively deuterated methyl groups of rhodopsin in aligned membranes enables determination of the average ligand conformation within the binding pocket. The relaxation data suggest that the β-ionone ring is not expelled from its hydrophobic pocket in the transition from the pre-activated metarhodopsin I to the active metarhodopsin II state. Rather, the major structural changes of the retinal cofactor occur already at the metarhodopsin I state in the activation process. The metarhodopsin I to metarhodopsin II transition involves mainly conformational changes of the protein within the membrane lipid bilayer rather than the ligand. The dynamics of the retinylidene methyl groups upon isomerization are explained by an activation mechanism involving cooperative rearrangements of extracellular loop E2 together with transmembrane helices H5 and H6. These activating movements are triggered by steric clashes of the isomerized all-trans retinal with the β4 strand of the E2 loop and the side chains of Glu¹²² and Trp²⁶⁵ within the binding pocket. The solid-state ²H NMR data are discussed with regard to the pathway of the energy flow in the receptor activation mechanism.