Distributed all-atom simulations of intrinsically unstructured proteins

The Computational Biophysics Group at the Universitat Pompeu Fabra (GRIB-UPF) hosts two unique computational resources dedicated to the execution of large scale molecular dynamics (MD) simulations: (a) the ACMD molecular-dynamics software, used on standard personal computers with graphical processing units (GPUs); and (b) the GPUGRID. net computing network, supported by users distributed worldwide that volunteer GPUs for biomedical research. We leveraged these resources and developed studies, protocols and open-source software to elucidate energetics and pathways of a number of biomolecular systems, with a special focus on flexible proteins with many degrees of freedom. First, we characterized ion permeation through the bactericidal model protein Gramicidin A conducting one of the largest studies to date with the steered MD biasing methodology. Next, we addressed an open problem in structural biology, the determination of drug-protein association kinetics; we reconstructed the binding free energy, association, and dissaciociation rates of a drug like model system through a spatial decomposition and a Makov-chain analysis. The work was published in the Proceedings of the National Academy of Sciences and become one of the few landmark papers elucidating a ligand-binding pathway. Furthermore, we investigated the unstructured Kinase Inducible Domain (KID), a 28-peptide central to signalling and transcriptional response; the kinetics of this challenging system was modelled with a Markovian approach in collaboration with Frank Noe’s group at the Freie University of Berlin. The impact of the funding includes three peer-reviewed publication on high-impact journals; three more papers under review; four MD analysis components, released as open-source software; MD protocols; didactic material, and code for the hosting group.

SelenoDB 1.0 : a database of selenoprotein genes, proteins and SECIS elements

Selenoproteins are a diverse group of proteinsusually misidentified and misannotated in sequencedatabases. The presence of an in-frame UGA (stop)codon in the coding sequence of selenoproteingenes precludes their identification and correctannotation. The in-frame UGA codons are recodedto cotranslationally incorporate selenocysteine,a rare selenium-containing amino acid. The developmentof ad hoc experimental and, more recently,computational approaches have allowed the efficientidentification and characterization of theselenoproteomes of a growing number of species.Today, dozens of selenoprotein families have beendescribed and more are being discovered in recentlysequenced species, but the correct genomic annotationis not available for the majority of thesegenes. SelenoDB is a long-term project that aims toprovide, through the collaborative effort of experimentaland computational researchers, automaticand manually curated annotations of selenoproteingenes, proteins and SECIS elements. Version 1.0 ofthe database includes an initial set of eukaryoticgenomic annotations, with special emphasis on thehuman selenoproteome, for immediate inspectionby selenium researchers or incorporation into moregeneral databases. SelenoDB is freely available athttp://www.selenodb.org.

Emerging functions of myelin associated proteins during development neuronal plasticity and and neurpdegeneration.

Adult mammalian central nervous system (CNS) axons have a limited regrowth capacity following injury. Myelin-associated inhibitors (MAIs) limit axonal outgrowth and their blockage improves the regeneration of damaged fiber tracts. Three of these proteins, Nogo-A, MAG and OMgp, share two common neuronal receptors: NgR1, together with its co-receptors (p75(NTR), TROY and Lingo-1), and the recently described paired immunoglobulin-like receptor B (PirB). These proteins impair neuronal regeneration by limiting axonal sprouting. Some of the elements involved in the myelin inhibitory pathways may still be unknown, but the discovery that blocking both PirB and NgR1 activities leads to near-complete release from myelin inhibition, sheds light on one of the most competitive and intense fields of neuroregeneration study during in recent decades. In parallel with the identification and characterization of the roles and functions of these inhibitory molecules in axonal regeneration, data gathered in the field strongly suggest that most of these proteins have roles other than axonal growth inhibition. The discovery of a new group of interacting partners for myelin-associated receptors and ligands, as well as functional studies within or outside the CNS environment, highlights the potential new physiological roles for these proteins in processes such as development, neuronal homeostasis, plasticity and neurodegeneration.

egr-4, a target of EGFR signaling, is required for the formation of the brian primordial and head regeneration in planarians

