Insights into oncogenic mutations of plexin-B1 based on the solution structure of the Rho GTPase binding domain

Local Structures in Proteins from Microsecond Molecular Dynamics Simulations: 2. The Role of Symmetry in GTPase Binding and Dimer Formation

The Rho GTPase binding domain of Plexin-B1 (RBD) prevails in solution as dimer. Under appropriate circumstances, it binds the small GTPase Rac1 to yield the complex RBD-Rac1. Here, we study RBD dimerization and complex formation from a symmetry-based perspective using data derived from 1 μs long MD simulations. The quantities investigated are the local potentials, u(MD), prevailing at the N-H sites of the protein. These potentials are statistical in character providing an empirical description of the local structure. To establish more methodical description, a method for approximating them by explicit functions, u(simulated), was developed in the preceding article in this journal issue. These functions are combinations of analytical Wigner functions, DL,K, belonging to the D2h point group. The D2h subgroups Ag and B2u are found to dominate u(simulated); the B1u subgroup contributes in some cases. The Ag (B2u) functions have axial or rhombic symmetry. For the first time, local potentials in proteins can be quantitatively characterized in terms of their strength (rhombicity) evaluated by axial Ag (rhombic Ag and B2u) contributions. Until now, the chain-segment [β3-L3-β4] and to some extent the α2-helix have been associated with GTPase binding. Here, we find that this process causes an increase (decrease) in the potential strength of β3 and β4 (the preceding L2 loop and the remote chain-segment [(α2-helix)-(α2/β5-turn)-(β5-strand)]), suggesting effects of counterbalancing and allostery. There is evidence for the L2 loop being associated with RBD-GTPase binding. Until now only the L4 loop has been associated with RBD dimerization. The latter process is found to cause an increase (decrease) in the potential strength and rhombicity of the L4 loop (the adjacent chain-segment [(α2-helix)-(α2/β5-turn)-(β5-strand)]), suggesting counterbalancing activity. On average, the RBD dimer features stronger local potentials than RBD-Rac1. The novel information inherent in these findings is mesoscopic in character. Prospects of interest include exploring relation to atomistic force-field parameters.

Plexin-mediated neuronal development and neuroinflammatory responses in the nervous system

Plexins are a large family of single-pass transmembrane proteins that mediate semaphorin signaling in multiple systems. Plexins were originally characterized for their role modulating cytoskeletal activity to regulate axon guidance during nervous system development. Thereafter, different semaphorin-plexin complexes were identified in the nervous system that have diverse functions in neurons, astrocytes, glia, oligodendrocytes, and brain derived-tumor cells, providing unexpected but meaningful insights into the biological activities of this protein family. Here, we review the overall structure and relevant downstream signaling cascades of plexins. We consider the current knowledge regarding the function of semaphorin-plexin cascades in the nervous system, including the most recent data regarding their roles in neuronal development, neuroinflammation, and glioma.

Microsecond MD Simulations of the Plexin-B1 RBD: N-H Probability Density as Descriptor of Structural Dynamics, Dimerization-Related Conformational Entropy, and Transient Dimer Asymmetry

Amide-bond equilibrium probability density, P_eq = exp(-u) (u, local potential), and associated conformational entropy, S_k = -∫P_eq (ln P_eq) dΩ ─ln ∫dΩ, are derived for the Rho GTPase binding domain of Plexin-B1 (RBD) as monomer and dimer from 1 μs MD simulations. The objective is to elucidate the effect of dimerization on the dynamic structure of the RBD. Dispersed (peaked) P_eq functions indicate "flexibility" ("rigidity"; the respective concepts are used below in this context). The L1 and L3 loops are throughout highly flexible, the L2 loop and the secondary structure elements are generally rigid, and the L4 loop is flexible in the monomer and rigid in the dimer. Overall, many residues are more flexible in the dimer. These features, and their implications, are discussed. Unexpectedly, we find that monomer unit 1 of the dimer (in short, d1) is unusually flexible, whereas monomer unit 2 (in short, d2) is as rigid as the RBD monomer. This is revealed due to their engagement in slow-to-intermediate conformational exchange detected previously by ¹⁵N relaxation experiments. Such motions occur with rates on the order of 10³-10⁴ s^-1; hence, they cannot be completely sampled over the course of 1 μs simulation. However, the extent to which rigid d2 is affected is small enough to enable physically relevant analysis. The entropy difference between d2 and the monomer yields an entropic contribution of -7 ± 0.7 kJ/mol to the free energy of RBD dimerization. In previous work aimed at similar objectives we used 50-100 ns MD simulations. Those results and the present result differ considerably. In summary, bond-vector P_eq functions derived directly from long MD simulations are useful descriptors of protein structural dynamics and provide accurate conformational entropy. Within the scope of slow conformational exchange, they can be useful, even in the presence of incomplete sampling.

Microsecond MD Simulations of the Plexin-B1 RBD: 2. N-H Probability Densities and Conformational Entropy in Ligand-Free, Rac1-Bound, and Dimer RBD

