George Khelashvili

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,6F 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 β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. Proteins 2009.

Combined Monte Carlo and molecular dynamics simulation of hydrated 18:0 sphingomyelin-cholesterol lipid bilayers

We have carried out atomic level molecular dynamics and Monte Carlo simulations of hydrated 18:0 sphingomyelin (SM)-cholesterol (CHOL) bilayers at temperatures of 20 and 50 degrees C. The simulated systems each contained 266 SM, 122 CHOL, and 11861 water molecules. Each simulation was run for 10 ns under semi-isotropic pressure boundary conditions. The particle-mesh Ewald method was used for long-range electrostatic interactions. Properties of the systems were calculated over the final 3 ns. We compare the properties of 20 and 50 degrees C bilayer systems with each other, with experimental data, and with experimental and simulated properties of pure SM bilayers and dipalmitoyl phospatidyl choline (DPPC)-CHOL bilayers. The simulations reveal an overall similarity of both systems, despite the 30 degrees C temperature difference which brackets the pure SM main phase transition. The area per molecule, lipid chain order parameter profiles, atom distributions, and electron density profiles are all very similar for the two simulated systems. Consistent with simulations from our lab and others, we find strong intramolecular hydrogen bonding in SM molecules between the phosphate ester oxygen and the hydroxyl hydrogen atoms. We also find that cholesterol hydroxyl groups tend to form hydrogen bonds primarily with SM carbonyl, methyl, and amide moieties and to a lesser extent methyl and hydroxyl oxygens.

Transport domain unlocking sets the uptake rate of an aspartate transporter

Glutamate transporters terminate neurotransmission by clearing synaptically released glutamate from the extracellular space, allowing repeated rounds of signalling and preventing glutamate-mediated excitotoxicity. Crystallographic studies of a glutamate transporter homologue from the archaeon Pyrococcus horikoshii, GltPh, showed that distinct transport domains translocate substrates into the cytoplasm by moving across the membrane within a central trimerization scaffold. Here we report direct observations of these ‘elevator-like’ transport domain motions in the context of reconstituted proteoliposomes and physiological ion gradients using single-molecule fluorescence resonance energy transfer (smFRET) imaging. We show that GltPh bearing two mutations introduced to impart characteristics of the human transporter exhibits markedly increased transport domain dynamics, which parallels an increased rate of substrate transport, thereby establishing a direct temporal relationship between transport domain motion and substrate uptake. Crystallographic and computational investigations corroborated these findings by revealing that the ‘humanizing’ mutations favour structurally ‘unlocked’ intermediate states in the transport cycle exhibiting increased solvent occupancy at the interface between the transport domain and the trimeric scaffold.

Quantitative Modeling of Membrane Deformations by Multihelical Membrane Proteins: Application to G-Protein Coupled Receptors

The interpretation of experimental observations of the dependence of membrane protein function on the properties of the lipid membrane environment calls for a consideration of the energy cost of protein-bilayer interactions, including the protein-bilayer hydrophobic mismatch. We present a novel (to our knowledge) multiscale computational approach for quantifying the hydrophobic mismatch-driven remodeling of membrane bilayers by multihelical membrane proteins. The method accounts for both the membrane remodeling energy and the energy contribution from any partial (incomplete) alleviation of the hydrophobic mismatch by membrane remodeling. Overcoming previous limitations, it allows for radially asymmetric bilayer deformations produced by multihelical proteins, and takes into account the irregular membrane-protein boundaries. The approach is illustrated by application to two G-protein coupled receptors: rhodopsin in bilayers of different thickness, and the serotonin 5-HT(2A) receptor bound to pharmacologically different ligands. Analysis of the results identifies the residual exposure that is not alleviated by bilayer adaptation, and its quantification at specific transmembrane segments is shown to predict favorable contact interfaces in oligomeric arrays. In addition, our results suggest how distinct ligand-induced conformations of G-protein coupled receptors may elicit different functional responses through differential effects on the membrane environment.

Ligand-dependent conformations and dynamics of the serotonin 5-HT(2A) receptor determine its activation and membrane-driven oligomerization properties.

From computational simulations of a serotonin 2A receptor (5-HT2AR) model complexed with pharmacologically and structurally diverse ligands we identify different conformational states and dynamics adopted by the receptor bound to the full agonist 5-HT, the partial agonist LSD, and the inverse agonist Ketanserin. The results from the unbiased all-atom molecular dynamics (MD) simulations show that the three ligands affect differently the known GPCR activation elements including the toggle switch at W6.48, the changes in the ionic lock between E6.30 and R3.50 of the DRY motif in TM3, and the dynamics of the NPxxY motif in TM7. The computational results uncover a sequence of steps connecting these experimentally-identified elements of GPCR activation. The differences among the properties of the receptor molecule interacting with the ligands correlate with their distinct pharmacological properties. Combining these results with quantitative analysis of membrane deformation obtained with our new method (Mondal et al, Biophysical Journal 2011), we show that distinct conformational rearrangements produced by the three ligands also elicit different responses in the surrounding membrane. The differential reorganization of the receptor environment is reflected in (i)-the involvement of cholesterol in the activation of the 5-HT2AR, and (ii)-different extents and patterns of membrane deformations. These findings are discussed in the context of their likely functional consequences and a predicted mechanism of ligand-specific GPCR oligomerization.

Cholesterol orientation and tilt modulus in DMPC bilayers.