During the regeneration of freshwater planarians, polarity and patterning programs play essential roles in determining whether a head or a tail regenerates at anterior or posterior-facing wounds. This decision is made very soon after amputation. The pivotal role of the Wnt/β-catenin and Hh signaling pathways in re-establishing anterior-posterior (AP) polarity has been well documented. However, the mechanisms that control the growth and differentiation of the blastema in accordance with its AP identity are less well understood. Previous studies have described a role of Smed-egfr-3, a planarian epidermal growth factor receptor, in blastema growth and differentiation. Here, we identify Smed-egr-4, a zinc-finger transcription factor belonging to the early growth response gene family, as a putative downstream target of Smed-egfr-3. Smed-egr-4 is mainly expressed in the central nervous system and its silencing inhibits anterior regeneration without affecting the regeneration of posterior regions. Single and combinatorial RNA interference to target different elements of the Wnt/β-catenin pathway, together with expression analysis of brain- and anterior-specific markers, revealed that Smed-egr-4: (1) is expressed in two phases - an early Smed-egfr-3-independent phase and a late Smed-egfr-3-dependent phase; (2) is necessary for the differentiation of the brain primordia in the early stages of regeneration; and (3) that it appears to antagonize the activity of the Wnt/β-catenin pathway to allow head regeneration. These results suggest that a conserved EGFR/egr pathway plays an important role in cell differentiation during planarian regeneration and indicate an association between early brain differentiation and the proper progression of head regeneration.

Insertion of a malE B-Galactosidase fusion protein into the envelope of Escherichia coli disrupts biogenesis of outer membrane proteins and processing of inner membrane proteins

The synthesis of a membrane-bound MalE ,B-galactosidase hybrid protein, when induced by growth of Escherichia coli on maltose, leads to inhibition of cell division and eventually a reduced rate of mass increase. In addition, the relative rate of synthesis of outer membrane proteins, but not that of inner membrane proteins, was reduced by about 50%o. Kinetic experiments demonstrated that this reduction coincided with the period of maximum synthesis of the hybrid protein (and another maltose-inducible protein, LamB). The accumulation of this abnormal protein in the envelope therefore appeared specifically to inhibit the synthesis, the assembly of outer membrane proteins, or both, indicating that the hybrid protein blocks some export site or causes the sequestration of some limiting factor(s) involved in the export process. Since the MalE protein is normally located in the periplasm, the results also suggest that the synthesis of periplasmic and outer membrane proteins may involve some steps in common. The reduced rate of synthesis of outer membrane proteins was also accompanied by the accumulation in the envelope of at least one outer membrane protein and at least two inner membrane proteins as higher-molecular-weight forms, indicating that processing (removal of the N-terminal signal sequence) was also disrupted by the presence of the hybrid protein. These results may indicate that the assembly of these membrane proteins is blocked at a relatively late step rather than at the level of primary recognition of some site by the signal sequence. In addition, the results suggest that some step common to the biogenesis of quite different kinds of envelope protein is blocked by the presence of the hybrid protein.

Competition between SOCS36E and Drk modulates Sevenless receptor tyrosine kinase activity

Modulation of signalling pathways can trigger different cellular responses, including differences in cell fate. This modulation can be achieved by controlling the pathway activity with great precision to ensure robustness and reproducibility of the specification of cell fate. The development of the photoreceptor R7 in the Drosophila melanogasterretina has become a model in which to investigate the control of cell signalling. During R7 specification, a burst of Ras small GTPase (Ras) and mitogen-activated protein kinase (MAPK) controlled by Sevenless receptor tyrosine kinase (Sev) is required. Several cells in each ommatidium express sev. However, the spatiotemporal expression of the boss ligand and the action of negative regulators of the Sev pathway will restrict the R7 fate to a single cell. The Drosophila suppressor of cytokine signalling 36E (SOCS36E) protein contains an SH2 domain and acts as a Sev signalling attenuator. By contrast, downstream of receptor kinase (Drk), the fly homolog of the mammalian Grb2 adaptor protein, which also contains an SH2 domain, acts as a positive activator of the pathway. Here, we apply the Förster resonance energy transfer (FRET) assay to transfected Drosophila S2 cells and demonstrate that Sev binds directly to either the suppressor protein SOCS36E or the adaptor protein Drk. We propose a mechanistic model in which the competition between these two proteins for binding to the same docking site results in either attenuation of the Sev transduction in cells that should not develop R7 photoreceptors or amplification of the Ras-MAPK signal only in the R7 precursor.

Unique features of odorant-binding proteins of the parasitoid wasp Nasonia vitripennis revealed by genome annotation and comparative analyses

