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

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.

Progress in elucidating the structural and dynamic character of G Protein-Coupled Receptor oligomers for use in drug discovery

G Protein-Coupled Receptors (GPCRs) are the most targeted group of proteins for the development of therapeutic drugs. Until the last decade, structural information about this family of membrane proteins was relatively scarce, and their mechanisms of ligand binding and signal transduction were modeled on the assumption that GPCRs existed and functioned as monomeric entities. New crystal structures of native and engineered GPCRs, together with important biochemical and biophysical data that reveal structural details of the activation mechanism(s) of this receptor family, provide a valuable framework to improve dynamic molecular models of GPCRs with the ultimate goal of elucidating their allostery and functional selectivity. Since the dynamic movements of single GPCR protomers are likely to be affected by the presence of neighboring interacting subunits, oligomeric arrangements should be taken into account to improve the predictive ability of computer-assisted structural models of GPCRs for effective use in drug design.

Membrane cholesterol in the function and organization of G-protein coupled receptors

Cholesterol is an essential component of higher eukaryotic membranes and plays a crucial role in membrane organization, dynamics and function. The 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 a modulatory role in the function of a number of GPCRs. Two possible mechanisms have been previously suggested by which membrane cholesterol could influence the structure and function of GPCRs (i) through a direct/specific interaction with GPCRs, or (ii) through an indirect way by altering membrane physical properties in which the receptor is embedded, or due to a combination of both. Recently reported crystal structures of GPCRs have shown structural evidence of cholesterol binding sites. Against this backdrop, we recently proposed a novel mechanism by which membrane cholesterol could affect structure and function of GPCRs. According to our hypothesis, cholesterol binding sites in GPCRs could represent 'nonannular' binding sites. Interestingly, previous work from our laboratory has demonstrated that membrane cholesterol is required for the function of the serotonin(1A) receptor (a representative GPCR), which could be due to specific interaction of the receptor with cholesterol. Based on these results, we envisage that there could be specific/nonannular cholesterol binding site(s) in the serotonin(1A) receptor. We have analyzed putative cholesterol binding sites from protein databases in the serotonin(1A) receptor. Our analysis shows that cholesterol binding sites are inherent characteristic features of serotonin(1A) receptors and are conserved through natural evolution. Progress in deciphering molecular details of the GPCR-cholesterol interaction in the membrane would lead to better insight into our overall understanding of GPCR function in health and disease, thereby enhancing our ability to design better therapeutic strategies to combat diseases related to malfunctioning of 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.

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 (2)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 (2)H labels have been introduced into the methyl groups of retinal and solid-state (2)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 (2)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 beta-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 beta4 strand of the E2 loop and the side chains of Glu(122) and Trp(265) within the binding pocket. The solid-state (2)H NMR data are discussed with regard to the pathway of the energy flow in the receptor activation mechanism.

Lipid protein interactions couple protonation to conformation in a conserved cytosolic domain of G protein-coupled receptors

The visual photoreceptor rhodopsin is a prototypical class I (rhodopsin-like) G protein-coupled receptor. Photoisomerization of the covalently bound ligand 11-cis-retinal leads to restructuring of the cytosolic face of rhodopsin. The ensuing protonation of Glu-134 in the class-conserved D(E)RY motif at the C-terminal end of transmembrane helix-3 promotes the formation of the G protein-activating state. Using transmembrane segments derived from helix-3 of bovine rhodopsin, we show that lipid protein interactions play a key role in this cytosolic "proton switch." Infrared and fluorescence spectroscopic pK(a) determinations reveal that the D(E)RY motif is an autonomous functional module coupling side chain neutralization to conformation and helix positioning as evidenced by side chain to lipid headgroup Foerster resonance energy transfer. The free enthalpies of helix stabilization and hydrophobic burial of the neutral carboxyl shift the side chain pK(a) into the range typical of Glu-134 in photoactivated rhodopsin. The lipid-mediated coupling mechanism is independent of interhelical contacts allowing its conservation without interference with the diversity of ligand-specific interactions in class I G protein-coupled receptors.

Structural and dynamic effects of cholesterol at preferred sites of interaction with rhodopsin identified from microsecond length molecular dynamics simulations

An unresolved question about GPCR function is the role of membrane components in receptor stability and activation. In particular, cholesterol is known to affect the function of membrane proteins, but the details of its effect on GPCRs are still elusive. Here, we describe how cholesterol modulates the behavior of the TM1-TM2-TM7-helix 8(H8) functional network that comprises the highly conserved NPxxY(x)(5,6)F motif, through specific interactions with the receptor. The inferences are based on the analysis of microsecond length molecular dynamics (MD) simulations of rhodopsin in an explicit membrane environment. Three regions on the rhodopsin exhibit the highest cholesterol density throughout the trajectory: the extracellular end of TM7, a location resembling the high-density sterol area from the electron microscopy data; the intracellular parts of TM1, TM2, and TM4, a region suggested as the cholesterol binding site in the recent X-ray crystallography data on beta(2)-adrenergic GPCR; and the intracellular ends of TM2-TM3, a location that was categorized as the high cholesterol density area in multiple independent 100 ns MD simulations of the same system. We found that cholesterol primarily affects specific local perturbations of the helical TM domains such as the kinks in TM1, TM2, and TM7. These local distortions, in turn, relate to rigid-body motions of the TMs in the TM1-TM2-TM7-H8 bundle. The specificity of the effects stems from the nonuniform distribution of cholesterol around the protein. Through correlation analysis we connect local effects of cholesterol on structural perturbations with a regulatory role of cholesterol in the structural rearrangements involved in GPCR function.