We performed molecular dynamics (MD) simulations of hydrated bilayers containing mixtures of dimyristoylphosphatidylcholine (DMPC) and cholesterol at various ratios, to study the effect of cholesterol concentration on its orientation, and to characterize the link between cholesterol tilt and overall phospholipid membrane organization. The simulations show a substantial probability for cholesterol molecules to transiently orient perpendicular to the bilayer normal, and suggest that cholesterol tilt may be an important factor for inducing membrane ordering. In particular, we find that as cholesterol concentration increases (1-40% cholesterol) the average cholesterol orientation changes in a manner strongly (anti)correlated with the variation in membrane thickness. Furthermore, cholesterol orientation is found to be determined by the aligning force exerted by other cholesterol molecules. To quantify this aligning field, we analyzed cholesterol orientation using, to our knowledge, the first estimates of the cholesterol tilt modulus chi from MD simulations. Our calculations suggest that the aligning field that determines chi is indeed strongly linked to sterol composition. This empirical parameter (chi) should therefore become a useful quantitative measure to describe cholesterol interaction with other lipids in bilayers, particularly in various coarse-grained force fields.

GPCR-OKB

SUMMARY Rapid expansion of available data about G Protein Coupled Receptor (GPCR) dimers/oligomers over the past few years requires an effective system to organize this information electronically. Based on an ontology derived from a community dialog involving colleagues using experimental and computational methodologies, we developed the GPCR-Oligomerization Knowledge Base (GPCR-OKB). GPCR-OKB is a system that supports browsing and searching for GPCR oligomer data. Such data were manually derived from the literature. While focused on GPCR oligomers, GPCR-OKB is seamlessly connected to GPCRDB, facilitating the correlation of information about GPCR protomers and oligomers. AVAILABILITY AND IMPLEMENTATION The GPCR-OKB web application is freely available at http://www.gpcr-okb.org

Membrane Driven Spatial Organization of GPCRs

Spatial organization of G-protein coupled receptors (GPCRs) into dimers and higher order oligomers has been demonstrated in vitro and in vivo. The pharmacological readout was shown to depend on the specific interfaces, but why particular regions of the GPCR structure are involved, and how ligand-determined states change them remains unknown. Here we show why protein-membrane hydrophobic matching is attained upon oligomerization at specific interfaces from an analysis of coarse-grained molecular dynamics simulations of the spontaneous diffusion-interaction of the prototypical beta2-adrenergic (β2AR) receptors in a POPC lipid bilayer. The energy penalty from mismatch is significantly reduced in the spontaneously emerging oligomeric arrays, making the spatial organization of the GPCRs dependent on the pattern of mismatch in the monomer. This mismatch pattern is very different for β2AR compared to the highly homologous and structurally similar β1AR, consonant with experimentally observed oligomerization patterns of β2AR and β1AR. The results provide a mechanistic understanding of the structural context of oligomerization.

PIP2 regulates psychostimulant behaviors through its interaction with a membrane protein.

Phosphatidylinositol (4,5)-bisphosphate (PIP2) regulates the function of ion channels and transporters. Here, we demonstrate that PIP2 directly binds the human dopamine (DA) transporter (hDAT), a key regulator of DA homeostasis and a target of the psychostimulant amphetamine (AMPH). This binding occurs through electrostatic interactions with positively charged hDAT N-terminal residues and is shown to facilitate AMPH-induced, DAT-mediated DA efflux and the psychomotor properties of AMPH. Substitution of these residues with uncharged amino acids reduces hDAT-PIP2 interactions and AMPH-induced DA efflux, without altering the hDAT physiological function of DA uptake. We evaluated, for the first time, the significance of this interaction in vivo using locomotion as a behavioral assay in Drosophila melanogaster. Expression of mutated hDAT with reduced PIP2 interaction in Drosophila DA neurons impairs AMPH-induced locomotion without altering basal locomotion. We present the first demonstration of how PIP2 interactions with a membrane protein can regulate the behaviors of complex organisms.

Modeling Membrane Deformations and Lipid Demixing upon Protein-Membrane Interaction: The BAR Dimer Adsorption

We use a self-consistent mean-field theory, designed to investigate membrane reshaping and lipid demixing upon interaction with proteins, to explore BAR domains interacting with large patches of lipid membranes of heterogeneous compositions. The computational model includes contributions to the system free energy from electrostatic interactions and elastic energies of the membrane, as well as salt and lipid mixing entropies. The results from our simulation of a single adsorbing Amphiphysin BAR dimer indicate that it is capable of stabilizing a significantly curved membrane. However, we predict that such deformations will occur only for membrane patches that have the inherent propensity for high curvature, reflected in the tendency to create local distortions that closely match the curvature of the BAR dimer itself. Such favorable preconditioning for BAR-membrane interaction may be the result of perturbations such as local lipid demixing induced by the interaction, or of a prior insertion of the BAR domains amphiphatic N-helix. From our simulations it appears that local segregation of charged lipids under the influence of the BAR dimer cannot produce high enough asymmetry between bilayer leaflets to induce significant bending. In the absence of additional energy contributions that favor membrane asymmetry, the membrane will remain nearly flat upon single BAR dimer adsorption, relative to the undulation expected from thermal fluctuations. Thus, we conclude that the N-helix insertions have a critical mechanistic role in the local perturbation and curving of the membrane, which is then stabilized by the electrostatic interaction with the BAR dimer. We discuss how these results can be used to estimate the tendency of BARs to bend membranes in terms of a spatially nonisotropic spontaneous curvature.