Expression of Robo/Slit and Semaphorin/Plexin/Neuropilin family members in the developing hypothalamic paraventricular and supraoptic nuclei

The hypothalamic paraventricular nucleus (PVN) and supraoptic nucleus (SON) contain neuroendocrine cells that modulate pituitary secretion to maintain homeostasis. These two nuclei have a common developmental origin but they eventually form at locations distant from each other. Little is known about the molecular cues that direct the segregation of PVN and SON. As a means to identify potential factors, we have documented expression patterns of genes with known guidance roles in neural migration. Here, we focus on two groups of ligand/receptor families classified to mediate chemo-repulsion of neurons and their axons: the Slit/Robo and the Semaphorin/Plexin/Neuropilin families. Their dynamic expression patterns within and around the common PVN/SON progenitor as well as the mature PVN and SON may provide a framework for understanding the formation of these two important nuclei.     

Semaphorin 3E Initiates Antiangiogenic Signaling through Plexin D1 by Regulating Arf6 and R-Ras

Recent studies revealed that a class III semaphorin, semaphorin 3E (Sema3E), acts through a single-pass transmembrane receptor, plexin D1, to provide a repulsive cue for plexin D1-expressing endothelial cells, thus providing a highly conserved and developmentally regulated signaling system guiding the growth of blood vessels. We show here that Sema3E acts as a potent inhibitor of adult and tumor-induced angiogenesis. Activation of plexin D1 by Sema3E causes the rapid disassembly of integrin-mediated adhesive structures, thereby inhibiting endothelial cell adhesion to the extracellular matrix (ECM) and causing the retraction of filopodia in endothelial tip cells. Sema3E acts on plexin D1 to initiate a two-pronged mechanism involving R-Ras inactivation and Arf6 stimulation, which affect the status of activation of integrins and their intracellular trafficking, respectively. Ultimately, our present study provides a molecular framework for antiangiogenesis signaling, thus impinging on a myriad of pathological conditions that are characterized by aberrant increase in neovessel formation, including cancer.Pathological angiogenesis characterizes numerous human diseases, ranging from chronic inflammation, atherosclerosis, diabetic retinopathy, and age-related macular degeneration to cancer (5, 11, 30). Thus, elucidating the mechanisms underlying normal and aberrant blood vessel growth may provide new therapeutic options for many highly prevalent disease conditions. Ultimately, normal angiogenesis results from a precise balance between pro- and antiangiogenic mediators. Among the former, the family of vascular endothelial growth factors (VEGFs), basic fibroblastic growth factor (bFGF), sphingosine-1-phosphate (S1P), and the chemokines interleukin-8/CXCL8 and SDF-1/CXCL12 and their receptors are some of the most widely investigated (reviewed in references 3, 5, 8, and 17). The best-known angiogenesis inhibitors are proteolytic cleavage products of extracellular matrix (ECM) or serum components, such as endostatin, angiostatin, arresten, and tumstatin (reviewed in references 11 and 20). Antiangiogenic cytokines have also been described, including interferons and certain interleukins, which appear to act indirectly by limiting the expression of proangiogenic mediators or inducing antiangiogenic molecules (reviewed in references 11 and 20). In contrast, there are few known developmentally regulated, naturally occurring antiangiogenic molecules, which include platelet factor 4 (18), thrombospondin 1 (14), and pigment epithelium-derived factor (PEDF) (9). Their precise mechanism of action is not fully understood, thus limiting the ability to design new molecularly based antiangiogenic strategies.Emerging evidence suggests that proteins involved in transmitting axonal guidance cues, including members of the netrin, slit, eph, and semaphorin families, also play a critical role in blood vessel guidance during physiological and pathological blood vessel development (6). For example, multiple secreted class III semaphorins, which regulate developmental axonal growth (23, 27), are now known to act through their receptors, the A family plexins (plexins A1, A2, and A3), and their coreceptors, neuropilin 1 and neuropilin 2, to initiate pro- and antiangiogenic responses (reviewed in references 6 and 19). However, neuropilins also act as coreceptors for multiple angiogenic factors, such as VEGF, thus limiting our ability to distinguish the downstream signaling events initiated by semaphorins from those resulting from their interplay with endothelial growth and motility factors (19). In this regard, recent studies revealed that a class III semaphorin, semaphorin 3E (Sema3E), acts through a single-pass transmembrane receptor, plexin D1, independently of neuropilins to control endothelial cell (EC) positioning and patterning of the developing vasculature (13, 15). These findings prompted us to explore whether Sema3E acts as a natural antiangiogenic molecule and, if so, to investigate the underlying molecular mechanism. Indeed, we show here that Sema3E is a potent inhibitor of adult and tumor-induced angiogenesis. Sema3E causes filopodial retraction in endothelial tip cells and inhibits endothelial cell adhesion by disrupting integrin-mediated adhesive structures. At the molecular level, this process involves the stimulation of plexin D1 by Sema3E, which in turn interferes with R-Ras function and leads to the rapid activation of Arf6, thus revealing a novel physiological antiangiogenic signaling route.     

