Membrane driven spatial organization of GPCRs

Dynamic Dimerization of Chemokine Receptors and Potential Inhibitory Role of Their Truncated Isoforms Revealed through Combinatorial Prediction

Chemokine receptors play crucial roles in fundamental biological processes. Their malfunction may result in many diseases, including cancer, autoimmune diseases, and HIV. The oligomerization of chemokine receptors holds significant functional implications that directly affect their signaling patterns and pharmacological responses. However, the oligomerization patterns of many chemokine receptors remain poorly understood. Furthermore, several chemokine receptors have highly truncated isoforms whose functional role is not yet clear. Here, we computationally show homo- and heterodimerization patterns of four human chemokine receptors, namely CXCR2, CXCR7, CCR2, and CCR7, along with their interaction patterns with their respective truncated isoforms. By combining the neural network-based AlphaFold2 and physics-based protein-protein docking tool ClusPro, we predicted 15 groups of complex structures and assessed the binding affinities in the context of atomistic molecular dynamics simulations. Our results are in agreement with previous experimental observations and support the dynamic and diverse nature of chemokine receptor dimerization, suggesting possible patterns of higher-order oligomerization. Additionally, we uncover the strong potential of truncated isoforms to block homo- and heterodimerization of chemokine receptors, also in a dynamic manner. Our study provides insights into the dimerization patterns of chemokine receptors and the functional significance of their truncated isoforms.

How Neuromembrane Lipids Modulate Membrane Proteins: Insights from G-Protein-Coupled Receptors (GPCRs) and Receptor Tyrosine Kinases (RTKs)

Lipids play a diverse and critical role in cellular processes in all tissues. The unique lipid composition of nerve membranes is particularly interesting because it contains, among other things, polyunsaturated lipids, such as docosahexaenoic acid, which the body only gets through the diet. The crucial role of lipids in neurological processes, especially in receptor-mediated cell signaling, is emphasized by the fact that in many neuropathological diseases there are significant deviations in the lipid composition of nerve membranes compared to healthy individuals. The lipid composition of neuromembranes can significantly affect the function of receptors by regulating the physical properties of the membrane or by affecting specific interactions between receptors and lipids. In addition, it is worth noting that the ligand-binding pocket of many receptors is located inside the cell membrane, due to which lipids can even modulate the binding of ligands to their receptors. These mechanisms highlight the importance of lipids in the regulation of membrane receptor activation and function. In this article, we focus on two major protein families: G-protein-coupled receptors (GPCRs) and receptor tyrosine kinases (RTKs) and discuss how lipids affect their function in neuronal membranes, elucidating the basic mechanisms underlying neuronal function and dysfunction.

The World of GPCR dimers - Mapping dopamine receptor D2 homodimers in different activation states and configuration arrangements

G protein-coupled receptors (GPCRs) are known to dimerize, but the molecular and structural basis of GPCR dimers is not well understood. In this study, we developed a computational framework to generate models of symmetric and asymmetric GPCR dimers using different monomer activation states and identified their most likely interfaces with molecular details. We chose the dopamine receptor D2 (D2R) homodimer as a case study because of its biological relevance and the availability of structural information. Our results showed that transmembrane domains 4 and 5 (TM4 and TM5) are mostly found at the dimer interface of the D2R dimer and that these interfaces have a subset of key residues that are mostly nonpolar from TM4 and TM5, which was in line with experimental studies. In addition, TM2 and TM3 appear to be relevant for D2R dimers. In some cases, the inactive configuration is unaffected by the partnered protomer, whereas in others, the active protomer adopts the properties of an inactive receptor. Additionally, the β-arrestin configuration displayed the properties of an active receptor in the absence of an agonist, suggesting that a switch to another meta-state during dimerization occurred. Our findings are consistent with the experimental data, and this method can be adapted to study heterodimers and potentially extended to include additional proteins such as G proteins or β-arrestins. In summary, this approach provides insight into the impact of the conformational status of partnered protomers on the overall quaternary GPCR macromolecular structure and dynamics.

Cyclodextrins: Only Pharmaceutical Excipients or Full-Fledged Drug Candidates?

Cyclodextrins, representing a versatile family of cyclic oligosaccharides, have extensive pharmaceutical applications due to their unique truncated cone-shaped structure with a hydrophilic outer surface and a hydrophobic cavity, which enables them to form non-covalent host-guest inclusion complexes in pharmaceutical formulations to enhance the solubility, stability and bioavailability of numerous drug molecules. As a result, cyclodextrins are mostly considered as inert carriers during their medical application, while their ability to interact not only with small molecules but also with lipids and proteins is largely neglected. By forming inclusion complexes with cholesterol, cyclodextrins deplete cholesterol from cellular membranes and thereby influence protein function indirectly through alterations in biophysical properties and lateral heterogeneity of bilayers. In this review, we summarize the general chemical principles of direct cyclodextrin-protein interactions and highlight, through relevant examples, how these interactions can modify protein functions in vivo, which, despite their huge potential, have been completely unexploited in therapy so far. Finally, we give a brief overview of disorders such as Niemann-Pick type C disease, atherosclerosis, Alzheimer's and Parkinson's disease, in which cyclodextrins already have or could have the potential to be active therapeutic agents due to their cholesterol-complexing or direct protein-targeting properties.

Computationally designed GPCR quaternary structures bias signaling pathway activation

Communication across membranes controls critical cellular processes and is achieved by receptors translating extracellular signals into selective cytoplasmic responses. While receptor tertiary structures can be readily characterized, receptor associations into quaternary structures are challenging to study and their implications in signal transduction remain poorly understood. Here, we report a computational approach for predicting receptor self-associations, and designing receptor oligomers with various quaternary structures and signaling properties. Using this approach, we designed chemokine receptor CXCR4 dimers with reprogrammed binding interactions, conformations, and abilities to activate distinct intracellular signaling proteins. In agreement with our predictions, the designed CXCR4s dimerized through distinct conformations and displayed different quaternary structural changes upon activation. Consistent with the active state models, all engineered CXCR4 oligomers activated the G protein Gi, but only specific dimer structures also recruited β-arrestins. Overall, we demonstrate that quaternary structures represent an important unforeseen mechanism of receptor biased signaling and reveal the existence of a bias switch at the dimer interface of several G protein-coupled receptors including CXCR4, mu-Opioid and type-2 Vasopressin receptors that selectively control the activation of G proteins vs β-arrestin-mediated pathways. The approach should prove useful for predicting and designing receptor associations to uncover and reprogram selective cellular signaling functions.