Orientational probability densities, P_eq = exp(-u) (u, local potential), of bond-vectors in proteins provide information on structural flexibility. The related conformational entropy, S_k = -∫P_eq(ln P_eq)dΩ - ln ∫dΩ, provides the entropic contribution to the free energy of the physical/biological process studied. We have developed a new method for deriving P_eq and S_k from MD simulations, using the N-H bond as probe. Recently we used it to study the dimerization of the Rho GTPase binding domain of Plexin-B1 (RBD). Here we use it to study RBD binding to the small GTPase Rac1. In both cases 1 μs MD simulations have been employed. The RBD has the ubiquitin fold with four mostly long loops. L3 is associated with GTPase binding, L4 with RBD dimerization, L2 participates in interdomain interactions, and L1 has not been associated with function. We find that RBD-Rac1 binding renders L1, L3, and L4 more rigid and the turns β₂/α₁ and α₂/β₅ more flexible. By comparison, RBD dimerization renders L4 more rigid, and the α-helices, the β-strands, and L2 more flexible. The rigidity of L1 in RBDRAC is consistent with L1-L3 contacts seen in previous MD simulations. The analysis of the L3-loop reveals two states of distinct flexibility which we associate with involvement in slow conformational exchange processes differing in their rates. Overall, the N-H bonds make an unfavorable entropic contribution of (5.9 ± 0.9) kJ/mol to the free energy of RBD-Rac1 binding; they were found to make a favorably contribution of (-7.0 ± 0.7) kJ/mol to the free energy of RBD dimerization. In summary, the present study provides a new perspective on the impact of Rac1 binding and dimerization on the flexibility characteristics of the RBD. Further studies are stimulated by the results of this work.

Early-stage structure-based drug discovery for small GTPases by NMR spectroscopy

Small GTPase or Ras superfamily, including Ras, Rho, Rab, Ran and Arf, are fundamental in regulating a wide range of cellular processes such as growth, differentiation, migration and apoptosis. They share structural and functional similarities for binding guanine nucleotides and hydrolyzing GTP. Dysregulations of Ras proteins are involved in the pathophysiology of multiple human diseases, however there is still a stringent need for effective treatments targeting these proteins. For decades, small GTPases were recognized as 'undruggable' targets due to their complex regulatory mechanisms and lack of deep pockets for ligand binding. NMR has been critical in deciphering the structural and dynamic properties of the switch regions that are underpinning molecular switch functions of small GTPases, which pave the way for developing new effective inhibitors. The recent progress of drug or lead molecule development made for small GTPases profoundly delineated how modern NMR techniques reshape the field of drug discovery. In this review, we will summarize the progress of structural and dynamic studies of small GTPases, the NMR techniques developed for structure-based drug screening and their applications in early-stage drug discovery for small GTPases.

Conformational Entropy from Mobile Bond Vectors in Proteins: A Viewpoint that Unifies NMR Relaxation Theory and Molecular Dynamics Simulation Approaches

A new method for determining conformational entropy in proteins is reported. Proteins prevail as conformational ensembles, p ∝ exp(-u). By selecting a bond vector (e.g., N-H) as a conformation representative, molecular dynamics simulations can provide (relative to a reference structure) p as approximate Boltzmann probability density and u as N-H potential of mean force (POMF). The latter is as accurate as implied by the force field but statistical in character; this limits the insights it can provide and its utilization. Conformational entropy is given exclusively by u. Deriving it from POMFs renders it accurate but statistical in character. Previously, we devised explicit (i.e., analytical but not exact) potentials made of Wigner functions, D_0K^L, with L ≤ 4, which closely resemble the corresponding POMFs in form; hence, they also approach the latter in accuracy. Such potentials can be beneficially characterized/compared in terms of composition, symmetry, and associated order parameters. In this study, we develop a method for deriving conformational entropy from them, which also features the benefits specified above. The method developed is applied to the dimerization of the Rho GTPase-binding domain of plexin-B1. Insights into local ordering, entropy compensation, and features of allostery are gained. In previous work, we developed the slowly relaxing local structure (SRLS) approach for the analysis of NMR relaxation from restricted bond vector motion in proteins. SRLS comprises explicit (restricting) potentials of the kind developed here. It also comprises diffusion tensors describing the local motion and related features of local geometry. The complete model fits experimental data. In future work, the explicit potentials developed here will be inserted unchanged in SRLS-based data fitting, thereby improving the picture of structural dynamics. Given that SRLS is unique in featuring potentials that can closely approach the corresponding POMFs in accuracy, the present study is an important step toward generally improving protein dynamics by NMR relaxation.

Local Ordering at the N-H Sites of the Rho GTPase Binding Domain of Plexin-B1: Impact of Dimerization

We have developed a new molecular dynamics (MD) based method for describing analytically local potentials at mobile N-H sites in proteins. Here we apply it to the monomer and dimer of the Rho GTPase binding domain (RBD) of the transmembrane receptor plexin-B1 to gain insight into dimerization, which can compete with Rho GTPase binding. In our method, the local potential is given by linear combinations, u^(D_L,K), of the real combinations of the Wigner rotation matrix elements, D_L,K, with L = 1-4 and appropriate symmetry. The combination that "fits best" the corresponding MD potential of mean force, u^(MD), is the potential we are seeking, u^{(D_L,K - BEST)}. For practical reasons the fitting process involves probability distributions, P_eq ∝ exp(-u), instead of potentials, u. The symmetry of the potential, u^(D_L,K), may be related to the irreducible representations of the D_2h point group. The monomer (dimer) potentials have mostly A_g and B_2u (B_1u and B_2u) symmetry. For the monomer, the associated probability distributions are generally dispersed in space, shallow, and centered at the "reference N-H orientation" (defined in section 3.1. below); for the dimer many are more concentrated, deep and centered away from the "reference N-H orientation". The u^(D_L,K) functions provide a consistent description of the potential energy landscape at protein N-H sites. The L1-loop of the plexin-B1 RBD is not seen in the crystal structure, and many resonances of the L4 loop are missing in the NMR ¹⁵N-¹H HSQC spectrum of the dimer; we suggest reasons for these features. An allosteric signal transmission pathway was reported previously for the monomer. We find that it has shallow N-H potentials at its ends, which become deeper as one proceeds toward the middle, complementing structurally the previously derived dynamic picture. Prospects of this study include correlating u^{(D_L,K - BEST)} with MD force-fields, and using them without further adjustment in NMR relaxation analysis schemes.

