Membrane model for the G-protein-coupled receptor rhodopsin: hydrophobic interface and dynamical structure

Exploring GPCR conformational dynamics using single-molecule fluorescence

G protein-coupled receptors (GPCRs) are membrane proteins that transmit specific external stimuli into cells by changing their conformation. This conformational change allows them to couple and activate G-proteins to initiate signal transduction. A critical challenge in studying and inferring these structural dynamics arises from the complexity of the cellular environment, including the presence of various endogenous factors. Due to the recent advances in cell-expression systems, membrane-protein purification techniques, and labeling approaches, it is now possible to study the structural dynamics of GPCRs at a single-molecule level both in vitro and in live cells. In this review, we discuss state-of-the-art techniques and strategies for expressing, purifying, and labeling GPCRs in the context of single-molecule research. We also highlight four recent studies that demonstrate the applications of single-molecule microscopy in revealing the dynamics of GPCRs. These techniques are also useful as complementary methods to verify the results obtained from other structural biology tools like cryo-electron microscopy and x-ray crystallography.

Interdisciplinary biophysical studies of membrane proteins bacteriorhodopsin and rhodopsin

The centenary of the birth of H. Gobind Khorana provides an auspicious opportunity to review the origins and evolution of parallel advances in biophysical methodology and molecular genetics technology used to study membrane proteins. Interdisciplinary work in the Khorana laboratory in the late 1970s and for the next three decades led to productive collaborations and fostered three subsequent scientific generations whose biophysical work on membrane proteins has led to detailed elucidation of the molecular mechanisms of energy transduction by the light-driven proton pump bacteriorhodopsin (bR) and signal transduction by the G protein-coupled receptor (GPCR) rhodopsin. This review will highlight the origins and advances of biophysical studies of membrane proteins made possible by the application of molecular genetics approaches to engineer site-specific alterations of membrane protein structures.

Regulation of membrane protein structure and function by their lipid nano-environment

Transmembrane proteins comprise ~30% of the mammalian proteome, mediating metabolism, signalling, transport and many other functions required for cellular life. The microenvironment of integral membrane proteins (IMPs) is intrinsically different from that of cytoplasmic proteins, with IMPs solvated by a compositionally and biophysically complex lipid matrix. These solvating lipids affect protein structure and function in a variety of ways, from stereospecific, high-affinity protein-lipid interactions to modulation by bulk membrane properties. Specific examples of functional modulation of IMPs by their solvating membranes have been reported for various transporters, channels and signal receptors; however, generalizable mechanistic principles governing IMP regulation by lipid environments are neither widely appreciated nor completely understood. Here, we review recent insights into the inter-relationships between complex lipidomes of mammalian membranes, the membrane physicochemical properties resulting from such lipid collectives, and the regulation of IMPs by either or both. The recent proliferation of high-resolution methods to study such lipid-protein interactions has led to generalizable insights, which we synthesize into a general framework termed the 'functional paralipidome' to understand the mutual regulation between membrane proteins and their surrounding lipid microenvironments.

Phospholipid headgroups govern area per lipid and emergent elastic properties of bilayers

Phospholipid bilayers are liquid-crystalline materials whose intermolecular interactions at mesoscopic length scales have key roles in the emergence of membrane physical properties. Here we investigated the combined effects of phospholipid polar headgroups and acyl chains on biophysical functions of membranes with solid-state 2H NMR spectroscopy. We compared the structural and dynamic properties of phosphatidylethanolamine and phosphatidylcholine with perdeuterated acyl chains in the solid-ordered (so) and liquid-disordered (ld) phases. Our analysis of spectral lineshapes of 1,2-diperdeuteriopalmitoyl-sn-glycero-3-phosphoethanolamine (DPPE-d62) and 1,2-diperdeuteriopalmitoyl-sn-glycero-3-phosphocholine (DPPC-d62) in the so (gel) phase indicated an all-trans rotating chain structure for both lipids. Greater segmental order parameters (SCD) were observed in the ld (liquid-crystalline) phase for DPPE-d62 than for DPPC-d62 membranes, while their mixtures had intermediate values irrespective of the deuterated lipid type. Our results suggest the SCD profiles of the acyl chains are governed by methylation of the headgroups and are averaged over the entire system. Variations in the acyl chain molecular dynamics were further investigated by spin-lattice (R1Z) and quadrupolar-order relaxation (R1Q) measurements. The two acyl-perdeuterated lipids showed distinct differences in relaxation behavior as a function of the order parameter. The R1Z rates had a square-law dependence on SCD, implying collective mesoscopic dynamics, with a higher bending rigidity for DPPE-d62 than for DPPC-d62 lipids. Remodeling of lipid average and dynamic properties by methylation of the headgroups thus provides a mechanism to control the actions of peptides and proteins in biomembranes.

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.

Insights into Membrane Curvature Sensing and Membrane Remodeling by Intrinsically Disordered Proteins and Protein Regions

Cellular membranes are highly dynamic in shape. They can rapidly and precisely regulate their shape to perform various cellular functions. The protein’s ability to sense membrane curvature is essential in various biological events such as cell signaling and membrane trafficking. As they are bound, these curvature-sensing proteins may also change the local membrane shape by one or more curvature driving mechanisms. Established curvature-sensing/driving mechanisms rely on proteins with specific structural features such as amphipathic helices and intrinsically curved shapes. However, the recent discovery and characterization of many proteins have shattered the protein structure–function paradigm, believing that the protein functions require a unique structural feature. Typically, such structure-independent functions are carried either entirely by intrinsically disordered proteins or hybrid proteins containing disordered regions and structured domains. It is becoming more apparent that disordered proteins and regions can be potent sensors/inducers of membrane curvatures. In this article, we outline the basic features of disordered proteins and regions, the motifs in such proteins that encode the function, membrane remodeling by disordered proteins and regions, and assays that may be employed to investigate curvature sensing and generation by ordered/disordered proteins.

