Conformational rearrangements of peroxysome proliferator activated receptor (PPARgamma) ligand-binding domain (LBD) that accompany the release and binding of ligands are not well understood. To determine the major events associated with the escape of the partial agonist GW0072, molecular dynamic (MD) simulations were performed using two different methods: reversed targeted molecular dynamics (TMD(-1)) and time-dependent distance restraints (TDR) using as restraints either the root mean square deviation from a reference structure (TMD(-1)) or the distance between the geometrical centers of the binding pocket and of the ligand (TDR). Both methods do not assume any a priori route for ligand extraction. To avoid artifacts, different initial simulation conditions were used and particular attention was paid for giving time to the protein to relax during the extraction process by running 10-12 ns simulations within explicit water. Two distinct exit gates A and B were found, independently of initial conditions and method. During the exit process no interaction between GW0072 and the transactivation AF-2 helix was observed. Our results suggest that the ligand uses the intrinsic flexibility of the protein to move within the receptor. Paths A and B are very similar to those found for other nuclear receptors, suggesting that these routes are a common characteristics of nuclear receptors that are used by different kinds of ligands. Finally, the knowledge of entry/exit pathways of a receptor should be very useful in discriminating between different ligands that could have been favorably docked in the binding pocket by introducing docking along these pathways into computational drug design protocols.

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.