Neuropilin 1 Directly Interacts with Fer Kinase to Mediate Semaphorin 3A-induced Death of Cortical Neurons

Neuropilins (NRPs) are receptors for the major chemorepulsive axonal guidance cue semaphorins (Sema). The interaction of Sema3A/NRP1 during development leads to the collapse of growth cones. Here we show that Sema3A also induces death of cultured cortical neurons through NRP1. A specific NRP1 inhibitory peptide ameliorated Sema3A-evoked cortical axonal retraction and neuronal death. Moreover, Sema3A was also involved in cerebral ischemia-induced neuronal death. Expression levels of Sema3A and NRP1, but not NRP2, were significantly increased early during brain reperfusion following transient focal cerebral ischemia. NRP1 inhibitory peptide delivered to the ischemic brain was potently neuroprotective and prevented the loss of motor functions in mice. The integrity of the injected NRP1 inhibitory peptide into the brain remained unchanged, and the intact peptide permeated the ischemic hemisphere of the brain as determined using MALDI-MS-based imaging. Mechanistically, NRP1-mediated axonal collapse and neuronal death is through direct and selective interaction with the cytoplasmic tyrosine kinase Fer. Fer RNA interference effectively attenuated Sema3A-induced neurite retraction and neuronal death in cortical neurons. More importantly, down-regulation of Fer expression using Fer-specific RNA interference attenuated cerebral ischemia-induced brain damage. Together, these studies revealed a previously unknown function of NRP1 in signaling Sema3A-evoked neuronal death through Fer in cortical neurons.     

Coarse-Grained Simulations of the HIV-1 Matrix Protein Anchoring: Revisiting Its Assembly on Membrane Domains

In the accepted model for human immunodeficiency virus preassembly in infected host cells, the anchoring to the intracellular leaflet of the membrane of the matrix domain (MA) that lies at the N-terminus of the viral Gag protein precursor appears to be one of the crucial steps for particle assembly. In this study, we simulated the membrane anchoring of human immunodeficiency virus-1 myristoylated MA protein using a coarse-grained representation of both the protein and the membrane. Our calculations first suggest that the myristoyl group could spontaneously release from its initial hydrophobic pocket before MA protein interacts with the lipid membrane. All-atom simulations confirmed this possibility with a related energy cost estimated to be ∼5 kcal.mol⁻¹. The phosphatidylinositol (4,5) bisphosphate (PI(4,5)P₂) head binds preferentially to the MA highly basic region as described in available NMR data, but interestingly without flipping of its 2′ acyl chain into the MA protein. Moreover, MA was able to confine PI(4,5)P₂ lipids all around its molecular surface after having found a stable orientation at the membrane surface. Our results suggest that this orientation is dependent on Myr anchoring and that this confinement induces a lateral segregation of PI(4,5)P₂ in domains. This is consistent with a PI(4,5)P₂ enrichment of the virus envelope as compared to the host cell membrane.     

Coarse-Grained Simulations of the HIV-1 Matrix Protein Anchoring: Revisiting Its Assembly on Membrane Domains

