Homomeric protein complexes: evolution and assembly

Heligeom: A web resource to generate, analyze, and visualize filament architectures based on pairwise association geometries of biological macromolecules

At the subcellular level, macromolecules self-assemble to form molecular machinery in which the assembly modes play critical roles: the structural integrity of cell walls that allows mechanical growth, the maintenance and repair of the genetic material, membrane flow control, protein chaperoning, and ATP production, to cite just a few examples. As molecular modeling expands its scope to such systems, structural biologists are faced with the difficulty of understanding the structure and dynamics of these supramolecular assemblies. We present Heligeom, a webserver that offers a simple and efficient means for analyzing and constructing oligomeric assemblies based on user-provided structures of two interacting units. The input 3D coordinates may result from structure determination, simulations, docking trials, or deep-learning tools such as AlphaFold. For a given interface, Heligeom outputs the mathematical helical parameters of the corresponding oligomeric form, including axis, pitch, handedness, number of monomers per turn, etc. The server also allows leveraging these parameters to construct oligomers of specified size, visualizing them interactively, and downloading them as PDB files. For subunits (protomers) having multiple interaction geometries, the different interfaces and their global assembly geometry can be visualized and compared. Heligeom thus allows explicitly linking protomer-protomer interfaces to the oligomeric architecture, illuminating possible sources of plasticity in protein filaments such as mutations or thermal, mechanical, or chemical perturbations. Heligeom thus constitutes an intuitive tool to accompany integrative modeling of oligomeric filamentous assemblies. Examples of its application at different structural levels are presented. Heligeom webserver can be accessed at https://heligeom.galaxy.ibpc.fr.

One for All, All for One: The Peculiar Dynamics of TNF-Receptor-Associated Factor (TRAF2) Subunits

TNF Receptor-Associated Factor 2 (TRAF2) is a homo-trimer belonging to the TNF-receptor-associated factor family (TRAFs). The TRAF2 oligomeric state is crucial for receptor binding, the interaction with other proteins (involved in the TNFR signaling), and the interaction with biological membranes. In this study, we present a computational analysis of the Molecular Dynamics of TRAF2-C (a truncated and soluble TRAF2 form) to identify patterns in the interactions between the three chains. We have performed a canonical analysis of the motion applied to molecular dynamics starting from the available crystal structure to identify correlated motions in TRAF2 dynamics. We have computed the displacement matrix, providing a frame-by-frame displacement for each residue in the dynamic. We provide the results in terms of the correlation matrix, which represents a detailed map of the correlated motions of residues. Eventually, we computed the so-called dynamical clusters, based on the Principal Component Analysis (PCA) of the motion (displacement) and the k means application on the first two principal components space. The results clearly indicate that, most of the time, two chains move in a strongly correlated motion, while the third chain follows a freer motion. A detailed analysis of the correlation matrix also shows that a few specific interface residues characterize the interaction of the more independent subunit with the other two. These findings suggest that the equilibrium between the trimer and the dissociated species (dimers and monomers) might be finely tuned by controlling a few critical residues in the protein quaternary structure, probably facilitating the regulation of oligomerization and dissociation in vivo.

Targeting protein self-association in drug design

Protein self-association is a universal phenomenon essential for stability and molecular recognition. Disrupting constitutive homomers constitutes an original and emerging strategy in drug design. Inhibition of homomeric proteins can be achieved through direct complex disruption, subunit intercalation, or by promoting inactive oligomeric states. Targeting self-interaction grants several advantages over active site inhibition because of the stimulation of protein degradation, the enhancement of selectivity, substoichiometric inhibition, and by-pass of compensatory mechanisms. This new landscape in protein inhibition is driven by the development of biophysical and biochemical tools suited for the study of homomeric proteins, such as differential scanning fluorimetry (DSF), native mass spectrometry (MS), Förster resonance energy transfer (FRET) spectroscopy, 2D nuclear magnetic resonance (NMR), and X-ray crystallography. In this review, we discuss the different aspects of this new paradigm in drug design.

Hexamerization and thermostability emerged very early during geranylgeranylglyceryl phosphate synthase evolution

