A systems approach to investigate GPCR-mediated Ras signaling network in chemoattractant sensing

A GPCR-mediated signaling network enables a chemotactic cell to generate adaptative Ras signaling in response to a large range of concentrations of a chemoattractant. To explore potential regulatory mechanisms of GPCR-controlled Ras signaling in chemosensing, we applied a software package, Simmune, to construct detailed spatiotemporal models simulating responses of the cAR1-mediated Ras signaling network. We first determined the dynamics of G-protein activation and Ras signaling in Dictyostelium cells in response to cAMP stimulations using live-cell imaging and then constructed computation models by incorporating potential mechanisms. Using simulations, we validated the dynamics of signaling events and predicted the dynamic profiles of those events in the cAR1-mediated Ras signaling networks with defective Ras inhibitory mechanisms, such as without RasGAP, with RasGAP overexpression, or with RasGAP hyperactivation. We describe a method of using Simmune to construct spatiotemporal models of a signaling network and run computational simulations without writing mathematical equations. This approach will help biologists to develop and analyze computational models that parallel live-cell experiments.     

Computational modeling of cellular signaling processes embedded into dynamic spatial contexts

Cellular signaling processes depend on spatiotemporal distributions of molecular components. Multicolor, high-resolution microscopy permits detailed assessment of such distributions, providing input for fine-grained computational models that explore mechanisms governing dynamic assembly of multimolecular complexes and their role in shaping cellular behavior. However, it is challenging to incorporate into such models both complex molecular reaction cascades and the spatial localization of signaling components in dynamic cellular morphologies. Here we introduce an approach to address these challenges by automatically generating computational representations of complex reaction networks based on simple bimolecular interaction rules embedded into detailed, adaptive models of cellular morphology. Using examples of receptor-mediated cellular adhesion and signal-induced localized mitogen-activated protein kinase (MAPK) activation in yeast, we illustrate the capacity of this simulation technique to provide insights into cell biological processes. The modeling algorithms, implemented in a new version of the Simmune toolset, are accessible through intuitive graphical interfaces and programming libraries.     

Computing the Local Field Potential (LFP) from Integrate-and-Fire Network Models

Leaky integrate-and-fire (LIF) network models are commonly used to study how the spiking dynamics of neural networks changes with stimuli, tasks or dynamic network states. However, neurophysiological studies in vivo often rather measure the mass activity of neuronal microcircuits with the local field potential (LFP). Given that LFPs are generated by spatially separated currents across the neuronal membrane, they cannot be computed directly from quantities defined in models of point-like LIF neurons. Here, we explore the best approximation for predicting the LFP based on standard output from point-neuron LIF networks. To search for this best “LFP proxy”, we compared LFP predictions from candidate proxies based on LIF network output (e.g, firing rates, membrane potentials, synaptic currents) with “ground-truth” LFP obtained when the LIF network synaptic input currents were injected into an analogous three-dimensional (3D) network model of multi-compartmental neurons with realistic morphology, spatial distributions of somata and synapses. We found that a specific fixed linear combination of the LIF synaptic currents provided an accurate LFP proxy, accounting for most of the variance of the LFP time course observed in the 3D network for all recording locations. This proxy performed well over a broad set of conditions, including substantial variations of the neuronal morphologies. Our results provide a simple formula for estimating the time course of the LFP from LIF network simulations in cases where a single pyramidal population dominates the LFP generation, and thereby facilitate quantitative comparison between computational models and experimental LFP recordings in vivo.     

Experimental and Computational Analysis of a Large Protein Network That Controls Fat Storage Reveals the Design Principles of a Signaling Network

An approach combining genetic, proteomic, computational, and physiological analysis was used to define a protein network that regulates fat storage in budding yeast (Saccharomyces cerevisiae). A computational analysis of this network shows that it is not scale-free, and is best approximated by the Watts-Strogatz model, which generates “small-world” networks with high clustering and short path lengths. The network is also modular, containing energy level sensing proteins that connect to four output processes: autophagy, fatty acid synthesis, mRNA processing, and MAP kinase signaling. The importance of each protein to network function is dependent on its Katz centrality score, which is related both to the protein’s position within a module and to the module’s relationship to the network as a whole. The network is also divisible into subnetworks that span modular boundaries and regulate different aspects of fat metabolism. We used a combination of genetics and pharmacology to simultaneously block output from multiple network nodes. The phenotypic results of this blockage define patterns of communication among distant network nodes, and these patterns are consistent with the Watts-Strogatz model.     