Insects are the most diverse group of animals on the planet, comprising over 90% of all metazoan life forms, and have adapted to a wide diversity of ecosystems in nearly all environments. They have evolved highly sensitive chemical senses that are central to their interaction with their environment and to communication between individuals. Understanding the molecular bases of insect olfaction is therefore of great importance from both a basic and applied perspective. Odorant binding proteins (OBPs) are some of most abundant proteins found in insect olfactory organs, where they are the first component of the olfactory transduction cascade, carrying odorant molecules to the olfactory receptors. We carried out a search for OBPs in the genome of the parasitoid wasp Nasonia vitripennis and identified 90 sequences encoding putative OBPs. This is the largest OBP family so far reported in insects. We report unique features of the N. vitripennis OBPs, including the presence and evolutionary origin of a new subfamily of double-domain OBPs (consisting of two concatenated OBP domains), the loss of conserved cysteine residues and the expression of pseudogenes. This study also demonstrates the extremely dynamic evolution of the insect OBP family: (i) the number of different OBPs can vary greatly between species; (ii) the sequences are highly diverse, sometimes as a result of positive selection pressure with even the canonical cysteines being lost; (iii) new lineage specific domain arrangements can arise, such as the double domain OBP subfamily of wasps and mosquitoes.

Contribution of S6K1/MAPK signaling pathways in the response to oxidative stress: activation of RSK and MSK by hydrogen peroxide

Cells respond to different kind of stress through the coordinated activation of signaling pathways such as MAPK or p53. To find which molecular mechanisms are involved, we need to understand their cell adaptation. The ribosomal protein, S6 kinase 1 (S6K1), is a common downstream target of signaling by hormonal or nutritional stress. Here, we investigated the initial contribution of S6K1/MAPK signaling pathways in the cell response to oxidative stress produced by hydrogen peroxide (H2O2). To analyze S6K1 activation, we used the commercial anti-phospho-Thr389-S6K1 antibody most frequently mentioned in the bibliography. We found that this antibody detected an 80-90 kDa protein that was rapidly phosphorylated in response to H2O2 in several human cells. Unexpectedly, this phosphorylation was insensitive to both mTOR and PI3K inhibitors, and knock-down experiments showed that this protein was not S6K1. RSK and MSK proteins were candidate targets of this phosphorylation. We demonstrated that H2O2 stimulated phosphorylation of RSK and MSK kinases at residues that are homologous to Thr389 in S6K1. This phosphorylation required the activity of either p38 or ERK MAP kinases. Kinase assays showed activation of RSK and MSK by H2O2. Experiments with mouse embryonic fibroblasts from p38 animals" knockout confirmed these observations. Altogether, these findings show that the S6K1 signaling pathway is not activated under these conditions, clarify previous observations probably misinterpreted by non-specific detection of proteins RSK and MSK by the anti-phospho-Thr389-S6K1 antibody, and demonstrate the specific activation of MAPK signaling pathways through ERK/p38/RSK/MSK by H2O2.

Id Gene regulation and function in the prosensory domains of the chicken inner ear: a link between Bmp Signaling and Atoh1

Bone morphogenetic proteins (Bmps) regulate the expression of the proneural gene Atoh1 and the generation of hair cells in the developing inner ear. The present work explored the role of Inhibitor of Differentiation genes (Id1-3) in this process. The results show that Id genes are expressed in the prosensory domains of the otic vesicle, along with Bmp4 and Bmp7. Those domains exhibit high levels of the phosphorylated form of Bmp-responding R-Smads (P-Smad1,5,8), and of Bmp-dependent Smad transcriptional activity as shown by the BRE-tk-EGFP reporter. Increased Bmp signaling induces the expression of Id1-3 along with the inhibition of Atoh1. Conversely, the Bmp antagonist Noggin or the Bmp-receptor inhibitor Dorsomorphin elicit opposite effects, indicating that Bmp signaling is necessary for Id expression and Atoh1 regulation in the otocyst. The forced expression of Id3 is sufficient to reduce Atoh1 expression and to prevent the expression of hair cell differentiation markers. Together, these results suggest that Ids are part of the machinery that mediates the regulation of hair cell differentiation exerted by Bmps. In agreement with that, during hair cell differentiation Bmp4 expression, P-Smad1,5,8 levels and Id expression are downregulated from hair cells. However, Ids are also downregulated from the supporting cells which contrarily to hair cells exhibit high levels of Bmp4 expression, P-Smad1,5,8, and BRE-tk-EGFP activity, suggesting that in these cells Ids escape from Bmp/Smad signaling. The differential regulation of Ids in time and space may underlie the multiple functions of Bmp signaling during sensory organ development.