A large number of archaea live in hyperthermophilic environments. In consequence, their proteins need to adopt to these harsh conditions, including the enzymes that catalyze the synthesis of their membrane ether lipids. The enzyme that catalyzes the formation of the first ether bond in these lipids, geranylgeranylglyceryl phosphate synthase (GGGPS), exists as a hexamer in many hyperthermophilic archaea, and a recent study suggested that hexamerization serves for a fine‐tuning of the flexibility – stability trade‐off under hyperthermophilic conditions. We have recently reconstructed the sequences of ancestral group II GGGPS enzymes and now present a detailed biochemical characterization of nine of these predecessors, which allowed us to trace back the evolution of hexameric GGGPS and to draw conclusions about the properties of extant GGGPS branches that were not accessible to experiments up to now. Almost all ancestral GGGPS proteins formed hexamers, which demonstrates that hexamerization is even more widespread among the GGGPS family than previously assumed. Furthermore, all experimentally studied ancestral proteins showed high thermostability. Our results indicate that the hexameric oligomerization state and thermostability were present very early during the evolution of group II GGGPS, while the fine tuning of the flexibility – stability trade‐off developed very late, independent of the emergence of hexamerization.

Quaternary variations in the structural assembly of N-acetylglucosamine-6-phosphate deacetylase from Pasteurella multocida

N-acetylglucosamine 6-phosphate deacetylase (NagA) catalyzes the conversion of N-acetylglucosamine-6-phosphate to glucosamine-6-phosphate in amino sugar catabolism. This conversion is an essential step in the catabolism of sialic acid in several pathogenic bacteria, including Pasteurella multocida, and thus NagA is identified as a potential drug target. Here, we report the unique structural features of NagA from P. multocida (PmNagA) resolved to 1.95 Å. PmNagA displays an altered quaternary architecture with unique interface interactions compared to its close homolog, the Escherichia coli NagA (EcNagA). We confirmed that the altered quaternary structure is not a crystallographic artifact using single particle electron cryo-microscopy. Analysis of the determined crystal structure reveals a set of hot-spot residues involved in novel interactions at the dimer-dimer interface. PmNagA binds to one Zn2+ ion in the active site and demonstrates kinetic parameters comparable to other bacterial homologs. Kinetic studies reveal that at high substrate concentrations (~10-fold the KM ), the tetrameric PmNagA displays hysteresis similar to its distant neighbor, the dimeric Staphylococcus aureus NagA (SaNagA). Our findings provide key information on structural and functional properties of NagA in P. multocida that could be utilized to design novel antibacterials.

Interplay of self-association and conformational flexibility in regulating protein function

Many functional roles have been attributed to homodimers, the most common mode of protein self-association, notably in the regulation of enzymes, ion channels, transporters and transcription factors. Here we review findings that offer new insights into the different roles conformational flexibility plays in regulating homodimer function. Intertwined homodimers of two-domain proteins and their related family members display significant conformational flexibility, which translates into concerted motion between structural domains. This flexibility enables the corresponding proteins to regulate function across family members by modulating the spatial positions of key recognition surfaces of individual domains, to either maintain subunit interfaces, alter them or break them altogether, leading to a variety of functional consequences. Many proteins may exist as monomers but carry out their biological function as homodimers or higher-order oligomers. We present early evidence that in such systems homodimer formation primes the protein for its functional role. It does so by inducing elevated mobility in protein regions corresponding to the binding epitopes of functionally important ligands. In some systems this process acts as an allosteric response elicited by the self-association reaction itself. Our analysis furthermore suggests that the induced extra mobility likely facilitates ligand binding through the mechanism of conformational selection.This article is part of a discussion meeting issue 'Allostery and molecular machines'.

Self-adaptive multiscaling algorithm for efficient simulations of many-protein systems in crowded conditions

A method is described for the efficient simulation of multiprotein systems in crowded environments. It is based on an adaptive, reversible structural coarsening algorithm that preserves relevant physical features of the proteins across scales. Water is treated implicitly whereas all the other components of the aqueous solution, such as ions, cosolutes, or osmolytes, are treated in atomic detail. The focus is on the analytical adaptation of the solvent model to different levels of molecular resolutions, which allows continuous, on-the-fly transitions between scales. This permits the analytical calculation of forces during dynamics and preserves detailed balance in Monte Carlo simulations. A major computational speedup can be achieved in systems containing hundreds of proteins without cutting off the long-range interactions. The method can be combined with a self-adaptive configurational-bias sampling technique described previously, designed to detect strong, weak, or ultra-weak protein associations and shown to improve sampling efficiency and convergence. The implementation aims to simulate early stages of multimeric complexation, aggregation, or self-assembly. The method can be adopted as the basis for a more general algorithm to identify vertices, edges, and hubs in protein interaction networks or to predict critical steps in signal transduction pathways.