Single Turnover Autophosphorylation Cycle of the PKA RIIβ Holoenzyme

To provide tight spatiotemporal signaling control, the cyclic adenosine monophosphate (cAMP)-dependent protein kinase (PKA) holoenzyme typically nucleates a macromolecular complex or a “PKA signalosome.” Using the RIIβ holoenzyme as a prototype, we show how autophosphorylation/dephosphorylation of the RIIβ subunit, as well as cAMP and metal ions, contribute to the dynamics of PKA signaling. While we showed previously that the RIIβ holoenzyme could undergo a single turnover autophosphorylation with adenosine triphosphate and magnesium (MgATP) and trap both products in the crystal lattice, we asked here whether calcium could trap an ATP:RIIβ holoenzyme since the RIIβ holoenzyme is located close to ion channels. The 2.8Å structure of an RIIβ^p ₂:C₂:(Ca₂ADP)₂ holoenzyme, supported by biochemical and biophysical data, reveals a trapped single phosphorylation event similar to MgATP. Thus, calcium can mediate a single turnover event with either ATP or adenosine-5''-(β,γ-imido)triphosphate (AMP-PNP), even though it cannot support steady-state catalysis efficiently. The holoenzyme serves as a “product trap” because of the slow off-rate of the pRIIβ subunit, which is controlled by cAMP, not by phosphorylation of the inhibitor site. By quantitatively defining the RIIβ signaling cycle, we show that release of pRIIβ in the presence of cAMP is reduced by calcium, whereas autophosphorylation at the phosphorylation site (P-site) inhibits holoenzyme reassociation with the catalytic subunit. Adding a single phosphoryl group to the preformed RIIβ holoenzyme thus creates a signaling cycle in which phosphatases become an essential partner. This previously unappreciated molecular mechanism is an integral part of PKA signaling for type II holoenzymes.     

Deciphering the Structural Basis of Eukaryotic Protein Kinase Regulation

Eukaryotic protein kinases (EPKs) regulate numerous signaling processes by phosphorylating targeted substrates through the highly conserved catalytic domain. Our previous computational studies proposed a model stating that a properly assembled nonlinear motif termed the Regulatory (R) spine is essential for catalytic activity of EPKs. Here we define the required intramolecular interactions and biochemical properties of the R-spine and the newly identified “Shell” that surrounds the R-spine using site-directed mutagenesis and various in vitro phosphoryl transfer assays using cyclic AMP-dependent protein kinase as a representative of the entire kinome. Analysis of the 172 available Apo EPK structures in the protein data bank (PDB) revealed four unique structural conformations of the R-spine that correspond with catalytic inactivation of various EPKs. Elucidating the molecular entities required for the catalytic activation of EPKs and the identification of these inactive conformations opens new avenues for the design of efficient therapeutic EPK inhibitors.     

Dynamic analysis of integrated signaling, metabolic, and regulatory networks

Extracellular cues affect signaling, metabolic, and regulatory processes to elicit cellular responses. Although intracellular signaling, metabolic, and regulatory networks are highly integrated, previous analyses have largely focused on independent processes (e.g., metabolism) without considering the interplay that exists among them. However, there is evidence that many diseases arise from multifunctional components with roles throughout signaling, metabolic, and regulatory networks. Therefore, in this study, we propose a flux balance analysis (FBA)–based strategy, referred to as integrated dynamic FBA (idFBA), that dynamically simulates cellular phenotypes arising from integrated networks. The idFBA framework requires an integrated stoichiometric reconstruction of signaling, metabolic, and regulatory processes. It assumes quasi-steady-state conditions for “fast” reactions and incorporates “slow” reactions into the stoichiometric formalism in a time-delayed manner. To assess the efficacy of idFBA, we developed a prototypic integrated system comprising signaling, metabolic, and regulatory processes with network features characteristic of actual systems and incorporating kinetic parameters based on typical time scales observed in literature. idFBA was applied to the prototypic system, which was evaluated for different environments and gene regulatory rules. In addition, we applied the idFBA framework in a similar manner to a representative module of the single-cell eukaryotic organism Saccharomyces cerevisiae. Ultimately, idFBA facilitated quantitative, dynamic analysis of systemic effects of extracellular cues on cellular phenotypes and generated comparable time-course predictions when contrasted with an equivalent kinetic model. Since idFBA solves a linear programming problem and does not require an exhaustive list of detailed kinetic parameters, it may be efficiently scaled to integrated intracellular systems that incorporate signaling, metabolic, and regulatory processes at the genome scale, such as the S. cerevisiae system presented here.     