Filling of a water-free void explains the allosteric regulation of the β1-adrenergic receptor by cholesterol

Recent high-pressure NMR results indicate that the preactive conformation of the β₁-adrenergic receptor (β₁AR) harbours completely empty cavities of ~100 ų volume, which disappear in the active conformation of the receptor. Here we have localized these cavities using X-ray crystallography of xenon-derivatized β₁AR crystals. One of the cavities is in direct contact with the cholesterol-binding pocket. Solution NMR shows that addition of the cholesterol analogue cholesteryl hemisuccinate impedes the formation of the active conformation of detergent-solubilized β₁AR by blocking conserved G protein-coupled receptor microswitches, concomitant with an affinity reduction of both isoprenaline and G protein-mimicking nanobody Nb80 for β₁AR detected by isothermal titration calorimetry. This wedge-like action explains the function of cholesterol as a negative allosteric modulator of β₁AR. A detailed understanding of G protein-coupled receptor regulation by cholesterol by filling of a dry void and the easy scouting for such voids by xenon may provide new routes for the development of allosteric drugs.

Molecular Simulations and Drug Discovery of Adenosine Receptors

G protein-coupled receptors (GPCRs) represent the largest family of human membrane proteins. Four subtypes of adenosine receptors (ARs), the A1AR, A2AAR, A2BAR and A3AR, each with a unique pharmacological profile and distribution within the tissues in the human body, mediate many physiological functions and serve as critical drug targets for treating numerous human diseases including cancer, neuropathic pain, cardiac ischemia, stroke and diabetes. The A1AR and A3AR preferentially couple to the Gi/o proteins, while the A2AAR and A2BAR prefer coupling to the Gs proteins. Adenosine receptors were the first subclass of GPCRs that had experimental structures determined in complex with distinct G proteins. Here, we will review recent studies in molecular simulations and computer-aided drug discovery of the adenosine receptors and also highlight their future research opportunities.

Assembly and Analysis of Cell-Scale Membrane Envelopes

The march toward exascale computing will enable routine molecular simulation of larger and more complex systems, for example, simulation of entire viral particles, on the scale of approximately billions of atoms─a simulation size commensurate with a small bacterial cell. Anticipating the future hardware capabilities that will enable this type of research and paralleling advances in experimental structural biology, efforts are currently underway to develop software tools, procedures, and workflows for constructing cell-scale structures. Herein, we describe our efforts in developing and implementing an efficient and robust workflow for construction of cell-scale membrane envelopes and embedding membrane proteins into them. A new approach for construction of massive membrane structures that are stable during the simulations is built on implementing a subtractive assembly technique coupled with the development of a structure concatenation tool (fastmerge), which eliminates overlapping elements based on volumetric criteria rather than adding successive molecules to the simulation system. Using this approach, we have constructed two "protocells" consisting of MARTINI coarse-grained beads to represent cellular membranes, one the size of a cellular organelle and another the size of a small bacterial cell. The membrane envelopes constructed here remain whole during the molecular dynamics simulations performed and exhibit water flux only through specific proteins, demonstrating the success of our methodology in creating tight cell-like membrane compartments. Extended simulations of these cell-scale structures highlight the propensity for nonspecific interactions between adjacent membrane proteins leading to the formation of protein microclusters on the cell surface, an insight uniquely enabled by the scale of the simulations. We anticipate that the experiences and best practices presented here will form the basis for the next generation of cell-scale models, which will begin to address the addition of soluble proteins, nucleic acids, and small molecules essential to the function of a cell.

Insights into the Role of Membrane Lipids in the Structure, Function and Regulation of Integral Membrane Proteins

Membrane proteins exist within the highly hydrophobic membranes surrounding cells and organelles, playing key roles in cellular function. It is becoming increasingly clear that the membrane does not just act as an appropriate environment for these proteins, but that the lipids that make up these membranes are essential for membrane protein structure and function. Recent technological advances in cryogenic electron microscopy and in advanced mass spectrometry methods, as well as the development of alternative membrane mimetic systems, have allowed experimental study of membrane protein–lipid complexes. These have been complemented by computational approaches, exploiting the ability of Molecular Dynamics simulations to allow exploration of membrane protein conformational changes in membranes with a defined lipid content. These studies have revealed the importance of lipids in stabilising the oligomeric forms of membrane proteins, mediating protein–protein interactions, maintaining a specific conformational state of a membrane protein and activity. Here we review some of the key recent advances in the field of membrane protein–lipid studies, with major emphasis on respiratory complexes, transporters, channels and G-protein coupled receptors.

Quantifying Membrane Curvature Sensing of Peripheral Proteins by Simulated Buckling and Umbrella Sampling

Membrane curvature plays an essential role in the organization and trafficking of membrane associated proteins. Comparison or prediction of the experimentally resolved protein concentrations adopted at different membrane curvatures requires direct quantification of the relative partitioning free energy. Here, we present a highly efficient and simple to implement a free-energy calculation method which is able to directly resolve the relative partitioning free energy of proteins as a direct function of membrane curvature, i.e., a curvature sensing profile, within (coarse-grained) molecular dynamics simulations. We demonstrate its utility by resolving these profiles for two known curvature sensing peptides, namely ALPS and α-synuclein, for a membrane curvature ranging from -1/6.5 to +1/6.5 nm^-1. We illustrate that the difference in relative partitioning (binding) free energy between these two extrema is only about 13 k_BT for both peptides, illustrating that the driving force of curvature sensing is subtle. Furthermore, we illustrate that ALPS and α-synuclein sense curvature via a contrasting mechanism, which is differentially affected by membrane composition. In addition, we demonstrate that the intrinsic spontaneous curvature of both of these peptides lies beyond the range of membrane curvature accessible in micropipette aspiration experiments, being about 1/7 nm ^-1. Our approach offers an efficient and simple to implement in silico tool for exploring and screening the membrane curvature sensing mechanisms of proteins.