In the accepted model for human immunodeficiency virus preassembly in infected host cells, the anchoring to the intracellular leaflet of the membrane of the matrix domain (MA) that lies at the N-terminus of the viral Gag protein precursor appears to be one of the crucial steps for particle assembly. In this study, we simulated the membrane anchoring of human immunodeficiency virus-1 myristoylated MA protein using a coarse-grained representation of both the protein and the membrane. Our calculations first suggest that the myristoyl group could spontaneously release from its initial hydrophobic pocket before MA protein interacts with the lipid membrane. All-atom simulations confirmed this possibility with a related energy cost estimated to be ∼5 kcal.mol⁻¹. The phosphatidylinositol (4,5) bisphosphate (PI(4,5)P₂) head binds preferentially to the MA highly basic region as described in available NMR data, but interestingly without flipping of its 2′ acyl chain into the MA protein. Moreover, MA was able to confine PI(4,5)P₂ lipids all around its molecular surface after having found a stable orientation at the membrane surface. Our results suggest that this orientation is dependent on Myr anchoring and that this confinement induces a lateral segregation of PI(4,5)P₂ in domains. This is consistent with a PI(4,5)P₂ enrichment of the virus envelope as compared to the host cell membrane.     

Peptide Nanopores and Lipid Bilayers: Interactions by Coarse-Grained Molecular-Dynamics Simulations

A set of 49 protein nanopore-lipid bilayer systems was explored by means of coarse-grained molecular-dynamics simulations to study the interactions between nanopores and the lipid bilayers in which they are embedded. The seven nanopore species investigated represent the two main structural classes of membrane proteins (α-helical and β-barrel), and the seven different bilayer systems range in thickness from ∼28 to ∼43 Å. The study focuses on the local effects of hydrophobic mismatch between the nanopore and the lipid bilayer. The effects of nanopore insertion on lipid bilayer thickness, the dependence between hydrophobic thickness and the observed nanopore tilt angle, and the local distribution of lipid types around a nanopore in mixed-lipid bilayers are all analyzed. Different behavior for nanopores of similar hydrophobic length but different geometry is observed. The local lipid bilayer perturbation caused by the inserted nanopores suggests possible mechanisms for both lipid bilayer-induced protein sorting and protein-induced lipid sorting. A correlation between smaller lipid bilayer thickness (larger hydrophobic mismatch) and larger nanopore tilt angle is observed and, in the case of larger hydrophobic mismatches, the simulated tilt angle distribution seems to broaden. Furthermore, both nanopore size and key residue types (e.g., tryptophan) seem to influence the level of protein tilt, emphasizing the reciprocal nature of nanopore-lipid bilayer interactions.     

Coarse-Grained Brownian Dynamics Simulations of the 10-23 DNAzyme

Deoxyribozymes (DNAzymes) are single-stranded DNA that catalyze nucleic acid biochemistry. Although a number of DNAzymes have been discovered by in vitro selection, the relationship between their tertiary structure and function remains unknown. We focus here on the well-studied 10-23 DNAzyme, which cleaves mRNA with a catalytic efficiency approaching that of RNase A. Using coarse-grained Brownian dynamics simulations, we find that the DNAzyme bends its substrate away from the cleavage point, exposing the reactive site and buckling the DNAzyme catalytic core. This hypothesized transition state provides microscopic insights into experimental observations concerning the size of the DNAzyme/substrate complex, the impact of the recognition arm length, and the sensitivity of the enzymatic activity to point mutations of the catalytic core. Upon cleaving the pertinent backbone bond in the substrate, we find that the catalytic core of the DNAzyme unwinds and the overall complex rapidly extends, in agreement with experiments on the related 8-17 DNAzyme. The results presented here provide a starting point for interpreting experimental data on DNAzyme kinetics, as well as developing more detailed simulation models. The results also demonstrate the limitations of using a simple physical model to understand the role of point mutations.     

Neuropilin 1: function and therapeutic potential in cancer

