# Transmembrane Recognition of the Semaphorin Co-Receptors Neuropilin 1 and Plexin A1: Coarse-Grained Simulations

Samia Aci-Sèche, Paul Sawma, Pierre Hubert, James N. Sturgis, Dominique Bagnard, Laurent Jacob, Monique Genest, Norbert Garnier

► To cite this version:

Samia Aci-Sèche, Paul Sawma, Pierre Hubert, James N. Sturgis, Dominique Bagnard, et al. Trans-membrane Recognition of the Semaphorin Co-Receptors Neuropilin 1 and Plexin A1: Coarse-Grained Simulations. PLoS ONE, 2014, 9 (5), pp.e97779. 10.1371/journal.pone.0097779. hal-01175479

HAL Id: hal-01175479

https://hal.science/hal-01175479v1

Submitted on 20 Sep 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Transmembrane Recognition of the Semaphorin Co-Receptors Neuropilin 1 and Plexin A1: Coarse-Grained Simulations

Samia Aci-Sèche¹*, Paul Sawma², Pierre Hubert², James N. Sturgis², Dominique Bagnard³, Laurent Jacob³, Monique Genest¹*, Norbert Garnier¹*

1 Centre de Biophysique Moléculaire, Centre National de la Recherche Scientifique UPR, Université d'Orléans, Orléans, France, 2 Laboratoire d'Ingénierie des Systèmes Macromoléculaires UMR, Université d'Aix-Marseille, Marseille, France, 3 Institut National de la Santé et de la Recherche Médicale, Labex Medalis, Université de Strasbourg, Strasbourg, France

<figure><img src="image_2.png" /></figure>

## Abstract

The cancer associated class 3 semaphorins require direct binding to neuropilins and association to plexins to trigger cell signaling. Here, we address the role of the transmembrane domains of neuropilin 1 and plexin A1 for the dimerization of the two receptors by characterizing the assembly in lipid bilayers using coarse-grained molecular dynamics simulations. From experimental evidence using a two-hybrid system showing the biochemical association of the two receptors transmembrane domains, we performed molecular simulations in DOPC and POPC demonstrating spontaneously assembly to form homodimers and heterodimers with a very high propensity for right-handed packing of the helices. Inversely, left-handed packing was observed with a very low propensity. This mode of packing was observed uniquely when the plexin A1 transmembrane domain was involved in association. Potential of mean force calculations were used to predict a hierarchy of self-association for the monomers: the two neuropilin 1 transmembrane domains strongly associated, neuropilin 1 and plexin A1 transmembrane domains associated less and the two plexin A1 transmembrane domains weakly but significantly associated. We demonstrated that homodimerization and heterodimerization are driven by GxxxG motifs, and that the sequence context modulates the packing mode of the plexin A1 transmembrane domains. This work presents major advances towards our understanding of membrane signaling platforms assembly through membrane domains and provides exquisite information for the design of antagonist drugs defining a novel class of therapeutic agents.

**Citation:** Aci-Sèche S, Sawma P, Hubert P, Sturgis JN, Bagnard D, et al. (2014) Transmembrane Recognition of the Semaphorin Co-Receptors Neuropilin 1 and Plexin A1: Coarse-Grained Simulations. PLoS ONE 9(5): e97779. doi:10.1371/journal.pone.0097779

**Editor:** Claudio M. Soares, Instituto de Tecnologia Química e Biológica, UNL, Portugal

**Received** February 6, 2014; **Accepted** April 23, 2014; **Published** May 23, 2014

**Copyright:** © 2014 Aci-Sèche et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This work was supported by a grant from the French Agence Nationale de la Recherche (InterferenceTM project; ANR-10-BLAN-1507). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* E-mail: norbert.garnier@cnrs-orleans.fr (NG); monique.genest@cnrs-orleans.fr (MG)

‡ Current address: Institut de Chimie Organique et Analytique UMR, Université d'Orléans, Orléans, France

## Introduction

Most of cell activities are controlled by membrane receptors, which in response to extracellular ligands trigger intracellular signaling reactions. A hallmark of these signaling events is the reversible formation of multicomponent complexes, which associate membrane and cytosolic proteins through specialized protein binding modules [1]. The first step in the assembly of these multiprotein complexes is usually the dimerization/oligomerization of the receptors themselves, or structural rearrangements of preformed oligomers. Examples of such receptor complexes include receptor tyrosine kinases such as the ErbB family or Ephrin receptors [2], the T-cell receptor complex which is part of the immunological synapse [3], integrins [4], G-protein coupled receptors [5] among others.

In many cases, the assembly of these signaling platforms involves interactions between the transmembrane domains of the receptors. Different studies have reported the functional importance of intramembrane interactions for membrane protein signaling [6–12]. The role of the transmembrane (TM) domains

depends on their propensity to associate and on the existence of specific dimerization motifs in their sequence. The core of the dimerization motif of TM domains is often a GxxxG sequence (where G is glycine or another small residue, and x represents any amino-acid). This motif was first recognized in the case of glycophorin A and was since found in many other examples of interacting TM helices [13–18]. A detailed view on how such TM interactions are involved in the mechanism of activation of a receptor has very recently been provided for the EGF receptor, a receptor tyrosine kinase. Using NMR, molecular dynamics simulations and biochemical assays, J. Kuriyan's group has shown that ligand binding causes a rearrangement of TM helices in an EGFR dimer, so that the juxtamembrane and kinase domains can form the activated asymmetric dimer [19,20].

Semaphorins form a large family of membrane-associated and secreted proteins playing many roles in a variety of cellular processes. They are bifunctional signaling molecules capable of growth promoting or growth inhibitory effects [21]. The diversity of functions encoded by semaphorins is related to the formation of specific receptor complexes. Together with their receptors, the
neuropilins and the plexins, semaphorins are the constituents of a complex regulatory system responsible for axon guidance during the development of the central nervous system [22,23]. Neuropilins are also involved in vascular development and tumorigenesis [24]. The secreted semaphorins of subclass 3 (Sema3A-Sema3G) bind directly to neuropilins with different affinities for the two members of the family, neuropilin 1 and neuropilin 2 [25,26] but do not bind directly to plexins [27–29]. The functional receptor for secreted semaphorins is a complex including neuropilins and plexins. Thus, neuropilins form signaling complexes by associating with type A plexins to activate downstream signal transduction cascades. In these complexes, neuropilins are not able to transduce semaphorin signals owing to their short intracellular domains but act as binding receptors while plexins act as signal-transducing components [30].

Various studies have revealed that the extracellular domain of plexin A1 (PLXA1) interacts with neuropilin 1 (NRP1), allowing a Sema3 signal directly to the cytoplasm via the intracellular domain of plexin A1 [31–33]. It appears that the sema domain at the extracellular region of semaphorins and plexins is a common structural element able to mediate a large variety of protein–protein interactions, particularly semaphorin-plexin binding and semaphorin-neuropilin binding. The sema domain of plexins has been suggested to mediate the formation of ligand independent homodimers. Moreover, with their ligands, homodimers represent the active forms of Sema3A and likely of all semaphorins [34,35]. Structural studies have partially elucidated the binding mode between plexins and semaphorins [36–40]. Nevertheless, the data do not clarify the mechanism by which semaphorin-plexin A1 interactions, together with neuropilin binding, lead to the activation of the cytoplasmic region of plexins.

Several hypotheses have been proposed [32,40–42]. Recent findings clearly show that semaphorin dimers are needed for signaling [36,39] and that the core of the mechanism of semaphorin-mediated plexin dimerization is central to the Sema3 function. The requirement of neuropilin as co-receptor to stabilize this complex is believed to induce neuropilin dimerization/oligomerization or receptor clustering and may possibly induce clusters of the plexin intracellular region [36,39,43,44]. The architecture of the Sema3-plexin A complexes and the role of the extracellular domains in dimerization have been revealed by crystal structures [36,38,39]. Recently, the central role of the extracellular domain of neuropilins in Sema3-plexin A signaling has been demonstrated [43]. Thus, the extracellular domains are required for formation of the complex, but very probably, other binding sites are necessary to trigger the function of the semaphorin-neuropilin 1-plexin A1 complex. Identifying these sites represents an important issue in understanding the dynamic behavior of these receptors.

One such site could be the TM domain of neuropilin 1 (NRP1 TM) which contains a conserved double GxxxG motif. Mutations of this motif confirmed its biological importance for Sema3A signaling [45]. Moreover, this same study showed that this domain exhibits a strong dimerization capacity. Furthermore, the TM domain of plexin A1 (PLXA1 TM) contains six Gly residues constituting a GxGGGG (or GxG5) motif. This Gly rich motif could thus be expected to be involved in helix-helix interactions. Unfortunately, experimental data on the role of this TM domain in association processes are not available to date. These motifs have been widely described as mediators for TM dimerization and oligomerization. For neuropilin 1, the ability of the TM domain to form dimers was initially demonstrated by a two-hybrid assay performed in bacterial membranes and by Förster Resonance Energy Transfer (FRET) analysis [45]. A major functional role for

