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.     

Solid-State H NMR spectroscopy of retinal proteins in aligned membranes

Solid-state ²H 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. ²H 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-C²H₃ 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 ²H 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 β-ionone ring within its hydrophobic pocket. For rhodopsin, the strain energy and dynamics of retinal as established by ²H 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 ²H NMR is thus illuminating in terms of their different biological functions.     

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

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

Deuterium NMR structure of retinal in the ground state of rhodopsin

The conformation of retinal bound to the G protein-coupled receptor rhodopsin is intimately linked to its photochemistry, which initiates the visual process. Site-directed deuterium ((2)H) NMR spectroscopy was used to investigate the structure of retinal within the binding pocket of bovine rhodopsin. Aligned recombinant membranes were studied containing rhodopsin that was regenerated with retinal (2)H-labeled at the C(5), C(9), or C(13) methyl groups by total synthesis. Studies were conducted at temperatures below the gel to liquid-crystalline phase transition of the membrane lipid bilayer, where rotational and translational diffusion of rhodopsin is effectively quenched. The experimental tilt series of (2)H NMR spectra were fit to a theoretical line shape analysis [Nevzorov, A. A., Moltke, S., Heyn, M. P., and Brown, M. F. (1999) J. Am. Chem. Soc. 121, 7636-7643] giving the retinylidene bond orientations with respect to the membrane normal in the dark state. Moreover, the relative orientations of pairs of methyl groups were used to calculate effective torsional angles between different planes of unsaturation of the retinal chromophore. Our results are consistent with significant conformational distortion of retinal, and they have important implications for quantum mechanical calculations of its electronic spectral properties. In particular, we find that the beta-ionone ring has a twisted 6-s-cis conformation, whereas the polyene chain is twisted 12-s-trans. The conformational strain of retinal as revealed by solid-state (2)H NMR is significant for explaining the quantum yields and mechanism of its ultrafast photoisomerization in visual pigments. This work provides a consensus view of the retinal conformation in rhodopsin as seen by X-ray diffraction, solid-state NMR spectroscopy, and quantum chemical calculations.     

Retinal Conformation Changes Rhodopsin’s Dynamic Ensemble

G protein-coupled receptors are vital membrane proteins that allosterically transduce biomolecular signals across the cell membrane. However, the process by which ligand binding induces protein conformation changes is not well understood biophysically. Rhodopsin, the mammalian dim-light receptor, is a unique test case for understanding these processes because of its switch-like activity; the ligand, retinal, is bound throughout the activation cycle, switching from inverse agonist to agonist after absorbing a photon. By contrast, the ligand-free opsin is outside the activation cycle and may behave differently. We find that retinal influences rhodopsin dynamics using an ensemble of all-atom molecular dynamics simulations that in aggregate contain 100 μs of sampling. Active retinal destabilizes the inactive state of the receptor, whereas the active ensemble was more structurally homogenous. By contrast, simulations of an active-like receptor without retinal present were much more heterogeneous than those containing retinal. These results suggest allosteric processes are more complicated than a ligand inducing protein conformational changes or simply capturing a shifted ensemble as outlined in classic models of allostery.     

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.     

Relative orientation between the beta-ionone ring and the polyene chain for the chromophore of rhodopsin in native membranes

Rotational resonance solid state nuclear magnetic resonance has been used to determine the relative orientation of the beta-ionone ring and the polyene chain of the chromophore 11-Z-retinylidene of rhodopsin in rod outer segment membranes from bovine retina. The bleached protein was regenerated with either 11-Z-[8,18-(13)C(2)]retinal or 11-Z-[8,16/17(13)C(2)]retinal, the latter having only one (13)C label at either of the chemically equivalent positions 16 and 17. Observation of (13)C selectively enriched in the ring methyl groups, C16/17, revealed alternative conformational states for the ring. Minor spectral components comprised around 26% of the chromophore. The major conformation (approximately 74%) has the chemical shift resolution required for measuring internuclear distances to (13)C in the retinal chain (C8) separately from each of these methyl groups. The resulting distance constraints, C8 to C16 and C17 (4.05 +/- 0.25 A) and from C8 to C18 (2.95 +/- 0.15 A), show that the major portion of retinylidene in rhodopsin has a twisted 6-s-cis conformation. The more precise distance measurement made here between C8 and C18 (2.95 A) predicts that the chain is twisted out-of-plane with respect to the ring by a modest amount (C5-C6-C7-C8 torsion angle = -28 +/- 7 degrees ).     

Retinal Flip in Rhodopsin Activation?

