1H, 15N and 13C Resonance assignments and secondary structure determination reveal that the minimal Rac1 GTPase binding domain of plexin-B1 has a ubiquitin fold

Balancing act of small GTPases downstream of plexin-A4 signaling motifs promotes dendrite elaboration in mammalian cortical neurons

The precise development of neuronal morphologies is crucial to the establishment of synaptic circuits and, ultimately, proper brain function. Signaling by the axon guidance cue semaphorin 3A (Sema3A) and its receptor complex of neuropilin-1 and plexin-A4 has multifunctional outcomes in neuronal morphogenesis. Downstream activation of the RhoGEF FARP2 through interaction with the lysine-arginine-lysine motif of plexin-A4 and consequent activation of the small GTPase Rac1 promotes dendrite arborization, but this pathway is dispensable for axon repulsion. Here, we investigated the interplay of small GTPase signaling mechanisms underlying Sema3A-mediated dendritic elaboration in mouse layer V cortical neurons in vitro and in vivo. Sema3A promoted the binding of the small GTPase Rnd1 to the amino acid motif lysine-valine-serine (LVS) in the cytoplasmic domain of plexin-A4. Rnd1 inhibited the activity of the small GTPase RhoA and the kinase ROCK, thus supporting the activity of the GTPase Rac1, which permitted the growth and branching of dendrites. Overexpression of a dominant-negative RhoA, a constitutively active Rac1, or the pharmacological inhibition of ROCK activity rescued defects in dendritic elaboration in neurons expressing a plexin-A4 mutant lacking the LVS motif. Our findings provide insights into the previously unappreciated balancing act between Rho and Rac signaling downstream of specific motifs in plexin-A4 to mediate Sema3A-dependent dendritic elaboration in mammalian cortical neuron development.

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.

Insights into the regulatory role of Plexin D1 signalling in cardiovascular development and diseases

Plexin D1 (PLXND1), which was previously thought to mediate semaphorin signalling, belongs to the Plexin family of transmembrane proteins. PLXND1 cooperates mostly with the coreceptor neuropilin and participates in many aspects of axonal guidance. PLXND1 can also act as both a tumour promoter and a tumour suppressor. Emerging evidence suggests that mutations in PLXND1 or Semaphorin 3E, the canonical ligand of PLXND1, can lead to serious cardiovascular diseases, such as congenital heart defects, CHARGE syndrome and systemic sclerosis. Upon ligand binding, PLXND1 can act as a GTPase‐activating protein (GAP) and modulate integrin‐mediated cell adhesion, cytoskeletal dynamics and cell migration. These effects may play regulatory roles in the development of the cardiovascular system and disease. The cardiovascular effects of PLXND1 signalling have gradually been elucidated. PLXND1 was recently shown to detect physical forces and translate them into intracellular biochemical signals in the context of atherosclerosis. Therefore, the role of PLXND1 in cardiovascular development and diseases is gaining research interest because of its potential as a biomarker and therapeutic target. In this review, we describe the cardiac effects, vascular effects and possible molecular mechanisms of PLXND1 signalling.

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.

Semaphorin signaling in cardiovascular development

Semaphorins were originally identified as neuronal guidance molecules mediating their attractive or repulsive signals by forming complexes with plexin and neuropilin receptors. Subsequent research has identified functions for semaphorin signaling in many organs and tissues outside of the nervous system. Vital roles for semaphorin signaling in vascular patterning and cardiac morphogenesis have been demonstrated, and impaired semaphorin signaling has been associated with various human cardiovascular disorders, including persistent truncus arteriosus, sinus bradycardia and anomalous pulmonary venous connections. Here, we review the functions of semaphorins and their receptors in cardiovascular development and disease and highlight important recent discoveries in the field.

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.

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.

An activity-regulated microRNA, miR-188, controls dendritic plasticity and synaptic transmission by downregulating neuropilin-2