A Structural Model of the Urease Activation Complex Derived from Ion Mobility-Mass Spectrometry and Integrative Modeling

The synthesis of active Klebsiella aerogenes urease via an 18-subunit enzyme apoprotein-accessory protein pre-activation complex has been well studied biochemically, but thus far this complex has remained refractory to direct structural characterization. Using ion mobility-mass spectrometry, we characterized several protein complexes between the core urease apoprotein and its accessory proteins, including the 610-kDa (UreABC)₃(UreDFG)₃ complex. Using our recently developed computational modeling workflow, we generated ensembles of putative (UreABC)₃(UreDFG)₃ species consistent with experimental restraints and characterized the structural ambiguity present in these models. By integrating structural information from previous studies, we increased the resolution of the ion mobility-mass spectrometry-derived models substantially, and we observe a discrete population of structures consistent with all of the available data for this complex.

Template-based quaternary structure prediction of proteins using enhanced profile-profile alignments

Proteins often exist as their multimeric forms when they function as so‐called biological assemblies consisting of the specific number and arrangement of protein subunits. Consequently, elucidating biological assemblies is necessary to improve understanding of protein function. Template‐Based Modeling (TBM), based on known protein structures, has been used widely for protein structure prediction. Actually, TBM has become an increasingly useful approach in recent years because of the increased amounts of information related to protein amino acid sequences and three‐dimensional structures. An apparently similar situation exists for biological assembly structure prediction as protein complex structures in the PDB increase, although the inference of biological assemblies is not a trivial task. Many methods using TBM, including ours, have been developed for protein structure prediction. Using enhanced profile–profile alignments, we participated in the 12th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP12), as the FONT team (Group # 480). Herein, we present experimental procedures and results of retrospective analyses using our approach for the Quaternary Structure Prediction category of CASP12. We performed profile–profile alignments of several types, based on FORTE, our profile–profile alignment algorithm, to identify suitable templates. Results show that these alignment results enable us to find templates in almost all possible cases. Moreover, we have come to understand the necessity of developing a model selection method that provides improved accuracy. Results also demonstrate that, to some extent, finding templates of protein complexes is useful even for MEDIUM and HARD assembly prediction.

The enigma of the near-symmetry of proteins: Domain swapping

The majority of proteins form oligomers which have rotational symmetry. Literature has suggested many functional advantages that the symmetric packing offers. Yet, despite these advantages, the vast majority of protein oligomers are only nearly symmetric. A key question in the field of proteins structure is therefore, if symmetry is so advantageous, why do oligomers settle for aggregates that do not maximize that structural property? The answer to that question is apparently multi-parametric, and involves distortions at the interaction zones of the monomer units of the oligomer in order to minimize the free energy, the dynamics of the protein, the effects of surroundings parameters, and the mechanism of oligomerization. The study of this problem is in its infancy: Only the first parameter has been explored so far. Here we focus on the last parameter-the mechanism of formation. To test this effect we have selected to focus on the domain swapping mechanism of oligomerization, by which oligomers form in a mechanism that swaps identical portions of monomeric units, resulting in an interwoven oligomer. We are using continuous symmetry measures to analyze in detail the oligomer formed by this mechanism, and found, that without exception, in all analyzed cases, perfect symmetry is given away, and we are able to identify that the main burden of distortion lies in the hinge regions that connect the swapped portions. We show that the continuous symmetry analysis method clearly identifies the hinge region of swapped domain proteins-considered to be a non-trivial task. We corroborate our conclusion about the central role of the hinge region in affecting the symmetry of the oligomers, by a special probability analysis developed particularly for that purpose.

Structural Symmetry in Membrane Proteins

Symmetry is a common feature among natural systems, including protein structures. A strong propensity toward symmetric architectures has long been recognized for water-soluble proteins, and this propensity has been rationalized from an evolutionary standpoint. Proteins residing in cellular membranes, however, have traditionally been less amenable to structural studies, and thus the prevalence and significance of symmetry in this important class of molecules is not as well understood. In the past two decades, researchers have made great strides in this area, and these advances have provided exciting insights into the range of architectures adopted by membrane proteins. These structural studies have revealed a similarly strong bias toward symmetric arrangements, which were often unexpected and which occurred despite the restrictions imposed by the membrane environment on the possible symmetry groups. Moreover, membrane proteins disproportionately contain internal structural repeats resulting from duplication and fusion of smaller segments. This article discusses the types and origins of symmetry in membrane proteins and the implications of symmetry for protein function.

How the Sequence of a Gene Specifies Structural Symmetry in Proteins

Internal symmetry is commonly observed in the majority of fundamental protein folds. Meanwhile, sufficient evidence suggests that nascent polypeptide chains of proteins have the potential to start the co-translational folding process and this process allows mRNA to contain additional information on protein structure. In this paper, we study the relationship between gene sequences and protein structures from the viewpoint of symmetry to explore how gene sequences code for structural symmetry in proteins. We found that, for a set of two-fold symmetric proteins from left-handed beta-helix fold, intragenic symmetry always exists in their corresponding gene sequences. Meanwhile, codon usage bias and local mRNA structure might be involved in modulating translation speed for the formation of structural symmetry: a major decrease of local codon usage bias in the middle of the codon sequence can be identified as a common feature; and major or consecutive decreases in local mRNA folding energy near the boundaries of the symmetric substructures can also be observed. The results suggest that gene duplication and fusion may be an evolutionarily conserved process for this protein fold. In addition, the usage of rare codons and the formation of higher order of secondary structure near the boundaries of symmetric substructures might have coevolved as conserved mechanisms to slow down translation elongation and to facilitate effective folding of symmetric substructures. These findings provide valuable insights into our understanding of the mechanisms of translation and its evolution, as well as the design of proteins via symmetric modules.

Paralogue Interference Affects the Dynamics after Gene Duplication

Proteins tend to form homomeric complexes of identical subunits, which act as functional units. By definition, the subunits are encoded from a single genetic locus. When such a gene is duplicated, the gene products are suggested initially to cross-interact when coexpressed, thus resulting in the phenomenon of paralogue interference. In this opinion article, we explore how paralogue interference can shape the fate of a duplicated gene. One important outcome is a prolonged time window in which both copies remain under selection increasing the chance to accumulate mutations and to develop new properties. Thereby, paralogue interference can mediate the coevolution of duplicates and here we illustrate the potential of this phenomenon in light of recent new studies.

Detection and characterization of nonspecific, sparsely populated binding modes in the early stages of complexation

A method is proposed to study protein-ligand binding in a system governed by specific and nonspecific interactions. Strong associations lead to narrow distributions in the proteins configuration space; weak and ultraweak associations lead instead to broader distributions, a manifestation of nonspecific, sparsely populated binding modes with multiple interfaces. The method is based on the notion that a discrete set of preferential first-encounter modes are metastable states from which stable (prerelaxation) complexes at equilibrium evolve. The method can be used to explore alternative pathways of complexation with statistical significance and can be integrated into a general algorithm to study protein interaction networks. The method is applied to a peptide-protein complex. The peptide adopts several low-population conformers and binds in a variety of modes with a broad range of affinities. The system is thus well suited to analyze general features of binding, including conformational selection, multiplicity of binding modes, and nonspecific interactions, and to illustrate how the method can be applied to study these problems systematically. The equilibrium distributions can be used to generate biasing functions for simulations of multiprotein systems from which bulk thermodynamic quantities can be calculated.

Structure, dynamics, assembly, and evolution of protein complexes

The assembly of individual proteins into functional complexes is fundamental to nearly all biological processes. In recent decades, many thousands of homomeric and heteromeric protein complex structures have been determined, greatly improving our understanding of the fundamental principles that control symmetric and asymmetric quaternary structure organization. Furthermore, our conception of protein complexes has moved beyond static representations to include dynamic aspects of quaternary structure, including conformational changes upon binding, multistep ordered assembly pathways, and structural fluctuations occurring within fully assembled complexes. Finally, major advances have been made in our understanding of protein complex evolution, both in reconstructing evolutionary histories of specific complexes and in elucidating general mechanisms that explain how quaternary structure tends to evolve. The evolution of quaternary structure occurs via changes in self-assembly state or through the gain or loss of protein subunits, and these processes can be driven by both adaptive and nonadaptive influences.

Life at the border: adaptation of proteins to anisotropic membrane environment

This review discusses main features of transmembrane (TM) proteins which distinguish them from water-soluble proteins and allow their adaptation to the anisotropic membrane environment. We overview the structural limitations on membrane protein architecture, spatial arrangement of proteins in membranes and their intrinsic hydrophobic thickness, co-translational and post-translational folding and insertion into lipid bilayers, topogenesis, high propensity to form oligomers, and large-scale conformational transitions during membrane insertion and transport function. Special attention is paid to the polarity of TM protein surfaces described by profiles of dipolarity/polarizability and hydrogen-bonding capacity parameters that match polarity of the lipid environment. Analysis of distributions of Trp resides on surfaces of TM proteins from different biological membranes indicates that interfacial membrane regions with preferential accumulation of Trp indole rings correspond to the outer part of the lipid acyl chain region-between double bonds and carbonyl groups of lipids. These "midpolar" regions are not always symmetric in proteins from natural membranes. We also examined the hydrophobic effect that drives insertion of proteins into lipid bilayer and different free energy contributions to TM protein stability, including attractive van der Waals forces and hydrogen bonds, side-chain conformational entropy, the hydrophobic mismatch, membrane deformations, and specific protein-lipid binding.

Structural templates for modeling homodimers

Oligomeric proteins are more abundant in nature than monomeric proteins, and involved in all biological processes. In the absence of an experimental structure, their subunits can be modeled from their sequence like monomeric proteins, but reliable procedures to build the oligomeric assembly are scarce. Template-based methods, which start from known protein structures, are commonly applied to model subunits. We present a method to model homodimers that relies on a structural alignment of the subunits, and test it on a set of 511 target structures recently released by the Protein Data Bank, taking as templates the earlier released structures of 3108 homodimeric proteins (H-set), and 2691 monomeric proteins that form dimer-like assemblies in crystals (M-set). The structural alignment identifies a H-set template for 97% of the targets, and in half of the cases, it yields a correct model of the dimer geometry and residue-residue contacts in the target. It also identifies a M-set template for most of the targets, and some of the crystal dimers are very similar to the target homodimers. The procedure efficiently detects homology at low levels of sequence identities, and points to erroneous quaternary structures in the Protein Data Bank. The high coverage of the target set suggests that the content of the Protein Data Bank already approaches the structural diversity of protein assemblies in nature, and that template-based methods should become the choice method for modeling oligomeric as well as monomeric proteins.

Different 3D domain-swapped oligomeric cyanovirin-N structures suggest trapped folding intermediates

Although it has long been established that the amino acid sequence encodes the fold of a protein, how individual proteins arrive at their final conformation is still difficult to predict, especially for oligomeric structures. Here, we present a comprehensive characterization of oligomeric species of cyanovirin-N that all are formed by a polypeptide chain with the identical amino acid sequence. Structures of the oligomers were determined by X-ray crystallography, and each one exhibits 3D domain swapping. One unique 3D domain-swapped structure is observed for the trimer, while for both dimer and tetramer, two different 3D domain-swapped structures were obtained. In addition to the previously identified hinge-loop region of the 3D domain-swapped dimer, which resides between strands β5 and β6 in the middle of the polypeptide sequence, another hinge-loop region is observed between strands β7 and β8 in the structures. Plasticity in these two regions allows for variability in dihedral angles and concomitant differences in chain conformation that results in the differently 3D domain-swapped multimers. Based on all of the different structures, we propose possible folding pathways for this protein. Altogether, our results illuminate the amazing ability of cyanovirin-N to proceed down different folding paths and provide general insights into oligomer formation via 3D domain swapping.

Influence of homology and node age on the growth of protein-protein interaction networks

Proteins participating in a protein-protein interaction network can be grouped into homology classes following their common ancestry. Proteins added to the network correspond to genes added to the classes, so the dynamics of the two objects are intrinsically linked. Here we first introduce a statistical model describing the joint growth of the network and the partitioning of nodes into classes, which is studied through a combined mean-field and simulation approach. We then employ this unified framework to address the specific issue of the age dependence of protein interactions through the definition of three different node wiring or divergence schemes. A comparison with empirical data indicates that an age-dependent divergence move is necessary in order to reproduce the basic topological observables together with the age correlation between interacting nodes visible in empirical data. We also discuss the possibility of nontrivial joint partition and topology observables.

Conservation of functionally important global motions in an enzyme superfamily across varying quaternary structures

The α-d-phosphohexomutase superfamily comprises enzymes involved in carbohydrate metabolism that are found in all kingdoms of life. Recent biophysical studies have shown for the first time that several of these enzymes exist as dimers in solution, prompting an examination of the oligomeric state of all proteins of known structure in the superfamily (11 different proteins; 31 crystal structures) via computational and experimental analyses. We find that these proteins range in quaternary structure from monomers to tetramers, with 6 of the 11 known structures being likely oligomers. The oligomeric state of these proteins not only is associated in some cases with enzyme subgroup (i.e., substrate specificity) but also appears to depend on domain of life, with the two archaeal proteins existing as higher-order oligomers. Within the oligomers, three distinct interfaces are observed, one of which is found in both archaeal and bacterial proteins. Normal mode analysis shows that the topological arrangement of the oligomers permits domain 4 of each protomer to move independently as required for catalysis. Our analysis suggests that the advantages associated with protein flexibility in this enzyme family are of sufficient importance to be maintained during the evolution of multiple independent oligomers. This study is one of the first showing that global motions may be conserved not only within protein families but also across members of a superfamily with varying oligomeric structures.

A novel method for protein-protein interaction site prediction using phylogenetic substitution models

Protein-protein binding events mediate many critical biological functions in the cell. Typically, functionally important sites in proteins can be well identified by considering sequence conservation. However, protein-protein interaction sites exhibit higher sequence variation than other functional regions, such as catalytic sites of enzymes. Consequently, the mutational behavior leading to weak sequence conservation poses significant challenges to the protein-protein interaction site prediction. Here, we present a phylogenetic framework to capture critical sequence variations that favor the selection of residues essential for protein-protein binding. Through the comprehensive analysis of diverse protein families, we show that protein binding interfaces exhibit distinct amino acid substitution as compared with other surface residues. On the basis of this analysis, we have developed a novel method, BindML, which utilizes the substitution models to predict protein-protein binding sites of protein with unknown interacting partners. BindML estimates the likelihood that a phylogenetic tree of a local surface region in a query protein structure follows the substitution patterns of protein binding interface and nonbinding surfaces. BindML is shown to perform well compared to alternative methods for protein binding interface prediction. The methodology developed in this study is very versatile in the sense that it can be generally applied for predicting other types of functional sites, such as DNA, RNA, and membrane binding sites in proteins.

Structural model of ligand-G protein-coupled receptor (GPCR) complex based on experimental double mutant cycle data: MT7 snake toxin bound to dimeric hM1 muscarinic receptor

The snake toxin MT7 is a potent and specific allosteric modulator of the human M1 muscarinic receptor (hM1). We previously characterized by mutagenesis experiments the functional determinants of the MT7-hM1 receptor interaction (Fruchart-Gaillard, C., Mourier, G., Marquer, C., Stura, E., Birdsall, N. J., and Servent, D. (2008) Mol. Pharmacol. 74, 1554-1563) and more recently collected evidence indicating that MT7 may bind to a dimeric form of hM1 (Marquer, C., Fruchart-Gaillard, C., Mourier, G., Grandjean, O., Girard, E., le Maire, M., Brown, S., and Servent, D. (2010) Biol. Cell 102, 409-420). To structurally characterize the MT7-hM1 complex, we adopted a strategy combining double mutant cycle experiments and molecular modeling calculations. First, thirty-three ligand-receptor proximities were identified from the analysis of sixty-one double mutant binding affinities. Several toxin residues that are more than 25 Å apart still contact the same residues on the receptor. As a consequence, attempts to satisfy all the restraints by docking the toxin onto a single receptor failed. The toxin was then positioned onto two receptors during five independent flexible docking simulations. The different possible ligand and receptor extracellular loop conformations were described by performing simulations in explicit solvent. All the docking calculations converged to the same conformation of the MT7-hM1 dimer complex, satisfying the experimental restraints and in which (i) the toxin interacts with the extracellular side of the receptor, (ii) the tips of MT7 loops II and III contact one hM1 protomer, whereas the tip of loop I binds to the other protomer, and (iii) the hM1 dimeric interface involves the transmembrane helices TM6 and TM7. These results structurally support the high affinity and selectivity of the MT7-hM1 interaction and highlight the atypical mode of interaction of this allosteric ligand on its G protein-coupled receptor target.