Unusual mode of dimerization of retinitis pigmentosa-associated F220C rhodopsin

Mutations in the G protein-coupled receptor (GPCR) rhodopsin are a common cause of autosomal dominant retinitis pigmentosa, a blinding disease. Rhodopsin self-associates in the membrane, and the purified monomeric apo-protein opsin dimerizes in vitro as it transitions from detergent micelles to reconstitute into a lipid bilayer. We previously reported that the retinitis pigmentosa-linked F220C opsin mutant fails to dimerize in vitro, reconstituting as a monomer. Using fluorescence-based assays and molecular dynamics simulations we now report that whereas wild-type and F220C opsin display distinct dimerization propensities in vitro as previously shown, they both dimerize in the plasma membrane of HEK293 cells. Unexpectedly, molecular dynamics simulations show that F220C opsin forms an energetically favored dimer in the membrane when compared with the wild-type protein. The conformation of the F220C dimer is unique, with transmembrane helices 5 and 6 splayed apart, promoting widening of the intracellular vestibule of each protomer and influx of water into the protein interior. FRET experiments with SNAP-tagged wild-type and F220C opsin expressed in HEK293 cells are consistent with this conformational difference. We speculate that the unusual mode of dimerization of F220C opsin in the membrane may have physiological consequences.

Addressing the Excessive Aggregation of Membrane Proteins in the MARTINI Model

The MARTINI model is a widely used coarse-grained force field popular for its capacity to represent a diverse array of complex biomolecules. However, efforts to simulate increasingly realistic models of membranes, involving complex lipid mixtures and multiple proteins, suggest that membrane protein aggregates are overstabilized by the MARTINI v2.2 force field. In this study, we address this shortcoming of the MARTINI model. We determined the free energy of dimerization of four transmembrane protein systems using the nonpolarizable MARTINI model. Comparison with experimental FRET-based estimates of the dimerization free energy was used to quantify the significant overstabilization of each protein homodimer studied. To improve the agreement between simulation and experiment, a single uniform scaling factor, α, was used to enhance the protein-lipid Lennard-Jones interaction. A value of α = 1.04-1.045 was found to provide the best fit to the dimerization free energies for the proteins studied while maintaining the specificity of contacts at the dimer interface. To further validate the modified force field, we performed a multiprotein simulation using both MARTINI v2.2 and the reparameterized MARTINI model. While the original MARTINI model predicts oligomerization of protein into a single aggregate, the reparameterized MARTINI model maintains a dynamic equilibrium between monomers and dimers as predicted by experimental studies. The proposed reparameterization is an alternative to the standard MARTINI model for use in simulations of realistic models of a biological membrane containing diverse lipids and proteins.

Membrane transporter dimerization driven by differential lipid solvation energetics of dissociated and associated states

Over two-thirds of integral membrane proteins of known structure assemble into oligomers. Yet, the forces that drive the association of these proteins remain to be delineated, as the lipid bilayer is a solvent environment that is both structurally and chemically complex. In this study, we reveal how the lipid solvent defines the dimerization equilibrium of the CLC-ec1 Cl^-/H⁺ antiporter. Integrating experimental and computational approaches, we show that monomers associate to avoid a thinned-membrane defect formed by hydrophobic mismatch at their exposed dimerization interfaces. In this defect, lipids are strongly tilted and less densely packed than in the bulk, with a larger degree of entanglement between opposing leaflets and greater water penetration into the bilayer interior. Dimerization restores the membrane to a near-native state and therefore, appears to be driven by the larger free-energy cost of lipid solvation of the dissociated protomers. Supporting this theory, we demonstrate that addition of short-chain lipids strongly shifts the dimerization equilibrium toward the monomeric state, and show that the cause of this effect is that these lipids preferentially solvate the defect. Importantly, we show that this shift requires only minimal quantities of short-chain lipids, with no measurable impact on either the macroscopic physical state of the membrane or the protein's biological function. Based on these observations, we posit that free-energy differentials for local lipid solvation define membrane-protein association equilibria. With this, we argue that preferential lipid solvation is a plausible cellular mechanism for lipid regulation of oligomerization processes, as it can occur at low concentrations and does not require global changes in membrane properties.

A cell’s outer membrane is made of molecules called lipids, which band together to form a flexible thin film, just two molecules thick. This membrane is dotted with proteins that transport materials in to and out of cells. Most of these membrane proteins join with other proteins to form structures known as oligomers. Except, how membrane-bound proteins assemble into oligomers – the physical forces driving these molecules to take shape – remains unclear.

This is partly because the structural, physical and chemical properties of fat-like lipid membranes are radically different to the cell’s watery interior. Consequently, the conditions under which membrane oligomers form are distinct from those surrounding proteins inside cells. Membrane proteins are also more difficult to study and characterize than water-soluble proteins inside the cell, and yet many therapeutic drugs such as antibiotics specifically target membrane proteins. Overall, our understanding of how the unique properties of lipid membranes affect the formation of protein structures embedded within, is lacking and warrants further investigation.

Now, Chadda, Bernhardt et al. focused on one membrane protein, known as CLC, which tends to exist in pairs – or dimers. To understand why these proteins form dimers (a process called dimerization) Chadda, Bernhardt et al. first used computer simulations, and then validated the findings in experimental tests.

These complementary approaches demonstrated that the main reason CLC proteins ‘dimerize’ lies in their interaction with the lipid membrane, and not the attraction of one protein to its partner. When CLC proteins are on their own, they deform the surrounding membrane and create structural defects that put the membrane under strain. But when two CLC proteins join as a dimer, this membrane strain disappears – making dimerization the more stable and energetically favorable option.

Chadda, Bernhardt et al. also showed that with the addition of a few certain lipids, specifically smaller lipids, cell membranes become more tolerant of protein-induced structural changes. This might explain how cells could use various lipids to fine-tune the activity of membrane proteins by controlling how oligomers form. However, the theory needs to be examined further. Altogether, this work has provided fundamental insights into the physical forces shaping membrane-bound proteins, relevant to researchers studying cell biology and pharmacology alike.

Molecular dynamics of the histamine H3 membrane receptor reveals different mechanisms of GPCR signal transduction