Functional divergence of Plexin B structural motifs in distinct steps of Drosophila olfactory circuit assembly

Plexins exhibit multitudinous, evolutionarily conserved functions in neural development. How Plexins employ their diverse structural motifs in vivo to perform distinct roles is unclear. We previously reported that Plexin B (PlexB) controls multiple steps during the assembly of the Drosophila olfactory circuit (Li et al., 2018b). Here, we systematically mutagenized structural motifs of PlexB and examined the function of these variants in these multiple steps: axon fasciculation, trajectory choice, and synaptic partner selection. We found that the extracellular Sema domain is essential for all three steps, the catalytic site of the intracellular RapGAP is engaged in none, and the intracellular GTPase-binding motifs are essential for trajectory choice and synaptic partner selection, but are dispensable for fasciculation. Moreover, extracellular PlexB cleavage serves as a regulatory mechanism of PlexB signaling. Thus, the divergent roles of PlexB motifs in distinct steps of neural development contribute to its functional versatility in neural circuit assembly.

Molecular Dynamics Simulations Reveal Isoform Specific Contact Dynamics between the Plexin Rho GTPase Binding Domain (RBD) and Small Rho GTPases Rac1 and Rnd1

The Plexin family of transmembrane receptors are unique in that their intracellular region interacts directly with small GTPases of the Rho family. The Rho GTPase binding domain of plexin (RBD)-which is responsible for these interactions-can bind with Rac1 as well as Rnd1 GTPases. GTPase complexes have been crystallized with the RBDs of plexinA1, -A2, and -B1. The protein association is thought to elicit different functional responses in a GTPase and plexin isoform specific manner, but the origin of this is unknown. In this project, we investigated complexes between several RBD and Rac1/Rnd1 GTPases using multimicrosecond length all atom molecular dynamics simulations, also with reference to the free forms of the RBDs and GTPases. In accord with the crystallographic data, the RBDs experience more structural changes than Rho-GTPases upon complex formation. Changes in protein dynamics and networks of correlated motions are revealed by analyzing dihedral angle fluctuations in the proteins. The extent of these changes differs between the different RBDs and also between the Rac1 and Rnd1 GTPases. While the RBDs in the free and bound states have similar-if not decreased-correlations, correlations within the GTPases are increased upon binding. Mapping highly correlated residues to the structures, it is found that the plexinA1, -B1, and -A2 RBDs all have similar communication pathways within the ubiquitin fold, but that different residues are involved. Dynamic network analyses indicate that plexinA1 and -B1 RBDs interact with small GTPases in a similar manner, whereas complexes with the plexinA2 RBD display different features. Importantly complexes with Rnd1 have a considerable number of dynamic correlations and network connections between the proteins, whereas such features are missing in the RBD-Rac1 complexes. Overall, the simulations suggest mechanisms that are consistent with the experimental data on plexinB1 and indicate RBD and GTPase isoform specific changes in protein dynamics upon complex formation.

Characterizing Plexin GTPase Interactions Using Gel Filtration, Surface Plasmon Resonance Spectrometry, and Isothermal Titration Calorimetry

Plexins are unique, as they are the first example of a transmembrane receptor that interacts directly with small GTPases, a family of proteins that are essential for cell motility and proliferation/survival. We and other laboratories have determined the structure of the Rho GTPase-binding domain (RBD) of several plexins and also of the entire intracellular region of plexin-B1. Structures of plexin complexes with Rho GTPases, Rac1 and Rnd1, and a structure with a Ras GTPase, Rap1b, have also been solved. The relationship between plexin-Rho and plexin-Ras interactions is still unclear and in vitro biophysical experiments that characterize the protein interactions of purified components play an important role in advancing our understanding of the molecular mechanisms that underlie the function of plexin. This chapter describes the use of gel filtration (also known as size-exclusion chromatography or SEC), surface plasmon resonance (SPR), and isothermal titration calorimetry (ITC) in studies of plexin-small GTPase interactions with plexin-B1:Rac1 as an example. Together with other assays and manipulations (e.g., by mutagenesis or protein domain truncation/deletion), these in vitro measurements provide an important reference for the role and extent of the interactions.

Modeling transmembrane domain dimers/trimers of plexin receptors: implications for mechanisms of signal transmission across the membrane