Exploration of the dynamic interplay between lipids and membrane proteins by hydrostatic pressure

Cell membranes represent a complex and variable medium in time and space of lipids and proteins. Their physico-chemical properties are determined by lipid components which can in turn influence the biological function of membranes. Here, we used hydrostatic pressure to study the close dynamic relationships between lipids and membrane proteins. Experiments on the β–barrel OmpX and the α–helical BLT2 G Protein-Coupled Receptor in nanodiscs of different lipid compositions reveal conformational landscapes intimately linked to pressure and lipids. Pressure can modify the conformational landscape of the membrane protein per se, but also increases the gelation of lipids, both being monitored simultaneously at high atomic resolution by NMR. Our study also clearly shows that a membrane protein can modulate, at least locally, the fluidity of the bilayer. The strategy proposed herein opens new perspectives to scrutinize the dynamic interplay between membrane proteins and their surrounding lipids.

Direct information on the dynamic interplay between membrane proteins and lipids is scarce. Here the authors report a detailed description of these close relationships by combining lipid nanodiscs and high-pressure NMR. They report the link between pressure and lipid compositions to the conformational landscape of the β-barrel OmpX and the α-helical BLT2 G Protein-Coupled Receptor in nanodiscs.

Dimeric Rhodopsin R135L Mutant-Transducin-like Complex Sheds Light on Retinitis Pigmentosa Misfunctions

Rhodopsin (RHO) is a light-sensitive pigment in the retina and the main prototypical protein of the G-protein-coupled receptor (GCPR) family. After receiving a light stimulus, RHO and its cofactor retinylidene undergo a series of structural changes that initiate an intricate transduction mechanism. Along with RHO, other partner proteins play key roles in the signaling pathway. These include transducin, a GTPase, kinases that phosphorylate RHO, and arrestin (Arr), which ultimately stops the signaling process and promotes RHO regeneration. A large number of RHO genetic mutations may lead to very severe retinal dysfunction and eventually to impaired dark adaptation disease called autosomal dominant retinitis pigmentosa (adRP). In this study, we used molecular dynamics (MD) simulations to evaluate the different behaviors of the dimeric form of wild-type RHO (WT dRHO) and its mutant at position 135 of arginine to leucine (dR135L), both in the free (noncomplexed) and in complex with the transducin-like protein (Gtl). Gtl is a heterotrimeric model composed of a mixture of human and bovine G proteins. Our calculations allow us to explain how the mutation causes structural changes in the RHO dimer and how this can affect the signal that transducin generates when it is bound to RHO. Moreover, the structural modifications induced by the R135L mutation can also account for other misfunctions observed in the up- and downstream signaling pathways. The mechanism of these dysfunctions, together with the transducin activity reduction, provides structure-based explanations of the impairment of some key processes that lead to adRP.

Recent developments in membrane curvature sensing and induction by proteins

BACKGROUND

Membrane-bound intracellular organelles have characteristic shapes attributed to different local membrane curvatures, and these attributes are conserved across species. Over the past decade, it has been confirmed that specific proteins control the large curvatures of the membrane, whereas many others due to their specific structural features can sense the curvatures and bind to the specific geometrical cues. Elucidating the interplay between sensing and induction is indispensable to understand the mechanisms behind various biological processes such as vesicular trafficking and budding.

SCOPE OF REVIEW

We provide an overview of major classes of membrane proteins and the mechanisms of curvature sensing and induction. We then discuss the importance of membrane elastic characteristics to induce the membrane shapes similar to intracellular organelles. Finally, we survey recently available assays developed for studying the curvature sensing and induction by many proteins.

MAJOR CONCLUSIONS

Recent theoretical/computational modeling along with experimental studies have uncovered fascinating connections between lipid membrane and protein interactions. However, the phenomena of protein localization and synchronization to generate spatiotemporal dynamics in membrane morphology are yet to be fully understood.

GENERAL SIGNIFICANCE

The understanding of protein-membrane interactions is essential to shed light on various biological processes. This further enables the technological applications of many natural proteins/peptides in therapeutic treatments. The studies of membrane dynamic shapes help to understand the fundamental functions of membranes, while the medicinal roles of various macromolecules (such as proteins, peptides, etc.) are being increasingly investigated.

Activation of G-protein-coupled receptors is thermodynamically linked to lipid solvation

Preferential lipid solvation of the G-protein-coupled A_2A adenosine receptor (A_2AR) is evaluated from 35 μs of all-atom molecular dynamics simulation. A coarse-grained transition matrix algorithm is developed to overcome slow equilibration of the first solvation shell, obtaining estimates of the free energy of solvation by different lipids for the receptor in different activation states. Results indicate preference for solvation by unsaturated chains, which favors the active receptor. A model for lipid-dependent G-protein-coupled receptor activity is proposed in which the chemical potential of lipids in the bulk membrane modulates receptor activity. The entropies associated with moving saturated and unsaturated lipids from bulk to A_2AR's first solvation shell are evaluated. Overall, the acyl chains are more disordered (i.e., obtain a favorable entropic contribution) when partitioning to the receptor surface, and this effect is augmented for the saturated chains, which are relatively more ordered in bulk.

Activation of the G-Protein-Coupled Receptor Rhodopsin by Water

Visual rhodopsin is an important archetype for G-protein-coupled receptors, which are membrane proteins implicated in cellular signal transduction. Herein, we show experimentally that approximately 80 water molecules flood rhodopsin upon light absorption to form a solvent-swollen active state. An influx of mobile water is necessary for activating the photoreceptor, and this finding is supported by molecular dynamics (MD) simulations. Combined force-based measurements involving osmotic and hydrostatic pressure indicate the expansion occurs by changes in cavity volumes, together with greater hydration in the active metarhodopsin-II state. Moreover, we discovered that binding and release of the C-terminal helix of transducin is coupled to hydration changes as may occur in visual signal amplification. Hydration-dehydration explains signaling by a dynamic allosteric mechanism, in which the soft membrane matter (lipids and water) has a pivotal role in the catalytic G-protein cycle.