Retinal conformation and dynamics in activation of rhodopsin illuminated by solid-state H NMR spectroscopy

Solid-state NMR spectroscopy gives a powerful avenue for investigating G protein-coupled receptors and other integral membrane proteins in a native-like environment. This article reviews the use of solid-state (2)H NMR to study the retinal cofactor of rhodopsin in the dark state as well as the meta I and meta II photointermediates. Site-specific (2)H NMR labels have been introduced into three regions (methyl groups) of retinal that are crucially important for the photochemical function of rhodopsin. Despite its phenomenal stability (2)H NMR spectroscopy indicates retinal undergoes rapid fluctuations within the protein binding cavity. The spectral lineshapes reveal the methyl groups spin rapidly about their three-fold (C(3)) axes with an order parameter for the off-axial motion of SC(3) approximately 0.9. For the dark state, the (2)H NMR structure of 11-cis-retinal manifests torsional twisting of both the polyene chain and the beta-ionone ring due to steric interactions of the ligand and the protein. Retinal is accommodated within the rhodopsin binding pocket with a negative pretwist about the C11=C12 double bond. Conformational distortion explains its rapid photochemistry and reveals the trajectory of the 11-cis to trans isomerization. In addition, (2)H NMR has been applied to study the retinylidene dynamics in the dark and light-activated states. Upon isomerization there are drastic changes in the mobility of all three methyl groups. The relaxation data support an activation mechanism whereby the beta-ionone ring of retinal stays in nearly the same environment, without a large displacement of the ligand. Interactions of the beta-ionone ring and the retinylidene Schiff base with the protein transmit the force of the retinal isomerization. Solid-state (2)H NMR thus provides information about the flow of energy that triggers changes in hydrogen-bonding networks and helix movements in the activation mechanism of the photoreceptor.

Superimposition of protein structures with dynamically weighted RMSD

In protein modeling, one often needs to superimpose a group of structures for a protein. A common way to do this is to translate and rotate the structures so that the square root of the sum of squares of coordinate differences of the atoms in the structures, called the root-mean-square deviation (RMSD) of the structures, is minimized. While it has provided a general way of aligning a group of structures, this approach has not taken into account the fact that different atoms may have different properties and they should be compared differently. For this reason, when superimposed with RMSD, the coordinate differences of different atoms should be evaluated with different weights. The resulting RMSD is called the weighted RMSD (wRMSD). Here we investigate the use of a special wRMSD for superimposing a group of structures with weights assigned to the atoms according to certain thermal motions of the atoms. We call such an RMSD the dynamically weighted RMSD (dRMSD). We show that the thermal motions of the atoms can be obtained from several sources such as the mean-square fluctuations that can be estimated by Gaussian network model analysis. We show that the superimposition of structures with dRMSD can successfully identify protein domains and protein motions, and that it has important implications in practice, e.g., in aligning the ensemble of structures determined by nuclear magnetic resonance.

FTIR analysis of GPCR activation using azido probes

We demonstrate the site-directed incorporation of an IR-active amino acid, p-azido-L-phenylalanine (azidoF, 1), into the G protein-coupled receptor rhodopsin using amber codon suppression technology. The antisymmetric stretch vibration of the azido group absorbs at approximately 2,100 cm(-1) in a clear spectral window and is sensitive to its electrostatic environment. We used FTIR difference spectroscopy to monitor the azido probe and show that the electrostatic environments of specific interhelical networks change during receptor activation.

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.

Atomistic insights into rhodopsin activation from a dynamic model

Rhodopsin, the light sensitive receptor responsible for blue-green vision, serves as a prototypical G protein-coupled receptor (GPCR). Upon light absorption, it undergoes a series of conformational changes that lead to the active form, metarhodopsin II (META II), initiating a signaling cascade through binding to the G protein transducin (G(t)). Here, we first develop a structural model of META II by applying experimental distance restraints to the structure of lumi-rhodopsin (LUMI), an earlier intermediate. The restraints are imposed by using a combination of biased molecular dynamics simulations and perturbations to an elastic network model. We characterize the motions of the transmembrane helices in the LUMI-to-META II transition and the rearrangement of interhelical hydrogen bonds. We then simulate rhodopsin activation in a dynamic model to study the path leading from LUMI to our META II model for wild-type rhodopsin and a series of mutants. The simulations show a strong correlation between the transition dynamics and the pharmacological phenotypes of the mutants. These results help identify the molecular mechanisms of activation in both wild type and mutant rhodopsin. While static models can provide insights into the mechanisms of ligand recognition and predict ligand affinity, a dynamic model of activation could be applicable to study the pharmacology of other GPCRs and their ligands, offering a key to predictions of basal activity and ligand efficacy.

Four-scale description of membrane sculpting by BAR domains