MicroRNAs (miRNAs) have recently come to be viewed as critical players that modulate a number of cellular features in various biological systems including the mature CNS by exerting regulatory control over the stability and translation of mRNAs. Despite considerable evidence for the regulatory functions of miRNAs, the identities of the miRNA species that are involved in the regulation of synaptic transmission and plasticity and the mechanisms by which these miRNAs exert functional roles remain largely unknown. In the present study, the expression of microRNA-188 (miR-188) was found to be upregulated by the induction of long-term potentiation (LTP). The protein level of neuropilin-2 (Nrp-2), one of the possible molecular targets for miR-188, was decreased during LTP induction. We also confirmed that the luciferase activity of the 3'-UTR of Nrp-2 was diminished by treatment with a miR-188 oligonucleotide but not with a scrambled miRNA oligonucleotide. Nrp-2 serves as a receptor for semaphorin 3F, which is a negative regulator of spine development and synaptic structure. In addition, miR-188 specifically rescued the reduction in dendritic spine density induced by Nrp-2 expression in hippocampal neurons from rat primary culture. Furthermore, miR-188 counteracted the decrease in the miniature EPSC frequency induced by Nrp-2 expression in hippocampal neurons from rat primary culture. These findings suggest that miR-188 serves to fine-tune synaptic plasticity by regulating Nrp-2 expression.

Diverse functions for the semaphorin receptor PlexinD1 in development and disease

Plexins are a family of single-pass transmembrane proteins that serve as cell surface receptors for Semaphorins during the embryonic development of animals. Semaphorin-Plexin signaling is critical for many cellular aspects of organogenesis, including cell migration, proliferation and survival. Until recently, little was known about the function of PlexinD1, the sole member of the vertebrate-specific PlexinD (PlxnD1) subfamily. Here we review novel findings about PlxnD1's roles in the development of the cardiovascular, nervous and immune systems and salivary gland branching morphogenesis and discuss new insights concerning the molecular mechanisms of PlxnD1 activity.

Thermodynamic characterization of two homologous protein complexes: associations of the semaphorin receptor plexin-B1 RhoGTPase binding domain with Rnd1 and active Rac1

Plexin receptors function in response to semaphorin guidance cues in a variety of developmental processes involving cell motility. Interactions with Rho, as well as Ras family small GTPases are critical events in the cell signaling mechanism. We have recently determined the structure of a cytoplasmic domain (RBD) of plexin-B1 and mapped its binding interface with several Rho-GTPases, Rac1, Rnd1, and RhoD. All three GTPases associate with a similar region of this plexin domain, but show different functional behavior in cells. To understand whether thermodynamic properties of the GTPase-RBD interaction contribute to such different behavior, we have examined the interaction at different temperatures, buffer, and pH conditions. Although the binding affinity of both Rnd1 and Rac1 with the plexin-B1 RBD is similar, the detailed thermodynamic properties of the interactions are considerably different. These data suggest that on Rac1 binding to the plexin-B1 RBD, the proteins become more rigid in the complex. By contrast, Rnd1 binding is consistent with unchanged or slightly increased flexibility in one or both proteins. Both GTPases show an appreciable reduction in affinity for the dimeric plexin-B1 RBD indicating that GTPase binding is not cooperative with dimer formation, but that a partial steric hindrance destabilizes the dimer. However, a reduced affinity binding mode to a disulphide stabilized model for the dimeric RBD is also possible. Consistent with cellular studies, the interaction thermodynamics imply that further levels of regulation involving additional binding partners and/or regions outside of the RhoGTPase binding domain are required for receptor activation.

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.

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

The plexin family of transmembrane receptors are important for axon guidance, angiogenesis, but also in cancer. Recently, plexin-B1 somatic missense mutations were found in both primary tumors and metastases of breast and prostate cancers, with several mutations mapping to the Rho GTPase binding domain (RBD) in the cytoplasmic region of the receptor. Here we present the NMR solution structure of this domain, confirming that the protein has both a ubiquitin-like fold and surface features. Oncogenic mutations T1795A and T1802A are located in a loop region, perturb the average structure locally, and have no effect on Rho GTPase binding affinity. Mutations L1815F and L1815P are located at the Rho GTPase binding site and are associated with a complete loss of binding for Rac1 and Rnd1. Both are found to disturb the conformation of the beta3-beta4 sheet and the orientation of surrounding side chains. Our study suggests that the oncogenic behavior of the mutants can be rationalized with reference to the structure of the RBD of plexin-B1.

Binding of Rac1, Rnd1, and RhoD to a novel Rho GTPase interaction motif destabilizes dimerization of the plexin-B1 effector domain