these TM interactions was evidenced by the potent antagonistic effect of a synthetic peptide mimicking the NRP1 TM domain on neuropilin signaling, both *in vitro* and *in vivo* [21,46]. Mutations of the GxxxG motifs rendered the peptide and the receptor completely inactive, thus confirming their functional importance. Moreover peptides targeting and inhibiting the TM of NRP1 turn out to be very potent anti-cancer drug [46] contributing to the identification of NRP1 as a major therapeutic target in cancer [47].

Insights into activation mechanisms require information about the conformational dynamics of TM interactions, not yet available experimentally. Consequently, to explore the relationship between semaphorin binding proteins, we sought to characterize the molecular interactions between these TM domains and determine how the putative interaction motifs contribute to the stabilization/destabilization of the associated TM complexes.

In this study, we first used a two-hybrid system to demonstrate homotypic and heterotypic interactions of the TM sequences of NRP1 and PLXA1. Then, we employed coarse-grained molecular dynamics (CG-MD) simulations to explore the self-assembly of the NRP1 and PLXA1 TM domains within lipid bilayers [48–50] composed of pure DOPC (1,2-di-oleoyl-sn-glycero-3-phosphocholine) and pure POPC (palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine). The success of these techniques in understanding the behavior of lipid membranes and membrane proteins [51] opens the route to investigating the association properties of these domains in homotypic and heterotypic forms. The potential of mean force (PMF) has been calculated and differences in free energies of association are compared. To the best of our knowledge, this work is the first devoted to the investigation of the association of the TM domains of the semaphorin co-receptors in membrane bilayers using CG-MD simulations.

## Results

### Dimerization Propensity of NRP1 and PLXA1 TM Sequences in BACTH

To assess first the propensities of the TM domains of NRP1 and PLXA1, we employed a bacterial adenylate cyclase two-hybrid assay (BACTH). We used the TM sequence of human Glycophorin A (GpA) as both a positive control for homotypic interactions and a negative control for heterotypic ones. Figure 1 depicts the results of BACTH functional complementation between the TM domains of GpA, NRP1 and PLXA1 (Table 1) determined by growing transformed cells on an indicator plate (panel A) and by measuring β-galactosidase activity (panel B). Both tests indicate clearly that NRP1 and PLXA1 TM domains can self-associate, but also form heterodimers. Expression levels of constructs was checked by Western blotting and found to be similar for all proteins (Sawma et al., manuscript submitted). The interactions order is as follows: NRP1-NRP1> NRP1-PLXA1> PLXA1-PLXA1 The NRP1-NRP1 interaction propensity is about 50% of that of GpA-GpA, and specificity of the assay is shown by the absence of interactions of the GpA sequence with the other two.

### Self-assembly of TM Domains by Molecular Dynamics

All the simulations performed in DOPC and POPC bilayers led to the association of the TM domains after their diffusion within the membrane. The time required for the formation of the two homodimers and the heterodimer was variable, from 1 µs to 10 µs, depending on the simulations. The overall simulations indicate that once the TM domains associate, helix-helix contacts were maintained for the rest of the simulation. The time during
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Bacterial adenyl cyclase two-hybrid (BACTH) assay of NRP1 and PLXA1 TM sequences. *E. coli* BTH101 host cells were cotransformed with a pair of plasmids encoding adenyl cyclase domains T18 and T25 fused to the indicated TM sequences. Cells were directly plated on MacConkey indicator plates (panel A). The fusions with the TM domain of glycophorin A (GpA) are a positive control for homotypic interactions. Red coloration of the colonies indicates a cya+ phenotype and thus interaction of the two TM domains. In panel B, the efficiencies of functional complementation between the indicated TM domains were quantified by measuring β-galactosidase activities in *E. coli* BTH101 harboring the corresponding plasmids as described in Materials and Methods. Values are normalized relatively to GpA taken as 100%. doi:10.1371/journal.pone.0097779.g001</figcaption></figure>

which the TM domains remained associated is significant compared to the total time of each simulation. No dissociation event was seen with the exception of only one simulation conducted for the two PLXA1 monomers embedded in a POPC bilayer.

All RMSD maps and time-courses of the crossing angles were calculated. Three of these maps and their corresponding crossing angles (Figure 2) were chosen to illustrate the formation of each of the three dimers and their behavior during a complete simulation performed in the POPC environment. Self-assembly of the TM domains may happen very early as illustrated here for the NRP1 and PLXA1 monomers (around 1 µs) or later as seen for the two NRP1 monomers (around 10 µs). RMSD maps reveal successive

periods for which RMSD values are lower than 0.6 nm, evidencing discrete dimer states very close to each other. These successive periods of low RMSD reflect the dynamics of the two associated monomers within the bilayer.