Structure, dynamics and lipid interactions of serotonin receptors: excitements and challenges

Serotonin (5-hydroxytryptamine, 5-HT) is an intrinsically fluorescent neurotransmitter found in organisms spanning a wide evolutionary range. Serotonin exerts its diverse actions by binding to distinct cell membrane receptors which are classified into many groups. Serotonin receptors are involved in regulating a diverse array of physiological signaling pathways and belong to the family of either G protein-coupled receptors (GPCRs) or ligand-gated ion channels. Serotonergic signaling appears to play a key role in the generation and modulation of various cognitive and behavioral functions such as sleep, mood, pain, anxiety, depression, aggression, and learning. Serotonin receptors act as drug targets for a number of diseases, particularly neuropsychiatric disorders. The signaling mechanism and efficiency of serotonin receptors depend on their amazing ability to rapidly access multiple conformational states. This conformational plasticity, necessary for the wide variety of functions displayed by serotonin receptors, is regulated by binding to various ligands. In this review, we provide a succinct overview of recent developments in generating and analyzing high-resolution structures of serotonin receptors obtained using crystallography and cryo-electron microscopy. Capturing structures of distinct conformational states is crucial for understanding the mechanism of action of these receptors, which could provide important insight for rational drug design targeting serotonin receptors. We further provide emerging information and insight from studies on interactions of membrane lipids (such as cholesterol) with serotonin receptors. We envision that a judicious combination of analysis of high-resolution structures and receptor-lipid interaction would allow a comprehensive understanding of GPCR structure, function and dynamics, thereby leading to efficient drug discovery.

Quantum Mechanical and Molecular Mechanics Modeling of Membrane-Embedded Rhodopsins

Computational chemistry provides versatile methods for studying the properties and functioning of biological systems at different levels of precision and at different time scales. The aim of this article is to review the computational methodologies that are applicable to rhodopsins as archetypes for photoactive membrane proteins that are of great importance both in nature and in modern technologies. For each class of computational techniques, from methods that use quantum mechanics for simulating rhodopsin photophysics to less-accurate coarse-grained methodologies used for long-scale protein dynamics, we consider possible applications and the main directions for improvement.

Flexible lipid nanomaterials studied by NMR spectroscopy

Advances in solid-state nuclear magnetic resonance spectroscopy inform the emergence of material properties from atomistic-level interactions in membrane lipid nanostructures.

Cholesterol Effects on the Physical Properties of Lipid Membranes Viewed by Solid-state NMR Spectroscopy

In this chapter, we review the physical properties of lipid/cholesterol mixtures involving studies of model membranes using solid-state NMR spectroscopy. The approach allows one to quantify the average membrane structure, fluctuations, and elastic deformation upon cholesterol interaction. Emphasis is placed on understanding the membrane structural deformation and emergent fluctuations at an atomistic level. Lineshape measurements using solid-state NMR spectroscopy give equilibrium structural properties, while relaxation time measurements study the molecular dynamics over a wide timescale range. The equilibrium properties of glycerophospholipids, sphingolipids, and their binary and tertiary mixtures with cholesterol are accessible. Nonideal mixing of cholesterol with other lipids explains the occurrence of liquid-ordered domains. The entropic loss upon addition of cholesterol to sphingolipids is less than for glycerophospholipids, and may drive formation of lipid rafts. The functional dependence of ²H NMR spin-lattice relaxation (R _1Z) rates on segmental order parameters (S _CD) for lipid membranes is indicative of emergent viscoelastic properties. Addition of cholesterol shows stiffening of the bilayer relative to the pure lipids and this effect is diminished for lanosterol. Opposite influences of cholesterol and detergents on collective dynamics and elasticity at an atomistic scale can potentially affect lipid raft formation in cellular membranes.

A Critical Analysis of Molecular Mechanisms Underlying Membrane Cholesterol Sensitivity of GPCRs

G protein-coupled receptors (GPCRs) are the largest and a diverse family of proteins involved in signal transduction across biological membranes. GPCRs mediate a wide range of physiological processes and have emerged as major targets for the development of novel drug candidates in all clinical areas. Since GPCRs are integral membrane proteins, regulation of their organization, dynamics, and function by membrane lipids, in particular membrane cholesterol, has emerged as an exciting area of research. Cholesterol sensitivity of GPCRs could be due to direct interaction of cholesterol with the receptor (specific effect). Alternately, GPCR function could be influenced by the effect of cholesterol on membrane physical properties (general effect). In this review, we critically analyze the specific and general mechanisms of the modulation of GPCR function by membrane cholesterol, taking examples from representative GPCRs. While evidence for both the proposed mechanisms exists, there appears to be no clear-cut distinction between these two mechanisms, and a combination of these mechanisms cannot be ruled out in many cases. We conclude that classifying the mechanism underlying cholesterol sensitivity of GPCR function merely into these two mutually exclusive classes could be somewhat arbitrary. A more holistic approach could be suitable for analyzing GPCR-cholesterol interaction.

GPCRs: What Can We Learn from Molecular Dynamics Simulations?

Advances in the structural biology of G-protein Coupled Receptors have resulted in a significant step forward in our understanding of how this important class of drug targets function at the molecular level. However, it has also become apparent that they are very dynamic molecules, and moreover, that the underlying dynamics is crucial in shaping the response to different ligands. Molecular dynamics simulations can provide unique insight into the dynamic properties of GPCRs in a way that is complementary to many experimental approaches. In this chapter, we describe progress in three distinct areas that are particularly difficult to study with other techniques: atomic level investigation of the conformational changes that occur when moving between the various states that GPCRs can exist in, the pathways that ligands adopt during binding/unbinding events and finally, the influence of lipids on the conformational dynamics of GPCRs.