Identification of potential small molecule binding pockets on Rho family GTPases

Rho GTPases are conformational switches that control a wide variety of signaling pathways critical for eukaryotic cell development and proliferation. They represent attractive targets for drug design as their aberrant function and deregulated activity is associated with many human diseases including cancer. Extensive high-resolution structures (.100) and recent mutagenesis studies have laid the foundation for the design of new structure-based chemotherapeutic strategies. Although the inhibition of Rho signaling with drug-like compounds is an active area of current research, very little attention has been devoted to directly inhibiting Rho by targeting potential allosteric non-nucleotide binding sites. By avoiding the nucleotide binding site, compounds may minimize the potential for undesirable off-target interactions with other ubiquitous GTP and ATP binding proteins. Here we describe the application of molecular dynamics simulations, principal component analysis, sequence conservation analysis, and ensemble small-molecule fragment mapping to provide an extensive mapping of potential small-molecule binding pockets on Rho family members. Characterized sites include novel pockets in the vicinity of the conformationaly responsive switch regions as well as distal sites that appear to be related to the conformations of the nucleotide binding region. Furthermore the use of accelerated molecular dynamics simulation, an advanced sampling method that extends the accessible time-scale of conventional simulations, is found to enhance the characterization of novel binding sites when conformational changes are important for the protein mechanism.

Calmodulin Regulates Intracellular Trafficking of Epidermal Growth Factor Receptor and the MAPK Signaling Pathway

The epidermal growth factor receptor (EGFR) is a member of the tyrosine kinase receptor family involved in signal transduction and the regulation of cellular proliferation and differentiation. It is also a calmodulin-binding protein. To examine the role of calmodulin in the regulation of EGFR, the effect of calmodulin antagonist, W-13, on the intracellular trafficking of EGFR and the MAPK signaling pathway was analyzed. W-13 did not alter the internalization of EGFR but inhibited its recycling and degradation, thus causing the accumulation of EGF and EGFR in enlarged early endosomal structures. In addition, we demonstrated that W-13 stimulated the tyrosine phosphorylation of EGFR and consequent recruitment of Shc adaptor protein with EGFR, presumably through inhibition of the calmodulin-dependent protein kinase II (CaM kinase II). W-13¿mediated EGFR phosphorylation was blocked by metalloprotease inhibitor, BB94, indicating a possible involvement of shedding in this process. However, MAPK activity was decreased by W-13; dissection of this signaling pathway showed that W-13 specifically interferes with Raf-1 activity. These data are consistent with the regulation of EGFR by calmodulin at several steps of the receptor signaling and trafficking pathways.

A system-level, molecular evolutionary analysis of mammalian phototransduction

Visual perception is initiated in the photoreceptor cells of the retina via the phototransduction system.This system has shown marked evolution during mammalian divergence in such complex attributes as activation time and recovery time. We have performed a molecular evolutionary analysis of proteins involved in mammalianphototransduction in order to unravel how the action of natural selection has been distributed throughout thesystem to evolve such traits. We found selective pressures to be non-randomly distributed according to both a simple protein classification scheme and a protein-interaction network representation of the signaling pathway. Proteins which are topologically central in the signaling pathway, such as the G proteins, as well as retinoid cycle chaperones and proteins involved in photoreceptor cell-type determination, were found to be more constrained in their evolution. Proteins peripheral to the pathway, such as ion channels and exchangers, as well as the retinoid cycle enzymes, have experienced a relaxation of selective pressures. Furthermore, signals of positive selection were detected in two genes: the short-wave (blue) opsin (OPN1SW) in hominids and the rod-specific Na+/Ca2+,K+ ion exchanger (SLC24A1) in rodents. The functions of the proteins involved in phototransduction and the topology of the interactions between them have imposed non-random constraints on their evolution. Thus, in shaping or conserving system-level phototransduction traits, natural selection has targeted the underlying proteins in a concerted manner.

Extrasynaptic neurotransmission in the modulation of brain function. Focus on the striatal neuronal glial networks