The RMSD map and the crossing angle representing the time-
course of the self-assembly of the PLXA1 TM domains summarize
all the characteristics observed for this dimer (Figure 2B). The two
periods between 9 µs and 14 µs and between 22 µs and 35 µs
reveal two different states, one with the monomers associated in a
left-handed manner and the second with the monomers associated
in a right-handed manner respectively. These two states are
separated by about 5 µs during which the two monomers are fully
<table><caption>**Table 1.** TM domain sequences of human Neuropilin-1 (NRP1) and Plexin-A1 (PLXA1), protein names, identifiers, entry names taken from the Uniprot database (http://www.uniprot.org/) and TM domain sequences used in CG-MD simulations.</caption><thead><tr><th>TM domain</th><th>Sequence</th><th>Protein residues</th></tr></thead><tbody><tr><td><strong>Neuropilin-1</strong><br>(O14786, NRP1_HUMAN)</td><td>ILITIIAMSALGVLLGAVCGVVL</td><td>857-879 (23 aa)</td></tr><tr><td>NRP1 TM domain (MD)</td><td><sup>1</sup>TLDPILITIIAMSALGVLLGAVCGVVLY<sup>28</sup></td><td>853-880 (28 aa)</td></tr><tr><td><strong>Plexin-A1</strong><br>(Q9UIW2, PLXA1_HUMAN)</td><td>AIVGIGGGGGLLLLVIVAVLI</td><td>1245-1265 (21 aa)</td></tr><tr><td>PLXA1 TM domain (MD)</td><td><sup>1</sup>LLTLPAIVGIGGGGGLLLLVIVAVLIAY<sup>28</sup></td><td>1240-1267 (28 aa)</td></tr></tbody></table>

For the sake of clarity, the TM sequences have been re-numbered 1 to 28 throughout this paper.
doi:10.1371/journal.pone.0097779.t001

dissociated. Among all the simulations, this is the only one that reveals a dissociation event.

The complete results of self-simulations of the three dimers embedded in a DOPC bilayer and in a POPC bilayer are given in Table 2. The main feature is that the right-handed association is highly preponderant for the three dimers in the two lipid environments. The crossing angles are similar in the two milieus between -33° and -37° on average with relatively small fluctuations of about 5°-7°. Besides this major property of the three dimers, each has its own particularities.

The very high propensity of the NRP1 monomers to dimerize in right-handed interactions is clearly evidenced in both types of

bilayers. This right-handed packing induces short inter-helix contact distances of 0.47 nm with very low fluctuations. Also, inter-helix distances are short, slightly shorter in POPC than in DOPC by 0.1 nm, indicating a weak effect of the nature of the lipids on the TM domains assembly.

The behavior of the PLXA1 TM domains is quite different. First, the two monomers associate spontaneously either in right-handed or in left-handed interactions. Then, for the rest of the simulation, the two monomers cross successively with a positive or a negative angle (or inversely) without dissociation of the monomers, except for the dissociation event mentioned above. For each association mode the mean values of the crossing angle

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Time-course of the self-assembly of the dimers embedded in a POPC bilayer. RMSD maps and crossing angles illustrate the right-handed association of the TM domains NRP1-NRP1 (A) and NRP1-PLXA1 (C) (negative value of the crossing angle). The two TM domains of PLXA1 (B) associate alternatively in a right-handed or in a left-handed manner along the trajectories. The arrow indicates a transient dissociation event for ~5 µs seen in one simulation from the left to the right-handed structure.
doi:10.1371/journal.pone.0097779.g002</figcaption></figure>
<table><caption>**Table 2.** Summary of the self-assembly simulations of the dimers embedded in a DOPC bilayer and in a POPC bilayer.</caption><thead><tr><th>TM Dimers</th><th>Lipid bilayer</th><th>Total simulation time (is)</th><th>Total association time (is)</th><th>Structure</th><th>% association time (L, R)</th><th>Crossing angle (°)</th><th>Inter-helix COM distance (nm)</th><th>Inter-helix contact distance (nm)</th></tr></thead><tbody><tr><td rowspan="2">NRP1 homodimer</td><td>DOPC</td><td>100 (4)</td><td>75</td><td>L</td><td>0</td><td>100</td><td>0.63 (0.12)</td><td>0.47 (0.01)</td></tr><tr><td>POPC</td><td>110 (4)</td><td>93</td><td>L</td><td>0</td><td>-34 (5)</td><td></td><td></td></tr><tr><td rowspan="3">PLXA1 homodimer</td><td>DOPC</td><td>80 (3)</td><td>64</td><td>L</td><td>100</td><td>-37 (5)</td><td>0.54 (0.12)</td><td>0.47 (0.01)</td></tr><tr><td>POPC</td><td>85 (3)</td><td>57</td><td>L</td><td>70</td><td>-33 (7)</td><td>0.59 (0.16)</td><td>0.48 (0.03)</td></tr><tr><td>POPC</td><td>118 (6)</td><td>83</td><td>L</td><td>13</td><td>15 (11)</td><td>0.89 (0.15)</td><td>0.48 (0.07)</td></tr><tr><td rowspan="3">NRP1-PLXA1 heterodimer</td><td>DOPC</td><td>120 (6)</td><td>101</td><td>L</td><td>87</td><td>-35 (6)</td><td>0.59 (0.15)</td><td>0.48 (0.01)</td></tr><tr><td>POPC</td><td>120 (6)</td><td>101</td><td>L</td><td>91</td><td>-33 (5)</td><td>0.57 (0.12)</td><td>0.47 (0.01)</td></tr><tr><td>POPC</td><td>120 (6)</td><td>101</td><td>L</td><td>1</td><td>20 (8)</td><td>0.95 (0.15)</td><td>0.48 (0.02)</td></tr><tr><td></td><td></td><td></td><td></td><td>R</td><td>99</td><td>-34 (5)</td><td>0.57 (0.12)</td><td>0.47 (0.01)</td></tr></tbody></table>

The total simulation time is the sum of the simulation time of the separate MDs. The number in parentheses represents the number of independent simulations. The total association time is the sum of the time of association of the TM domains observed along the separate MDs. The percentage of the association time of the right-handed structures (R) and the left-handed (L) structures is calculated over the total association time. The values of the crossing angle and the inter-helix distances (COM distances and contact distances) are an average calculated over the last 10 µs of each of the simulations in the cases of the NRP1 homodimer and the heterodimer. For the PLXA1 homodimer the values are an average calculated over the total time of association in right-handed and left-handed interactions.
doi:10.1371/journal.pone.0097779.t002
are practically identical in DOPC and POPC in the limit of the fluctuations. These fluctuations are almost twofold greater for a positive angle than for a negative angle, and their magnitude is in the same order as the mean values of the crossing angle. These data reveal the instability of the left-handed structure for this homodimer. The strong dynamics of the monomers crossed in left-handed interactions does not affect their close approach. The inter-helix contact distances are 0.48 nm on average, as found for the right-handed structures in both types of lipids. It should be pointed out that the inter-helix distances depend on the manner in which the two TM domains are crossed, as attested by a difference of about 0.3 nm between the two modes of association.

The propensity of the PLXA1 TM domains to pack in a right-handed or in a left-handed structure is modulated by the nature of the lipids. Left-handed structures are observed for 30% of the total time of association in the DOPC bilayer while in the POPC bilayer they are observed for 13% of this total time. This finding is of interest because it suggests the possibility of equilibrium between the two modes of association.

The self-assembly of the NRP1 and PLXA1 TM domains shows an intermediate behavior between what is observed for the two TM domains NRP1 and PLXA1. The simulations demonstrate a net propensity for the formation of a right-handed structure. Interestingly, left-handed structures are formed but this mode of association is much less common than for the PLXA1 homodimer. These structures are observed for a very low percentage of time in the DOPC bilayer (8%) and for a negligible time in the POPC bilayer (≤1%). We note that for this left-handed association, the two NRP1 and PLXA1 monomers are slightly more crossed than the two PLXA1 monomers by about 7° in DOPC and 5° in POPC on average. Inter-helix contact distances are equal on average to those reported for the two homodimers. Inter-helix distances are in the same range as those measured for the PLXA1 homodimer. As mentioned for the PLXA1 homodimer, inter-helix distances in left-handed structures are about 0.4 nm larger than in right-handed structures.

## Helix-helix Interface

Helix-helix contacts were evaluated by calculating the average distance matrices. The contact maps exhibit the residue pairs having the smallest distances between the backbone beads on average over a time range for which the dimer is completely stabilized. All the distance matrices corresponding to the two homodimers and the heterodimer were calculated and the resulting contact maps were compared. An overview of all these contact maps indicates that the number of residue pairs having the smallest distance between their backbone beads is greater in POPC than in DOPC. The maps in figure 3 display the major characteristics of the helix-helix interfaces of each dimer. Representative structures in the lipid bilayer are shown in figure 4 with key interfacial residues highlighted.

NRP1 homodimer contact maps are all asymmetric, suggesting that the two helices are slightly shifted along the Z axis and/or slightly rotated about their own helical axis. The variability of the helix-helix interfaces is clearly evidenced but all exhibit the same key residues that confer dimer stability, namely residues Met12, Gly16 and Gly20, which define the interacting motif M12xxxG16xxxG20. Gly24, the last residue of the double motif is not systematically observed at the interface and the contact distances are generally greater than 0.6 nm. Gly24 does not seem to be a key residue for helix-helix stability. Other residues in proximity to these key residues or aligned on the same helix face are observed at the interface, such as the Thr8 residue seen in the POPC bilayer.

The right-handed and the left-handed structures of the PLXA1 homodimer induce significant differences in helix-helix interfaces. The contact maps associated to right-handed packing are practically symmetric whereas those associated to left-handed packing are asymmetric. The right-handed structures exhibit well conserved interfacing residues namely Gly11, Gly15, and Leu19 included in the motif G11xxxG15xxxL19. The presence of Leu19 at the interface is better represented in POPC than in DOPC. The surrounding residues Ile7 and Gly12 are frequently observed at the interface. The predominant interfacing residues of the left-handed structures are Ala6, Gly9, and Gly13. These three key residues form the interacting motif A6xxG9xxxG13 both in DOPC and in POPC. Ile10 is also frequently observed at the helix-helix interface.

The right-handed structure of the heterodimer shows that each of the two monomers exhibits the same interfacing residues identified in the two right-handed structures of the homodimers. Contact maps are generally asymmetric but all exhibit the same pattern with the key interfacial residues Met12, Gly16 and Gly20 for the NRP1 TM domain and Gly11, Gly15, and Leu19 for the PLXA1 TM domain. Residues adjacent to the key residues are also involved in interfacial contacts. Contact distances are generally larger than 0.6 nm but their presence is almost systematic.

## Free Energy of Dissociation

The free energies of dissociation of the TM domains were evaluated in the DOPC and POPC environments. For the two milieus, the PMF profiles were calculated by considering the right-handed structure of the homo and hetero dimers. In addition, the free energy of dissociation of the PLXA1 TM domains packed in a left-handed manner was evaluated only in DOPC where this packing mode is well represented (see Table 2). The very weak representation of this structure in POPC suggests a very unstable dimer making this state not pertinent for PMF calculations. The results are summarized in Table 3 and the PMF profiles are shown in figure 5.

PMF profiles calculated in DOPC show that the lowest global minimum is obtained for the NRP1 homodimer with a free energy of dissociation estimated at -68.7 kJ mol⁻¹. Compared to the values evaluated for the right-handed structures of the heterodimer and the PLXA1 homodimer, differences in free energies are about 10 kJ mol⁻¹ and 32 kJ mol⁻¹ respectively. Therefore, the separation of the two PLXA1 monomers crossed in right-handed interactions appears relatively easy compared to the two other dimers. However, the separation of these monomers crossed in left-handed interactions is even easier. The very large difference in free energy of dissociation of 49 kJ mol⁻¹, relative to the lowest minimum, evidences a very weak power of association in a left-handed mode. We note that these differences in free energies calculated from equation (2) are in the same order of magnitude.

PMF profiles differ by the deepness of the global minima but also by their position. The minimum is at 0.47 nm for the right-handed PLXA1 homodimer, which corresponds to the shortest inter-helix distance compared to 0.56 nm for the heterodimer and 0.65 nm for the NRP1 homodimer. The position of the minimum for the left-handed structure of the PLXA1 homodimer corresponds to the largest inter-helix distance, 0.32 nm larger than for the right-handed mode. As a result, the distance of approach of the monomers is not correlated to the strength of the interactions. In both cases the free energy increases more or less rapidly with different regimes during the separation of the monomers but without energy barriers. The full separated state is reached at ~1.6 nm for all the dimers.
<figure><img src="image_5.png" /><figcaption>**Figure 3. Contact maps of dimers embedded in a POPC bilayer.** Top: right-handed structure of the NRP1 homodimer and the NRP1-PLXA1 heterodimer. Bottom: left-handed structure (L) and right-handed structure (R) of the PLXA1 homodimer. Average interfaces highlight the key interfacing residues at distances 0.5 nm<d<0.8 nm. Squares in black, dark gray, and light gray denote the residue pairs at backbone bead distances <0.6 nm, <0.7 nm and <0.8 nm respectively.
doi:10.1371/journal.pone.0097779.g003</figcaption></figure>

PMFs calculated in the POPC environment exhibit the same features as those observed in DOPC (figure 5). The deepest global minimum is observed for the NRP1 homodimer (−73.4 kJ mol⁻¹) very close to the deep minimum observed for the heterodimer. The difference in free energies is around ~6 kJ mol⁻¹. More remarkable is the very weak free energy of dissociation of the two PLXA1 TM domains evaluated at −29.3 kJ mol⁻¹ which means a destabilization of ~44 kJ mol⁻¹ and 38 kJ mol⁻¹ compared to the

NRP1 homodimer and the heterodimer respectively. Here again, the scale of the free energies is in agreement with that calculated from equation (2). As observed previously, PMF profiles exhibit several regimes over the course of the separation of the monomers without energy barriers. The positions of the global minima are practically identical to those observed in DOPC, albeit slightly shifted towards larger inter-helix distances, at the most 0.02 nm. The shortest inter-helix distance is observed for the PLXA1
<figure><img src="image_6.png" /><figcaption>**Figure 4.** **Representation of the dimers in the POPC bilayer.** Key interfacing residues depicted on the contact maps (Figure 3) are highlighted in bead representations. A surface representation figures out the TM domain backbones. The NRP1 TM monomers are figured in blue or cyan. The PLXA1 TM monomers are figured in red or pink. The grey surfaces represent the polar heads of the POPC. **A:** NRP1 homodimer; **B:** NRP1-PLXA1 heterodimer; **C** and **D:** PLXA1 homodimer in left-handed structure (L) and in right handed structure (R) respectively. doi:10.1371/journal.pone.0097779.g004</figcaption></figure>

homodimer. The full separated states occur at inter-helix distances of 1.8 nm for the NRP1 homodimer and the heterodimer and 1.6 nm for the PLXA1 homodimer.

PMFs calculations demonstrate the effect of the environment on the dissociation of the monomers. This effect is weak for the NRP1 homodimer and the heterodimer. Free energies of dissociation are slightly greater in POPC than in DOPC with a gain of −4.7 kJ mol⁻¹ and −8.1 kJ mol⁻¹ respectively. These values indicate a weak stabilization effect of the assembly of the TM domains. Lipid

effects are more clearly evidenced for the PLXA1 homodimer. This dimer is unstable in the two milieus, but this instability is intensified in POPC. The differences in the apparent dissociation free energies ΔΔGdis reveal an increase in destabilization of the PLXA1 homodimer of over 15 kJ mol⁻¹ in POPC with respect to the NRP1 homodimer.

These calculations evidence the effect of the lipids on the dimerization power of TM domains. They clearly demonstrate the weakness for self-association of the PLXA1 TM domains both in

<figure><img src="image_7.png" /><figcaption>**Figure 5.** **Potential of mean force evaluated by umbrella sampling calculations at a sampling of 2 µs per window in DOPC and POPC of the NRP1 and PLXA1 homodimers and the NRP1-PLXA1 heterodimer.** The PMF profiles correspond to the right-handed structure of the dimers in both bilayers. The PMF profile evaluated for the left-handed structure of the PLXA1 homodimer seen in the DOPC bilayer is shown. doi:10.1371/journal.pone.0097779.g005</figcaption></figure>
<table><caption>**Table 3.** Characteristics of the PMF profiles calculated for the dimers in DOPC and POPC bilayers.</caption><thead><tr><th rowspan="2">Dimer</th><th colspan="2">DOPC</th><th colspan="2">POPC</th><th colspan="2">DOPC</th><th colspan="2">POPC</th></tr><tr><th>E<sub>min</sub></th><th>D<sub>min</sub></th><th>E<sub>min</sub></th><th>D<sub>min</sub></th><th>ΔΔG<sup>ds</sup><sub>(I)-(NRP1)</sub></th><th>D<sub>min</sub></th><th>ΔΔG<sup>ds</sup><sub>(I)-(NRP1)</sub></th><th>D<sub>min</sub></th></tr></thead><tbody><tr><td>PLXA1 homodimer</td><td>-36.6 (1.1)</td><td>0.47</td><td>-29.3 (0.9)</td><td>0.49</td><td>28.9</td><td>0.49</td><td>44.4</td><td>0.49</td></tr><tr><td>NRP1 homodimer</td><td>-19.4 (0.8)</td><td>0.79</td><td>-73.4 (2.2)</td><td>0.64</td><td>47.7</td><td>0.64</td><td>0</td><td>0.64</td></tr><tr><td>NRP1-PLXA1 heterodimer</td><td>-59.2 (1.7)</td><td>0.56</td><td>-67.3 (1.4)</td><td>0.55</td><td>8.0</td><td>0.55</td><td>5.1</td><td>0.55</td></tr></tbody></table>

The PMFs correspond to the separation of the monomers associated in right-handed interactions (R). The PMF profile was calculated in the case of the separation of the PLXA1 monomers associated in left-handed (L) within a DOPC bilayer. E<sub>min</sub> is the difference in free energy between the associated state and the full separation of the monomers. D<sub>min</sub> is the distance separating the center of mass of the backbone beads of the two helices at the minimum of the PMF profile. ΔΔG<sup>ds</sup><sub>(I)-(NRP1)</sub> = ΔΔG<sup>ds</sup><sub>(I)-(NRP1)</sub> is the difference between the apparent dissociation free energy of the PLXA1 homodimer or the heterodimer (I) and the apparent dissociation free energy of the NRP1 homodimer (see equation 2). Energies are given in kJ mol<sup>-1</sup> and distances are in nm.
doi:10.1371/journal.pone.0097779.t003

right-handed and in left-handed interactions. This finding could suggest the possibility of equilibrium between the two forms of packing.

Thus, the PMF profiles calculated for the homodimers and the heterodimer embedded within a DOPC and a POPC bilayer reveal the existence of a hierarchy of interaction strengths between the TM domains. The results indicate that the two NRP1 TM domains exhibit a very strong power of association, slightly stronger than that of the NRP1-PLXA1 TM domains. Conversely, the strength of association of the two PLXA1 TM domains is considerably lower. The scale of interaction energies of the intramembrane domains can be summarized as follows: NRP1-NRP1≥ NRP1-PLXA1>>PLXA1-PLXA1, a hierarchy in agreement with the results of the BACTH two-hybrid assay.

## Discussion

This work focuses on the understanding of the Sema3 receptor complexes involving plexins A and neuropilins. The multiple functions of these secreted semaphorins (Sema3s) depend on signaling through class A plexins. Alone, plexins cannot function as semaphorin receptors and neuropilins are required as co-receptors. Consequently, for signaling through plexin A, class 3 Semaphorins have to interact with neuropilins that in turn have to interact with plexins.

In this paper, we have addressed the contribution of the transmembrane domains of neuropilin 1 and plexin A1 in receptor association using CG-MD simulations in lipid bilayers. Our focus on TM interactions comes from the intriguing presence of single or double GxxxG motifs in almost all known neuropilin partners. Here our observations showing that the PLXA1 TM sequence is also able to dimerize with NRP1 TM suggest that blocking the interaction of NRP1 TM with PLXA1 TM may also be an efficient therapeutic strategy. The simulations presented here provide a unique source of information contributing to the understanding of the action mechanism of such peptidic antagonists. They demonstrate the existence of dynamic TM-TM interactions which determine the hierarchy of the receptor complex formation.

## CG-MD Simulations Predict a Hierarchy of TM Interactions

CG-MD simulations in DOPC and POPC bilayers demonstrated that the NRP1 and PLXA1 TM domains spontaneously assemble to form homodimers and heterodimers. Our results point to two major findings. One is the very high propensity for a right-handed packing of the helices and the second is the low propensity for a left-handed packing mode when at least one PLXA1 TM monomer is involved in association (Figure S1).

We show that the PLXA1 homodimer is allowed to visit the two binding modes at a significant level in the two lipid environments, contrarily to the heterodimer for which the left-handed association is observed at low level. The propensity of the TM helices to associate in right or left handed manner is modulated by the nature of the environment. Also, the strength of helix-helix interactions differ in the two milieus as deduced from PMF calculations. Similar effects of lipids on the self-assembly of peptides have been reported in several other CG-MD studies [52–54]. We have to note here that only one starting structure was used to evaluate these effects. In addition, because of the high computational cost, we did not test the dependence on starting conditions nor the sensitivity of the PMFs on the initial configurations used in each window.
The differences in energy minima of the PMF profiles calculated in DOPC and POPC are around 8 kJ mol⁻¹ or less as evaluated for the NRP1 homodimer. Lipids effects on the power of TM domain association exist but are weak. According to various studies, these differences in energies may not be significant because of statistical and systematic errors that could influence free energies [52,55].

Intriguing is the large destabilization of the right-handed binding mode of the PLXA1 homodimer in POPC (more than 15 kJ mol⁻¹) with respect to the more stable mode of the NRP1 homodimer. This result may suggest that the PLXA1 TM domains are more sensitive to the nature of the environment during the association than the two other dimers. We also remark that there is no relation between the preferred binding mode and the strength of TM interactions. Left-handed interactions are more often observed in DOPC while right-handed interactions are more often observed in POPC, but despite this preference of binding mode, the strength of the interactions is weaker than in DOPC. We speculate that these results are both related to the differences in the sequence of the TMs and to the nature of the lipids. An effect of the bilayer thickness and its fluidity can be hypothesized. In this study we did not focus on the role of the lipids in modulating association of TM domains.

The left-handed packing of the PLXA1 homodimer leads to the weakest free energy of dissociation as estimated from PMF calculations, revealing a very unstable dimerization state. Weak interactions facilitate the complete dissociation of the monomers, as seen in POPC. Within the bilayer the monomers can move separately for more than 5 µs before they re-associate in right-handed interactions for the rest of the simulation. This behavior reveals the existence of a dynamic equilibrium between monomers and dimers. Monomer-dimer equilibrium has been reported in CG simulations of the GpA dimer [56]. In addition, the alternate succession of helices crossed with a positive and a negative angle (or inversely) without dissociation demonstrated the existence of a rapid exchange between the two modes of association.

Association/dissociation events have been also described in the literature, such as for example on WALP peptides and GpA wild and mutant peptides [55,57-59]. Together, the existence of a bimodal packing of the PLXA1 TM helices and of an equilibrium between monomers and dimers suggests the possibility of dynamic exchanges between several other intramembrane domains of the various plexin partners including neuropilins. This is consistent with the demonstration that a peptide mimicking the TM of PlexA1 is as efficient as a peptide mimicking the TM of NRP1 to perturb receptor assembly and subsequent signaling pathway.

An important finding is the possibility of left-handed packing when both the NRP1 and PLXA1 TM domains are embedded in DOPC or in POPC bilayers. As evidenced, these interactions are less frequent than for the PLXA1 homodimer but, similarly to what was found for this homodimer, DOPC seems to favor this type of interaction.

Contrarily to these two dimers, NRP1-NRP1 TMs were never detected in a left-handed binding mode, neither in DOPC (Figure S1) nor in POPC (data not shown) during a total of 200 µs simulations. A unique mode of binding is evidenced. The presence of the double GxxxG motif in the TM sequence is in total coherence with this finding. However, taking in mind the problem of conformational sampling and the lack of binding/unbinding events within the time scale simulated, we cannot exclude other potentially relevant binding modes. Alternative approaches such as 2D PMF along the distance and rotation coordinates could be applied for the search of other binding modes [60]. In the same line, we have performed analysis of the 1D PMF runs in the two

environments along the distance coordinates. No other binding mode than the one described from our CG-MD simulations was distinguished. Once the two helices are far apart, helix-helix reorganization was not observed (Figure S2).

Thus, the present results tend to demonstrate that the PLXA1 TM domain encodes the properties for both left-handed and right-handed binding mode and the NRP1 TM domain encodes the properties for right-handed binding mode.

PMF analysis gives only approximate magnitudes on the strength of TM association but makes it possible to predict a hierarchy of self-association for the TM monomers. We estimate that the NRP1 TM domain displays a very strong power for homo dimerization and hetero dimerization with the PLXA1 TM domain. The strong power of interactions involving at least one NRP1 TM domain reflects specific interactions which could drive dimerization of the full length protein. Neuropilin 1 dimers preexist at the cell surface in the absence of ligand before the formation of complexes [26,45]. Contrarily, the PLXA1 TM domain displays a very weak power for homo association and this could be the result of unspecific interactions, explaining the modulator role of this receptor in the formation of neuropilin 1 receptor complexes.

## GxxxG Motifs Drive TM Interactions

Dimer structures are well defined and the same interacting helix faces were identified in DOPC and POPC bilayers. Each helical TM exhibits its own face for dimerization and we demonstrate that the TM sequences encode an intrinsic propensity to self-associate in the presence of the GxxxG motifs. For NRP1 the key interacting residues constitute the M₁₂xxxG₁₆xxxG₂₀ motif as a guide for homo and hetero dimerization. The C-terminal Gly₂₄ of the double motif (G₁₆xxxG₂₀xxxG₂₄) is not involved in tight packing albeit present at the interface. The Gly rich sequence G₉IGGGGG₁₅ of PLXA1 is also involved both in homo and hetero dimerization and, together with surrounding residues, gives rise to two distinct association properties. With the small Ala₆ residue at the N-terminus, the A₆xxG₉xxxG₁₃ motif confers the specificity of left-handed interactions for the homodimer. With Leu₁₉ at the C-terminus, the G₁₁xxxG₁₅xxxL₁₉ motif confers the specificity of right-handed interactions for the homodimer and the heterodimer.

The role of the sequence context in stabilizing helix-helix interactions mediated by a GxxxG motif is here very well illustrated. According to the general principles for transmembrane stability, the GxxxG motif is necessary but not sufficient to achieve strong TM helix-helix association and the sequence context modulates the stability of the TM helix-helix packing [61].

In the case of NRP1, Met₁₂ is aligned on the same helix face as G₁₆, G₂₀ and G₂₄ and with the flanking Ser residue contributes to dimer stabilization through polar and Van der Waals interactions. In the case of PLXA1, Ala₆ is a weak stabilizing residue, but with Ile₇ it contributes to Van der Waals stabilization of the left-handed structure. Leu₁₉, with the proximal residues, contributes to Van der Waals stabilization of the right-handed structure. These are the factors identified as involved in controlling the stability of a GxxxG stabilized TM helix dimer [8].

## "Versatile" PLXA1 TM Domain

Interacting key residues responsible for the specificity of left-handed and right-handed interactions of the PLXA1 TM homodimer are distributed on the two opposed faces of the intramembrane helix. This major finding is clearly illustrated in figure 6. This result is of considerable importance because the two helix faces can serve for distinct interactions leading to a large
diversity of partners. These properties are very probably useful in
explaining the versatile character of semaphorin - neuropilin -
plexin complexes known to be regulators of tumor progression and
tumor angiogenesis [22]. The dynamic equilibrium between the
two modes of helix packing governed by the presence of distinct
motifs could induce the possibility of switching from one mode of
packing to another one depending on the partner. This hypothesis
is supported by a similar behavior reported for the ErbB receptors
[62-65].

Within the TM domain of most human ErbB receptor tyrosine
kinases, two GxxxG like motifs are conserved and interactions
mediated by one motif guide the receptor in active state, whereas
interactions mediated by the second motif guide the receptor in
inactive state. Thus, the ErbB receptor function could depend on
the TM dimer alternating in a switch-like fashion between two
structures, stabilized by either of the two GxxxG-like motifs.
Moreover, some residues involved in ErbB2 homodimerization are
also involved in ErbB1/ErbB2 heterodimer formation. The core
of dimerization is conserved and the surrounding residues which
contribute to the stability of the structures of the ErbB2
homodimer or the ErbB1/ErbB2 heterodimer differ [66]. Thus,
the properties of the ErbB TM domains parallel those described
here for the PLXA1 TM domain. The core of dimerization
appears to be the G₁₁G₁₂G₁₃ sequence common to the two motifs
A₆ₓₓG₉ₓₓG₁₃ and G₁₁ₓₓG₁₅ₓₓL₁₉ which allows rapid switch-
ing between the two modes of association. The first motif is
responsible for homodimerization and guides left-handed interac-
tions and the second motif is responsible for homodimerization
and heterodimerization and guides right-handed interactions.

Conclusion

Toward the Design of Therapeutic Agents

As a conclusion, CG-MD simulations enable the description of
models for TM interactions of Sema3 co-receptors. To the best of
our knowledge, this work is the first devoted to investigations on
the stability and specificity of TM interactions of the Sema3 co-
receptor in membrane bilayers.

The success of coarse-grained models has been proven in many
studies on membrane protein systems [67] and thus gives us
confidence in the biological relevance of our results. Ongoing
research aims at exploring these complex systems in more
sophisticated bilayers including mixed lipids and cholesterol to
corroborate our models.

<figure><img src="image_8.png" /><figcaption>**Figure 6.** Position of the key residues involved in PLXA1 TM homodimerization. R denotes the helix face conferring a right-handed packing and L denotes the helix face conferring a left-handed packing.</figcaption></figure>

doi:10.1371/journal.pone.0097779.g006

The GxxxG motifs play critical roles in transmembrane
association and it can be expected that TM interactions will
influence the association of full-length Sema3 receptors. The
present data agree with the first experimental studies on
neuropilins 1 [21,45,46] and open the route for additional studies
to explore experimentally the role of the TM domains of class A
Plexins and neuropilin members that can now be extended to
receptor tyrosine kinases to have a full picture of the Semaphorin
signaling platform [68].

Many types of tumor cells express semaphorin receptors. These
transmembrane receptors represent major potential advances
toward the design of a new class of efficient therapeutic agents
for cancer treatment [10,21]. The therapeutic potential of the
NRP1 TM domain has been already demonstrated [45,46] and we
anticipate that the properties of association of the intramembrane
domains of plexin A1 and neuropilin 1 reported here will
contribute to therapeutic perspectives. The computational rational
design of hydrophobic peptides able to specifically bind to TM
domains of receptors, and alter their functions, has already been
successfully employed in the case of integrins [69]. Our present
work should pave the way for the design of such interfering
peptides, defining a novel class of therapeutic agents.

**Materials and Methods**

Estimation of TM Domains Dimerization Propensity
(BACTH Method)

A bacterial two-hybrid system based on the recombination of
adenylate cyclase CyaA from *Bordetella pertussis* was used to
measure both homo- and heterodimerization propensities of
transmembrane domains of interest. This system relies on
reconstitution of the catalytic domain of CyaA which can be
separated into two complementary fragments [70]. When each
fragment is fused to a protein of interest, a functional adenylate
cyclase can be reassembled upon interaction of the two proteins,
which is followed by the production of cyclic AMP in an *E. coli*
strain lacking its own adenylate cyclase. The synthesis of cyclic
AMP is proportional to the association propensity of the two
fragments, and can be measured in a number of ways [70]. This
system has been used widely and is readily amenable for study of
membrane proteins interactions. Very briefly, we have modified
the BACTH pKTN25 and pUT18 plasmids so that they encode
for hybrid proteins containing an OmpA signal sequence followed
by the different TM domains of interest and the T25 and T18
fragments of adenylate cyclase (Sawma et al., manuscript submit-
ted). The following transmembrane sequences were used:
ILITIIAMSALGVLLGAVCGVVLYRKR for NRP1,
LLTLPAIVGIGGGGLLLLLIVAVLYAYRRK for PLXA1
and IEITLIIFGVMAGVIGTILLISYGIRRK for the GpA con-
trol. These sequences contain a C-terminal tribasic sequence to
ensure their correct insertion, as was previously done for studies
with synthetic peptides [45,46]. Double transformation in
BTH101 (cya') *E. coli* cells, cell growth and induction, colorimetric
assay on MacConkey reporter plates and beta-Galactosidase assay
in 96-well arrays on a TECAN machine were performed as
described [70].

**Molecular Dynamics System Setup**

CG-MD simulations were performed to investigate the homo
and hetero association of NRP1 and PLXA1 TM domains using
the GROMACS software package (version 4.5.4) [71,72]. The
TM peptides, the lipids and the water CG particles were described
with the MARTINI force-field version 2.1 [49,50] extensively used
for the study of various peptide-lipid membrane systems [54].
Images were created using the VMD software [73] (http://www ks.uiuc.edu/Research/vmd/).

## TM Peptide Sequences

The two TM sequences used in the simulations correspond to the native potential sequences as found in the Uniprot database (http://www.uniprot.org/) (Table 1). We adjusted the length of the TM domains by adding a few residues on each side, in order to design a 28 amino acid sequence. For the sake of clarity, the TM sequences have been re-numbered 1 to 28 throughout this paper.

The TM domain of NRP1 is characterized by two successive GxxxG motifs $G_{16}xxxG_{20}xxxG_{24}$. The PLXA1 TM sequence is characterized by the presence of 6 Gly residues, 5 being consecutive in the sequence constituting the $G_{11}GGGG_{15}$ motif. Another GxxxG like motif, $G_9IGGG_{13}$, is observed overlapping the previous one. Considering the α helical structure of the TM domains, the 6 Gly residues of the PLXA1 TM domain are distributed over all the helix faces, covering almost two helix turns, while the 3 Gly residues in the NRP1 TM domain are aligned on only one helix face.

## Bilayer Preparation

We prepared large CG DOPC (1,2-di-oleoyl-*sn*-glycero-3-phosphocholine) and POPC (palmitoyl-2-oleoyl-*sn*-glycero-3-phosphocholine) bilayers for TM domain self-assembly studies. The pre-equilibrated DOPC bilayer taken from the MARTINI web site (http://md.chem.rug.nl/cgmartini/) was replicated onto a 2×2 grid in the bilayer plane (X, Y directions). The POPC bilayer was built from the original pre-equilibrated POPE bilayer. The POPE bilayer was replicated onto a 2×2 grid in the bilayer plane as for the DOPC bilayer. After minimization, 1 µs of equilibration was performed using the conditions given in the following. The POPC bilayer was obtained by replacing the first CG particle (Qd NH3 in POPE was replaced by Qo NC3 in POPC).

A pure DOPC bilayer containing 492 lipids and 6000 water particles, and a pure POPC bilayer composed of 451 lipids and 5844 water particles were obtained. After minimization, each bilayer was equilibrated for 1 µs using the following conditions. The temperature was coupled to a thermostat at 300 K using a Berendsen algorithm [74] with a time constant of 1 ps applied separately to the lipids and the water groups. The pressure was coupled using a semi-isotropic scheme (coupling time 0.5 ps, compressibility $4.5 \times 10^{-5} \text{ bar}^{-1}$) in which the lateral and perpendicular pressures are coupled independently at 1 bar. A time step of 0.025 ps was used. After equilibration the X, Y, Z dimensions of the bilayer were 13.61, 12.1, 8.73 nm and 11.96, 11.94, 9.22 nm for the DOPC bilayer and the POPC bilayer respectively. The Z axis represents the normal to the bilayer plane (X, Y).

## System Preparation

The CG structures of the NRP1 and the PLXA1 TM domains were modeled as ideal α helices from the atomistic structure initially built using AMBER [75,76]. The C-terminus was capped with a negative charge, and the N-terminus was capped with a positive charge as given in the MARTINI force field. The aspartic acid residue included in the NRP1 TM sequence was considered to be uncharged. Thus, the two TM domains were electrically neutral. Both the NRP1 and PLXA1 monomers were aligned along the helical axis before applying the insertion procedure.

Two monomeric peptides were inserted in a parallel manner in the equilibrated bilayer (DOPC/POPC) at a distance of 6 nm from each other. Prior to insertion, the TM peptides were aligned along the normal to the bilayer and were randomly rotated about the helix axis. This leads to a random orientation of GxxxG motifs

of one monomer relative to the second. The resulting initial helix-helix interface is thus totally arbitrary. A small number of lipids and water molecules (a few tens depending on the systems) were removed from the sites of peptide insertion.

The system was energy minimized for 10000 steps using the steep integrator without peptide particle restraint. MD simulations were then performed at 300 K with a coupling time of 1 ps applied separately to the protein, the lipids (DOPC/POPC) and the water groups using a v-rescale algorithm. The pressure was coupled using a semi-isotropic scheme, in which the lateral and perpendicular pressures were coupled independently at 1 bar with a coupling time of 5 ps and a compressibility of $3 \times 10^{-4} \text{ bar}^{-1}$. A time step of 0.040 ps was used for the production step. For all the simulations (pure bilayers and peptide-lipid bilayers) Lennard Jones and Coulombic interactions were shifted to zero between 0.9 and 1.2 nm and 0.0 and 1.2 nm respectively.

Six systems were studied. Three systems contained two peptides embedded in a DOPC bilayer and three systems contained the same two peptides embedded in a POPC bilayer, the two peptides being identical (homodimer) or different (heterodimer). The simulations were repeated several times, varying the initial orientation of the two monomers relative to each other. Four distinct simulations were performed for self-association studies of the NRP1 TM domains, three in the case of the PLXA1 TM domains, and six in the case of the heterodimer, for each type of bilayer. Simulations of 15 µs or longer were performed depending upon the dimer type, but sufficiently long to observe the assembly of the two peptides and the stability of association for at least 10 µs. A total of 13 simulations were performed in each lipid bilayer, resulting in a total simulation time of 613 µs (Table 2).

## System Analysis

Tools provided in the Gromacs package were used for the analyses. The associated state of the dimer was defined from the analysis of inter-helix distance (distance between the center of mass of the helix backbones, COM) and of inter-helix contact distance (minimal distance between the helix backbones). The two helices were considered associated when the closest contact distance was less than 0.5 nm.

The RMSD (Root Mean Square Deviation) matrix, which compares all the conformations generated during the simulation, was calculated to identify the periods of the simulation during which the dimer structure is stabilized. A dimer is stable when the RMSD value (calculated over the backbone beads) is less than 0.3 nm. The longest period of the dimer stabilization was chosen to calculate the contact matrix that represents the smallest distances between the backbone beads of the residue pairs, on average. The helix-helix interface is then extracted from this contact matrix.

For the NRP1 homodimer and the NRP1-PLXA1 heterodimer, the simulations led generally to only one major state for the association of the helices. For the PLXA1 homodimer two distinct states were observed along each of the six simulations. As detailed in the following, the two helices were successively crossed with a positive or a negative crossing angle (see results). Then, for the analysis, the trajectory portions corresponding to the positive mode of association were merged and similarly, the trajectory portions corresponding to the negative mode of association were merged.

## PMF via Umbrella Sampling

To compare the strength of interaction of the NRP1 homodimer, the PLXA1 homodimer and the NRP1-PLXA1 heterodimer,
we calculated the potential of mean force (PMF) following the procedure described by Marrink and coll [54].

PMFs were calculated by carrying out a series of umbrella sampling simulations, before unbiasing with the weighted histogram analysis method (WHAM) [77]. To conduct umbrella sampling, we generated a series of configurations along the reaction coordinate, $ζ$, chosen as the distance between the centers of mass of the backbone beads of the two helices. A pulling simulation was performed over the course of 420 ns in the X and Y dimensions (membrane plane). Starting from two associated helices, taken from the self-assembly simulations, the two helices were dissociated by restraining the backbone beads of one helix with a force constant of 1000 kJ mol$^{-1}$ nm$^{-2}$ and applying constant velocity pulling (pull rate = 0.01 nm.ps$^{-1}$). For each system, 22 configurations, corresponding to a 0.1 nm shift of the monomer, were chosen as starting configurations for the umbrella sampling windows that were run independently. Each configuration was equilibrated for 200 ns followed by 2 µs production. The distances between the centers of mass of the backbone beads of the two helices were sampled from 0.5 nm to 2.6 nm for the NRP1 homodimer and the heterodimer, and from 0.4 nm to 2.5 nm for the PLXA1 homodimer. Seven PMF profiles were calculated corresponding to 310 µs of total simulation time.

PMF extraction using WHAM yielded the $ΔG$ for the binding/unbinding process. By integrating the PMF profile, we calculated the association constant $K_A$ which is defined in cylindrical coordinates by the following expression [52,54].

$$K_A = C \int_{0}^{\zeta_{\max}} 2\pi\zeta e^{-\beta G(\zeta)} d\zeta \quad (1)$$

$ζ_{max}$ represents the cylindrical radius separating the associated and dissociated states of the two TM helices defined to be 2.5 nm. $C = 1$ for the heterodimer and $C = 1/2$ for the two homodimers. From the association constant, we determine the apparent dissociation free energy given by $ΔG^{dis} = -RT\ln K_A$. Since a direct comparison between the calculated value and quantitative experimental data is difficult, we determined $ΔΔG^{dis}$ given by the following expression:

$$\Delta\Delta G^{dis}(i)-(NRP1) = (\Delta G^{dis}(i) - \Delta G^{dis}(NRP1)) = -RT \ln \left( \frac{K_{A(i)}}{K_{A(NRP1)}} \right) \quad (2)$$

With $RT = 2.4943$ kJ mol$^{-1}$

## References

1. Cebecauer M, Spitaler M, Serge A, Magee AI (2010) Signalling complexes and clusters: functional advantages and methodological hurdles. Journal of Cell Science 123: 309-320.

2. Lemmon MA, Schlessinger J (2010) Cell Signaling by Receptor Tyrosine Kinases. Cell 141: 1117-1134.

3. Call ME, Wucherpfennig KW (2007) Common themes in the assembly and architecture of activating immune receptors. Nature Reviews Immunology 7: 841-850.

4. Campbell ID, Humphries MJ (2011) Integrin Structure, Activation, and Interactions. Cold Spring Harbor Perspectives in Biology 3.

5. Audet M, Bouvier M (2012) Restructuring G-Protein-Coupled Receptor Activation. Cell 151: 14-23.

6. Beevers AJ, Dixon AM (2010) Helical membrane peptides to modulate cell function. Chem Soc Rev 39: 2146-2157.

7. Bordag N, Keller S (2010) Alpha-helical transmembrane peptides: a “divide and conquer” approach to membrane proteins. Chem Phys Lipids 163: 1-26.

8. Cymer F, Veerappan A, Schneider D (2011) Transmembrane helix-helix interactions are modulated by the sequence context and by lipid bilayer properties. Biochim Biophys Acta 1818: 963-973.

$ΔG^{dis}(i)$ is the apparent dissociation free energy of the PLXA1 homodimer or the NRP1-PLXA1 heterodimer.

$ΔG^{dis}(NRP1)$ is the apparent dissociation free energy of the NRP1 homodimer taken as reference.

## Supporting Information

**Figure S1 Crossing angle as a function of simulation time.** Time evolution of the crossing angle of the helices for the three dimers NRP1-NRP1, NRP1-PLXA1, PLXA1-PLXA1 in DOPC. Several independent CG-MD simulations were performed for each dimer. In the case of NRP1-NRP1 only the right-handed binding mode is observed (negative value of the crossing angle). The multiple independent simulations arrived at the same dimer structure. The right-handed binding mode is preferred for the heterodimer, but left-handed transitions are observed. The PLXA1-PLXA1 homodimer evidences frequent right and left-handed transitions.

(TIF)

**Figure S2 The motions of one helix relatively to the other.** Analysis of the 1D PMF runs for the dimers NRP1-NRP1, NRP1-PLXA1 in the right-handed binding mode and PLXA1-PLXA1 in right and left-handed binding modes. The graphs represent the X Y coordinates of the center of mass of each helix during the Umbrella Sampling calculations in DOPC (22 windows). Helix H1 is restrained at the center of the lipid box and H2 is allowed to move in the XY plane for each value of the reaction coordinate $ζ$. The NRP1-NRP1 and NRP1-PLXA1 dimers dissociate without reorganization. For the right-handed PLXA1 homodimer, a reorganization of the helices is suggested at the first steps of the dissociation process.

(TIF)

## Acknowledgments

The technical assistance of A. Boyer is greatly appreciated. Many thanks to Dr. E. Rowley-Jolivet for providing language help.

## Author Contributions

Conceived and designed the experiments: PH JS DB MG. Performed the experiments: SAS PS IJ MG NG. Wrote the paper: SAS MG NG.

9. Fink A, Sal-Man N, Gerber D, Shai Y (2011) Transmembrane domains interactions within the membrane milieu: Principles, advances and challenges. Biochim Biophys Acta 1818: 974-983.

10. Hubert P, Sawma P, Duneau JP, Khao J, Henin J, et al. (2010) Single-spanning transmembrane domains in cell growth and cell-cell interactions: More than meets the eye? Cell Adh Migr 4: 313-324.

11. Li E, Wimley WC, Hristova K (2012) Transmembrane helix dimerization: beyond the search for sequence motifs. Biochim Biophys Acta 1818: 183-193.

12. Matthews EE, Zoonens M, Engelman DM (2006) Dynamic helix interactions in transmembrane signaling. Cell 127: 447-450.

13. Lemmon MA, Flanagan JM, Hunt JF, Adair BD, Bormann BJ, et al. (1992) Glycophorin A dimerization is driven by specific interactions between transmembrane alpha-helices. J Biol Chem 267: 7683-7689.

14. Mackenzie KR (2006) Folding and stability of alpha-helical integral membrane proteins. Chem Rev 106: 1931-1977.

15. MacKenzie KR, Engelman DM (1998) Structure-based prediction of the stability of transmembrane helix-helix interactions: the sequence dependence of glycophorin A dimerization. Proc Natl Acad Sci U S A 95: 3583-3590.

16. MacKenzie KR, Prestegard JH, Engelman DM (1997) A transmembrane helix dimer: structure and implications. Science 276: 131-133.
17. Moore DT, Berger BW, DeGrado WF (2008) Protein-protein interactions in the membrane: sequence, structural, and biological motifs. Structure 16: 991-1001.

18. Russ WP, Engelman DM (2000) The GxxxG motif: a framework for transmembrane helix-helix association. J Mol Biol 296: 911-919.

19. Arkhipov A, Shan YB, Das R, Endres NF, Eastwood MP, et al. (2013) Architecture and Membrane Interactions of the EGF Receptor. Cell 152: 557-569.

20. Endres NF, Das R, Smith AW, Arkhipov A, Kovacs E, et al. (2013) Conformational Coupling across the Plasma Membrane in Activation of the EGF Receptor. Cell 152: 543-556.

21. Roth L, Koncina E, Satkauskas S, Cremel G, Aunis D, et al. (2009) The many faces of semaphorins: from development to pathology. Cell Mol Life Sci 66: 649-666.

22. Neufeld G, Kessler O (2008) The semaphorins: versatile regulators of tumour progression and tumour angiogenesis. Nat Rev Cancer 8: 632-645.

23. Neufeld G, Sabag AD, Rabinovicz N, Kessler O (2012) Semaphorins in angiogenesis and tumor progression. Cold Spring Harb Perspect Med 2: a006718.

24. Soker S, Takashima S, Miao HQ, Neufeld G, Klagsbrun M (1998) Neuropilin-1 is expressed by endothelial and tumor cells as an isoform-specific receptor for vascular endothelial growth factor. Cell 92: 735-745.

25. Chen H, Chedotal A, He Z, Goodman CS, Tessier-Lavigne M (1997) Neuropilin-2, a novel member of the neuropilin family, is a high affinity receptor for the semaphorins Sema E and Sema IV but not Sema III. Neuron 19: 547-559.

26. Takahashi T, Nakamura F, Jin Z, Kalb RG, Strittmatter SM (1998) Semaphorins A and E act as antagonists of neuropilin-1 and agonists of neuropilin-2 receptors. Nat Neurosci 1: 487-493.

27. Feiner L, Koppel AM, Kobayashi H, Raper JA (1997) Secreted chick semaphorins bind recombinant neuropilin with similar affinities but bind different subsets of neurons in situ. Neuron 19: 539-545.

28. Kolodkin AL, Levengood DV, Rowe EG, Tai YT, Giger RJ, et al. (1997) Neuropilin is a semaphorin III receptor. Cell 90: 753-762.

29. Nakamura F, Tanaka M, Takahashi T, Kalb RG, Strittmatter SM (1998) Neuropilin-1 extracellular domains mediate semaphorin D/III-induced growth cone collapse. Neuron 21: 1093-1100.

30. Sharma A, Verhaagen J, Harvey AR (2012) Receptor complexes for each of the Class 3 Semaphorins. Front Cell Neurosci 6: 28.

31. Rohm B, Ottemeyer A, Lohrum M, Puschel AW (2000) Plexin/neuropilin complexes mediate repulsion by the axonal guidance signal semaphorin 3A. Mech Dev 93: 95-104.

32. Takahashi T, Fournier A, Nakamura F, Wang LH, Murakami Y, et al. (1999) Plexin-neuropilin-1 complexes form functional semaphorin-3A receptors. Cell 99: 59-69.

33. Tamagnone L, Artigiani S, Chen H, He Z, Ming GI, et al. (1999) Plexins are a large family of receptors for transmembrane, secreted, and GPI-anchored semaphorins in vertebrates. Cell 99: 71-80.

34. Antipenko A, Himanen JP, van Leyen K, Nardi-Dei V, Lesniak J, et al. (2003) Structure of the semaphorin-3A receptor binding module. Neuron 39: 589-598.

35. Love CA, Harlos K, Mavaddat N, Davis SJ, Stuart DI, et al. (2003) The ligand-binding face of the semaphorins revealed by the high-resolution crystal structure of SEMA4D. Nat Struct Biol 10: 843-848.

36. Janssen BJ, Robinson RA, Perez-Branguli F, Bell CH, Mitchell KJ, et al. (2010) Structural basis of semaphorin-plexin signalling. Nature 467: 1118-1122.

37. Klostermann A, Lohrum M, Adams RH, Puschel AW (1998) The chemorepulsive activity of the axonal guidance signal semaphorin D requires dimerization. J Biol Chem 273: 7326-7331.

38. Liu H, Juo ZS, Shim AH, Focia PJ, Chen X, et al. (2010) Structural basis of semaphorin-plexin recognition and viral mimicry from Sema7A and A39R complexes with PlexinC1. Cell 142: 749-761.

39. Nogi T, Yasui N, Mihara E, Matsunaga Y, Noda M, et al. (2010) Structural basis for semaphorin signalling through the plexin receptor. Nature 467: 1123-1127.

40. Tong Y, Chugha P, Hota PK, Alviani RS, Li M, et al. (2007) Binding of Rac1, Rnd1, and RhoD to a novel Rho GTPase interaction motif destabilizes dimerization of the plexin-B1 effector domain. J Biol Chem 282: 37215-37224.

41. Oinuma I, Katoh H, Negishi M (2004) Molecular dissection of the semaphorin 4D receptor plexin-B1-stimulated R-Ras GTPase-activating protein activity and neurite remodeling in hippocampal neurons. J Neurosci 24: 11473-11480.

42. Takahashi T, Strittmatter SM (2001) PlexinA1 autoinhibition by the plexin sema domain. Neuron 29: 429-439.

43. Janssen BJ, Malinauskas T, Weir GA, Cader MZ, Siebold C, et al. (2012) Neuropilins lock secreted semaphorins onto plexins in a ternary signaling complex. Nat Struct Mol Biol 19: 1293-1299.

44. Tamagnone L, Comoglio PM (2000) Signalling by semaphorin receptors: cell guidance and beyond. Trends Cell Biol 10: 377-383.

45. Roth L, Nasarre C, Dirrig-Grosch S, Aunis D, Cremel G, et al. (2008) Transmembrane domain interactions control biological functions of neuropilin-1. Mol Biol Cell 19: 646-654.

46. Nasarre C, Roth M, Jacob L, Roth L, Koncina E, et al. (2010) Peptide-based interference of the transmembrane domain of neuropilin-1 inhibits glioma growth in vivo. Oncogene 29: 2381-2392.

47. Chaudhary B, Khaled YS, Ammori BJ, Elkord E (2013) Neuropilin 1: function and therapeutic potential in cancer. Cancer Immunol Immunother.

48. Marrink SJ, de Vries AH, Mark AE (2004) Coarse grained model for semi-quantitative lipid simulations.. 108: 750-760.

49. Marrink SJ, Risselada HJ, Yefimov S, Tieleman DP, de Vries AH (2007) The MARTINI force field: coarse grained model for biomolecular simulations. J Phys Chem B 111: 7812-7824.

50. Monticelli L, Kandasamy SK, Periole X, Larson RG, Tieleman DP, et al. (2008) The MARTINI coarse grained forcefield: extension to proteins. J Chem Theory and Comput 4: 819-834.

51. Hall BA, Chetwynd AP, Sansom MS (2011) Exploring peptide-membrane interactions with coarse-grained MD simulations. Biophys J 100: 1940-1948.

52. Henin J, Pohorille A, Chipot C (2005) Insights into the recognition and association of transmembrane alpha-helices. The free energy of alpha-helix dimerization in glycophorin A. J Am Chem Soc 127: 8478-8484.

53. Janosi L, Prakash A, Doxastakis M (2010) Lipid-modulated sequence-specific association of glycophorin A in membranes. Biophys J 99: 284-292.

54. Sengupta D, Marrink SJ (2010) Lipid-mediated interactions tune the association of glycophorin A helix and its disruptive mutants in membranes. Phys Chem Chem Phys 12: 12987-12996.

55. Psachoulia E, Fowler PW, Bond PJ, Sansom MS (2008) Helix-helix interactions in membrane proteins: coarse-grained simulations of glycophorin A a helix dimerization. Biochemistry 47: 10503-10512.

56. Bond PJ, Sansom MS (2006) Insertion and assembly of membrane proteins via simulation. J Am Chem Soc 128: 2697-2704.

57. Domanski J, Marrink SJ, Schafer LV (2012) Transmembrane helices can induce domain formation in crowded model membranes. Biochim Biophys Acta 1818: 984-994.

58. Psachoulia E, Marshall DP, Sansom MS (2010) Molecular dynamics simulations of the dimerization of transmembrane alpha-helices. Acc Chem Res 43: 388-396.

59. Schafer LV, de Jong DH, Holt A, Rzepiela AJ, de Vries AH, et al. (2011) Lipid packing drives the segregation of transmembrane helices into disordered lipid domains in model membranes. Proc Natl Acad Sci U S A 108: 1343-1348.

60. Periole X, Knepp AM, Sakmar TP, Marrink SJ, Huber T (2012) Structural Determinants of the Supramolecular Organization of G Protein-Coupled Receptors in Bilayers. Journal of the American Chemical Society 134: 10959-10965.

61. Doura AK, Kobus FJ, Dubrovsky L, Hibbard E, Fleming KG (2004) Sequence context modulates the stability of a GxxxG-mediated transmembrane helix-helix dimer. J Mol Biol 341: 991-998.

62. Cymer F, Schneider D (2010) Transmembrane helix-helix interactions involved in ErbB receptor signaling. Cell Adh Migr 4: 299-312.

63. Escher C, Cymer F, Schneider D (2009) Two GxxxG-like motifs facilitate promiscuous interactions of the human ErbB transmembrane domains. J Mol Biol 389: 10-16.

64. Fleishman SJ, Schlessinger J, Ben-Tal N (2002) A putative molecular-activation switch in the transmembrane domain of erbB2. Proc Natl Acad Sci U S A 99: 15937-15940.

65. Mendrola JM, Berger MB, King MC, Lemmon MA (2002) The single transmembrane domains of ErbB receptors self-associate in cell membranes. J Biol Chem 277: 4704-4712.

66. Mineev KS, Bocharov EV, Pustovalova YE, Bocharova OV, Chupin VV, et al. (2010) Spatial structure of the transmembrane domain heterodimer of ErbB1 and ErbB2 receptor tyrosine kinases. J Mol Biol 400: 231-243.

67. Periole X, Marrink SJ (2013) The Martini coarse-grained force field. Methods Mol Biol 924: 533-565.

68. Zachary IC (2011) How neuropilin-1 regulates receptor tyrosine kinase signalling: the knowns and known unknowns. Biochemical Society Transactions 39: 1583-1591.

69. Yin H, Slusky JS, Berger BW, Walters RS, Vilaire G, et al. (2007) Computational design of peptides that target transmembrane helices. Science 315: 1817-1822.

70. Battesti A, Bouveret E (2012) The bacterial two-hybrid system based on adenylate cyclase reconstitution in Escherichia coli. Methods 58: 325-334.

71. Hess B, Kutzner C, van der Spoel D, Lindahl E (2008) GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. J Chem Theory Comp 4: 435-447.

72. Van der Spoel D, Lindahl B, Hess G, Groenhof G, Mark AE, et al. (2005) GROMACS: fast, flexible, and free. J Comput Chem 26: 1701-1718.

73. Humphrey W, Dalke A, Schulten K (1996) VMD - Visual Molecular Dynamics. J Molec Graphics 14: 33-38.

74. Berendsen HJC, Postma JPM, van Gunsteren WF, DiNola A, Haak JR (1984) Molecular dynamics with coupling to an external bath. J Chem Phys 81: 3684-3690.

75. Case DA, Cheatham TE 3rd, Darden T, Gohlke H, Luo R, et al. (2005) The Amber biomolecular simulation programs. J Comput Chem 26: 1668-1688.

76. Ponder JW, Case DA (2003) Force fields for protein simulations. Adv Protein Chem 66: 27-85.

77. Kumar S, Rosenberg JM, Bouzida D, Swendsen RH, Kollman PA (1992) THE weighted histogram analysis method for free-energy calculations on biomolecules. I. The method. J Comput Chem 13: 1011-1021.