Concepts and Methods of Solid-State NMR Spectroscopy Applied to Biomembranes

Concepts of solid-state NMR spectroscopy and applications to fluid membranes are reviewed in this paper. Membrane lipids with ²H-labeled acyl chains or polar head groups are studied using ²H NMR to yield knowledge of their atomistic structures in relation to equilibrium properties. This review demonstrates the principles and applications of solid-state NMR by unifying dipolar and quadrupolar interactions and highlights the unique features offered by solid-state ²H NMR with experimental illustrations. For randomly oriented multilamellar lipids or aligned membranes, solid-state ²H NMR enables direct measurement of residual quadrupolar couplings (RQCs) due to individual C-²H-labeled segments. The distribution of RQC values gives nearly complete profiles of the segmental order parameters S_CD⁽ⁱ⁾ as a function of acyl segment position (i). Alternatively, one can measure residual dipolar couplings (RDCs) for natural abundance lipid samples to obtain segmental S_CH order parameters. A theoretical mean-torque model provides acyl-packing profiles representing the cumulative chain extension along the normal to the aqueous interface. Equilibrium structural properties of fluid bilayers and various thermodynamic quantities can then be calculated, which describe the interactions with cholesterol, detergents, peptides, and integral membrane proteins and formation of lipid rafts. One can also obtain direct information for membrane-bound peptides or proteins by measuring RDCs using magic-angle spinning (MAS) in combination with dipolar recoupling methods. Solid-state NMR methods have been extensively applied to characterize model membranes and membrane-bound peptides and proteins, giving unique information on their conformations, orientations, and interactions in the natural liquid-crystalline state.

Interplay of G Protein-Coupled Receptors with the Membrane: Insights from Supra-Atomic Coarse Grain Molecular Dynamics Simulations

G protein-coupled receptors (GPCRs) are central to many fundamental cellular signaling pathways. They transduce signals from the outside to the inside of cells in physiological processes ranging from vision to immune response. It is extremely challenging to look at them individually using conventional experimental techniques. Recently, a pseudo atomistic molecular model has emerged as a valuable tool to access information on GPCRs, more specifically on their interactions with their environment in their native cell membrane and the consequences on their supramolecular organization. This approach uses the Martini coarse grain (CG) model to describe the receptors, lipids, and solvent in molecular dynamics (MD) simulations and in enough detail to allow conserving the chemical specificity of the different molecules. The elimination of unnecessary degrees of freedom has opened up large-scale simulations of the lipid-mediated supramolecular organization of GPCRs. Here, after introducing the Martini CGMD method, we review these studies carried out on various members of the GPCR family, including rhodopsin (visual receptor), opioid receptors, adrenergic receptors, adenosine receptors, dopamine receptor, and sphingosine 1-phosphate receptor. These studies have brought to light an interesting set of novel biophysical principles. The insights range from revealing localized and heterogeneous deformations of the membrane bilayer at the surface of the protein, specific interactions of lipid molecules with individual GPCRs, to the effect of the membrane matrix on global GPCR self-assembly. The review ends with an overview of the lessons learned from the use of the CGMD method, the biophysical-chemical findings on lipid-protein interplay.

Cholesterol-induced suppression of membrane elastic fluctuations at the atomistic level

Applications of solid-state NMR spectroscopy for investigating the influences of lipid-cholesterol interactions on membrane fluctuations are reviewed in this paper. Emphasis is placed on understanding the energy landscapes and fluctuations at an emergent atomistic level. Solid-state (2)H NMR spectroscopy directly measures residual quadrupolar couplings (RQCs) due to individual C-(2)H labeled segments of the lipid molecules. Moreover, residual dipolar couplings (RDCs) of (13)C-(1)H bonds are obtained in separated local-field NMR spectroscopy. The distributions of RQC or RDC values give nearly complete profiles of the order parameters as a function of acyl segment position. Measured equilibrium properties of glycerophospholipids and sphingolipids including their binary and tertiary mixtures with cholesterol show unequal mixing associated with liquid-ordered domains. The entropic loss upon addition of cholesterol to sphingolipids is less than for glycerophospholipids and may drive the formation of lipid rafts. In addition relaxation time measurements enable one to study the molecular dynamics over a wide time-scale range. For (2)H NMR the experimental spin-lattice (R1Z) relaxation rates follow a theoretical square-law dependence on segmental order parameters (SCD) due to collective slow dynamics over mesoscopic length scales. The functional dependence for the liquid-crystalline lipid membranes is indicative of viscoelastic properties as they emerge from atomistic-level interactions. A striking decrease in square-law slope upon addition of cholesterol denotes stiffening relative to the pure lipid bilayers that is diminished in the case of lanosterol. Measured equilibrium properties and relaxation rates infer opposite influences of cholesterol and detergents on collective dynamics and elasticity at an atomistic scale that potentially affects lipid raft formation in cellular membranes.

Micelle-Enhanced Bioorthogonal Labeling of Genetically Encoded Azido Groups on the Lipid-Embedded Surface of a GPCR

Genetically encoded p-azido-phenylalanine (azF) residues in G protein-coupled receptors (GPCRs) can be targeted with dibenzocyclooctyne-modified (DIBO-modified) fluorescent probes by means of strain-promoted [3+2] azide-alkyne cycloaddition (SpAAC). Here we show that azF residues situated on the transmembrane surfaces of detergent-solubilized receptors exhibit up to 1000-fold rate enhancement relative to azF residues on water-exposed surfaces. We show that the amphipathic moment of the labeling reagent, consisting of hydrophobic DIBO coupled to hydrophilic Alexa dye, results in strong partitioning of the DIBO group into the hydrocarbon core of the detergent micelle and consequently high local reactant concentrations. The observed rate constant for the micelleenhanced SpAAC is comparable with those of the fastest bioorthogonal labeling reactions known. Targeting hydrophobic regions of membrane proteins by use of micelle-enhanced SpAAC should expand the utility of bioorthogonal labeling strategies.