Extrasynaptic neurotransmission is an important short distance form of volume transmission (VT) and describes the extracellular diffusion of transmitters and modulators after synaptic spillover or extrasynaptic release in the local circuit regions binding to and activating mainly extrasynaptic neuronal and glial receptors in the neuroglial networks of the brain. Receptor-receptor interactions in G protein-coupled receptor (GPCR) heteromers play a major role, on dendritic spines and nerve terminals including glutamate synapses, in the integrative processes of the extrasynaptic signaling. Heteromeric complexes between GPCR and ion-channel receptors play a special role in the integration of the synaptic and extrasynaptic signals. Changes in extracellular concentrations of the classical synaptic neurotransmitters glutamate and GABA found with microdialysis is likely an expression of the activity of the neuron-astrocyte unit of the brain and can be used as an index of VT-mediated actions of these two neurotransmitters in the brain. Thus, the activity of neurons may be functionally linked to the activity of astrocytes, which may release glutamate and GABA to the extracellular space where extrasynaptic glutamate and GABA receptors do exist. Wiring transmission (WT) and VT are fundamental properties of all neurons of the CNS but the balance between WT and VT varies from one nerve cell population to the other. The focus is on the striatal cellular networks, and the WT and VT and their integration via receptor heteromers are described in the GABA projection neurons, the glutamate, dopamine, 5-hydroxytryptamine (5-HT) and histamine striatal afferents, the cholinergic interneurons, and different types of GABA interneurons. In addition, the role in these networks of VT signaling of the energy-dependent modulator adenosine and of endocannabinoids mainly formed in the striatal projection neurons will be underlined to understand the communication in the striatal cellular networks

Chemokines induce axon outgrowth downstream of Hepatocyte Growth Factor and TCF/β-catenin signaling

Axon morphogenesis is a complex process regulated by a variety of secreted molecules, including morphogens and growth factors, resulting in the establishment of the neuronal circuitry. Our previous work demonstrated that growth factors [Neurotrophins (NT) and Hepatocyte Growth Factor (HGF)] signal through β-catenin during axon morphogenesis. HGF signaling promotes axon outgrowth and branching by inducing β-catenin phosphorylation at Y142 and transcriptional regulation of T-Cell Factor (TCF) target genes. Here, we asked which genes are regulated by HGF signaling during axon morphogenesis. An array screening indicated that HGF signaling elevates the expression of chemokines of the CC and CXC families. In line with this, CCL7, CCL20, and CXCL2 significantly increase axon outgrowth in hippocampal neurons. Experiments using blocking antibodies and chemokine receptor antagonists demonstrate that chemokines act downstream of HGF signaling during axon morphogenesis. In addition, qPCR data demonstrates that CXCL2 and CCL5 expression is stimulated by HGF through Met/b-catenin/TCF pathway. These results identify CC family members and CXCL2 chemokines as novel regulators of axon morphogenesis downstream of HGF signaling.

The SUD1 Gene Encodes a Putative E3 Ubiquitin Ligase and Is a Positive Regulator of 3-Hydroxy-3-Methylglutaryl Coenzyme A Reductase Activity in Arabidopsis

The 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR) enzyme catalyzes the major rate-limiting step of the mevalonic acid (MVA) pathway from which sterols and other isoprenoids are synthesized. In contrast with our extensive knowledge of the regulation of HMGR in yeast and animals, little is known about this process in plants. To identify regulatory components of the MVA pathway in plants, we performed a genetic screen for second-site suppressor mutations of the Arabidopsis thaliana highly drought-sensitive drought hypersensitive2 (dry2) mutant that shows decreased squalene epoxidase activity. We show that mutations in SUPPRESSOR OF DRY2 DEFECTS1 (SUD1) gene recover most developmental defects in dry2 through changes in HMGR activity. SUD1 encodes a putative E3 ubiquitin ligase that shows sequence and structural similarity to yeast Degradation of a factor (Doa10) and human TEB4, components of the endoplasmic reticulum-associated degradation C (ERAD-C) pathway. While in yeast and animals, the alternative ERAD-L/ERAD-M pathway regulates HMGR activity by controlling protein stability, SUD1 regulates HMGR activity without apparent changes in protein content. These results highlight similarities, as well as important mechanistic differences, among the components involved in HMGR regulation in plants, yeast, and animals.