Neuropilin 1 (NRP1) is a transmembrane glycoprotein that acts as a co-receptor for a number of extracellular ligands including class III/IV semaphorins, certain isoforms of vascular endothelial growth factor and transforming growth factor beta. An exact understanding of the role of NRP1 in the immune system has been obscured by the differences in NRP1 expression observed between mice and humans. In mice, NRP1 is selectively expressed on thymic-derived Tregs and greatly enhances immunosuppressive function. In humans, NRP1 is expressed on plasmacytoid dendritic cells (pDCs) where it aids in priming immune responses and on a subset of T regulatory cells (Tregs) isolated from secondary lymph nodes. Preliminary studies that show NRP1 expression on T cells confers enhanced immunosuppressive activity. However, the mechanism by which this activity is mediated remains unclear. NRP1 expression has also been identified on activated T cells and Tregs isolated from inflammatory microenvironments, suggesting NRP1 might represent a novel T cell activation marker. Of clinical interest, NRP1 may enhance Treg tumour infiltration and a decrease in NRP1+ Tregs correlates with successful chemotherapy, suggesting a specific role for NRP1 in cancer pathology. As a therapeutic target, NRP1 allows simultaneous targeting of NRP1-expressing tumour vasculature, NRP1+ Tregs and pDCs. With the development of anti-NRP1 monoclonal antibodies and cell-penetrating peptides, NRP1 represents a promising new target for cancer therapies. This paper reviews current knowledge on the role and function of NRP1 in Tregs and pDCs, both in physiological and cancer settings, as well as its potential as a therapeutic target in cancer.     

Molecular Dynamics Simulations Reveal the HIV-1 Vpu Transmembrane Protein to Form Stable Pentamers

The human immunodeficiency virus type I (HIV-1) Vpu protein is 81 residues long and has two cytoplasmic and one transmembrane (TM) helical domains. The TM domain oligomerizes to form a monovalent cation selective ion channel and facilitates viral release from host cells. Exactly how many TM domains oligomerize to form the pore is still not understood, with experimental studies indicating the existence of a variety of oligomerization states. In this study, molecular dynamics (MD) simulations were performed to investigate the propensity of the Vpu TM domain to exist in tetrameric, pentameric, and hexameric forms. Starting with an idealized α-helical representation of the TM domain, a thorough search for the possible orientations of the monomer units within each oligomeric form was carried out using replica-exchange MD simulations in an implicit membrane environment. Extensive simulations in a fully hydrated lipid bilayer environment on representative structures obtained from the above approach showed the pentamer to be the most stable oligomeric state, with interhelical van der Waals interactions being critical for stability of the pentamer. Atomic details of the factors responsible for stable pentamer structures are presented. The structural features of the pentamer models are consistent with existing experimental information on the ion channel activity, existence of a kink around the Ile17, and the location of tetherin binding residues. Ser23 is proposed to play an important role in ion channel activity of Vpu and possibly in virus propagation.     

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.     

Molecular Dynamics Simulations of Micelle Formation around Dimeric Glycophorin A Transmembrane Helices

Insertion and formation of membrane proteins involves the interaction of protein helices with one another in lipid environments. Researchers have studied glycophorin A (GpA) transmembrane helices embedded in sodium dodecyl sulfate (SDS) micelles to identify contacts significant for helix dimerization. However, a detailed picture of the conformation and dynamics of the GpA-SDS system cannot be obtained solely through experiment. Molecular dynamics simulations of SDS and a GpA dimer can provide an atomic-level picture of SDS aggregation and helix association. We report 2.5-ns simulations of GpA wild-type and mutants in a preformed micelle as well as a 32-ns simulation showing the formation of a complete micelle around wild-type GpA from an initially random placement of SDS molecules in an aqueous environment. In the latter case, an initial instability of GpA helices in water is reversed after the helices become surrounded by SDS. The properties of the spontaneously formed micelle surrounding the GpA are indistinguishable from those of the preformed micelle surrounding the GpA dimer.     

Probing the Disordered Domain of the Nuclear Pore Complex through Coarse-Grained Molecular Dynamics Simulations