BAR domains are proteins that sense and sculpt curved membranes in cells, furnishing a relatively well-studied example of mechanisms employed in cellular morphogenesis. We report a computational study of membrane bending by BAR domains at four levels of resolution, described by 1), all-atom molecular dynamics; 2), residue-based coarse-graining (resolving single amino acids and lipid molecules); 3), shape-based coarse-graining (resolving overall protein and membrane shapes); and 4), a continuum elastic membrane model. Membrane sculpting performed by BAR domains collectively is observed in agreement with experiments. Different arrangements of BAR domains on the membrane surface are found to lead to distinct membrane curvatures and bending dynamics.

Regulation of membrane proteins by dietary lipids: effects of cholesterol and docosahexaenoic acid acyl chain-containing phospholipids on rhodopsin stability and function

Purified bovine rhodopsin was reconstituted into vesicles consisting of 1-stearoyl-2-oleoyl phosphatidylcholine or 1-stearoyl-2-docosahexaenoyl phosphatidylcholine with and without 30 mol % cholesterol. Rhodopsin stability was examined using differential scanning calorimetry (DSC). The thermal unfolding transition temperature (T(m)) of rhodopsin was scan rate-dependent, demonstrating the presence of a rate-limited component of denaturation. The activation energy of this kinetically controlled process (E(a)) was determined from DSC thermograms by four separate methods. Both T(m) and E(a) varied with bilayer composition. Cholesterol increased the T(m) both the presence and absence of docosahexaenoic acid acyl chains (DHA). In contrast, cholesterol lowered E(a) in the absence of DHA, but raised E(a) in the presence of 20 mol % DHA-containing phospholipid. The relative acyl chain packing order was determined from measurements of diphenylhexatriene fluorescence anisotropy decay. The T(m) for thermal unfolding was inversely related to acyl chain packing order. Rhodopsin kinetic stability (E(a)) was reduced in highly ordered or disordered membranes. Maximal kinetic stability was found within the range of acyl chain order found in native bovine rod outer segment disk membranes. The results demonstrate that membrane composition has distinct effects on the thermal versus kinetic stabilities of membrane proteins, and suggests that a balance between membrane constituents with opposite effects on acyl chain packing, such as DHA and cholesterol, may be required for maximum protein stability.

Mechanism of signal propagation upon retinal isomerization: insights from molecular dynamics simulations of rhodopsin restrained by normal modes

As one of the best studied members of the pharmaceutically relevant family of G-protein-coupled receptors, rhodopsin serves as a prototype for understanding the mechanism of G-protein-coupled receptor activation. Here, we aim at exploring functionally relevant conformational changes and signal transmission mechanisms involved in its photoactivation brought about through a cis-trans photoisomerization of retinal. For this exploration, we propose a molecular dynamics simulation protocol that utilizes normal modes derived from the anisotropic network model for proteins. Deformations along multiple low-frequency modes of motion are used to efficiently sample collective conformational changes in the presence of explicit membrane and water environment, consistent with interresidue interactions. We identify two highly stable regions in rhodopsin, one clustered near the chromophore, the other near the cytoplasmic ends of transmembrane helices H1, H2, and H7. Due to redistribution of interactions in the neighborhood of retinal upon stabilization of the trans form, local structural rearrangements in the adjoining H3-H6 residues are efficiently propagated to the cytoplasmic end of these particular helices. In the structures obtained by our simulations, all-trans retinal interacts with Cys(167) on H4 and Phe(203) on H5, which were not accessible in the dark state, and exhibits stronger interactions with H5, while some of the contacts made (in the cis form) with H6 are lost.

Structure and dynamics of polyunsaturated hydrocarbon chains in lipid bilayers-significance for GPCR function

This review summarizes results of our recent solid-state NMR investigations on polyunsaturated 18:0-22:6n3-PC/PE/PS and 18:0-22:5n6-PC bilayers. Data on structure and dynamics of the polyunsaturated docosahexaenoyl (DHAn3, 22:6n3) and docosapentaenoyl chains (DPAn6, 22:5n6), investigated at physiological conditions, are reported. Lipid-protein interaction was studied on reconstituted bilayers containing the G-protein coupled membrane receptor (GPCR) rhodopsin as well as on rod outer segment (ROS) disk membranes prepared from bovine retinas. Results reveal surprisingly rapid conformational transitions of polyunsaturated chains and existence of weakly specific interactions of DHAn3 with spatially distinct sites on rhodopsin.

Advances in the Development and Application of Computational Methodologies for Structural Modeling of G-Protein Coupled Receptors

BACKGROUND: Despite the large amount of experimental data accumulated in the past decade on G-protein coupled receptor (GPCR) structure and function, understanding of the molecular mechanisms underlying GPCR signaling is still far from being complete, thus impairing the design of effective and selective pharmaceuticals. OBJECTIVE: Understanding of GPCR function has been challenged even further by more recent experimental evidence that several of these receptors are organized in the cell membrane as homo- or hetero-oligomers, and that they may exhibit unique pharmacological properties. Given the complexity of these new signaling systems, researcher's efforts are turning increasingly to molecular modeling, bioinformatics and computational simulations for mechanistic insights of GPCR functional plasticity. METHODS: We review here current advances in the development and application of computational approaches to improve prediction of GPCR structure and dynamics, thus enhancing current understanding of GPCR signaling. RESULTS/CONCLUSIONS: Models resulting from use of these computational approaches further supported by experiments are expected to help elucidate the complex allosterism that propagates through GPCR complexes, ultimately aiming at successful structure-based rational drug design.