All-atom molecular dynamics simulation of a photosystem i/detergent complex

All-atom molecular dynamics (MD) simulation was used to investigate the solution structure and dynamics of the photosynthetic pigment-protein complex photosystem I (PSI) from Thermosynechococcus elongatus embedded in a toroidal belt of n-dodecyl-β-d-maltoside (DDM) detergent. Evaluation of root-mean-square deviations (RMSDs) relative to the known crystal structure show that the protein complex surrounded by DDM molecules is stable during the 200 ns simulation time, and root-mean-square fluctuation (RMSF) analysis indicates that regions of high local mobility correspond to solvent-exposed regions such as turns in the transmembrane α-helices and flexible loops on the stromal and lumenal faces. Comparing the protein-detergent complex to a pure detergent micelle, the detergent surrounding the PSI trimer is found to be less densely packed but with more ordered detergent tails, contrary to what is seen in most lipid bilayer models. We also investigated any functional implications for the observed conformational dynamics and protein-detergent interactions, discovering interesting structural changes in the psaL subunits associated with maintaining the trimeric structure of the protein. Importantly, we find that the docking of soluble electron mediators such as cytochrome c6 and ferredoxin to PSI is not significantly impacted by the solubilization of PSI in detergent.

Structural dynamics of the cell wall precursor lipid II in the presence and absence of the lantibiotic nisin

Representing a physiological "Achilles' heel", the cell wall precursor lipid II (LII) is a prime target for various classes of antibiotics. Over the years LII-binding agents have been recognized as promising candidates and templates in the search for new antibacterial compounds to complement or replace existing drugs. To elucidate the molecular structural basis underlying LII functional mechanism and to better understand if and how lantibiotic binding alters the molecular behavior of LII, we performed molecular dynamics (MD) simulations of phospholipid membrane-embedded LII in the absence and presence of the LII-binding lantibiotic nisin. In a series of 2×4 independent, unbiased 100ns MD simulations we sampled the conformational dynamics of nine LII as well as nine LII-nisin complexes embedded in an aqueous 150mM NaCl/POPC phospholipid membrane environment. We found that nisin binding to LII induces a reduction of LII mobility and flexibility, an outward shift of the LII pentapeptide, an inward movement of the LII disaccharide section, and an overall deeper insertion of the LII tail group into the membrane. The latter effect might indicate an initial step in adopting a stabilizing, scaffold-like structure in the process of nisin-induced membrane leakage. At the same time nisin conformation and LII interaction remain similar to the 1WCO LII-nisin NMR solution structure.

Molecular dynamics simulations of the adenosine A2a receptor in POPC and POPE lipid bilayers: effects of membrane on protein behavior

Analysis of 300 ns (ns) molecular dynamics (MD) simulations of an adenosine A2a receptor (A2a AR) model, conducted in triplicate, in 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) and 1-palmitoyl-2-oleoylphosphatidylethanolamine (POPE) bilayers reveals significantly different protein dynamical behavior. Principal component analysis (PCA) shows that the dissimilarities stem from interhelical rather than intrahelical motions. The difference in the hydrophobic thicknesses of these simulated lipid bilayers is potentially a significant reason for the observed difference in results. The distinct lipid headgroups might also lead to different molecular interactions and hence different protein loop motions. Overall, the A2a AR shows higher mobility and flexibility in POPC as compared to POPE.

Retinal ligand mobility explains internal hydration and reconciles active rhodopsin structures

Rhodopsin, the mammalian dim-light receptor, is one of the best-characterized G-protein-coupled receptors, a pharmaceutically important class of membrane proteins that has garnered a great deal of attention because of the recent availability of structural information. Yet the mechanism of rhodopsin activation is not fully understood. Here, we use microsecond-scale all-atom molecular dynamics simulations, validated by solid-state (2)H nuclear magnetic resonance spectroscopy, to understand the transition between the dark and metarhodopsin I (Meta I) states. Our analysis of these simulations reveals striking differences in ligand flexibility between the two states. Retinal is much more dynamic in Meta I, adopting an elongated conformation similar to that seen in the recent activelike crystal structures. Surprisingly, this elongation corresponds to both a dramatic influx of bulk water into the hydrophobic core of the protein and a concerted transition in the highly conserved Trp265(6.48) residue. In addition, enhanced ligand flexibility upon light activation provides an explanation for the different retinal orientations observed in X-ray crystal structures of active rhodopsin.

Coarse-grained molecular dynamics provides insight into the interactions of lipids and cholesterol with rhodopsin

Protein function is a complicated interplay between structure and dynamics, which can be heavily influenced by environmental factors and conditions. This is particularly true in the case of membrane proteins, such as the visual receptor rhodopsin. It has been well documented that lipid headgroups, polyunsaturated tails, and the concentration of cholesterol in membranes all play a role in the function of rhodopsin. Recently, we used all-atom simulations to demonstrate that different lipid species have preferential interactions and possible binding sites on rhodopsin's surface, consistent with experiment. However, the limited timescales of the simulations meant that the statistical uncertainty of these results was substantial. Accordingly, we present here 32 independent 1.6 μs coarse-grained simulations exploring lipids and cholesterols surrounding rhodopsin and opsin, in lipid bilayers mimicking those found naturally. Our results agree with those found experimentally and in previous simulations, but with far better statistical certainty. The results demonstrate the value of combining all-atom and coarse-grained models with experiment to provide a well-rounded view of lipid-protein interactions.

Solid-state NMR spectroscopy to study protein-lipid interactions