Single-pass transmembrane (TM) receptors transmit signals across lipid bilayers by helix association or by configurational changes within preformed dimers. The structure determination for such TM regions is challenging and has mostly been accomplished by NMR spectroscopy. Recently, the computational prediction of TM dimer structures is becoming recognized for providing models, including alternate conformational states, which are important for receptor regulation. Here we pursued a strategy to predict helix oligomers that is based on packing considerations (using the PREDDIMER webserver) and is followed by a refinement of structures, utilizing microsecond all-atom molecular dynamics simulations. We applied this method to plexin TM receptors, a family of 9 human proteins, involved in the regulation of cell guidance and motility. The predicted models show that, overall, the preferences identified by PREDDIMER are preserved in the unrestrained simulations and that TM structures are likely to be diverse across the plexin family. Plexin-B1 and -B3 TM helices are regular and tend to associate, whereas plexin-A1, -A2, -A3, -A4, -C1 and -D1 contain sequence elements, such as poly-Glycine or aromatic residues that distort helix conformation and association. Plexin-B2 does not form stable dimers due to the presence of TM prolines. No experimental structural information on the TM region is available for these proteins, except for plexin-C1 dimeric and plexin-B1 - trimeric structures inferred from X-ray crystal structures of the intracellular regions. Plexin-B1 TM trimers utilize Ser and Thr sidechains for interhelical contacts. We also modeled the juxta-membrane (JM) region of plexin-C1 and plexin-B1 and show that it synergizes with the TM structures. The structure and dynamics of the JM region and TM-JM junction provide determinants for the distance and distribution of the intracellular domains, and for their binding partners relative to the membrane. The structures suggest experimental tests and will be useful for the interpretation of future studies.

Loss of plexin-B3 in hepatocellular carcinoma

Plexins are the primary receptors of semaphorins, and participate in the majority of intracellular pathways triggered by semaphorins, including the regulation of cell adhesion and the motility of numerous cell types. Recently, several studies have reported that plexins can significantly affect different aspects of cancer cell biology, and the aberrant expression of plexins has been observed in a wide variety of tumor types. However, the expression and role of plexin-B3 in hepatocellular carcinoma (HCC) is yet to be investigated. In the present study, plexin-B3 expression was measured in 14 paired HCC samples and the corresponding adjacent non-cancerous tissue by quantitative polymerase chain reaction and western blot analysis. The results indicated that the mRNA and protein expression levels of plexin-B3 were downregulated in HCC samples when compared with the corresponding adjacent non-cancerous tissue. In order to elucidate the correlation between clinicopathological data and the expression of plexin-B3 in patients with HCC, 84 HCC archived specimens were analyzed by immunohistochemistry (IHC). The IHC results revealed that the protein expression level of plexin-B3 was lower in the HCC samples compared with the corresponding adjacent non-cancerous tissue, and plexin-B3 underexpression was correlated with the patient gender and tumor size. In conclusion, these results indicated that loss of plexin-B3 in HCC may be of predictive value for the occurrence and progression of HCC. Thus, plexin-B3 may be a promising biomarker for the diagnosis and treatment of tumors in the future.

The Rho GTPase Rnd1 suppresses mammary tumorigenesis and EMT by restraining Ras-MAPK signalling

We identified the Rho GTPase Rnd1 as a candidate metastasis suppressor in basal-like and triple-negative breast cancer through bioinformatics analysis. Depletion of Rnd1 disrupted epithelial adhesion and polarity, induced epithelial-to-mesenchymal transition, and cooperated with deregulated expression of c-Myc or loss of p53 to cause neoplastic conversion. Mechanistic studies revealed that Rnd1 suppresses Ras signalling by activating the GAP domain of Plexin B1, which inhibits Rap1. Rap1 inhibition in turn led to derepression of p120 Ras-GAP, which was able to inhibit Ras. Inactivation of Rnd1 in mammary epithelial cells induced highly undifferentiated and invasive tumours in mice. Conversely, Rnd1 expression inhibited spontaneous and experimental lung colonization in mouse models of metastasis. Genomic studies indicated that gene deletion in combination with epigenetic silencing or, more rarely, point mutation inactivates RND1 in human breast cancer. These results reveal a previously unappreciated mechanism through which Rnd1 restrains activation of Ras-MAPK signalling and breast tumour initiation and progression.

Structure and dynamics analysis on plexin-B1 Rho GTPase binding domain as a monomer and dimer

Plexin-B1 is a single-pass transmembrane receptor. Its Rho GTPase binding domain (RBD) can associate with small Rho GTPases and can also self-bind to form a dimer. In total, more than 400 ns of NAMD molecular dynamics simulations were performed on RBD monomer and dimer. Different analysis methods, such as root mean squared fluctuation (RMSF), order parameters (S(2)), dihedral angle correlation, transfer entropy, principal component analysis, and dynamical network analysis, were carried out to characterize the motions seen in the trajectories. RMSF results show that after binding, the L4 loop becomes more rigid, but the L2 loop and a number of residues in other regions become slightly more flexible. Calculating order parameters (S(2)) for CH, NH, and CO bonds on both backbone and side chain shows that the L4 loop becomes essentially rigid after binding, but part of the L1 loop becomes slightly more flexible. Backbone dihedral angle cross-correlation results show that loop regions such as the L1 loop including residues Q25 and G26, the L2 loop including residue R61, and the L4 loop including residues L89-R91, are highly correlated compared to other regions in the monomer form. Analysis of the correlated motions at these residues, such as Q25 and R61, indicate two signal pathways. Transfer entropy calculations on the RBD monomer and dimer forms suggest that the binding process should be driven by the L4 loop and C-terminal. However, after binding, the L4 loop functions as the motion responder. The signal pathways in RBD were predicted based on a dynamical network analysis method using the pathways predicted from the dihedral angle cross-correlation calculations as input. It is found that the shortest pathways predicted from both inputs can overlap, but signal pathway 2 (from F90 to R61) is more dominant and overlaps all of the routes of pathway 1 (from F90 to P111). This project confirms the allosteric mechanism in signal transmission inside the RBD network, which was in part proposed in the previous experimental study.