Rhodopsin is a well-characterized structural model of a G protein-coupled receptor. Photoisomerization of the covalently bound retinal triggers activation. Surprisingly, the x-ray crystal structure of the active Meta-II state has a 180° rotation about the long-axis of the retinal polyene chain. Unbiased microsecond-timescale all-atom molecular dynamics simulations show that the retinal cofactor can flip back to the orientation observed in the inactive state of rhodopsin under conditions favoring the Meta-I state. Our results provide, to our knowledge, the first evidence from molecular dynamics simulations showing how rotation of the retinal ligand within its binding pocket can occur in the activation mechanism of rhodopsin.     

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.     

Agonist-induced conformational changes in bovine rhodopsin: insight into activation of G-protein-coupled receptors

Activation of G-protein-coupled receptors (GPCRs) is initiated by conformational changes in the transmembrane (TM) helices and the intra- and extracellular loops induced by ligand binding. Understanding the conformational changes in GPCRs leading to activation is imperative in deciphering the role of these receptors in the pathology of diseases. Since the crystal structures of activated GPCRs are not yet available, computational methods and biophysical techniques have been used to predict the structures of GPCR active states. We have recently applied the computational method LITiCon to understand the ligand-induced conformational changes in β₂-adrenergic receptor by ligands of varied efficacies. Here we report a study of the conformational changes associated with the activation of bovine rhodopsin for which the crystal structure of the inactive state is known. Starting from the inactive (dark) state, we have predicted the TM conformational changes that are induced by the isomerization of 11-cis retinal to all-trans retinal leading to the fully activated state, metarhodopsin II. The predicted active state of rhodopsin satisfies all of the 30 known experimental distance constraints. The predicted model also correlates well with the experimentally observed conformational switches in rhodopsin and other class A GPCRs, namely, the breaking of the ionic lock between R135^3.50 at the intracellular end of TM3 (part of the DRY motif) and E247^6.30 on TM6, and the rotamer toggle switch on W265^6.48 on TM6. We observe that the toggling of the W265^6.48 rotamer modulates the bend angle of TM6 around the conserved proline. The rotamer toggling is facilitated by the formation of a water wire connecting S298^7.45, W265^6.48 and H211^5.46. As a result, the intracellular ends of TMs 5 and 6 move outward from the protein core, causing large conformational changes at the cytoplasmic interface. The predicted outward movements of TM5 and TM6 are in agreement with the recently published crystal structure of opsin, which is proposed to be close to the active-state structure. In the predicted active state, several residues in the intracellular loops, such as R69, V139^3.54, T229, Q237, Q239, S240, T243 and V250^6.33, become more water exposed compared to the inactive state. These residues may be involved in mediating the conformational signal from the receptor to the G protein. From mutagenesis studies, some of these residues, such as V139^3.54, T229 and V250^6.33, are already implicated in G-protein activation. The predicted active state also leads to the formation of new stabilizing interhelical hydrogen-bond contacts, such as those between W265^6.48 and H211^5.46 and E122^3.37 and C167^4.56. These hydrogen-bond contacts serve as potential conformational switches offering new opportunities for future experimental investigations. The calculated retinal binding energy surface shows that binding of an agonist makes the receptor dynamic and flexible and accessible to many conformations, while binding of an inverse agonist traps the receptor in the inactive state and makes the other conformations inaccessible.     

Molecular dynamics simulations of retinal in rhodopsin: from the dark-adapted state towards lumirhodopsin

The formation of photointermediates and conformational changes observed in the retinal chromophore of bilayer-embedded rhodopsin during the early steps of the protein activation have been studied by molecular dynamics (MD) simulation. In particular, the lysine-bound retinal has been examined, focusing on its conformation in the dark-adapted state (10 ns) and on the early steps after the isomerization of the 11-cis bond to trans (up to 10 ns). The parametrization for the chromophore is based on a recent quantum study [Sugihara, M., Buss, V., Entel, P., Elstner, M., and Frauenheim, T. (2002) Biochemistry 41, 15259-15266] and shows good conformational agreement with recent experimental results. The isomerization, induced by switching the function governing the dihedral angle for the C11=C12 bond, was repeated with several different starting conformations. From the repeated simulations, it is shown that the retinal model exhibits a conserved activation pattern. The conformational changes are sequential and propagate outward from the C11=C12 bond, starting with isomerization of the C11=C12 bond, then a rotation of methyl group C20, and followed by increased fluctuations at the beta-ionone ring. The dynamics of these changes suggest that they are linked with photointermediates observed by spectroscopy. The exact moment when these events occur after the isomerization is modulated by the starting conformation, suggesting that retinal isomerizes through multiple pathways that are slightly different. The amplitudes of the structural fluctuations observed for the protein in the dark-adapted state and after isomerization of the retinal are similar, suggesting a subtle mechanism for the transmission of information from the chromophore to the protein.     