The appropriate lipid environment is crucial for the proper function of membrane proteins. There is a tremendous variety of lipid molecules in the membrane and so far it is often unclear which component of the lipid matrix is essential for the function of a respective protein. Lipid molecules and proteins mutually influence each other; parameters such as acyl chain order, membrane thickness, membrane elasticity, permeability, lipid-domain and annulus formation are strongly modulated by proteins. More recent data also indicates that the influence of proteins goes beyond a single annulus of next-neighbor boundary lipids. Therefore, a mesoscopic approach to membrane lipid-protein interactions in terms of elastic membrane deformations has been developed. Solid-state NMR has greatly contributed to the understanding of lipid-protein interactions and the modern view of biological membranes. Methods that detect the influence of proteins on the membrane as well as direct lipid-protein interactions have been developed and are reviewed here. Examples for solid-state NMR studies on the interaction of Ras proteins, the antimicrobial peptide protegrin-1, the G protein-coupled receptor rhodopsin, and the K(+) channel KcsA are discussed. This article is part of a Special Issue entitled Tools to study lipid functions.

Genetically-encoded molecular probes to study G protein-coupled receptors

To facilitate structural and dynamic studies of G protein-coupled receptor (GPCR) signaling complexes, new approaches are required to introduce informative probes or labels into expressed receptors that do not perturb receptor function. We used amber codon suppression technology to genetically-encode the unnatural amino acid, p-azido-L-phenylalanine (azF) at various targeted positions in GPCRs heterologously expressed in mammalian cells. The versatility of the azido group is illustrated here in different applications to study GPCRs in their native cellular environment or under detergent solubilized conditions. First, we demonstrate a cell-based targeted photocrosslinking technology to identify the residues in the ligand-binding pocket of GPCR where a tritium-labeled small-molecule ligand is crosslinked to a genetically-encoded azido amino acid. We then demonstrate site-specific modification of GPCRs by the bioorthogonal Staudinger-Bertozzi ligation reaction that targets the azido group using phosphine derivatives. We discuss a general strategy for targeted peptide-epitope tagging of expressed membrane proteins in-culture and its detection using a whole-cell-based ELISA approach. Finally, we show that azF-GPCRs can be selectively tagged with fluorescent probes. The methodologies discussed are general, in that they can in principle be applied to any amino acid position in any expressed GPCR to interrogate active signaling complexes.

Surface plasmon resonance applied to G protein-coupled receptors

G protein-coupled receptors (GPCR) are integral membrane proteins that transmit signals from external stimuli to the cell interior via activation of GTP-binding proteins (G proteins) thereby mediating key sensorial, hormonal, metabolic, immunological, and neurotransmission processes. Elucidating their structure and mechanism of interaction with extracellular and intracellular binding partners is of fundamental importance and highly relevant to rational design of new effective drugs. Surface plasmon resonance (SPR) has become a method of choice for studying biomolecular interactions at interfaces because measurements take place in real-time and do not require labeling of any of the interactants. However, due to the particular challenges imposed by the high hydrophobicity of membrane proteins and the great diversity of receptor-stimulating ligands, the application of this technique to characterize interactions of GPCR is still in the developmental phase. Here we give an overview of the principle of SPR and analyze current approaches for the preparation of the sensor chip surface, capture and stabilization of GPCR, and experimental design to characterize their interaction with ligands, G proteins and specific antibodies.

Unnatural amino acid mutagenesis of GPCRs using amber codon suppression and bioorthogonal labeling

To advance dynamic, temporal, and kinetic studies of the G protein-coupled receptor (GPCR) signalosome, new approaches are required to introduce non- or minimally perturbing labels or probes into expressed receptors. We report here a series of methods that are based on unnatural amino acid mutagenesis of GPCRs using amber codon suppression technology. We show that labeling reactions at genetically introduced keto moieties (p-acetyl-L-Phe/AcF and p-benzoyl-L-Phe/BzF) are not completely bioorthogonal due to protein oxidation, which causes high background. However, labeling reactions that target p-azido-L-Phe (azF) using the Staudinger-Bertozzi ligation and the strain-promoted alkyne-azide cycloaddition are bioorthogonal and are satisfactory for introducing labels or probes at near quantitative efficiency under mild labeling conditions. To our knowledge, this is the first report of a site-specific modification of an azF residue with a dibenzocyclooctyne-derivatized fluorophore. The methodologies we discuss are general, in that they can be applied in principle to any amino acid position in any expressed GPCR.

Retinal conformation governs pKa of protonated Schiff base in rhodopsin activation

We have explored the relationship between conformational energetics and the protonation state of the Schiff base in retinal, the covalently bound ligand responsible for activating the G protein-coupled receptor rhodopsin, using quantum chemical calculations. Guided by experimental structural determinations and large-scale molecular simulations on this system, we examined rotation about each bond in the retinal polyene chain, for both the protonated and deprotonated states that represent the dark and photoactivated states, respectively. Particular attention was paid to the torsional degrees of freedom that determine the shape of the molecule, and hence its interactions with the protein binding pocket. While most torsional degrees of freedom in retinal are characterized by large energetic barriers that minimize structural fluctuations under physiological temperatures, the C6-C7 dihedral defining the relative orientation of the β-ionone ring to the polyene chain has both modest barrier heights and a torsional energy surface that changes dramatically with protonation of the Schiff base. This surprising coupling between conformational degrees of freedom and protonation state is further quantified by calculations of the pKa as a function of the C6-C7 dihedral angle. Notably, pKa shifts of greater than two units arise from torsional fluctuations observed in molecular dynamics simulations of the full ligand-protein-membrane system. It follows that fluctuations in the protonation state of the Schiff base occur prior to forming the activated MII state. These new results shed light on important mechanistic aspects of retinal conformational changes that are involved in the activation of rhodopsin in the visual process.

Curvature forces in membrane lipid-protein interactions