Fluid mechanical matching of H+-ATP synthase subunit c-ring with lipid membranes revealed by 2H solid-state NMR

The F(1)F(o)-ATP synthase utilizes the transmembrane H(+) gradient for the synthesis of ATP. F(o) subunit c-ring plays a key role in transporting H(+) through F(o) in the membrane. We investigated the interactions of Escherichia coli subunit c with dimyristoylphosphatidylcholine (DMPC-d(54)) at lipid/protein ratios of 50:1 and 20:1 by means of (2)H-solid-state NMR. In the liquid-crystalline state of DMPC, the (2)H-NMR moment values and the order parameter (S(CD)) profile were little affected by the presence of subunit c, suggesting that the bilayer thickness in the liquid-crystalline state is matched to the transmembrane hydrophobic surface of subunit c. On the other hand, hydrophobic mismatch of subunit c with the lipid bilayer was observed in the gel state of DMPC. Moreover, the viscoelasticity represented by a square-law function of the (2)H-NMR relaxation was also little influenced by subunit c in the fluid phase, in contrast with flexible nonionic detergents or rigid additives. Thus, the hydrophobic matching of the lipid bilayer to subunit c involves at least two factors, the hydrophobic length and the fluid mechanical property. These findings may be important for the torque generation in the rotary catalytic mechanism of the F(1)F(o)-ATPse molecular motor.

Rapid incorporation of functional rhodopsin into nanoscale apolipoprotein bound bilayer (NABB) particles

Human apolipoprotein A-I (apo A-I) and its engineered constructs form discoidal lipid bilayers upon interaction with lipids in vitro. We now report the cloning, expression, and purification of apo A-I derived from zebrafish (Danio rerio), which combines with phospholipids to form similar discoidal bilayers and may prove to be superior to human apo A-I constructs for rapid reconstitution of seven-transmembrane helix receptors into nanoscale apolipoprotein bound bilayers (NABBs). We characterized NABBs by gel-filtration chromatography, native polyacrylamide gradient gel electrophoresis, UV-visible photobleaching difference spectroscopy, and fluorescence spectroscopy. We used electron microscopy to determine the stoichiometry and orientation of rhodopsin (rho)-containing NABBs prepared under various conditions and correlated stability and signaling efficiency of rho in NABBs with either one or two receptors. We discovered that the specific activity of G protein coupling for single rhos sequestered in individual NABBs was nearly identical with that of two rhos per NABB under conditions where stoichiometry and orientation could be inferred by electron microscopy imaging. Thermal stability of rho in NABBs was superior to that of rho in various commonly used detergents. We conclude that the NABB system using engineered zebrafish apo A-I is a native-like membrane mimetic system for G-protein-coupled receptors and discuss strategies for rapid incorporation of expressed membrane proteins into NABBs.

X-ray scattering and solid-state deuterium nuclear magnetic resonance probes of structural fluctuations in lipid membranes

Molecular fluctuations are a dominant feature of biomembranes. Cellular functions might rely on these properties in ways yet to be determined. This expectation is suggested by the fact that membrane deformation and rigidity, which govern molecular fluctuations, have been implicated in a number of cellular functions. However, fluctuations are more challenging to measure than average structures, which partially explain the small number of dedicated studies. Here, it is shown that two accessible laboratory methods, small-angle X-ray scattering and solid-state deuterium nuclear magnetic resonance (NMR), can be used as complementary probes of structural fluctuations in lipid membranes. In the case of X-ray scattering, membrane undulations give rise to logarithmically varying positional correlations that generate scattering peaks with long (power-law) tails. In the case of 2H NMR spectroscopy, fluctuations in the magnetic-coupling energies resulting from molecular motions cause relaxation among the various spin energy levels, and yield a powerful probe of orientational fluctuations of the lipid molecules. A unified interpretation of the combined scattering and 2H NMR data is provided by a liquid-crystalline membrane deformation model. The importance of this approach is that it is possible to utilize a microscopic model for positional and orientational observables to calculate bulk material properties of liquid-crystalline systems.

Solid-state 2H NMR spectroscopy of retinal proteins in aligned membranes