Plexins are the first known transmembrane receptors that interact directly with small GTPases. On binding to certain Rho family GTPases, the receptor regulates the remodeling of the actin cytoskeleton and alters cell movement in response to semaphorin guidance cues. In a joint solution NMR spectroscopy and x-ray crystallographic study, we characterize a 120-residue cytoplasmic independent folding domain of plexin-B1 that directly binds three Rho family GTPases, Rac1, Rnd1, and RhoD. The NMR data show that, surprisingly, the Cdc42/Rac interactive binding-like motif of plexin-B1 is not involved in this interaction. Instead, all three GTPases interact with the same region, beta-strands 3 and 4 and a short alpha-helical segment of the plexin domain. The 2.0 A resolution x-ray structure shows that these segments are brought together by the tertiary structure of the ubiquitin-like fold. In the crystal, the protein is dimerized with C2 symmetry through a four-stranded antiparallel beta-sheet that is formed outside the fold by a long loop between the monomers. This region is adjacent to the GTPase binding motifs identified by NMR. Destabilization of the dimer in solution by binding of any one of the three GTPases suggests a model for receptor regulation that involves bidirectional signaling. The model implies a multifunctional role for the GTPase-plexin interaction that includes conformational change and a localization of active receptors in the signaling mechanism.

Prediction of Ras-effector interactions using position energy matrices

MOTIVATION

One of the more challenging problems in biology is to determine the cellular protein interaction network. Progress has been made to predict protein-protein interactions based on structural information, assuming that structural similar proteins interact in a similar way. In a previous publication, we have determined a genome-wide Ras-effector interaction network based on homology models, with a high accuracy of predicting binding and non-binding domains. However, for a prediction on a genome-wide scale, homology modelling is a time-consuming process. Therefore, we here successfully developed a faster method using position energy matrices, where based on different Ras-effector X-ray template structures, all amino acids in the effector binding domain are sequentially mutated to all other amino acid residues and the effect on binding energy is calculated. Those pre-calculated matrices can then be used to score for binding any Ras or effector sequences.

RESULTS

Based on position energy matrices, the sequences of putative Ras-binding domains can be scanned quickly to calculate an energy sum value. By calibrating energy sum values using quantitative experimental binding data, thresholds can be defined and thus non-binding domains can be excluded quickly. Sequences which have energy sum values above this threshold are considered to be potential binding domains, and could be further analysed using homology modelling. This prediction method could be applied to other protein families sharing conserved interaction types, in order to determine in a fast way large scale cellular protein interaction networks. Thus, it could have an important impact on future in silico structural genomics approaches, in particular with regard to increasing structural proteomics efforts, aiming to determine all possible domain folds and interaction types.

AVAILABILITY

All matrices are deposited in the ADAN database (http://adan-embl.ibmc.umh.es/).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Semaphorins command cells to move

Semaphorins are secreted or transmembrane proteins that regulate cell motility and attachment in axon guidance, vascular growth, immune cell regulation and tumour progression. The main receptors for semaphorins are plexins, which have established roles in regulating Rho-family GTPases. Recent work shows that plexins can also influence R-Ras, which, in turn, can regulate integrins. Such regulation is probably a common feature of semaphorin signalling and contributes substantially to our understanding of semaphorin biology.

The ubiquitin domain superfold: structure-based sequence alignments and characterization of binding epitopes

Ubiquitin-like domains are present, apart from ubiquitin-like proteins themselves, in many multidomain proteins involved in different signal transduction processes. The sequence conservation for all ubiquitin superfold family members is rather poor, even between subfamily members, leading to mistakes in sequence alignments using conventional sequence alignment methods. However, a correct alignment is essential, especially for in silico methods that predict binding partners on the basis of sequence and structure. In this study, using 3D-structural information we have generated and manually corrected sequence alignments for proteins of the five ubiquitin superfold subfamilies. On the basis of this alignment, we suggest domains for which structural information will be useful to allow homology modelling. In addition, we have analysed the energetic and electrostatic properties of ubiquitin-like domains in complex with various functional binding proteins using the protein design algorithm FoldX. On the basis of an in silico alanine-scanning mutagenesis, we provide a detailed binding epitope mapping of the hotspots of the ubiquitin domain fold, involved in the interaction with different domains and proteins. Finally, we provide a consensus fingerprint sequence that identifies all sequences described to belong to the ubiquitin superfold family. It is possible that the method that we describe may be applied to other domain families sharing a similar fold but having low levels of sequence homology.

Ubiquitin and ubiquitin-like proteins as multifunctional signals

Protein ubiquitylation is a recognized signal for protein degradation. However, it is increasingly realized that ubiquitin conjugation to proteins can be used for many other purposes. Furthermore, there are many ubiquitin-like proteins that control the activities of proteins. The central structural element of these post-translational modifications is the ubiquitin superfold. A common ancestor based on this superfold has evolved to give various proteins that are involved in diverse activities in the cell.