In this work, we studied the mechanisms of classical activation and inactivation of signal transduction by the histamine H3 receptor, a 7-helix transmembrane bundle G-Protein Coupled Receptor through long-time-scale atomistic molecular dynamics simulations of the receptor embedded in a hydrated double layer of dipalmitoyl phosphatidyl choline, a zwitterionic polysaturated ordered lipid. Three systems were prepared: the apo receptor, representing the constitutively active receptor; and two holo-receptors-the receptor coupled to the antagonist/inverse agonist ciproxifan, representing the inactive state of the receptor, and the receptor coupled to the endogenous agonist histamine and representing the active state of the receptor. An extensive analysis of the simulation showed that the three states of H3R present significant structural and dynamical differences as well as a complex behavior given that the measured properties interact in multiple and interdependent ways. In addition, the simulations described an unexpected escape of histamine from the orthosteric binding site, in agreement with the experimental modest affinities and rapid off-rates of agonists.

Symmetry of membrane protein polyhedra with heterogeneous protein size

In experiments on membrane protein polyhedral nanoparticles (MPPNs) [Basta et al., Proc. Natl. Acad. Sci. USA 111, 670 (2014)PNASA60027-842410.1073/pnas.1321936111], it has been observed that membrane proteins and lipids can self-assemble into closed lipid bilayer vesicles with a polyhedral arrangement of membrane proteins. In particular, MPPNs formed from the mechanosensitive channel of small conductance (MscS) were found to have the symmetry of the snub cube-a chiral, Archimedean solid-with one MscS protein located at each one of the 24 vertices of the snub cube. It is currently unknown whether MPPNs with heterogeneous protein composition maintain a high degree of symmetry. Inspired by previous work on viral capsid symmetry, we employ here computational modeling to study the symmetry of MPPNs with heterogeneous protein size. We focus on MPPNs formed from MscS proteins, which can exist in closed or open conformational states with distinct sizes. We find that, as an increasing number of closed-state MscS proteins transitions to the open conformational state of MscS, the minimum-energy MscS arrangement in MPPNs follows a strikingly regular pattern, with the dominant MPPN symmetry always being provided by the snub cube. Our results suggest that MPPNs with heterogeneous protein size can be highly symmetric, with a well-defined polyhedral ordering of membrane proteins of different sizes.

Membrane lipids are both the substrates and a mechanistically responsive environment of TMEM16 scramblase proteins

Recent discoveries about functional mechanisms of proteins in the TMEM16 family of phospholipid scramblases have illuminated the dual role of the membrane as both the substrate and a mechanistically responsive environment in the wide range of physiological processes and genetic disorders in which they are implicated. This is highlighted in the review of recent findings from our collaborative investigations of molecular mechanisms of TMEM16 scramblases that emerged from iterative functional, structural, and computational experimentation. In the context of this review, we present new MD simulations and trajectory analyses motivated by the fact that new structural information about the TMEM16 scramblases is emerging from cryo-EM determinations in lipid nanodiscs. Because the functional environment of these proteins in in vivo and in in vitro is closer to flat membranes, we studied comparatively the responses of the membrane to the TMEM16 proteins in flat membranes and nanodiscs. We find that bilayer shapes in the nanodiscs are very different from those observed in the flat membrane systems, but the function-related slanting of the membrane observed at the nhTMEM16 boundary with the protein is similar in the nanodiscs and in the flat bilayers. This changes, however, in the bilayer composed of longer-tail lipids, which is thicker near the phospholipid translocation pathway, which may reflect an enhanced tendency of the long tails to penetrate the pathway and create, as shown previously, a nonconductive environment. These findings support the correspondence between the mechanistic involvement of the lipid environment in the flat membranes, and the nanodiscs. © 2019 Wiley Periodicals, Inc.

Structural instability and divergence from conserved residues underlie intracellular retention of mammalian odorant receptors

Mammalian odorant receptors are a diverse and rapidly evolving set of G protein-coupled receptors expressed in olfactory cilia membranes. Most odorant receptors show little to no cell surface expression in nonolfactory cells due to endoplasmic reticulum retention, which has slowed down biochemical studies. Here we provide evidence that structural instability and divergence from conserved residues of individual odorant receptors underlie intracellular retention using a combination of large-scale screening of odorant receptors cell surface expression in heterologous cells, point mutations, structural modeling, and machine learning techniques. We demonstrate the importance of conserved residues by synthesizing consensus odorant receptors that show high levels of cell surface expression similar to conventional G protein-coupled receptors. Furthermore, we associate in silico structural instability with poor cell surface expression using molecular dynamics simulations. We propose an enhanced evolutionary capacitance of olfactory sensory neurons that enable the functional expression of odorant receptors with cryptic mutations.

Prediction and targeting of GPCR oligomer interfaces

GPCR oligomerization has emerged as a hot topic in the GPCR field in the last years. Receptors that are part of these oligomers can influence each other's function, although it is not yet entirely understood how these interactions work. The existence of such a highly complex network of interactions between GPCRs generates the possibility of alternative targets for new therapeutic approaches. However, challenges still exist in the characterization of these complexes, especially at the interface level. Different experimental approaches, such as FRET or BRET, are usually combined to study GPCR oligomer interactions. Computational methods have been applied as a useful tool for retrieving information from GPCR sequences and the few X-ray-resolved oligomeric structures that are accessible, as well as for predicting new and trustworthy GPCR oligomeric interfaces. Machine-learning (ML) approaches have recently helped with some hindrances of other methods. By joining and evaluating multiple structure-, sequence- and co-evolution-based features on the same algorithm, it is possible to dilute the issues of particular structures and residues that arise from the experimental methodology into all-encompassing algorithms capable of accurately predict GPCR-GPCR interfaces. All these methods used as a single or a combined approach provide useful information about GPCR oligomerization and its role in GPCR function and dynamics. Altogether, we present experimental, computational and machine-learning methods used to study oligomers interfaces, as well as strategies that have been used to target these dynamic complexes.

Rearrangement of the transmembrane domain interfaces associated with the activation of a GPCR hetero-oligomer