Solid-state 2H NMR spectroscopy gives a powerful avenue to investigating the structures of ligands and cofactors bound to integral membrane proteins. For bacteriorhodopsin (bR) and rhodopsin, retinal was site-specifically labeled by deuteration of the methyl groups followed by regeneration of the apoprotein. 2H NMR studies of aligned membrane samples were conducted under conditions where rotational and translational diffusion of the protein were absent on the NMR time scale. The theoretical lineshape treatment involved a static axial distribution of rotating C-C2H3 groups about the local membrane frame, together with the static axial distribution of the local normal relative to the average normal. Simulation of solid-state 2H NMR lineshapes gave both the methyl group orientations and the alignment disorder (mosaic spread) of the membrane stack. The methyl bond orientations provided the angular restraints for structural analysis. In the case of bR the retinal chromophore is nearly planar in the dark- and all-trans light-adapted states, as well upon isomerization to 13-cis in the M state. The C13-methyl group at the "business end" of the chromophore changes its orientation to the membrane upon photon absorption, moving towards W182 and thus driving the proton pump in energy conservation. Moreover, rhodopsin was studied as a prototype for G protein-coupled receptors (GPCRs) implicated in many biological responses in humans. In contrast to bR, the retinal chromophore of rhodopsin has an 11-cis conformation and is highly twisted in the dark state. Three sites of interaction affect the torsional deformation of retinal, viz. the protonated Schiff base with its carboxylate counterion; the C9-methyl group of the polyene; and the beta-ionone ring within its hydrophobic pocket. For rhodopsin, the strain energy and dynamics of retinal as established by 2H NMR are implicated in substituent control of activation. Retinal is locked in a conformation that is twisted in the direction of the photoisomerization, which explains the dark stability of rhodopsin and allows for ultra-fast isomerization upon absorption of a photon. Torsional strain is relaxed in the meta I state that precedes subsequent receptor activation. Comparison of the two retinal proteins using solid-state 2H NMR is thus illuminating in terms of their different biological functions.

Role of lipid charge in organization of water/lipid bilayer interface: insights via computer simulations

Anionic unsaturated lipid bilayers represent suitable model systems that mimic real cell membranes: they are fluid and possess a negative surface charge. Understanding of detailed molecular organization of water-lipid interfaces in such systems may provide an important insight into the mechanisms of proteins' binding to membranes. Molecular dynamics (MD) of full-atom hydrated lipid bilayers is one of the most powerful tools to address this problem in silico. Unfortunately, wide application of computational methods for such systems is limited by serious technical problems. They are mainly related to correct treatment of long-range electrostatic effects. In this study a physically reliable model of an anionic unsaturated bilayer of 1,2-dioleoyl-sn-glycero-3-phosphoserine (DOPS) was elaborated and subjected to long-term MD simulations. Electrostatic interactions were treated with two different algorithms: spherical cutoff function and particle-mesh Ewald summation (PME). To understand the role of lipid charge in the system behavior, similar calculations were also carried out for zwitterionic bilayer composed of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). It was shown that, for the charged DOPS bilayer, the PME protocol performs much better than the cutoff scheme. In the last case a number of artifacts in the structural organization of the bilayer were observed. All of them were attributed to inadequate treatment of electrostatic interactions of lipid headgroups with counterions. Electrostatic properties, along with structural and dynamic parameters, of both lipid bilayers were investigated. Comparative analysis of the MD data reveals that the water-lipid interface of the DOPC bilayer is looser than that for DOPS. This makes possible deeper penetration of water molecules inside the zwitterionic (DOPC) bilayer, where they strongly interact with carbonyls of lipids. This can lead to thickening of the membrane interface in zwitterionic as compared to negatively charged bilayers.

The signaling pathway of rhodopsin

The signal-transduction mechanism of rhodopsin was studied by molecular dynamics (MD) simulations of the high-resolution, inactive structure in an explicit membrane environment. The simulations were employed to calculate equal-time correlations of the fluctuating interaction energy of residue pairs. The resulting interaction-correlation matrix was used to determine a network that couples retinal to the cytoplasmic interface, where transducin binds. Two highly conserved motifs, D(E)RY and NPxxY, were found to have strong interaction correlation with retinal. MD simulations with restraints on each transmembrane helix indicated that the major signal-transduction pathway involves the interdigitating side chains of helices VI and VII. The functional roles of specific residues were elucidated by the calculated effect of retinal isomerization from 11-cis to all-trans on the residue-residue interaction pattern. It is suggested that Glu134 may act as a "signal amplifier" and that Asp83 may introduce a threshold to prevent background noise from activating rhodopsin.

Dynamic structure of retinylidene ligand of rhodopsin probed by molecular simulations

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

Flexibility of ras lipid modifications studied by 2H solid-state NMR and molecular dynamics simulations

Human posttranslationally modified N-ras oncogenes are known to be implicated in numerous human cancers. Here, we applied a combination of experimental and computational techniques to determine structural and dynamical details of the lipid chain modifications of an N-ras heptapeptide in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) membranes. Experimentally, 2H NMR spectroscopy was used to study oriented membranes that incorporated ras heptapeptides with two covalently attached perdeuterated hexadecyl chains. Atomistic molecular dynamics simulations of the same system were carried out over 100 ns including 60 DMPC and 4 ras molecules. Several structural and dynamical experimental parameters could be directly compared to the simulation. Experimental and simulated 2H NMR order parameters for the methylene groups of the ras lipid chains exhibited a systematic difference attributable to the absence of collective motions in the simulation and to geometrical effects. In contrast, experimental 2H NMR spin-lattice relaxation rates for Zeeman order were well reproduced in the simulation. The lack of slower collective motions in the simulation did not appreciably influence the relaxation rates at a Larmor frequency of 115.1 MHz. The experimental angular dependence of the 2H NMR relaxation rates with respect to the external magnetic field was also relatively well simulated. These relaxation rates showed a weak angular dependence, suggesting that the lipid modifications of ras are very flexible and highly mobile in agreement with the low order parameters. To quantify these results, the angular dependence of the 2H relaxation rates was calculated by an analytical model considering both molecular and collective motions. Peptide dynamics in the membrane could be modeled by an anisotropic diffusion tensor with principal values of Dparallel=2.1x10(9) s(-1) and Dperpendicular=4.5x10(5) s(-1). A viscoelastic fitting parameter describing the membrane elasticity, viscosity, and temperature was found to be relatively similar for the ras peptide and the DMPC host matrix. Large motional amplitudes and relatively short correlation times facilitate mixing and dispersal with the lipid bilayer matrix, with implications for the role of the full-length ras protein in signal transduction and oncogenesis.