Interaction characteristics of Plexin-B1 with Rho family proteins

Plexin-B1 regulates various cellular processes interacting directly with several Rho proteins. Molecular details of these interactions are, however, not well understood. In this study, we examined in vitro and in silico the interaction of the Rho binding domain (B1RBD) of human Plexin-B1 with 11 different Rho proteins. We show that B1RBD binds in a GTP-dependent manner to Rac1, Rac2, Rac3, Rnd1, Rnd2, Rnd3, and RhoD, but not to RhoA, Cdc42, RhoG, or Rif. Interestingly, Rnd1 competitively displaces the Rac1 from B1RBD but not vice versa. Structure-function analysis revealed a negatively charged loop region, called B1L(31), which may facilitate a selective B1RBD interaction with Rho proteins.

Semaphorins in cancer: biological mechanisms and therapeutic approaches

The hallmarks of cancer include multiple alterations in the physiological processes occurring in normal tissues, such as cell proliferation, apoptosis, and restricted cell migration. These aberrant behaviors are due to genetic and epigenetic changes that affect signaling pathways controlling cancer cells, as well as the surrounding "normal" cells in the tumor microenvironment. Semaphorins and their receptors (mainly plexins and neuropilins) are aberrantly expressed in human tumors, and multiple family members are emerging as pivotal signals deregulated in cancer. Notably, different semaphorins can promote or inhibit tumor progression, depending on the implicated receptor complexes and responsive cell type. The important role of semaphorin signals in the regulation of tumor angiogenesis, invasion and metastasis has initiated multiple experimental approaches aimed at targeting these pathways to inhibit cancer.

Structural insights into semaphorins and their receptors

Ten years ago nothing was known of the three-dimensional structure of members of the semaphorin family of cell guidance cues, nor of their major receptors, the plexins. The structural biology of this cell surface ligand-receptor system has now come of age. Detailed atomic level information is available on the architecture of semaphorin and plexin ectodomains and their recognition complexes. Similarly the structure of the plexin cytoplasmic region, and its interactions with members of the Rho family of small GTPases have been unveiled. These structural analyses, in combination with biochemical, biophysical and cellular studies, have progressed our understanding of this signalling system into the realm of molecular mechanism.

Analysis of 15N-1H NMR relaxation in proteins by a combined experimental and molecular dynamics simulation approach: picosecond-nanosecond dynamics of the Rho GTPase binding domain of plexin-B1 in the dimeric state indicates allosteric pathways

We investigate picosecond–nanosecond dynamics of the Rho-GTPase Binding Domain (RBD) of plexin-B1, which plays a key role in plexin-mediated cell signaling. Backbone 15N relaxation data of the dimeric RBD are analyzed with the model-free (MF) method, and with the slowly relaxing local structure/molecular dynamics (SRLS-MD) approach. Independent analysis of the MD trajectories, based on the MF paradigm, is also carried out. MF is a widely popular and simple method, SRLS is a general approach, and SRLS-MD is an integrated approach we developed recently. Corresponding parameters from the RBD dimer, a previously studied RBD monomer mutant, and the previously studied complex of the latter with the GTPase Rac1, are compared. The L2, L3, and L4 loops of the plexin-B1 RBD are involved in interactions with other plexin domains, GTPase binding, and RBD dimerization, respectively. Peptide groups in the loops of both the monomeric and dimeric RBD are found to experience weak and moderately asymmetric local ordering centered approximately at the C(i–1)(α)–C(i)(α) axes, and nanosecond backbone motion. Peptide groups in the α-helices and the β-strands of the dimer (the β-strands of the monomer) experience strong and highly asymmetric local ordering centered approximately at the C(i–1)(α)–C(i)(α) axes (N–H bonds). N–H fluctuations occur on the picosecond time scale. An allosteric pathway for GTPase binding, providing new insights into plexin function, is delineated.

Plexin structures are coming: opportunities for multilevel investigations of semaphorin guidance receptors, their cell signaling mechanisms, and functions

Plexin transmembrane receptors and their semaphorin ligands, as well as their co-receptors (Neuropilin, Integrin, VEGFR2, ErbB2, and Met kinase) are emerging as key regulatory proteins in a wide variety of developmental, regenerative, but also pathological processes. The diverse arenas of plexin function are surveyed, including roles in the nervous, cardiovascular, bone and skeletal, and immune systems. Such different settings require considerable specificity among the plexin and semaphorin family members which in turn are accompanied by a variety of cell signaling networks. Underlying the latter are the mechanistic details of the interactions and catalytic events at the molecular level. Very recently, dramatic progress has been made in solving the structures of plexins and of their complexes with associated proteins. This molecular level information is now suggesting detailed mechanisms for the function of both the extracellular as well as the intracellular plexin regions. Specifically, several groups have solved structures for extracellular domains for plexin-A2, -B1, and -C1, many in complex with semaphorin ligands. On the intracellular side, the role of small Rho GTPases has been of particular interest. These directly associate with plexin and stimulate a GTPase activating (GAP) function in the plexin catalytic domain to downregulate Ras GTPases. Structures for the Rho GTPase binding domains have been presented for several plexins, some with Rnd1 bound. The entire intracellular domain structure of plexin-A1, -A3, and -B1 have also been solved alone and in complex with Rac1. However, key aspects of the interplay between GTPases and plexins remain far from clear. The structural information is helping the plexin field to focus on key questions at the protein structural, cellular, as well as organism level that collaboratoria of investigations are likely to answer.