Membrane biochemists are becoming increasingly aware of the role of lipid-protein interactions in diverse cellular functions. This review describes how conformational changes in membrane proteins, involving folding, stability, and membrane shape transitions, potentially involve elastic remodeling of the lipid bilayer. Evidence suggests that membrane lipids affect proteins through interactions of a relatively long-range nature, extending beyond a single annulus of next-neighbor boundary lipids. It is assumed the distance scale of the forces is large compared to the molecular range of action. Application of the theory of elasticity to flexible soft surfaces derives from classical physics and explains the polymorphism of both detergents and membrane phospholipids. A flexible surface model (FSM) describes the balance of curvature and hydrophobic forces in lipid-protein interactions. Chemically nonspecific properties of the lipid bilayer modulate the conformational energetics of membrane proteins. The new biomembrane model challenges the standard model (the fluid mosaic model) found in biochemistry texts. The idea of a curvature force field based on data first introduced for rhodopsin gives a bridge between theory and experiment. Influences of bilayer thickness, nonlamellar-forming lipids, detergents, and osmotic stress are all explained by the FSM. An increased awareness of curvature forces suggests that research will accelerate as structural biology becomes more closely entwined with the physical chemistry of lipids in explaining membrane structure and function.

Why GPCRs behave differently in cubic and lamellar lipidic mesophases

Recent successes in the crystallographic determination of structures of transmembrane proteins in the G protein-coupled receptor (GPCR) family have established the lipidic cubic phase (LCP) environment as the medium of choice for growing structure-grade crystals by the method termed "in meso". The understanding of in meso crystallogenesis is currently at a descriptive level. To enable an eventual quantitative, energy-based description of the nucleation and crystallization mechanism, we have examined the properties of the lipidic cubic phase system and the dynamics of the GPCR rhodopsin reconstituted into the LCP with coarse-grained molecular dynamics simulations with the Martini force-field. Quantifying the differences in the hydrophobic/hydrophilic exposure of the GPCR to lipids in the cubic and lamellar phases, we found that the highly curved geometry of the cubic phase provides more efficient shielding of the protein from unfavorable hydrophobic exposure, which leads to a lesser hydrophobic mismatch and less unfavorable hydrophobic-hydrophilic interactions between the protein and lipid-water interface in the LCP, compared to the lamellar phase. Since hydrophobic mismatch is considered a driving force for oligomerization, the differences in exposure mismatch energies between the LCP and the lamellar structures suggest that the latter provide a more favorable setting in which GPCRs can oligomerize as a prelude to nucleation and crystal growth. These new findings lay the foundation for future investigations of in meso crystallization mechanisms related to the transition from the LCP to the lamellar phase and studies aimed at an improved rational approach for generating structure-quality crystals of membrane proteins.

Structural determinants of the supramolecular organization of G protein-coupled receptors in bilayers

The G protein-coupled receptor (GPCR) rhodopsin self-assembles into supramolecular structures in native bilayers, but the structural determinants of receptor oligomerization are not known. We carried out multiple self-assembly coarse-grained molecular dynamics (CGMD) simulations of model membranes containing up to 64 molecules of the visual receptor rhodopsin over time scales reaching 100 μs. The simulations show strong preferential interaction modes between receptors. Two primary modes of receptor-receptor interactions are consistent with umbrella sampling/potential of mean force (PMF) calculations as a function of the distance between a pair of receptors. The preferential interfaces, involving helices (H) 1/8, 4/5 and 5, present no energy barrier to forming a very stable receptor dimer. Most notably, the PMFs show that the preferred rhodopsin dimer exists in a tail-to-tail conformation, with the interface comprising transmembrane H1/H2 and amphipathic H8 at the extracellular and cytoplasmic surfaces, respectively. This dimer orientation is in line with earlier electron microscopy, X-ray, and cross-linking experiments of rhodopsin and other GPCRs. Less stable interfaces, involving H4 and H6, have a free energy barrier for desolvation (delipidation) of the interfaces and appear to be designed to stabilize "lubricated" (i.e., lipid-coated) dimers. The overall CGMD strategy used here is general and can be applied to study the homo- and heterodimerization of GPCRs and other transmembrane proteins. Systematic extension of the work will deepen our understanding of the forces involved in the membrane organization of integral membrane proteins.

Characterization of adenosine receptor in its native environment: insights from molecular dynamics simulations of palmitoylated/glycosylated, membrane-integrated human A(2B) adenosine receptor

Selective A(2B) receptor antagonists and agonists may play a role in important pathologies such as gastrointestinal, neurological (i.e., Alzheimer disease and dementia) and hypersensitive disorders (i.e., asthma), diabetes, atherosclerosis, restenosis and cancer. Hence, it is regarded as a good target for the development of clinically useful agents. In this study, the effects of lipid bilayer, N-acetylglucosamine and S-palmitoyl on the dynamic behavior of A(2B)AR model is explored. Homology modeling, molecular docking and molecular dynamics simulations were performed to explore structural features of A(2B)AR in the presence of lipid bilayer. Twenty ns MD simulation was performed on the constructed model inserted in a hydrated lipid bilayer to examine stability of the best model. OSIP339391 as the most potent antagonist was docked in the active site of the model. Another MD simulation was performed on the ligand-protein complex to explore effects of the bilayer on this complex. A similar procedure was performed for the modified protein with N-acetylglucosamine and S-palmitoyl moieties in its structure. Phe173 and Glu174 located in EL2 were determined to be involved in ligand-receptor interactions through π-π stacking and hydrogen bonding. Asn254 was crucial to form hydrogen-bonding. The reliability of the model was assessed through docking using both commercial and synthetic antagonists and an r(2) of 0.70 was achieved. Our results show that molecular dynamics simulations of palmitoylated/glycosylated, membrane-integrated human A(2B)AR in its native environment is a possible approach and this model can be used for designing potent and selective A(2B)AR antagonists.

Quantifying the differential effects of DHA and DPA on the early events in visual signal transduction