Molecular Dynamics Simulations of Rhodopsin in Different One-Component Lipid Bilayers

Four 20 ns molecular dynamic simulations of rhodopsin embedded in different one-component lipid bilayers have been carried out to ascertain the importance of membrane lipids on the protein structure. Specifically, dimyristoyl phosphatidylcholine (DMPC), dipalmitoyl phosphatidylcholine (DPPC), palmitoyl oleoyl phosphatidylcholine (POPC), and palmitoyl linoleyl phosphatidylcholine (PLPC) lipid bilayers have been considered for the present work. The results reported here provide information on the hydrophobic matching between the protein and the bilayer and about the differential effects of the protein on the thickness of the different membranes. Furthermore, a careful analysis of the individual protein-lipid interactions permits the identification of residues that exhibit permanent interactions with atoms of the lipid environment that may putatively act as hooks of the protein to the membrane. The analysis of the trajectories also provides information about the effect of the bilayer on the protein structure, including secondary structural elements, salt bridges, and rigid-body motions.

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

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

Linkage Between the Intramembrane H-bond Network Around Aspartic Acid 83 and the Cytosolic Environment of Helix 8 in Photoactivated Rhodopsin

Understanding the coupling between conformational changes in the intramembrane domain and at the membrane-exposed surface of the bovine photoreceptor rhodopsin, a prototypical G protein-coupled receptor (GPCR), is crucial for the elucidation of molecular mechanisms in GPCR activation. Here, we have combined Fourier transform infrared (FTIR) and fluorescence spectroscopy to address the coupling between conformational changes in the intramembrane region around the retinal and the environment of helix 8, a putative cytosolic surface switch region in class-I GPCRs. Using FTIR/fluorescence cross-correlation we show specifically that surface alterations monitored by emission changes of fluorescein bound to Cys316 in helix 8 of rhodopsin are highly correlated with (i) H-bonding to Asp83 proximal of the retinal Schiff base but not to Glu122 close to the beta-ionone and (ii) with a metarhodopsin II (MII)-specific 1643 cm(-1) IR absorption change, indicative of a partial loss of secondary structure in helix 8 upon MII formation. These correlations are disrupted by limited C-terminal proteolysis but are maintained upon binding of a transducin alpha-subunit (G(talpha))-derived peptide, which stabilizes the MII state. Our results suggest that additional C-terminal cytosolic loop contacts monitored by an amide II absorption at 1557 cm(-1) play a functionally crucial role in keeping helix 8 in the position in which its environment is strongly coupled to the retinal-binding site near the Schiff base. In the intramembrane region, this coupling is mediated by the H-bonding network that connects Asp83 to the NPxxY(x)F motif preceding helix 8.

Predisposition of the dark state of rhodopsin to functional changes in structure

As the only member of the family of G-protein-coupled receptors for which atomic coordinates are available, rhodopsin is widely studied for insight into the molecular mechanism of G-protein-coupled receptor activation. The currently available structures refer to the inactive, dark state, of rhodopsin, rather than the light-activated metarhodopsin II (Meta II) state. A model for the Meta II state is proposed here by analyzing elastic network normal modes in conjunction with experimental data. Key mechanical features and interactions broken/formed in the proposed model are found to be consistent with the experimental data. The model is further tested by using a set of Meta II fluorescence decay rates measured to empirically characterize the deactivation of rhodopsin mutants. The model is found to correctly predict 93% of the experimentally observed effects in 119 rhodopsin mutants for which the decay rates and misfolding data have been measured, including a systematic analysis of Cys-->Ser replacements reported here. Based on the detailed comparison between model and experiments, a cooperative activation mechanism is deduced that couples retinal isomerization to concerted changes in conformation, facilitated by the intrinsic dynamics of rhodopsin. A global hinge site is identified near the retinal-binding pocket that ensures the efficient propagation of signals from the central transmembrane region to both cytoplasmic and extracellular ends. The predicted activation mechanism opens the transmembrane helices at the critical G-protein binding cytoplasmic domain. This model provides a detailed, mechanistic description of the activation process, extending experimental observations and yielding new insights for further tests.

How a small change in retinal leads to G-protein activation: initial events suggested by molecular dynamics calculations