G protein-coupled receptors (GPCRs) can integrate extracellular signals via allosteric interactions within dimers and higher-order oligomers. However, the structural bases of these interactions remain unclear. Here, we use the GABA_B receptor heterodimer as a model as it forms large complexes in the brain. It is subjected to genetic mutations mainly affecting transmembrane 6 (TM6) and involved in human diseases. By cross-linking, we identify the transmembrane interfaces involved in GABA_B1-GABA_B2, as well as GABA_B1-GABA_B1 interactions. Our data are consistent with an oligomer made of a row of GABA_B1. We bring evidence that agonist activation induces a concerted rearrangement of the various interfaces. While the GB1-GB2 interface is proposed to involve TM5 in the inactive state, cross-linking of TM6s lead to constitutive activity. These data bring insight for our understanding of the allosteric interaction between GPCRs within oligomers.

G protein-coupled receptors (GPCRs), such as GABA_B, can integrate extracellular signals via allosteric interactions within dimers and oligomers. Here authors use crosslinking and identify two transmembrane interfaces in GABA_B which undergo a concerted rearrangement upon agonist activation.

Supramolecular organization of membrane proteins with anisotropic hydrophobic thickness

Experiments have revealed that membrane proteins often self-assemble into locally ordered clusters. Such membrane protein lattices can play key roles in the functional organization of cell membranes. Membrane protein organization can be driven, at least in part, by bilayer-mediated elastic interactions between membrane proteins. For membrane proteins with anisotropic hydrophobic thickness, bilayer-mediated protein interactions are inherently directional. Here we establish general relations between anisotropy in membrane protein hydrophobic thickness and supramolecular membrane protein organization. We show that protein symmetry is distinctively reflected in the energy landscape of bilayer-mediated protein interactions, favoring characteristic lattice architectures of membrane protein clusters. We find that, in the presence of thermal fluctuations, anisotropy in protein hydrophobic thickness can induce membrane proteins to form mesh-like structures dividing the membrane into compartments. Our results help to elucidate the physical principles and mechanisms underlying the functional organization of cell membranes.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Emerging Diversity in Lipid-Protein Interactions

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid-protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid-protein interactions, and to link lipid-protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid-protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid-protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid-protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid-protein interactions.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

State-dependent Lipid Interactions with the A2a Receptor Revealed by MD Simulations Using In Vivo-Mimetic Membranes

Membranes are known to have modulatory effects on G protein-coupled receptors (GPCRs) via specific lipid interactions. However, the mechanisms of such modulations in physiological conditions and how they influence GPCR functions remain unclear. Here we report coarse-grained molecular dynamics simulations on the Adenosine A2a receptor in different conformational states embedded in an in vivo-mimetic membrane model. Nine lipid interaction sites were revealed. The strength of lipid interactions with these sites showed a degree of dependence on the conformational states of the receptor, suggesting that these lipids may regulate the conformational dynamics of the receptor. In particular, we revealed a dual role of PIP₂ on A2aR activation that involves both stabilization of the characteristic outward tilt of TM6 and enhancement of A2aR-mini-Gs association. Our results demonstrated that the bound lipids allosterically regulate the functional properties of GPCRs. These protein-lipid interactions provide a springboard for design of allosteric modulators of GPCRs.

Use multiscale simulation to explore the effects of the homodimerizations between different conformation states on the activation and allosteric pathway for the μ-opioid receptor

Recently, oligomers of G-protein coupled receptors (GPCRs) have been an important topic in the GPCR fields. However, knowledge about their structures and activation mechanisms is very limited due to the absence of crystal structures reported. In this work, we used multiscale simulations to study the effects of homodimerization between different conformation states on their activation, dynamic behaviors, and allosteric communication pathways for μ-OR. The results indicated that the dimerization of one inactive monomer with either one inactive monomer or one active one could enhance its constitutive activation. However, the conformation state of the other protomer (e.g., active or inactive) can influence the activated extent. The dimerization between the two inactive protomers leads to a negative cooperativity for their activation, which should contribute to the asymmetric activation of GPCR dimers observed in some experiments. On the other hand, for the active monomer, its dimerization with one inactive receptor could alleviate its deactivation, whereby negative and positive cooperativities can be observed between the two subunits of the dimer, depending on the different regions. Observations from protein structure network (PSN) analysis indicated that the dimerization of one inactive monomer with one active one would cause a significant drop in the number of main pathways from the ligand binding pocket to the G-protein coupled region for the inactive protomer, while the impact is minor for the active protomer. But, for the active monomer or the inactive one, its dimerization with one inactive monomer would significantly change the types of residues participating in the pathway with the highest frequency.

Atom-based 3D-QSAR, induced fit docking, and molecular dynamics simulations study of thieno[2,3-b]pyridines negative allosteric modulators of mGluR5

Atom-based three dimensional-quantitative structure-activity relationship (3D-QSAR) model was developed on the basis of 5-point pharmacophore hypothesis (AARRR) with two hydrogen bond acceptors (A) and three aromatic rings for the derivatives of thieno[2,3-b]pyridine, which modulates the activity to inhibit the mGluR5 receptor. Generation of a highly predictive 3D-QSAR model was performed using the alignment of predicted pharmacophore hypothesis for the training set (R² = 0.84, SD = 0.26, F = 45.8, N = 29) and test set (Q² = 0.74, RMSE = 0.235, Pearson-R = 0.94, N = 9). The best pharmacophore hypothesis AARRR was selected, and developed three dimensional-quantitative structure activity relationship (3D-QSAR) model also supported the outcome of this study by means of favorable and unfavorable electron withdrawing group and hydrophobic regions of most active compound 42d and least active compound 18b. Following, induced fit docking and binding free energy calculations reveals the reliable binding orientation of the compounds. Finally, molecular dynamics simulations for 100 ns were performed to depict the protein-ligand stability. We anticipate that the resulted outcome could be supportive to discover potent negative allosteric modulators for metabotropic glutamate receptor 5 (mGluR5).

Dynamic tuneable G protein-coupled receptor monomer-dimer populations

