Thomas Huber

Many integral membrane proteins assemble to form oligomeric structures in biological membranes. In particular, seven-transmembrane helical G protein-coupled receptors (GPCRs) appear to self-assemble constitutively in membranes, but the mechanism and physiological role of this assembly are unknown. We developed and employed coarse-grain molecular dynamics (CGMD) models to investigate the molecular basis of how the physicochemical properties of the phospholipid bilayer membrane affect self-assembly of visual rhodopsin, a prototypical GPCR. The CGMD method is a mesoscopic simulation technique in which groups of atoms are mapped to particles on the basis of a four-to-one rule. This systematic reduction of the degrees of freedom allows for computationally efficient calculation of the structure and dynamics of molecular assemblies for larger time and length scales than accessible to atomistic models, providing here an unprecedented view of spontaneous protein assembly in biomembranes. Systems with up to 16 rhodopsin molecules at a protein-to-lipid ratio of 1:100 were simulated for time scales of up to 8 micros. The results obtained for four different phospholipid environments showed that localized adaptation of the membrane bilayer to the presence of receptors is reproducibly most pronounced near transmembrane helices 2, 4, and 7. This local membrane deformation appears to be a key factor defining the rate, extent, and orientational preference of protein-protein association. The implications of our findings are discussed within a framework of a generalized mechanism of membrane protein self-assembly.

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