The distribution of disordered proteins (FG-nups) that line the transport channel of the nuclear pore complex (NPC) is investigated by means of coarse-grained molecular dynamics simulations. A one-bead-per-amino-acid model is presented that accounts for the hydrophobic/hydrophilic and electrostatic interactions between different amino acids, polarity of the solvent, and screening of free ions. The results indicate that the interaction of the FG-nups forms a high-density, doughnut-like distribution inside the NPC, which is rich in FG-repeats. We show that the obtained distribution is encoded in the amino-acid sequence of the FG-nups and is driven by both electrostatic and hydrophobic interactions. To explore the relation between structure and function, we have systematically removed different combinations of FG-nups from the pore to simulate inviable and viable NPCs that were previously studied experimentally. The obtained density distributions show that the maximum density of the FG-nups inside the pore does not exceed 185 mg/mL in the inviable NPCs, whereas for the wild-type and viable NPCs, this value increases to 300 mg/mL. Interestingly, this maximum density is not correlated to the total mass of the FG-nups, but depends sensitively on the specific combination of essential Nups located in the central plane of the NPC.     

Curvature Generation and Pressure Profile Modulation in Membrane by Lysolipids: Insights from Coarse-Grained Simulations

Although many membrane additives are known to modulate the activities of membrane proteins via perturbing the properties of lipid membrane, the underlying mechanism is often not precisely understood. In this study, we investigate the impact of asymmetrically incorporating single-tailed lysophosphatidylcholine (LPC) into a membrane bilayer using coarse-grained molecular dynamics simulations. Using a simple computational protocol designed to approximately mimic a micropipette setting, we show that asymmetric incorporation of LPC can lead to significant curvature in a bilayer. Detailed analysis of geometrical and mechanical properties (pressure profile) of the resulting mound structure indicates that the degree of pressure profile perturbation is determined not by the local curvature per se but by the packing of lipid headgroups (i.e., area-per-lipid). The findings help provide a concrete basis for understanding the activation mechanism of mechanosensitive channels by asymmetric incorporation of LPC into membrane patches in patch-clamp experiments. The calculated local pressure profiles are valuable to the construction of realistic membrane models for the analysis of mechanosensation in a continuum mechanics framework.     

Probing the Disordered Domain of the Nuclear Pore Complex through Coarse-Grained Molecular Dynamics Simulations

The distribution of disordered proteins (FG-nups) that line the transport channel of the nuclear pore complex (NPC) is investigated by means of coarse-grained molecular dynamics simulations. A one-bead-per-amino-acid model is presented that accounts for the hydrophobic/hydrophilic and electrostatic interactions between different amino acids, polarity of the solvent, and screening of free ions. The results indicate that the interaction of the FG-nups forms a high-density, doughnut-like distribution inside the NPC, which is rich in FG-repeats. We show that the obtained distribution is encoded in the amino-acid sequence of the FG-nups and is driven by both electrostatic and hydrophobic interactions. To explore the relation between structure and function, we have systematically removed different combinations of FG-nups from the pore to simulate inviable and viable NPCs that were previously studied experimentally. The obtained density distributions show that the maximum density of the FG-nups inside the pore does not exceed 185 mg/mL in the inviable NPCs, whereas for the wild-type and viable NPCs, this value increases to 300 mg/mL. Interestingly, this maximum density is not correlated to the total mass of the FG-nups, but depends sensitively on the specific combination of essential Nups located in the central plane of the NPC.     

Insight into Signal Transduction: Structural Alterations in Transmembrane Helices Probed by Multi-1 ns Molecular Dynamics Simulations

Abstract

The hypothesis of structural alteration in transmembrane helices for signal transduction process is viewed by molecular dynamics simulation techniques. For the c-erbB-2 transmembrane domain involved in oncogenicity, the occurrence of conformational changes has been previously described as transition from the α to π helix. This dynamical feature is thoroughly analyzed for the wild phenotype and oncogenic sequences from a series of 18 simulations carried out on one nanosecond time scale. We show that these structural events do not depend upon the conditions of simulations like force field or starting helix coordinates. We demonstrate that the oncogenic mutations Val659 Glu, Gin and Asp do not prevent the transition. Furthermore, we show that β branched residues, in conjunction with Gly residues in the c-erbB-2 sequence, act as destabilizers for the α helix structure, π deformations are tightly related to other local structural motifs found in soluble and membrane proteins. These structural alterations are discussed in term of structure-activity relationships for the c-erbB-2 activating mechanism mediated by transmembrane domain dimerization.