The Convallis Rule for Unsupervised Learning in Cortical Networks

The phenomenology and cellular mechanisms of cortical synaptic plasticity are becoming known in increasing detail, but the computational principles by which cortical plasticity enables the development of sensory representations are unclear. Here we describe a framework for cortical synaptic plasticity termed the “Convallis rule”, mathematically derived from a principle of unsupervised learning via constrained optimization. Implementation of the rule caused a recurrent cortex-like network of simulated spiking neurons to develop rate representations of real-world speech stimuli, enabling classification by a downstream linear decoder. Applied to spike patterns used in in vitro plasticity experiments, the rule reproduced multiple results including and beyond STDP. However STDP alone produced poorer learning performance. The mathematical form of the rule is consistent with a dual coincidence detector mechanism that has been suggested by experiments in several synaptic classes of juvenile neocortex. Based on this confluence of normative, phenomenological, and mechanistic evidence, we suggest that the rule may approximate a fundamental computational principle of the neocortex.     

The Oral Tolerance as a Complex Network Phenomenon

The phenomenon of oral tolerance refers to a local and systemic state of tolerance induced in the gut after its exposure to innocuous antigens. Recent findings have shown the interrelationship between cellular and molecular components of oral tolerance, but its representation through a network of interactions has not been investigated. Our work aims at identifying the causal relationship of each element in an oral tolerance network, and also to propose a phenomenological model that’s capable of predicting the stochastic behavior of this network when under manipulation. We compared the changes of a “healthy” network caused by “knock-outs” (KOs) in two approaches: an analytical approach by the Perron Frobenius theory; and a computational approach, which we describe within this work in order to find numerical results for the model. Both approaches have shown the most relevant immunological components for this phenomena, that happens to corroborate the empirical results from animal models. Besides explain in a intelligible fashion how the components interacts in a complex manner, we also managed to describe and quantify the importance of KOs that hasn’t been empirically tested.     

Protein molecular function prediction by Bayesian phylogenomics

We present a statistical graphical model to infer specific molecular function for unannotated protein sequences using homology. Based on phylogenomic principles, SIFTER (Statistical Inference of Function Through Evolutionary Relationships) accurately predicts molecular function for members of a protein family given a reconciled phylogeny and available function annotations, even when the data are sparse or noisy. Our method produced specific and consistent molecular function predictions across 100 Pfam families in comparison to the Gene Ontology annotation database, BLAST, GOtcha, and Orthostrapper. We performed a more detailed exploration of functional predictions on the adenosine-5′-monophosphate/adenosine deaminase family and the lactate/malate dehydrogenase family, in the former case comparing the predictions against a gold standard set of published functional characterizations. Given function annotations for 3% of the proteins in the deaminase family, SIFTER achieves 96% accuracy in predicting molecular function for experimentally characterized proteins as reported in the literature. The accuracy of SIFTER on this dataset is a significant improvement over other currently available methods such as BLAST (75%), GeneQuiz (64%), GOtcha (89%), and Orthostrapper (11%). We also experimentally characterized the adenosine deaminase from Plasmodium falciparum, confirming SIFTER's prediction. The results illustrate the predictive power of exploiting a statistical model of function evolution in phylogenomic problems. A software implementation of SIFTER is available from the authors.     

Design of multi-specificity in protein interfaces

Interactions in protein networks may place constraints on protein interface sequences to maintain correct and avoid unwanted interactions. Here we describe a “multi-constraint” protein design protocol to predict sequences optimized for multiple criteria, such as maintaining sets of interactions, and apply it to characterize the mechanism and extent to which 20 multi-specific proteins are constrained by binding to multiple partners. We find that multi-specific binding is accommodated by at least two distinct patterns. In the simplest case, all partners share key interactions, and sequences optimized for binding to either single or multiple partners recover only a subset of native amino acid residues as optimal. More interestingly, for signaling interfaces functioning as network “hubs,” we identify a different, “multi-faceted” mode, where each binding partner prefers its own subset of wild-type residues within the promiscuous binding site. Here, integration of preferences across all partners results in sequences much more “native-like” than seen in optimization for any single binding partner alone, suggesting these interfaces are substantially optimized for multi-specificity. The two strategies make distinct predictions for interface evolution and design. Shared interfaces may be better small molecule targets, whereas multi-faceted interactions may be more “designable” for altered specificity patterns. The computational methodology presented here is generalizable for examining how naturally occurring protein sequences have been selected to satisfy a variety of positive and negative constraints, as well as for rationally designing proteins to have desired patterns of altered specificity.     