G protein-coupled receptors (GPCRs) are the largest class of membrane receptors, playing a key role in the regulation of processes as varied as neurotransmission and immune response. Evidence for GPCR oligomerisation has been accumulating that challenges the idea that GPCRs function solely as monomeric receptors; however, GPCR oligomerisation remains controversial primarily due to the difficulties in comparing evidence from very different types of structural and dynamic data. Using a combination of single-molecule and ensemble FRET, double electron–electron resonance spectroscopy, and simulations, we show that dimerisation of the GPCR neurotensin receptor 1 is regulated by receptor density and is dynamically tuneable over the physiological range. We propose a “rolling dimer” interface model in which multiple dimer conformations co-exist and interconvert. These findings unite previous seemingly conflicting observations, provide a compelling mechanism for regulating receptor signalling, and act as a guide for future physiological studies.

Evidence suggests oligomerisation of G protein-coupled receptors in membranes, but this is controversial. Here, authors use single-molecule and ensemble FRET, and spectroscopy to show that the neurotensin receptor 1 forms multiple dimer conformations that interconvert - “rolling” interfaces.

Closely related, yet unique: Distinct homo- and heterodimerization patterns of G protein coupled chemokine receptors and their fine-tuning by cholesterol

Chemokine receptors, a subclass of G protein coupled receptors (GPCRs), play essential roles in the human immune system, they are involved in cancer metastasis as well as in HIV-infection. A plethora of studies show that homo- and heterodimers or even higher order oligomers of the chemokine receptors CXCR4, CCR5, and CCR2 modulate receptor function. In addition, membrane cholesterol affects chemokine receptor activity. However, structural information about homo- and heterodimers formed by chemokine receptors and their interplay with cholesterol is limited. Here, we report homo- and heterodimer configurations of the chemokine receptors CXCR4, CCR5, and CCR2 at atomistic detail, as obtained from thousands of molecular dynamics simulations. The observed homodimerization patterns were similar for the closely related CC chemokine receptors, yet they differed significantly between the CC receptors and CXCR4. Despite their high sequence identity, cholesterol modulated the CC homodimer interfaces in a subtype-specific manner. Chemokine receptor heterodimers display distinct dimerization patterns for CXCR4/CCR5 and CXCR4/CCR2. Furthermore, associations between CXCR4 and CCR5 reveal an increased cholesterol-sensitivity as compared to CXCR4/CCR2 heterodimerization patterns. This work provides a first comprehensive structural overview over the complex interaction network between chemokine receptors and indicates how heterodimerization and the interaction with the membrane environment diversifies the function of closely related GPCRs.

Membrane-Mediated Cooperativity of Proteins

Besides direct protein-protein interactions, indirect interactions mediated by membranes play an important role for the assembly and cooperative function of proteins in membrane shaping and adhesion. The intricate shapes of biological membranes are generated by proteins that locally induce membrane curvature. Indirect curvature-mediated interactions between these proteins arise because the proteins jointly affect the bending energy of the membranes. These curvature-mediated interactions are attractive for crescent-shaped proteins and are a driving force in the assembly of the proteins during membrane tubulation. Membrane adhesion results from the binding of receptor and ligand proteins that are anchored in the apposing membranes. The binding of these proteins strongly depends on nanoscale shape fluctuations of the membranes, leading to a fluctuation-mediated binding cooperativity. A length mismatch between receptor-ligand complexes in membrane adhesion zones causes repulsive curvature-mediated interactions that are a driving force for the length-based segregation of proteins during membrane adhesion.

G protein-coupled receptor-receptor interactions give integrative dynamics to intercellular communication

The proposal of receptor-receptor interactions (RRIs) in the early 1980s broadened the view on the role of G protein-coupled receptors (GPCR) in the dynamics of the intercellular communication. RRIs, indeed, allow GPCR to operate not only as monomers but also as receptor complexes, in which the integration of the incoming signals depends on the number, spatial arrangement, and order of activation of the protomers forming the complex. The main biochemical mechanisms controlling the functional interplay of GPCR in the receptor complexes are direct allosteric interactions between protomer domains. The formation of these macromolecular assemblies has several physiologic implications in terms of the modulation of the signaling pathways and interaction with other membrane proteins. It also impacts on the emerging field of connectomics, as it contributes to set and tune the synaptic strength. Furthermore, recent evidence suggests that the transfer of GPCR and GPCR complexes between cells via the exosome pathway could enable the target cells to recognize/decode transmitters and/or modulators for which they did not express the pertinent receptors. Thus, this process may also open the possibility of a new type of redeployment of neural circuits. The fundamental aspects of GPCR complex formation and function are the focus of the present review article.

Exploring oligomeric state of the serotonin1A receptor utilizing photobleaching image correlation spectroscopy: implications for receptor function

Photobleaching image correlation spectroscopy (pbICS) reveals that membrane cholesterol modulates the oligomeric state of the serotonin_1A receptor.

An engineered opsin monomer scrambles phospholipids

The G protein-coupled receptor opsin is a phospholipid scramblase that facilitates rapid transbilayer phospholipid exchange in liposomes. The mechanism by which opsin scrambles lipids is unknown. It has been proposed that lipid translocation may occur at protein-protein interfaces of opsin dimers. To test this possibility, we rationally engineered QUAD opsin by tryptophan substitution of four lipid-facing residues in transmembrane helix 4 (TM4) that is known to be important for dimerization. Atomistic molecular dynamics simulations of wild type and QUAD opsins combined with continuum modeling revealed that the tryptophan substitutions lower the energetically unfavorable residual hydrophobic mismatch between TM4 and the membrane, reducing the drive of QUAD opsin to dimerize. We purified thermostable wild type and QUAD opsins, with or without a SNAP tag for fluorescence labeling. Single molecule fluorescence measurements of purified SNAP-tagged constructs revealed that both proteins are monomers. Fluorescence-based activity assays indicated that QUAD opsin is a fully functional scramblase. However, unlike wild type opsin which dimerizes en route to insertion into phospholipid vesicles, QUAD opsin reconstitutes as a monomer. We conclude that an engineered opsin monomer can scramble phospholipids, and that the lipid-exposed face of TM4 is unlikely to contribute to transbilayer phospholipid exchange.

Excessive aggregation of membrane proteins in the Martini model