Effect of cancer-associated mutations in the PlexinB1 gene

Background

Semaphorins act as chemotactic cues for cell movement via their transmembrane receptors, plexins. Somatic missense mutations in the plexinB1 gene coupled with overexpression of the protein frequently occur in prostate tumours, indicating a role for plexinB1 in the pathogenesis of prostate cancer.

Results

Two specific mutations found in prostate cancer enhance RhoD binding and one other mutation results in loss of inhibition of Rac-dependent Pak1 phosphorylation and lamellipodia formation and in impairment of trafficking of plexinB1 to the membrane. None of the three characterised mutations affect PDZRhoGEF binding, RhoA activity, the interaction of plexinB1with the oncogenes ErbB2 or c-Met or ErbB2 phosphorylation. The mutations have the net effect of increasing cell motility by blocking plexinB1-mediated inhibition of Rac while enhancing the interaction with RhoD, an anti-migratory factor.

Conclusions

PlexinB1 mutations block plexinB1-mediated signalling pathways that inhibit cell motility.

Combining NMR and molecular dynamics studies for insights into the allostery of small GTPase-protein interactions

Combinations of experimentally derived data from nuclear magnetic resonance spectroscopy and analyses of molecular dynamics trajectories increasingly allow us to obtain a detailed description of the molecular mechanisms by which proteins function in signal transduction. This chapter provides an introduction into these two methodologies, illustrated by example of a small GTPase-effector interaction. It is increasingly becoming clear that new insights are provided by the combination of experimental and computational methods. Understanding the structural and protein dynamical contributions to allostery will be useful for the engineering of new binding interfaces and protein functions, as well as for the design/in silico screening of chemical agents that can manipulate the function of small GTPase-protein interactions in diseases such as cancer.

A direct coupling between global and internal motions in a single domain protein? MD investigation of extreme scenarios

Proteins are not rigid molecules, but exhibit internal motions on timescales ranging from femto- to milliseconds and beyond. In solution, proteins also experience global translational and rotational motions, sometimes on timescales comparable to those of the internal fluctuations. The possibility that internal and global motions may be directly coupled has intriguing implications, given that enzymes and cell signaling proteins typically associate with binding partners and cellular scaffolds. Such processes alter their global motion and may affect protein function. Here, we present molecular dynamics simulations of extreme case scenarios to examine whether a possible relationship exists. In our model protein, a ubiquitin-like RhoGTPase binding domain of plexin-B1, we removed either internal or global motions. Comparisons with unrestrained simulations show that internal and global motions are not appreciably coupled in this single-domain protein. This lack of coupling is consistent with the observation that the dynamics of water around the protein, which is thought to permit, if not stimulate, internal dynamics, is also largely independent of global motion. We discuss implications of these results for the structure and function of proteins.

A dual binding mode for RhoGTPases in plexin signalling

Plexins are cell surface receptors for the semaphorin family of cell guidance cues. The cytoplasmic region comprises a Ras GTPase-activating protein (GAP) domain and a RhoGTPase binding domain. Concomitant binding of extracellular semaphorin and intracellular RhoGTPase triggers GAP activity and signal transduction. The mechanism of this intricate regulation remains elusive. We present two crystal structures of the human Plexin-B1 cytoplasmic region in complex with a constitutively active RhoGTPase, Rac1. The structure of truncated Plexin-B1-Rac1 complex provides no mechanism for coupling RhoGTPase and Ras binding sites. On inclusion of the juxtamembrane helix, a trimeric structure of Plexin-B1-Rac1 complexes is stabilised by a second, novel, RhoGTPase binding site adjacent to the Ras site. Site-directed mutagenesis combined with cellular and biophysical assays demonstrate that this new binding site is essential for signalling. Our findings are consistent with a model in which extracellular and intracellular plexin clustering events combine into a single signalling output.

Structural basis of Rnd1 binding to plexin Rho GTPase binding domains (RBDs)

Plexin receptors regulate cell adhesion, migration, and guidance. The Rho GTPase binding domain (RBD) of plexin-A1 and -B1 can bind GTPases, including Rnd1. By contrast, plexin-C1 and -D1 reportedly bind Rnd2 but associate with Rnd1 only weakly. The structural basis of this differential Rnd1 GTPase binding to plexin RBDs remains unclear. Here, we solved the structure of the plexin-A2 RBD in complex with Rnd1 and the structures of the plexin-C1 and plexin-D1 RBDs alone, also compared with the previously determined plexin-B1 RBD.Rnd1 complex structure. The plexin-A2 RBD·Rnd1 complex is a heterodimer, whereas plexin-B1 and -A2 RBDs homodimerize at high concentration in solution, consistent with a proposed model for plexin activation. Plexin-C1 and -D1 RBDs are monomeric, consistent with major residue changes in the homodimerization loop. In plexin-A2 and -B1, the RBD β3-β4 loop adjusts its conformation to allow Rnd1 binding, whereas minimal structural changes occur in Rnd1. The plexin-C1 and -D1 RBDs lack several key non-polar residues at the corresponding GTPase binding surface and do not significantly interact with Rnd1. Isothermal titration calorimetry measurements on plexin-C1 and -D1 mutants reveal that the introduction of non-polar residues in this loop generates affinity for Rnd1. Structure and sequence comparisons suggest a similar mode of Rnd1 binding to the RBDs, whereas mutagenesis suggests that the interface with the highly homologous Rnd2 GTPase is different in detail. Our results confirm, from a structural perspective, that Rnd1 does not play a role in the activation of plexin-C1 and -D1. Plexin functions appear to be regulated by subfamily-specific mechanisms, some of which involve different Rho family GTPases.