Electron crystallography reveals the structure of metarhodopsin I

Rhodopsin is the prototypical G protein-coupled receptor, responsible for detection of dim light in vision. Upon absorption of a photon, rhodopsin undergoes structural changes, characterised by distinct photointermediates. Currently, only the ground-state structure has been described. We have determined a density map of a photostationary state highly enriched in metarhodopsin I, to a resolution of 5.5 A in the membrane plane, by electron crystallography. The map shows density for helix 8, the cytoplasmic loops, the extracellular plug, all tryptophan residues, an ordered cholesterol molecule and the beta-ionone ring. Comparison of this map with X-ray structures of the ground state reveals that metarhodopsin I formation does not involve large rigid-body movements of helices, but there is a rearrangement close to the bend of helix 6, at the level of the retinal chromophore. There is no gradual build-up of the large conformational change known to accompany metarhodopsin II formation. The protein remains in a conformation similar to that of the ground state until late in the photobleaching process.     

Lipid Gymnastics: Evidence of Complete Acyl Chain Reversal in Oxidized Phospholipids from Molecular Simulations

In oxidative environments, biomembranes contain oxidized lipids with short, polar acyl chains. Two stable lipid oxidation products are PoxnoPC and PazePC. PoxnoPC has a carbonyl group, and PazePC has an anionic carboxyl group pendant at the end of the short, oxidized acyl chain. We have used MD simulations to explore the possibility of complete chain reversal in OXPLs in POPC-OXPL mixtures. The polar AZ chain of PazePC undergoes chain reversal without compromising the lipid bilayer integrity at concentrations up to 25% OXPL, and the carboxyl group points into the aqueous phase. Counterintuitively, the perturbation of overall membrane structural and dynamic properties is stronger for PoxnoPC than for PazePC. This is because of the overall condensing and ordering effect of sodium ions bound strongly to the lipids in the PazePC simulations. The reorientation of AZ chain is similar for two different lipid force fields. This work provides the first molecular evidence of the “extended lipid conformation” in phospholipid membranes. The chain reversal of PazePC lipids decorates the membrane interface with reactive, negatively charged functional groups. Such chain reversal is likely to exert a profound influence on the structure and dynamics of biological membranes, and on membrane-associated biological processes.     

Studies on the conformational state of the chromophore group (11-cis-retinal) in rhodopsin by computer molecular simulation methods

The molecular dynamics of the rhodopsin chromophore (11-cis-retinal) has been followed over a 3-ns path, whereby 3 × 10⁶ discrete conformational states of the molecule were recorded. It is shown that within a short time, 0.3–0.4 ns from the start of simulation, the retinal β-ionone ring rotates about the C6–C7 bond through ～60° relative to the initial configuration, and the whole chromophore becomes twisted. The results of ab initio quantum chemical calculations indicate that for the final conformation of the chromophore center (t = 3 ns) the rhodopsin absorption maximum is shifted by 10 nm toward longer wavelengths as compared with the initial state (t = 0). In other words, the energy of transition of such a system into the excited singlet state S1 upon photon capture will be lower than that for the molecule where the β-ionone ring of the chromophore is coplanar to its polyene chain.     

Crystallographic Study of the LUMI Intermediate of Squid Rhodopsin

Upon absorption of light, the retinal chromophore in rhodopsin isomerizes from the 11-cis to the trans configuration, initiating a photoreaction cycle. The primary photoreaction state, bathorhodopsin (BATHO), relaxes thermally through lumirhodopsin (LUMI) into a photoactive state, metarhodopsin (META), which stimulates the conjugated G-protein. Previous crystallographic studies of squid and bovine rhodopsins have shown that the structural change in the primary photoreaction of squid rhodopsin is considerably different from that observed in bovine rhodopsin. It would be expected that there is a fundamental difference in the subsequent thermal relaxation process between vertebrate and invertebrate rhodopsins. In this work, we performed crystallographic analyses of the LUMI state of squid rhodopsin using the P62 crystal. When the crystal was illuminated at 100 K with blue light, a half fraction of the protein was converted into BATHO. This reaction state relaxed into LUMI when the illuminated crystal was warmed in the dark to 170 K. It was found that, whereas trans retinal is largely twisted in BATHO, it takes on a more planar configuration in LUMI. This relaxation of retinal is accompanied by reorientation of the Schiff base NH bond, the hydrogen-bonding partner of which is switched to Asn185 in LUMI. Unlike bovine rhodopsin, the BATHO-to-LUMI transition in squid rhodopsin was accompanied by no significant change in the position/orientation of the beta-ionone ring of retinal.     