Rhodopsin is the prototypical G-protein coupled receptor, coupling light activation with high efficiency to signaling molecules. The dark-state X-ray structures of the protein provide a starting point for consideration of the relaxation from initial light activation to conformational changes that may lead to signaling. In this study we create an energetically unstable retinal in the light activated state and then use molecular dynamics simulations to examine the types of compensation, relaxation, and conformational changes that occur following the cis-trans light activation. The results suggest that changes occur throughout the protein, with changes in the orientation of Helices 5 and 6, a closer interaction between Ala 169 on Helix 4 and retinal, and a shift in the Schiff base counterion that also reflects changes in sidechain interactions with the retinal. Taken together, the simulation is suggestive of the types of changes that lead from local conformational change to light-activated signaling in this prototypical system.

The physical chemistry of biological membranes

Physical chemistry explains the principles of self-organization of lipids into bilayers that form the matrix of biological membranes, and continuum theory of membrane energetics is successful in explaining many biological processes. With increasing sophistication of investigative tools, there is now a growing appreciation for lipid diversity and for the role of individual lipids and specific lipid-protein interactions in membrane structure and function.

Effect of different treatments of long-range interactions and sampling conditions in molecular dynamic simulations of rhodopsin embedded in a dipalmitoyl phosphatidylcholine bilayer

The present study analyzes the effect of the simulation conditions on the results of molecular dynamics simulations of G-protein coupled receptors (GPCRs) performed with an explicit lipid bilayer. Accordingly, the present work reports the analysis of different simulations of bovine rhodopsin embedded in a dipalmitoyl phosphatidylcholine (DPPC) lipid bilayer using two different sampling conditions and two different approaches for the treatment of long-range electrostatic interactions. Specifically, sampling was carried out either by using the statistical ensembles NVT or NPT (constant number of atoms, a pressure of 1 atm in all directions and fixed temperature), and the electrostatic interactions were treated either by using a twin-cutoff, or the particle mesh Ewald summation method (PME). The results of the present study suggest that the use of the NPT ensemble in combination with the PME method provide more realistic simulations. The use of NPT during the equilibration avoids the need of an a priori estimation of the box dimensions, giving the correct area per lipid. However, once the system is equilibrated, the simulations are irrespective of the sampling conditions used. The use of an electrostatic cutoff induces artifacts on both lipid thickness and the ion distribution, but has no direct effect on the protein and water molecules.

Curvature and hydrophobic forces drive oligomerization and modulate activity of rhodopsin in membranes

G protein-coupled receptors (GPCRs) are essential components of cellular signaling pathways. They are the targets of many current pharmaceuticals and are postulated to dimerize or oligomerize in cellular membranes in conjunction with their functional mechanisms. We demonstrate using fluorescence resonance energy transfer how association of rhodopsin occurs by long-range lipid-protein interactions due to geometrical forces, yielding greater receptor crowding. Constitutive association of rhodopsin is promoted by a reduction in membrane thickness (hydrophobic mismatch), but also by an increase in protein/lipid molar ratio, showing the importance of interactions extending well beyond a single annulus of boundary lipids. The fluorescence data correlate with the pK(a) for the MI-to-MII transition of rhodopsin, where deprotonation of the retinylidene Schiff base occurs in conjunction with helical movements leading to activation of the photoreceptor. A more dispersed membrane environment optimizes formation of the MII conformation that results in visual function. A flexible surface model explains both the dispersal and activation of rhodopsin in terms of bilayer curvature deformation (strain) and hydrophobic solvation energy. The bilayer stress is related to the lateral pressure profile in terms of the spontaneous curvature and associated bending rigidity. Transduction of the strain energy (frustration) of the bilayer drives protein oligomerization and conformational changes in a coupled manner. Our findings illuminate the physical principles of membrane protein association due to chemically nonspecific interactions in fluid lipid bilayers. Moreover, they yield a conceptual framework for understanding how the tightly regulated lipid compositions of cellular membranes influence their protein-mediated functions.

Structural dynamics of the mitochondrial ADP/ATP carrier revealed by molecular dynamics simulation studies

The mitochondrial adenosine diphosphate/adenosine triphosphate (ADP/ATP) carrier has been recently crystallized in complex with its specific inhibitor carboxyatractyloside (CATR). In the crystal structure, the six-transmembrane helix bundle that defines the nucleotide translocation pathway is closed on the matrix side due to sharp kinks in the odd-numbered helices. The closed conformation is further sealed by the loops protruding into the matrix that interact through an intricate network of charge-pairs. To gain insight into its structural dynamics we performed molecular dynamics (MD) simulation studies of the ADP/ATP carrier with and without its cocrystallized inhibitor. The two trajectories sampled a conformational space around two different configurations characterized by distinct salt-bridge networks with a significant shift from inter- to intrarepeat bonding on the matrix side in the absence of CATR. Analysis of the geometrical parameters defining the transmembrane helices showed that even-numbered helices can undergo a face rotation, whereas odd-numbered helices can undergo a change in the wobble angle with a conserved proline acting as molecular hinge. Our results provide new information on the dynamical properties of the ADP/ATP carrier and for the first time yield a detailed picture of a stable carrier conformation in absence of the inhibitor.

Dynamic models of G-protein coupled receptor dimers: indications of asymmetry in the rhodopsin dimer from molecular dynamics simulations in a POPC bilayer