The Structural Basis for Endotoxin-induced Allosteric Regulation of the Toll-like Receptor 4 (TLR4) Innate Immune Receptor

As part of the innate immune system, Toll-like receptor 4 (TLR4) recognizes bacterial cell surface lipopolysaccharide (LPS) by forming a complex with a lipid-binding co-receptor, MD-2. In the presence of agonist, TLR4·MD-2 dimerizes to form an active receptor complex, leading to initiation of intracellular inflammatory signals. TLR4 is of great biomedical interest, but its pharmacological manipulation is complicated because even subtle variations in the structure of LPS can profoundly impact the resultant immunological response. Here, we use atomically detailed molecular simulations to gain insights into the nature of the molecular signaling mechanism. We first demonstrate that MD-2 is extraordinarily flexible. The “clamshell-like” motions of its β-cup fold enable it to sensitively match the volume of its hydrophobic cavity to the size and shape of the bound lipid moiety. We show that MD-2 allosterically transmits this conformational plasticity, in a ligand-dependent manner, to a phenylalanine residue (Phe-126) at the cavity mouth previously implicated in TLR4 activation. Remarkably, within the receptor complex, we observe spontaneous transitions between active and inactive signaling states of Phe-126, and we confirm that Phe-126 is indeed the “molecular switch” in endotoxic signaling.     

Phyletic Profiling with Cliques of Orthologs Is Enhanced by Signatures of Paralogy Relationships

New microbial genomes are sequenced at a high pace, allowing insight into the genetics of not only cultured microbes, but a wide range of metagenomic collections such as the human microbiome. To understand the deluge of genomic data we face, computational approaches for gene functional annotation are invaluable. We introduce a novel model for computational annotation that refines two established concepts: annotation based on homology and annotation based on phyletic profiling. The phyletic profiling-based model that includes both inferred orthologs and paralogs—homologs separated by a speciation and a duplication event, respectively—provides more annotations at the same average Precision than the model that includes only inferred orthologs. For experimental validation, we selected 38 poorly annotated Escherichia coli genes for which the model assigned one of three GO terms with high confidence: involvement in DNA repair, protein translation, or cell wall synthesis. Results of antibiotic stress survival assays on E. coli knockout mutants showed high agreement with our model''s estimates of accuracy: out of 38 predictions obtained at the reported Precision of 60%, we confirmed 25 predictions, indicating that our confidence estimates can be used to make informed decisions on experimental validation. Our work will contribute to making experimental validation of computational predictions more approachable, both in cost and time. Our predictions for 998 prokaryotic genomes include ∼400000 specific annotations with the estimated Precision of 90%, ∼19000 of which are highly specific—e.g. “penicillin binding,” “tRNA aminoacylation for protein translation,” or “pathogenesis”—and are freely available at http://gorbi.irb.hr/.     

The Brain's Router: A Cortical Network Model of Serial Processing in the Primate Brain

The human brain efficiently solves certain operations such as object recognition and categorization through a massively parallel network of dedicated processors. However, human cognition also relies on the ability to perform an arbitrarily large set of tasks by flexibly recombining different processors into a novel chain. This flexibility comes at the cost of a severe slowing down and a seriality of operations (100–500 ms per step). A limit on parallel processing is demonstrated in experimental setups such as the psychological refractory period (PRP) and the attentional blink (AB) in which the processing of an element either significantly delays (PRP) or impedes conscious access (AB) of a second, rapidly presented element. Here we present a spiking-neuron implementation of a cognitive architecture where a large number of local parallel processors assemble together to produce goal-driven behavior. The precise mapping of incoming sensory stimuli onto motor representations relies on a “router” network capable of flexibly interconnecting processors and rapidly changing its configuration from one task to another. Simulations show that, when presented with dual-task stimuli, the network exhibits parallel processing at peripheral sensory levels, a memory buffer capable of keeping the result of sensory processing on hold, and a slow serial performance at the router stage, resulting in a performance bottleneck. The network captures the detailed dynamics of human behavior during dual-task-performance, including both mean RTs and RT distributions, and establishes concrete predictions on neuronal dynamics during dual-task experiments in humans and non-human primates.     

A Boolean Approach to Linear Prediction for Signaling Network Modeling