The coarse-grained Martini model is employed extensively to study membrane protein oligomerization. While this approach is exceptionally promising given its computational efficiency, it is alarming that a significant fraction of these studies demonstrate unrealistic protein clusters, whose formation is essentially an irreversible process. This suggests that the protein-protein interactions are exaggerated in the Martini model. If this held true, then it would limit the applicability of Martini to study multi-protein complexes, as the rapidly clustering proteins would not be able to properly sample the correct dimerization conformations. In this work we first demonstrate the excessive protein aggregation by comparing the dimerization free energies of helical transmembrane peptides obtained with the Martini model to those determined from FRET experiments. Second, we show that the predictions provided by the Martini model for the structures of transmembrane domain dimers are in poor agreement with the corresponding structures resolved using NMR. Next, we demonstrate that the first issue can be overcome by slightly scaling down the Martini protein-protein interactions in a manner, which does not interfere with the other Martini interaction parameters. By preventing excessive, irreversible, and non-selective aggregation of membrane proteins, this approach renders the consideration of lateral dynamics and protein-lipid interactions in crowded membranes by the Martini model more realistic. However, this adjusted model does not lead to an improvement in the predicted dimer structures. This implicates that the poor agreement between the Martini model and NMR structures cannot be cured by simply uniformly reducing the interactions between all protein beads. Instead, a careful amino-acid specific adjustment of the protein-protein interactions is likely required.

Membrane lipid order of sub-synaptic T cell vesicles correlates with their dynamics and function

During an immune response, T cells survey antigen presenting cells for antigenic peptides via the formation of an interface known as an immunological synapse. Among the complex and dynamic biophysical phenomena occurring at this interface is the trafficking of sub-synaptic vesicles carrying a variety of proximal signalling molecules. Here, we show that rather than being a homogeneous population, these vesicles display a diversity of membrane lipid order profiles, as measured using the environmentally sensitive dye di-4-ANEPPDHQ and multi-spectral TIRF microscopy. Using live-cell imaging, vesicle tracking and a variety of small molecule drugs to manipulate components of the actin and tubulin cytoskeleton, we show that the membrane lipid order of these vesicles correlate with their dynamics. Furthermore, we show that the key proximal signalling molecule Linker for Activation of T cells (LAT) is enriched in specific vesicle populations as defined by their higher membrane order. These results imply that vesicle lipid order may represent a novel regulatory mechanism for the sorting and trafficking of signalling molecules at the immunological synapse, and, potentially, other cellular structures.

Membrane cholesterol effect on the 5-HT2A receptor: Insights into the lipid-induced modulation of an antipsychotic drug target

The serotonin 5-hydroxytryptamine 2A (5-HT_2A ) receptor is a G-protein-coupled receptor (GPCR) relevant for the treatment of CNS disorders. In this regard, neuronal membrane composition in the brain plays a crucial role in the modulation of the receptor functioning. Since cholesterol is an essential component of neuronal membranes, we have studied its effect on the 5-HT_2A receptor dynamics through all-atom MD simulations. We find that the presence of cholesterol in the membrane increases receptor conformational variability in most receptor segments. Importantly, detailed structural analysis indicates that conformational variability goes along with the destabilization of hydrogen bonding networks not only within the receptor but also between receptor and lipids. In addition to increased conformational variability, we also find receptor segments with reduced variability. Our analysis suggests that this increased stabilization is the result of stabilizing effects of tightly bound cholesterol molecules to the receptor surface. Our finding contributes to a better understanding of membrane-induced alterations of receptor dynamics and points to cholesterol-induced stabilizing and destabilizing effects on the conformational variability of GPCRs.

An Ensemble-Based Protocol for the Computational Prediction of Helix-Helix Interactions in G Protein-Coupled Receptors using Coarse-Grained Molecular Dynamics

The accurate identification of the specific points of interaction between G protein-coupled receptor (GPCR) oligomers is essential for the design of receptor ligands targeting oligomeric receptor targets. A coarse-grained molecular dynamics computer simulation approach would provide a compelling means of identifying these specific protein–protein interactions and could be applied both for known oligomers of interest and as a high-throughput screen to identify novel oligomeric targets. However, to be effective, this in silico modeling must provide accurate, precise, and reproducible information. This has been achieved recently in numerous biological systems using an ensemble-based all-atom molecular dynamics approach. In this study, we describe an equivalent methodology for ensemble-based coarse-grained simulations. We report the performance of this method when applied to four different GPCRs known to oligomerize using error analysis to determine the ensemble size and individual replica simulation time required. Our measurements of distance between residues shown to be involved in oligomerization of the fifth transmembrane domain from the adenosine A_2A receptor are in very good agreement with the existing biophysical data and provide information about the nature of the contact interface that cannot be determined experimentally. Calculations of distance between rhodopsin, CXCR4, and β₁AR transmembrane domains reported to form contact points in homodimers correlate well with the corresponding measurements obtained from experimental structural data, providing an ability to predict contact interfaces computationally. Interestingly, error analysis enables identification of noninteracting regions. Our results confirm that GPCR interactions can be reliably predicted using this novel methodology.

Methods used to study the oligomeric structure of G-protein-coupled receptors

G-protein-coupled receptors (GPCRs), which constitute the largest family of cell surface receptors, were originally thought to function as monomers, but are now recognized as being able to act in a wide range of oligomeric states and indeed, it is known that the oligomerization state of a GPCR can modulate its pharmacology and function. A number of experimental techniques have been devised to study GPCR oligomerization including those based upon traditional biochemistry such as blue-native PAGE (BN-PAGE), co-immunoprecipitation (Co-IP) and protein-fragment complementation assays (PCAs), those based upon resonance energy transfer, FRET, time-resolved FRET (TR-FRET), FRET spectrometry and bioluminescence resonance energy transfer (BRET). Those based upon microscopy such as FRAP, total internal reflection fluorescence microscopy (TIRFM), spatial intensity distribution analysis (SpIDA) and various single molecule imaging techniques. Finally with the solution of a growing number of crystal structures, X-ray crystallography must be acknowledged as an important source of discovery in this field. A different, but in many ways complementary approach to the use of more traditional experimental techniques, are those involving computational methods that possess obvious merit in the study of the dynamics of oligomer formation and function. Here, we summarize the latest developments that have been made in the methods used to study GPCR oligomerization and give an overview of their application.

Impact of Lipid Composition and Receptor Conformation on the Spatio-temporal Organization of μ-Opioid Receptors in a Multi-component Plasma Membrane Model