Based on the growing evidence that G-protein coupled receptors (GPCRs) form homo- and hetero-oligomers, models of GPCR signaling are now considering macromolecular assemblies rather than monomers, with the homo-dimer regarded as the minimal oligomeric arrangement required for functional coupling to the G-protein. The dynamic mechanisms of such signaling assemblies are unknown. To gain some insight into properties of GPCR dimers that may be relevant to functional mechanisms, we study their current structural prototype, rhodopsin. We have carried out nanosecond time-scale molecular dynamics (MD) simulations of a rhodopsin dimer and compared the results to the monomer simulated in the same type of bilayer membrane model composed of an equilibrated unit cell of hydrated palmitoyl-oleoyl-phosphatidyl choline (POPC). The dynamic representation of the homo-dimer reveals the location of structural changes in several regions of the monomeric subunits. These changes appear to be more pronounced at the dimerization interface that had been shown to be involved in the activation process [Proc Natl Acad Sci USA 102:17495, 2005]. The results are consistent with a model of GPCR activation that involves allosteric modulation through a single GPCR subunit per dimer.

Ligand Binding and G-protein Coupling of the Serotonin1A Receptor in Cholesterol-enriched Hippocampal Membranes

The serotonin1A receptor is the most extensively studied member of the family of seven transmembrane domain G-protein coupled serotonin receptors. Since a large portion of such transmembrane receptors remains in contact with the membrane lipid environment, lipid–protein interactions assume importance in the structure-function analysis of such receptors. We have earlier reported the requirement of cholesterol for serotonin1A receptor function in native hippocampal membranes by specific depletion of cholesterol using methyl- β-cyclodextrin. In this paper, we monitored the serotonin1A receptor function in membranes that are enriched in cholesterol using a complex prepared from cholesterol and methyl-β-cyclodextrin. Our results indicate that ligand binding and receptor/G-protein interaction of the serotonin1A receptor do not exhibit significant difference in native and cholesterol-enriched hippocampal membranes indicating that further enrichment of cholesterol has little functional consequence on the serotonin1A receptor function. These results therefore provide new information on the effect of cholesterol enrichment on the hippocampal serotonin1A receptor function.

A role for direct interactions in the modulation of rhodopsin by omega-3 polyunsaturated lipids

Rhodopsin, the G protein-coupled receptor primarily responsible for sensing light, is found in an environment rich in polyunsaturated lipid chains and cholesterol. Biophysical experiments have shown that lipid unsaturation and cholesterol both have significant effects on rhodopsin's stability and function; omega-3 polyunsaturated chains, such as docosahexaenoic acid (DHA), destabilize rhodopsin and enhance the kinetics of the photocycle, whereas cholesterol has the opposite effect. Here, we use molecular dynamics simulations to investigate the possibility that polyunsaturated chains modulate rhodopsin stability and kinetics via specific direct interactions. By analyzing the results of 26 independent 100-ns simulations of dark-adapted rhodopsin, we found that DHA routinely forms tight associations with the protein in a small number of specific locations qualitatively different from the nonspecific interactions made by saturated chains and cholesterol. Furthermore, the presence of tightly packed DHA molecules tends to weaken the interhelical packing. These results are consistent with recent NMR work, which proposes that rhodopsin binds DHA, and they suggest a molecular rationale for DHA's effects on rhodopsin stability and kinetics.

Role of cholesterol in the function and organization of G-protein coupled receptors

Cholesterol is an essential component of eukaryotic membranes and plays a crucial role in membrane organization, dynamics and function. The modulatory role of cholesterol in the function of a number of membrane proteins is well established. This effect has been proposed to occur either due to a specific molecular interaction between cholesterol and membrane proteins or due to alterations in the membrane physical properties induced by the presence of cholesterol. The contemporary view regarding heterogeneity in cholesterol distribution in membrane domains that sequester certain types of membrane proteins while excluding others has further contributed to its significance in membrane protein function. The seven transmembrane domain G-protein coupled receptors (GPCRs) are among the largest protein families in mammals and represent approximately 2% of the total proteins coded by the human genome. Signal transduction events mediated by this class of proteins are the primary means by which cells communicate with and respond to their external environment. GPCRs therefore represent major targets for the development of novel drug candidates in all clinical areas. In view of their importance in cellular signaling, the interaction of cholesterol with such receptors represents an important determinant in functional studies of such receptors. This review focuses on the effect of cholesterol on the membrane organization and function of GPCRs from a variety of sources, with an emphasis on the more contemporary role of cholesterol in maintaining a domain-like organization of such receptors on the cell surface. Importantly, the recently reported role of cholesterol in the function and organization of the neuronal serotonin(1A) receptor, a representative of the GPCR family which is present endogenously in the hippocampal region of the brain, will be highlighted.

Molecular dynamics simulations of proteins in lipid bilayers

With recent advances in X-ray crystallography of membrane proteins promising many new high-resolution structures, molecular dynamics simulations will become increasingly valuable for understanding membrane protein function, as they can reveal the dynamic behavior concealed in the static structures. Dramatic increases in computational power, in synergy with more efficient computational methodologies, now allow us to carry out molecular dynamics simulations of any structurally known membrane protein in its native environment, covering timescales of up to 0.1 micros. At the frontiers of membrane protein simulations are ion channels, aquaporins, passive and active transporters, and bioenergetic proteins.