Novel modulation factor quantifies the role of water molecules in protein interactions

Water molecules decrease the potential of mean force of a hydrogen bond (H-bond), as well as modulate (de)solvation forces, but exactly how much has not been easy to determine. Crystallographic water molecules provide snapshots of optimal solutions for the role of solvent in protein interactions, information that is often ignored by implicit solvent models. Motivated by high-resolution crystal structures, we describe a simple quantitative approach to explicitly incorporate the role of molecular water in protein interactions. Applications to protein-DNA interactions show that the accuracy of binding free-energy estimates improves significantly if a distinction is made between H-bonds that are desolvated (or only contact crystal waters), solvated by mobile waters trapped at the binding interface, or partially solvated through connections to bulk water. These different environments are modeled by a unique "water" scaling factor that decreases or increases the strength of hydrogen bonds depending on whether water contacts the acceptor or donor atoms or the bond is fully desolvated, respectively. Our empirical energies are fully consistent with mobile water molecules having a strong polarization effect in direct intermolecular interactions.

Integrated computational approach to the analysis of NMR relaxation in proteins: application to ps-ns main chain 15N-1H and global dynamics of the Rho GTPase binding domain of plexin-B1

An integrated computational methodology for interpreting NMR spin relaxation in proteins has been developed. It combines a two-body coupled-rotator stochastic model with a hydrodynamics-based approach for protein diffusion, together with molecular dynamics based calculations for the evaluation of the coupling potential of mean force. The method is applied to ¹⁵N relaxation of N-H bonds in the Rho GTPase binding (RBD) domain of plexin-B1, which exhibits intricate internal mobility. Bond vector dynamics are characterized by a rhombic local ordering tensor, S, with principal values S₀² and S₂², and an axial local diffusion tensor, D₂, with principal values D(2,||) and D(2,⊥). For α-helices and β-sheets we find that S₀² ~ -0.5 (strong local ordering), -1.2 < S₂² < -0.8 (large S tensor anisotropy), D(2,⊥) ~ D₁ = 1.93 × 10⁷ s⁻¹ (D₁ is the global diffusion rate), and log(D(2,||)/D₁) ~ 4. For α-helices the z-axis of the local ordering frame is parallel to the C(α)-C(α) axis. For β-sheets the z-axes of the S and D₂ tensors are parallel to the N-H bond. For loops and terminal chain segments the local ordering is generally weaker and more isotropic. On average, D(2,⊥) ~ D₁ also, but log(D(2,||)/D₁) is on the order of 1-2. The tensor orientations are diversified. This study sets forth an integrated computational approach for treating NMR relaxation in proteins by combining stochastic modeling and molecular dynamics. The approach developed provides new insights by its application to a protein that experiences complex dynamics.

Roles of Sema4D and Plexin-B1 in tumor progression

Sema4D, also known as CD100, is a protein belonging to class IV semaphorin. Its physiologic roles in the immune and nervous systems have been extensively explored. However, the roles of Sema4D have extended beyond these traditionally studied territories. Via interaction with its high affinity receptor Plexin-B1, Sema4D-Plexin-B1 involvement in tumor progression is strongly implied. Here, we critically review and delineate the Sema4D-Plexin-B1 interaction in many facets of tumor progression: tumor angiogenesis, regulation of tumor-associated macrophages and control of invasive growth. We correlate the in vitro and in vivo experimental data with the clinical study outcomes, and present a molecular mechanistic basis accounting for the intriguingly contradicting results from these recent studies.

Semaphorin 5A and plexin-B3 inhibit human glioma cell motility through RhoGDIalpha-mediated inactivation of Rac1 GTPase

Semaphorins and plexins are implicated in the progression of various types of cancer, although the molecular basis has not been fully elucidated. Here, we report the expression of plexin-B3 in glioma cells, which upon stimulation by its ligand Sema5A results in significant inhibition of cell migration and invasion. A search for the underlying mechanism revealed direct interaction of plexin-B3 with RhoGDP dissociation inhibitor α (RhoGDIα), a negative regulator of RhoGTPases that blocks guanine nucleotide exchange and sequesters them away from the plasma membrane. Glioma cells challenged with Sema5A indeed showed a marked reduction in Rac1-GTP levels by 60%, with a concomitant disruption of lamellipodia. The inactivation of Rac1 was corroborated to contribute to the impediment of glioma cell invasion by Sema5A, as supported by the abolishment of effect upon forced expression of a constitutively active Rac1 mutant. Furthermore, silencing the endogenous expression of RhoGDIα in glioma cells was found to be sufficient in abrogating the down-regulation of Rac1-GTP and the ensuing suppression of glioma cell motility induced by Sema5A. Mechanistically, we provide evidence that Sema5A promotes Rac1 recruitment to RhoGDIα and reduces its membrane localization in a plexin-B3-dependent manner, thereby preventing Rac1 activation. This represents a novel signaling of semaphorin and plexin in the control of cell motility by indirect inactivation of Rac1 through RhoGDIα.

Structure and function of the intracellular region of the plexin-b1 transmembrane receptor