The lipid composition of cell membranes has increasingly been recognized as playing an important role in the function of various membrane proteins, including G Protein-Coupled Receptors (GPCRs). For instance, experimental and computational evidence has pointed to lipids influencing receptor oligomerization directly, by physically interacting with the receptor, and/or indirectly, by altering the bulk properties of the membrane. While the exact role of oligomerization in the function of class A GPCRs such as the μ-opioid receptor (MOR) is still unclear, insight as to how these receptors oligomerize and the relevance of the lipid environment to this phenomenon is crucial to our understanding of receptor function. To examine the effect of lipids and different MOR conformations on receptor oligomerization we carried out extensive coarse-grained molecular dynamics simulations of crystal structures of inactive and/or activated MOR embedded in an idealized mammalian plasma membrane composed of 63 lipid types asymmetrically distributed across the two leaflets. The results of these simulations point, for the first time, to specific direct and indirect effects of the lipids, as well as the receptor conformation, on the spatio-temporal organization of MOR in the plasma membrane. While sphingomyelin-rich, high-order lipid regions near certain transmembrane (TM) helices of MOR induce an effective long-range attractive force on individual protomers, both long-range lipid order and interface formation are found to be conformation dependent, with a larger number of different interfaces formed by inactive MOR compared to active MOR.

The μ-opioid receptor (MOR) is an important pharmaceutical target in the treatment of pain. In order to develop novel pain therapies, devoid of the serious side-effects of present opioid analgesics, we need to understand the fundamentals of how MOR works on the molecular level. While some studies suggest that oligomers of MOR could play a role in signaling, how MOR forms dimers, which interfaces form, and the exact role of oligomers in MOR function remain unclear. While research has shown that the membrane environment can affect membrane protein function, most previous computational work to study oligomerization has been performed in a very simple membrane. Here, we use molecular dynamics simulations of MOR in a heterogeneous plasma membrane model (comprising 63 lipid types) to investigate how the presence of the protein modulates its lipid environment, affecting species distribution and sculpting characteristic order and thickness profiles around the receptors. Such modulations, in turn, induce long-range interactions between the proteins and favor the formation of specific dimeric conformations.

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

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

Membrane-Mediated Oligomerization of G Protein Coupled Receptors and Its Implications for GPCR Function

The dimerization or even oligomerization of G protein coupled receptors (GPCRs) causes ongoing, controversial debates about its functional role and the coupled biophysical, biochemical or biomedical implications. A continously growing number of studies hints to a relation between oligomerization and function of GPCRs and strengthens the assumption that receptor assembly plays a key role in the regulation of protein function. Additionally, progress in the structural analysis of GPCR-G protein and GPCR-ligand interactions allows to distinguish between actively functional and non-signaling complexes. Recent findings further suggest that the surrounding membrane, i.e., its lipid composition may modulate the preferred dimerization interface and as a result the abundance of distinct dimeric conformations. In this review, the association of GPCRs and the role of the membrane in oligomerization will be discussed. An overview of the different reported oligomeric interfaces is provided and their capability for signaling discussed. The currently available data is summarized with regard to the formation of GPCR oligomers, their structures and dependency on the membrane microenvironment as well as the coupling of oligomerization to receptor function.

Insertion into lipid bilayer of truncated pHLIP®

The investigation of pH-dependent membrane-associated folding has both fundamental interest and practical applications for targeting of acidic tumors and specific delivery of therapeutic molecules across membrane of cancer cells. We and others investigated molecular mechanism and medical uses of class of water soluble membrane peptides, pH (Low) Insertion Peptides (pHLIP^® peptides). Here we employed optical spectroscopy methods to study interactions of the truncated pHLIP^® peptide (Short pHLIP^®) with lipid bilayer of membrane. Tryptophan fluorescence, CD and OCD data indicate on pH-triggered formation of transmembrane helical structure. Dual quenching and FRET assays demonstrated that Short pHLIP^® peptide spans lipid bilayer of membrane similar to Long pHLIP^® peptides. Truncated pHLIP^® peptides with multiple charged and protonatable residues in their sequences potentially can make these peptides to be less hydrophobic compared to Long pHLIP^® peptides, and might have utility in tumor imaging, and potentially, in pH-regulated cytoplasmic delivery of moderately hydrophobic drugs.

•

Protonation of single Asp triggers pH-dependent insertion into membrane.

•

Truncated pHLIP^® peptide adopts TM helical orientation at low pH.

•

Truncated pHLIP^® peptide spans lipid bilayer similar to full-length long peptide.

•

Potential applications in cytoplasmic delivery of small hydrophobic molecules.

Molecular dynamics simulations of membrane proteins and their interactions: from nanoscale to mesoscale

Molecular dynamics simulations provide a computational tool to probe membrane proteins and systems at length scales ranging from nanometers to close to a micrometer, and on microsecond timescales. All atom and coarse-grained simulations may be used to explore in detail the interactions of membrane proteins and specific lipids, yielding predictions of lipid binding sites in good agreement with available structural data. Building on the success of protein-lipid interaction simulations, larger scale simulations reveal crowding and clustering of proteins, resulting in slow and anomalous diffusional dynamics, within realistic models of cell membranes. Current methods allow near atomic resolution simulations of small membrane organelles, and of enveloped viruses to be performed, revealing key aspects of their structure and functionally important dynamics.

Cholesterol-dependent Conformational Plasticity in GPCR Dimers

The organization and function of the serotonin_1A receptor, an important member of the GPCR family, have been shown to be cholesterol-dependent, although the molecular mechanism is not clear. We performed a comprehensive structural and dynamic analysis of dimerization of the serotonin_1A receptor by coarse-grain molecular dynamics simulations totaling 3.6 ms to explore the molecular details of its cholesterol-dependent association. A major finding is that the plasticity and flexibility of the receptor dimers increase with increased cholesterol concentration. In particular, a dimer interface formed by transmembrane helices I-I was found to be sensitive to cholesterol. The modulation of dimer interface appears to arise from a combination of direct cholesterol occupancy and indirect membrane effects. Interestingly, the presence of cholesterol at the dimer interface is correlated with increased dimer plasticity and flexibility. These results represent an important step in characterizing the molecular interactions in GPCR organization with potential relevance to therapeutic interventions.