Component volumes of unsaturated phosphatidylcholines in fluid bilayers: a densitometric study

The specific volumes of six 1,2-diacylphosphatidylcholines with monounsaturated acyl chains (diCn:1PC, n=14-24 is the even number of acyl chain carbons) in fluid bilayers in multilamellar vesicles dispersed in H(2)O were determined by the vibrating tube densitometry as a function of temperature. From the data obtained with diCn:1PC (n=14-22) vesicles in combination with the densitometric data from Tristram-Nagle et al. [Tristram-Nagle, S., Petrache, H.I., Nagle, J.F., 1998. Structure and interactions of fully hydrated dioleoylphosphatidylcholine bilayers. Biophys. J. 75, 917-925.] and Koenig and Gawrisch [Koenig, B.W., Gawrisch, K., 2005. Specific volumes of unsaturated phosphatidylcholines in the liquid crystalline lamellar phase. Biochim. Biophys. Acta 1715, 65-70.], the component volumes of phosphatidylcholines in fully hydrated fluid bilayers at 30 degrees C were obtained. The volume of the acyl chain CH and CH(2) group is V(CH)=22.30 A(3) and V(CH2) =A(3), respectively. The volume of the headgroup including the glyceryl and acyl carbonyls, V(H), and the ratio of acyl chain methyl and methylene group volumes, r=V(CH3):V(CH2) are linearly interdependent: V(H)=a-br, where a=434.41 A(3) and b=-55.36 A(3) at 30 degrees C. From the temperature dependencies of component volumes, their isobaric thermal expansivities (alpha(X)=V(X)(-1)(partial differential V(X)/ partial differential T) where X=CH(2), CH, or H were calculated: alpha(CH2)=118.4x10(-5)K(-1), alpha(CH)=71.0x10(-5)K(-1), alpha(H)=7.9x10(-5)K(-1) (for r=2) and alpha(H)=9.6x10(-5)K(-1) (for r=1.9). The specific volume of diC24:1PC changes at the main gel-fluid phase transition temperature, t(m)=26.7 degrees C, by 0.0621 ml/g, its specific volume is 0.9561 and 1.02634 ml/g at 20 and 30 degrees C, respectively, and its isobaric thermal expansivity alpha=68.7x10(-5) and 109.2x10(-5)K(-1) below and above t(m), respectively. The component volumes and thermal expansivities obtained can be used for the interpretation of X-ray and neutron scattering and diffraction experiments and for the guiding and testing molecular dynamics simulations of phosphatidylcholine bilayers in the fluid state.     

Binding and interactions of L-BABP to lipid membranes studied by molecular dynamic simulations

Chicken liver bile acid-binding protein (L-BABP) is a member of the fatty acid-binding proteins super family. The common fold is a β-barrel of ten strands capped with a short helix-loop-helix motif called portal region, which is involved in the uptake and release of non-polar ligands. Using multiple-run molecular dynamics simulations we studied the interactions of L-BABP with lipid membranes of anionic and zwitterionic phospholipids. The simulations were in agreement with our experimental observations regarding the electrostatic nature of the binding and the conformational changes of the protein in the membrane. We observed that L-BABP migrated from the initial position in the aqueous bulk phase to the interface of anionic lipid membranes and established contacts with the head groups of phospholipids through the side of the barrel that is opposite to the portal region. The conformational changes in the protein occurred simultaneously with the binding to the membrane. Remarkably, these conformational changes were observed in the portal region which is opposite to the zone where the protein binds directly to the lipids. The protein was oriented with its macrodipole aligned in the configuration of lowest energy within the electric field of the anionic membrane, which indicates the importance of the electrostatic interactions to determine the preferred orientation of the protein. We also identified this electric field as the driving force for the conformational change. For all the members of the fatty acid-binding protein family, the interactions with lipid membranes is a relevant process closely related to the uptake, release and transfer of the ligand. The observations presented here suggest that the ligand transfer might not necessarily occur through the domain that directly interacts with the lipid membrane. The interactions with the membrane electric field that determine orientation and conformational changes described here can also be relevant for other peripheral proteins.     

The Photocycle and Proton Translocation Pathway in a Cyanobacterial Ion-Pumping Rhodopsin

The genome of thylakoidless cyanobacterium Gloeobacter violaceus encodes a fast-cycling rhodopsin capable of light-driven proton transport. We characterize the dark state, the photocycle, and the proton translocation pathway of GR spectroscopically. The dark state of GR contains predominantly all-trans-retinal and, similar to proteorhodopsin, does not show the light/dark adaptation. We found an unusually strong coupling between the conformation of the retinal and the site of Glu¹³², the homolog of Asp⁹⁶ of BR. Although the photocycle of GR is similar to that of proteorhodopsin in general, it differs in accumulating two intermediates typical for BR, the L-like and the N-like states. The latter state has a deprotonated cytoplasmic proton donor and is spectrally distinct from the strongly red-shifted N intermediate known for proteorhodopsin. The proton uptake precedes the release and occurs during the transition to the O intermediate. The proton translocation pathway of GR is similar to those of other proton-pumping rhodopsins, involving homologs of BR Schiff base proton acceptor and donor Asp⁸⁵ and Asp⁹⁶ (Asp¹²¹ and Glu¹³²). We assigned a pair of FTIR bands (positive at 1749 cm⁻¹ and negative at 1734 cm⁻¹) to the protonation and deprotonation, respectively, of these carboxylic acids.     

Structure and fluctuations of charged phosphatidylserine bilayers in the absence of salt

Using x-ray diffraction and NMR spectroscopy, we present structural and material properties of phosphatidylserine (PS) bilayers that may account for the well documented implications of PS headgroups in cell activity. At 30 degrees C, the 18-carbon monounsaturated DOPS in the fluid state has a cross-sectional area of 65.3 A(2) which is remarkably smaller than the area 72.5 A(2) of the DOPC analog, despite the extra electrostatic repulsion expected for charged PS headgroups. Similarly, at 20 degrees C, the 14-carbon disaturated DMPS in the gel phase has an area of 40.8 A(2) vs. 48.1 A(2) for DMPC. This condensation of area suggests an extra attractive interaction, perhaps hydrogen bonding, between PS headgroups. Unlike zwitterionic lipids, stacks of PS bilayers swell indefinitely as water is added. Data obtained for osmotic pressure versus interbilayer water spacing for fluid phase DOPS are well fit by electrostatic interactions calculated for the Gouy-Chapman regime. It is shown that the electrostatic interactions completely dominate the fluctuational pressure. Nevertheless, the x-ray data definitively exhibit the effects of fluctuations in fluid phase DOPS. From our measurements of fluctuations, we obtain the product of the bilayer bending modulus K(C) and the smectic compression modulus B. At the same interbilayer separation, the interbilayer fluctuations are smaller in DOPS than for DOPC, showing that B and/or K(C) are larger. Complementing the x-ray data, (31)P-chemical shift anisotropy measured by NMR suggest that the DOPS headgroups are less sensitive to osmotic pressure than DOPC headgroups, which is consistent with a larger K(C) in DOPS. Quadrupolar splittings for D(2)O decay less rapidly with increasing water content for DOPS than for DOPC, indicating greater perturbation of interlamellar water and suggesting a greater interlamellar hydration force in DOPS. Our comparisons between bilayers of PS and PC lipids with the same chains and the same temperature enable us to focus on the effects of these headgroups on bilayer properties.     

Acyl Transfer from Membrane Lipids to Peptides Is a Generic Process

The generality of acyl transfer from phospholipids to membrane-active peptides has been probed using liquid chromatography–mass spectrometry analysis of peptide–lipid mixtures. The peptides examined include melittin, magainin II, PGLa, LAK1, LAK3 and penetratin. Peptides were added to liposomes with membrane lipid compositions ranging from pure phosphatidylcholine (PC) to mixtures of PC with phosphatidylethanolamine, phosphatidylserine or phosphatidylglycerol. Experiments were typically conducted at pH 7.4 at modest salt concentrations (90 mM NaCl). In favorable cases, lipidated peptides were further characterized by tandem mass spectrometry methods to determine the sites of acylation. Melittin and magainin II were the most reactive peptides, with significant acyl transfer detected under all conditions and membrane compositions. Both peptides were lipidated at the N-terminus by transfer from PC, phosphatidylethanolamine, phosphatidylserine or phosphatidylglycerol, as well as at internal sites: lysine for melittin; serine and lysine for magainin II. Acyl transfer could be detected within 3 h of melittin addition to negatively charged membranes. The other peptides were less reactive, but for each peptide, acylation was found to occur in at least one of the conditions examined. The data demonstrate that acyl transfer is a generic process for peptides bound to membranes composed of diacylglycerophospholipids. Phospholipid membranes cannot therefore be considered as chemically inert toward peptides and by extension proteins.