Members of the plexin family are unique transmembrane receptors in that they interact directly with Rho family small GTPases; moreover, they contain a GTPase-activating protein (GAP) domain for R-Ras, which is crucial for plexin-mediated regulation of cell motility. However, the functional role and structural basis of the interactions between the different intracellular domains of plexins remained unclear. Here we present the 2.4 A crystal structure of the complete intracellular region of human plexin-B1. The structure is monomeric and reveals that the GAP domain is folded into one structure from two segments, separated by the Rho GTPase binding domain (RBD). The RBD is not dimerized, as observed previously. Instead, binding of a conserved loop region appears to compete with dimerization and anchors the RBD to the GAP domain. Cell-based assays on mutant proteins confirm the functional importance of this coupling loop. Molecular modeling based on structural homology to p120(GAP).H-Ras suggests that Ras GTPases can bind to the plexin GAP region. Experimentally, we show that the monomeric intracellular plexin-B1 binds R-Ras but not H-Ras. These findings suggest that the monomeric form of the intracellular region is primed for GAP activity and extend a model for plexin activation.

Crystal structure of the plexin A3 intracellular region reveals an autoinhibited conformation through active site sequestration

Plexin cell surface receptors bind to semaphorin ligands and transduce signals for regulating neuronal axon guidance. The intracellular region of plexins is essential for signaling and contains a R-Ras/M-Ras GTPase activating protein (GAP) domain that is divided into two segments by a Rho GTPase-binding domain (RBD). The regulation mechanisms for plexin remain elusive, although it is known that activation requires both binding of semaphorin to the extracellular region and a Rho-family GTPase (Rac1 or Rnd1) to the RBD. Here we report the crystal structure of the plexin A3 intracellular region. The structure shows that the N- and C-terminal portions of the GAP homologous regions together form a GAP domain with an overall fold similar to other Ras GAPs. However, the plexin GAP domain adopts a closed conformation and cannot accommodate R-Ras/M-Ras in its substrate-binding site, providing a structural basis for the autoinhibited state of plexins. A comparison with the plexin B1 RBD/Rnd1 complex structure suggests that Rnd1 binding alone does not induce a conformational change in plexin, explaining the requirement of both semaphorin and a Rho GTPase for activation. The structure also identifies an N-terminal segment that is important for regulation. Both the N-terminal segment and the RBD make extensive interactions with the GAP domain, suggesting the presence of an allosteric network connecting these three domains that integrates semaphorin and Rho GTPase signals to activate the GAP. The importance of these interactions in plexin signaling is shown by both cell-based and in vivo axon guidance assays.

Tyrosine phosphorylation in semaphorin signalling: shifting into overdrive

The semaphorins constitute a large family of molecular signals with regulatory functions in neuronal development, angiogenesis, cancer progression and immune responses. Accumulating data indicate that semaphorins might trigger multiple signalling pathways, and mediate different and sometimes opposing effects, depending on the cellular context and the particular plexin-associated subunits of the receptor complex, which can include receptor-type or cytoplasmic tyrosine kinases such as MET, ERBB2, VEGFR2, FYN, FES, PYK2 and SRC. It has also been shown that a specific plexin can alternatively associate with different tyrosine kinase receptors, eliciting divergent signalling pathways and functional outcomes. Tyrosine phosphorylation is a pivotal post-translational protein modification that regulates intracellular signalling. Therefore, phosphorylation of tyrosines in the intracellular domain of plexins could determine or modify their interactions with additional signal transducers. Here, we discuss the potential relevance of tyrosine phosphorylation in semaphorin-induced signalling, with an emphasis on its probable role in dictating the choice between multiple pathways and functional outcomes. The identification of implicated tyrosine kinases will pave the way to target individual semaphorin-mediated functions.

Compensatory and long-range changes in picosecond-nanosecond main-chain dynamics upon complex formation: 15N relaxation analysis of the free and bound states of the ubiquitin-like domain of human plexin-B1 and the small GTPase Rac1

The formation of a complex between Rac1 and the cytoplasmic domain of plexin-B1 is one of the first documented cases of a direct interaction between a small guanosine 5'-triphosphatase (GTPase) and a transmembrane receptor. Structural studies have begun to elucidate the role of this interaction for the signal transduction mechanism of plexins. Mapping of the Rac1 GTPase surface that contacts the Rho GTPase binding domain of plexin-B1 by solution NMR spectroscopy confirms the plexin domain as a GTPase effector protein. Regions neighboring the GTPase switch I and II regions are also involved in the interaction and there is considerable interest to examine the changes in protein dynamics that take place upon complex formation. Here we present main-chain nitrogen-15 relaxation measurements for the unbound proteins as well as for the Rho GTPase binding domain and Rac1 proteins in their complexed state. Derived order parameters, S2, show that considerable motions are maintained in the bound state of plexin. In fact, some of the changes in S2 on binding appear compensatory, exhibiting decreased as well as increased dynamics. Fluctuations in Rac1, already a largely rigid protein on the picosecond-nanosecond timescale, are overall diminished, but isomerization dynamics in the switch I and II regions of the GTPase are retained in the complex and appear to be propagated to the bound plexin domain. Remarkably, fluctuations in the GTPase are attenuated at sites, including helices alpha6 (the Rho-specific insert helix), alpha7 and alpha8, that are spatially distant from the interaction region with plexin. This effect of binding on long-range dynamics appears to be communicated by hinge sites and by subtle conformational changes in the protein. Similar to recent studies on other systems, we suggest that dynamical protein features are affected by allosteric mechanisms. Altered protein fluctuations are likely to prime the Rho GTPase-plexin complex for interactions with additional binding partners.