A range of evidence from animal, clinical and epidemiological studies indicates that highly polyunsaturated acyl chains play important roles in development, cognition, vision and other aspects of neurological function. In a number of these studies n3 polyunsaturated fatty acids (PUFAs) appear to be more efficacious than n6 PUFAs. In a previous study of retinal rod outer segments obtained from rats raised on either an n3 adequate or deficient diet, we demonstrated that the replacement of 22:6n3 by 22:5n6 in the n3 deficient rats led to functional deficits in each step in the visual signaling process (Niu et al., 2004). In this study, we examined rhodopsin and phosphodiesterase function and acyl chain packing properties in membranes consisting of phosphatidylcholines with sn-1=18:0, and sn-2=22:6n3, 22:5n6, or 22:5n3 in order to determine if differences in function are due to the loss of one double bond or due to differences in double bond location. At 37 °C the n6 lipid shifted the equilibrium between the active metarhodopsin II (MII) state and inactive metarhodopsin I (MI) state towards MI. In addition, 22:5n6 reduced the rates of MII formation and MII-transducin complex formation by 2- and 6-fold, respectively. At a physiologically relevant level of rhodopsin light stimulation, the activity of phosphodiesterase was reduced by 50% in the 22:5n6 membrane, relative to either of the n3 membranes. Activity levels in the two n3 membranes were essentially identical. Ensemble acyl chain order was assessed with time-resolved fluorescence measurements of the membrane probe diphenylhexatriene (DPH). Analysis in terms of the orientational distribution of DPH showed that acyl chain packing in the two n3 membranes is quite similar, while in the 22:5n6 membrane there was considerably less packing disorder in the bilayer midplane. These results demonstrate that the n3 bond configuration uniquely optimizes the early steps in signaling via a mechanism which may involve acyl chain packing deep in the bilayer.

Structure and activation of rhodopsin

Rhodopsin is the first G-protein-coupled receptor (GPCR) with its three-dimensional structure solved by X-ray crystallography. The crystal structure of rhodopsin has revealed the molecular mechanism of photoreception and signal transduction in the visual system. Although several other GPCR crystal structures have been reported over the past few years, the rhodopsin structure remains an important model for understanding the structural and functional characteristics of other GPCRs. This review summarizes the structural features, the photoactivation, and the G protein signal transduction of rhodopsin.

Molecular simulations and solid-state NMR investigate dynamical structure in rhodopsin activation

Rhodopsin has served as the primary model for studying G protein-coupled receptors (GPCRs)-the largest group in the human genome, and consequently a primary target for pharmaceutical development. Understanding the functions and activation mechanisms of GPCRs has proven to be extraordinarily difficult, as they are part of a complex signaling cascade and reside within the cell membrane. Although X-ray crystallography has recently solved several GPCR structures that may resemble the activated conformation, the dynamics and mechanism of rhodopsin activation continue to remain elusive. Notably solid-state ((2))H NMR spectroscopy provides key information pertinent to how local dynamics of the retinal ligand change during rhodopsin activation. When combined with molecular mechanics simulations of proteolipid membranes, a new paradigm for the rhodopsin activation process emerges. Experiment and simulation both suggest that retinal isomerization initiates the rhodopsin photocascade to yield not a single activated structure, but rather an ensemble of activated conformational states. This article is part of a Special Issue entitled: Membrane protein structure and function.

Solid-state 2H NMR relaxation illuminates functional dynamics of retinal cofactor in membrane activation of rhodopsin

Rhodopsin is a canonical member of the family of G protein-coupled receptors, which transmit signals across cellular membranes and are linked to many drug interventions in humans. Here we show that solid-state (2)H NMR relaxation allows investigation of light-induced changes in local ps-ns time scale motions of retinal bound to rhodopsin. Site-specific (2)H labels were introduced into methyl groups of the retinal ligand that are essential to the activation process. We conducted solid-state (2)H NMR relaxation (spin-lattice, T(1Z), and quadrupolar-order, T(1Q)) experiments in the dark, Meta I, and Meta II states of the photoreceptor. Surprisingly, we find the retinylidene methyl groups exhibit site-specific differences in dynamics that change upon light excitation--even more striking, the C9-methyl group is a dynamical hotspot that corresponds to a crucial functional hotspot of rhodopsin. Following 11-cis to trans isomerization, the (2)H NMR data suggest the β-ionone ring remains in its hydrophobic binding pocket in all three states of the protein. We propose a multiscale activation mechanism with a complex energy landscape, whereby the photonic energy is directed against the E2 loop by the C13-methyl group, and toward helices H3 and H5 by the C5-methyl of the β-ionone ring. Changes in retinal structure and dynamics initiate activating fluctuations of transmembrane helices H5 and H6 in the Meta I-Meta II equilibrium of rhodopsin. Our proposals challenge the Standard Model whereby a single light-activated receptor conformation yields the visual response--rather an ensemble of substates is present, due to the entropy gain produced by photolysis of the inhibitory retinal lock.

Identification of cholesterol recognition amino acid consensus (CRAC) motif in G-protein coupled receptors

G-protein coupled receptors (GPCRs) are the largest class of molecules involved in signal transduction across membranes, and represent major targets in the development of novel drug candidates in all clinical areas. Membrane cholesterol has been reported to have an important role in the function of a number of GPCRs. Several structural features of proteins, believed to result in preferential association with cholesterol, have been recognized. Cholesterol recognition/interaction amino acid consensus (CRAC) sequence represents such a motif. Many proteins that interact with cholesterol have been shown to contain the CRAC motif in their sequence. We report here the presence of CRAC motifs in three representative GPCRs, namely, rhodopsin, the β(2)-adrenergic receptor, and the serotonin(1A) receptor. Interestingly, the function of these GPCRs has been previously shown to be dependent on membrane cholesterol. The presence of CRAC motifs in GPCRs indicates that interaction of cholesterol with GPCRs could be specific in nature. Further analysis shows that CRAC motifs are inherent characteristic features of the serotonin(1A) receptor and are conserved over natural evolution. These results constitute the first report of the presence of CRAC motifs in GPCRs and provide novel insight in the molecular nature of GPCR-cholesterol interaction.