The task of the DREAM4 (Dialogue for Reverse Engineering Assessments and Methods) “Predictive signaling network modeling” challenge was to develop a method that, from single-stimulus/inhibitor data, reconstructs a cause-effect network to be used to predict the protein activity level in multi-stimulus/inhibitor experimental conditions. The method presented in this paper, one of the best performing in this challenge, consists of 3 steps: 1. Boolean tables are inferred from single-stimulus/inhibitor data to classify whether a particular combination of stimulus and inhibitor is affecting the protein. 2. A cause-effect network is reconstructed starting from these tables. 3. Training data are linearly combined according to rules inferred from the reconstructed network. This method, although simple, permits one to achieve a good performance providing reasonable predictions based on a reconstructed network compatible with knowledge from the literature. It can be potentially used to predict how signaling pathways are affected by different ligands and how this response is altered by diseases.     

Bringing Statistics Up to Speed with Data in Analysis of Lymphocyte Motility

Two-photon (2P) microscopy provides immunologists with 3D video of the movement of lymphocytes in vivo. Motility parameters extracted from these videos allow detailed analysis of lymphocyte motility in lymph nodes and peripheral tissues. However, standard parametric statistical analyses such as the Student’s t-test are often used incorrectly, and fail to take into account confounds introduced by the experimental methods, potentially leading to erroneous conclusions about T cell motility. Here, we compare the motility of WT T cell versus PKCθ-/-, CARMA1-/-, CCR7-/-, and PTX-treated T cells. We show that the fluorescent dyes used to label T cells have significant effects on T cell motility, and we demonstrate the use of factorial ANOVA as a statistical tool that can control for these effects. In addition, researchers often choose between the use of “cell-based” parameters by averaging multiple steps of a single cell over time (e.g. cell mean speed), or “step-based” parameters, in which all steps of a cell population (e.g. instantaneous speed) are grouped without regard for the cell track. Using mixed model ANOVA, we show that we can maintain cell-based analyses without losing the statistical power of step-based data. We find that as we use additional levels of statistical control, we can more accurately estimate the speed of T cells as they move in lymph nodes as well as measure the impact of individual signaling molecules on T cell motility. As there is increasing interest in using computational modeling to understand T cell behavior in in vivo, these quantitative measures not only give us a better determination of actual T cell movement, they may prove crucial for models to generate accurate predictions about T cell behavior.     

Bacterial Signal Transduction by Cyclic Di-GMP and Other Nucleotide Second Messengers

The first International Symposium on c-Di-GMP Signaling in Bacteria (22 to 25 March 2015, Harnack-Haus, Berlin, Germany) brought together 131 molecular microbiologists from 17 countries to discuss recent progress in our knowledge of bacterial nucleotide second messenger signaling. While the focus was on signal input, synthesis, degradation, and the striking diversity of the modes of action of the current second messenger paradigm, i.e., cyclic di-GMP (c-di-GMP), “classics” like cAMP and (p)ppGpp were also presented, in novel facets, and more recent “newcomers,” such as c-di-AMP and c-AMP-GMP, made an impressive appearance. A number of clear trends emerged during the 30 talks, on the 71 posters, and in the lively discussions, including (i) c-di-GMP control of the activities of various ATPases and phosphorylation cascades, (ii) extensive cross talk between c-di-GMP and other nucleotide second messenger signaling pathways, and (iii) a stunning number of novel effectors for nucleotide second messengers that surprisingly include some long-known master regulators of developmental pathways. Overall, the conference made it amply clear that second messenger signaling is currently one of the most dynamic fields within molecular microbiology, with major impacts in research fields ranging from human health to microbial ecology.     

Quorum-Sensing Dysbiotic Shifts in the HIV-Infected Oral Metabiome

We implemented a Systems Biology approach using Correlation Difference Probability Network (CDPN) analysis to provide insights into the statistically significant functional differences between HIV-infected patients and uninfected individuals. The analysis correlates bacterial microbiome (“bacteriome”), fungal microbiome (“mycobiome”), and metabolome data to model the underlying biological processes comprising the Human Oral Metabiome. CDPN highlights the taxa-metabolite-taxa differences between the cohorts that frequently capture quorum-sensing modifications that reflect communication disruptions in the dysbiotic HIV cohort. The results also highlight the significant role of cyclic mono and dipeptides as quorum-sensing (QS) mediators between oral bacteria and fungal genus. The developed CDPN approach allowed us to model the interactions of taxa and key metabolites, and hypothesize their possible contribution to the etiology of Oral Candidiasis (OC).