ALBERT

Description: by Johnatan Aljadeff, Ronen Segev, Michael J. Berry, Tatyana O. Sharpee Many biological systems perform computations on inputs that have very large dimensionality. Determining the relevant input combinations for a particular computation is often key to understanding its function. A common way to find the relevant input dimensions is to examine the difference in variance between the input distribution and the distribution of inputs associated with certain outputs. In systems neuroscience, the corresponding method is known as spike-triggered covariance (STC). This method has been highly successful in characterizing relevant input dimensions for neurons in a variety of sensory systems. So far, most studies used the STC method with weakly correlated Gaussian inputs. However, it is also important to use this method with inputs that have long range correlations typical of the natural sensory environment. In such cases, the stimulus covariance matrix has one (or more) outstanding eigenvalues that cannot be easily equalized because of sampling variability. Such outstanding modes interfere with analyses of statistical significance of candidate input dimensions that modulate neuronal outputs. In many cases, these modes obscure the significant dimensions. We show that the sensitivity of the STC method in the regime of strongly correlated inputs can be improved by an order of magnitude or more. This can be done by evaluating the significance of dimensions in the subspace orthogonal to the outstanding mode(s). Analyzing the responses of retinal ganglion cells probed with Gaussian noise, we find that taking into account outstanding modes is crucial for recovering relevant input dimensions for these neurons.

Description: by Alex Bateman, Janet Kelso, Daniel Mietchen, Geoff Macintyre, Tomás Di Domenico, Thomas Abeel, Darren W. Logan, Predrag Radivojac, Burkhard Rost

Description: by Emek Demir, Özgün Babur, Igor Rodchenkov, Bülent Arman Aksoy, Ken I. Fukuda, Benjamin Gross, Onur Selçuk Sümer, Gary D. Bader, Chris Sander A rapidly growing corpus of formal, computable pathway information can be used to answer important biological questions including finding non-trivial connections between cellular processes, identifying significantly altered portions of the cellular network in a disease state and building predictive models that can be used for precision medicine. Due to its complexity and fragmented nature, however, working with pathway data is still difficult. We present Paxtools, a Java library that contains algorithms, software components and converters for biological pathways represented in the standard BioPAX language. Paxtools allows scientists to focus on their scientific problem by removing technical barriers to access and analyse pathway information. Paxtools can run on any platform that has a Java Runtime Environment and was tested on most modern operating systems. Paxtools is open source and is available under the Lesser GNU public license (LGPL), which allows users to freely use the code in their software systems with a requirement for attribution. Source code for the current release (4.2.0) can be found in Software S1. A detailed manual for obtaining and using Paxtools can be found in Protocol S1. The latest sources and release bundles can be obtained from biopax.org/paxtools.

Description: by Noa Liscovitch, Gal Chechik The transcriptome of the brain changes during development, reflecting processes that determine functional specialization of brain regions. We analyzed gene expression, measured using in situ hybridization across the full developing mouse brain, to quantify functional specialization of brain regions. Surprisingly, we found that during the time that the brain becomes anatomically regionalized in early development, transcription specialization actually decreases reaching a low, “neurotypic”, point around birth. This decrease of specialization is brain-wide, and mainly due to biological processes involved in constructing brain circuitry. Regional specialization rises again during post-natal development. This effect is largely due to specialization of plasticity and neural activity processes. Post-natal specialization is particularly significant in the cerebellum, whose expression signature becomes increasingly different from other brain regions. When comparing mouse and human expression patterns, the cerebellar post-natal specialization is also observed in human, but the regionalization of expression in the human Thalamus and Cortex follows a strikingly different profile than in mouse.

Description: by Linhui Hao, Qiuling He, Zhishi Wang, Mark Craven, Michael A. Newton, Paul Ahlquist Systematic, genome-wide RNA interference (RNAi) analysis is a powerful approach to identify gene functions that support or modulate selected biological processes. An emerging challenge shared with some other genome-wide approaches is that independent RNAi studies often show limited agreement in their lists of implicated genes. To better understand this, we analyzed four genome-wide RNAi studies that identified host genes involved in influenza virus replication. These studies collectively identified and validated the roles of 614 cell genes, but pair-wise overlap among the four gene lists was only 3% to 15% (average 6.7%). However, a number of functional categories were overrepresented in multiple studies. The pair-wise overlap of these enriched-category lists was high, ∼19%, implying more agreement among studies than apparent at the gene level. Probing this further, we found that the gene lists implicated by independent studies were highly connected in interacting networks by independent functional measures such as protein-protein interactions, at rates significantly higher than predicted by chance. We also developed a general, model-based approach to gauge the effects of false-positive and false-negative factors and to estimate, from a limited number of studies, the total number of genes involved in a process. For influenza virus replication, this novel statistical approach estimates the total number of cell genes involved to be ∼2,800. This and multiple other aspects of our experimental and computational results imply that, when following good quality control practices, the low overlap between studies is primarily due to false negatives rather than false-positive gene identifications. These results and methods have implications for and applications to multiple forms of genome-wide analysis.

Description: by Hassan Masum, Aarthi Rao, Benjamin M. Good, Matthew H. Todd, Aled M. Edwards, Leslie Chan, Barry A. Bunin, Andrew I. Su, Zakir Thomas, Philip E. Bourne

Description: by Dilek Eren, Burak Alakent Signals created by local perturbations are known to propagate long distances through proteins via backbone connectivity and nonbonded interactions. In the current study, signal propagation from the flexible ligand binding loop to the rest of Protein Tyrosine Phosphatase 1B (PTP1B) was investigated using frequency response techniques. Using restrained Targeted Molecular Dynamics (TMD) potential on WPD and R loops, PTP1B was driven between its crystal structure conformations at different frequencies. Propagation of the local perturbation signal was manifested via peaks at the fundamental frequency and upper harmonics of 1/ f distributed spectral density of atomic variables, such as C α atoms, dihedral angles, or polar interaction distances. Frequency of perturbation was adjusted high enough (simulation length 〉∼10×period of a perturbation cycle) not to be clouded by random diffusional fluctuations, and low enough (

Description: by Dipak Barua, William S. Hlavacek In colorectal cancer cells, APC, a tumor suppressor protein, is commonly expressed in truncated form. Truncation of APC is believed to disrupt degradation of β—catenin, which is regulated by a multiprotein complex called the destruction complex. The destruction complex comprises APC, Axin, β—catenin, serine/threonine kinases, and other proteins. The kinases and , which are recruited by Axin, mediate phosphorylation of β—catenin, which initiates its ubiquitination and proteosomal degradation. The mechanism of regulation of β—catenin degradation by the destruction complex and the role of truncation of APC in colorectal cancer are not entirely understood. Through formulation and analysis of a rule-based computational model, we investigated the regulation of β—catenin phosphorylation and degradation by APC and the effect of APC truncation on function of the destruction complex. The model integrates available mechanistic knowledge about site-specific interactions and phosphorylation of destruction complex components and is consistent with an array of published data. We find that the phosphorylated truncated form of APC can outcompete Axin for binding to β—catenin, provided that Axin is limiting, and thereby sequester β—catenin away from Axin and the Axin-recruited kinases and . Full-length APC also competes with Axin for binding to β—catenin; however, full-length APC is able, through its SAMP repeats, which bind Axin and which are missing in truncated oncogenic forms of APC, to bring β—catenin into indirect association with Axin and Axin-recruited kinases. Because our model indicates that the positive effects of truncated APC on β—catenin levels depend on phosphorylation of APC, at the first 20-amino acid repeat, and because phosphorylation of this site is mediated by , we suggest that is a potential target for therapeutic intervention in colorectal cancer. Specific inhibition of is predicted to limit binding of β—catenin to truncated APC and thereby to reverse the effect of APC truncation.

Description: by Geoff Macintyre, Magali Michaut, Thomas Abeel The International Society for Computational Biology (ISCB) Student Council was launched in 2004 to facilitate interaction between young scientists in the fields of bioinformatics and computational biology. Since then, the Student Council has successfully run events and programs to promote the development of the next generation of computational biologists. However, in its early years, the Student Council faced a major challenge, in that students from different geographical regions had different needs; no single activity or event could address the needs of all students. To overcome this challenge, the Student Council created the Regional Student Group (RSG) program. The program consists of locally organised and run student groups that address the specific needs of students in their region. These groups usually encompass a given country, and, via affiliation with the international Student Council, are provided with financial support, organisational support, and the ability to share information with other RSGs. In the last five years, RSGs have been created all over the world and organised activities that have helped develop dynamic bioinformatics student communities. In this article series, we present common themes emerging from RSG initiatives, explain their goals, and highlight the challenges and rewards through specific examples. This article, the first in the series, introduces the Student Council and provides a high-level overview of RSG activities. Our hope is that the article series will be a valuable source of information and inspiration for initiating similar activities in other regions and scientific communities.

Description: by Ben Mitchinson, Tony J. Prescott Spatial attention is most often investigated in the visual modality through measurement of eye movements, with primates, including humans, a widely-studied model. Its study in laboratory rodents, such as mice and rats, requires different techniques, owing to the lack of a visual fovea and the particular ethological relevance of orienting movements of the snout and the whiskers in these animals. In recent years, several reliable relationships have been observed between environmental and behavioural variables and movements of the whiskers, but the function of these responses, as well as how they integrate, remains unclear. Here, we propose a unifying abstract model of whisker movement control that has as its key variable the region of space that is the animal's current focus of attention, and demonstrate, using computer-simulated behavioral experiments, that the model is consistent with a broad range of experimental observations. A core hypothesis is that the rat explicitly decodes the location in space of whisker contacts and that this representation is used to regulate whisker drive signals. This proposition stands in contrast to earlier proposals that the modulation of whisker movement during exploration is mediated primarily by reflex loops. We go on to argue that the superior colliculus is a candidate neural substrate for the siting of a head-centred map guiding whisker movement, in analogy to current models of visual attention. The proposed model has the potential to offer a more complete understanding of whisker control as well as to highlight the potential of the rodent and its whiskers as a tool for the study of mammalian attention.

Description: by Audrey Salles, Cyrille Billaudeau, Arnauld Sergé, Anne-Marie Bernard, Marie-Claire Phélipot, Nicolas Bertaux, Mathieu Fallet, Pierre Grenot, Didier Marguet, Hai-Tao He, Yannick Hamon We introduce a series of experimental procedures enabling sensitive calcium monitoring in T cell populations by confocal video-microscopy. Tracking and post-acquisition analysis was performed using Methods for Automated and Accurate Analysis of Cell Signals (MAAACS), a fully customized program that associates a high throughput tracking algorithm, an intuitive reconnection routine and a statistical platform to provide, at a glance, the calcium barcode of a population of individual T-cells. Combined with a sensitive calcium probe, this method allowed us to unravel the heterogeneity in shape and intensity of the calcium response in T cell populations and especially in naive T cells, which display intracellular calcium oscillations upon stimulation by antigen presenting cells.

Description: by Jianhui Song, Sheung Chun Ng, Peter Tompa, Kevin A. W. Lee, Hue Sun Chan Molecular recognition by intrinsically disordered proteins (IDPs) commonly involves specific localized contacts and target-induced disorder to order transitions. However, some IDPs remain disordered in the bound state, a phenomenon coined “fuzziness”, often characterized by IDP polyvalency, sequence-insensitivity and a dynamic ensemble of disordered bound-state conformations. Besides the above general features, specific biophysical models for fuzzy interactions are mostly lacking. The transcriptional activation domain of the Ewing's Sarcoma oncoprotein family (EAD) is an IDP that exhibits many features of fuzziness, with multiple EAD aromatic side chains driving molecular recognition. Considering the prevalent role of cation-π interactions at various protein-protein interfaces, we hypothesized that EAD-target binding involves polycation- π contacts between a disordered EAD and basic residues on the target. Herein we evaluated the polycation-π hypothesis via functional and theoretical interrogation of EAD variants. The experimental effects of a range of EAD sequence variations, including aromatic number, aromatic density and charge perturbations, all support the cation-π model. Moreover, the activity trends observed are well captured by a coarse-grained EAD chain model and a corresponding analytical model based on interaction between EAD aromatics and surface cations of a generic globular target. EAD-target binding, in the context of pathological Ewing's Sarcoma oncoproteins, is thus seen to be driven by a balance between EAD conformational entropy and favorable EAD-target cation-π contacts. Such a highly versatile mode of molecular recognition offers a general conceptual framework for promiscuous target recognition by polyvalent IDPs.

Description: by Danielle S. Bassett, Nicholas F. Wymbs, M. Puck Rombach, Mason A. Porter, Peter J. Mucha, Scott T. Grafton As a person learns a new skill, distinct synapses, brain regions, and circuits are engaged and change over time. In this paper, we develop methods to examine patterns of correlated activity across a large set of brain regions. Our goal is to identify properties that enable robust learning of a motor skill. We measure brain activity during motor sequencing and characterize network properties based on coherent activity between brain regions. Using recently developed algorithms to detect time-evolving communities, we find that the complex reconfiguration patterns of the brain's putative functional modules that control learning can be described parsimoniously by the combined presence of a relatively stiff temporal core that is composed primarily of sensorimotor and visual regions whose connectivity changes little in time and a flexible temporal periphery that is composed primarily of multimodal association regions whose connectivity changes frequently. The separation between temporal core and periphery changes over the course of training and, importantly, is a good predictor of individual differences in learning success. The core of dynamically stiff regions exhibits dense connectivity, which is consistent with notions of core-periphery organization established previously in social networks. Our results demonstrate that core-periphery organization provides an insightful way to understand how putative functional modules are linked. This, in turn, enables the prediction of fundamental human capacities, including the production of complex goal-directed behavior.

Description: by Stephen A. Smith, Joseph W. Brown, Cody E. Hinchliff Phylogenetic trees are used to analyze and visualize evolution. However, trees can be imperfect datatypes when summarizing multiple trees. This is especially problematic when accommodating for biological phenomena such as horizontal gene transfer, incomplete lineage sorting, and hybridization, as well as topological conflict between datasets. Additionally, researchers may want to combine information from sets of trees that have partially overlapping taxon sets. To address the problem of analyzing sets of trees with conflicting relationships and partially overlapping taxon sets, we introduce methods for aligning, synthesizing and analyzing rooted phylogenetic trees within a graph, called a tree alignment graph (TAG). The TAG can be queried and analyzed to explore uncertainty and conflict. It can also be synthesized to construct trees, presenting an alternative to supertrees approaches. We demonstrate these methods with two empirical datasets. In order to explore uncertainty, we constructed a TAG of the bootstrap trees from the Angiosperm Tree of Life project. Analysis of the resulting graph demonstrates that areas of the dataset that are unresolved in majority-rule consensus tree analyses can be understood in more detail within the context of a graph structure, using measures incorporating node degree and adjacency support. As an exercise in synthesis (i.e., summarization of a TAG constructed from the alignment trees), we also construct a TAG consisting of the taxonomy and source trees from a recent comprehensive bird study. We synthesized this graph into a tree that can be reconstructed in a repeatable fashion and where the underlying source information can be updated. The methods presented here are tractable for large scale analyses and serve as a basis for an alternative to consensus tree and supertree methods. Furthermore, the exploration of these graphs can expose structures and patterns within the dataset that are otherwise difficult to observe.

Description: by Javier Garcia-Bernardo, Mary J. Dunlop Cells live in uncertain, dynamic environments and have many mechanisms for sensing and responding to changes in their surroundings. However, sudden fluctuations in the environment can be catastrophic to a population if it relies solely on sensory responses, which have a delay associated with them. Cells can reconcile these effects by using a tunable stochastic response, where in the absence of a stressor they create phenotypic diversity within an isogenic population, but use a deterministic response when stressors are sensed. Here, we develop a stochastic model of the multiple antibiotic resistance network of Escherichia coli and show that it can produce tunable stochastic pulses in the activator MarA. In particular, we show that a combination of interlinked positive and negative feedback loops plays an important role in setting the dynamics of the stochastic pulses. Negative feedback produces a pulsatile response that is tunable, while positive feedback serves to amplify the effect. Our simulations show that the uninduced native network is in a parameter regime that is of low cost to the cell (taxing resistance mechanisms are expressed infrequently) and also elevated noise strength (phenotypic variability is high). The stochastic pulsing can be tuned by MarA induction such that variability is decreased once stresses are sensed, avoiding the detrimental effects of noise when an optimal MarA concentration is needed. We further show that variability in the expression of MarA can act as a bet hedging mechanism, allowing for survival in time-varying stress environments, however this effect is tunable to allow for a fully induced, deterministic response in the presence of a stressor.

Description: by Kate A. Stafford, Paul Robustelli, Arthur G. Palmer The relationship between inherent internal conformational processes and enzymatic activity or thermodynamic stability of proteins has proven difficult to characterize. The study of homologous proteins with differing thermostabilities offers an especially useful approach for understanding the functional aspects of conformational dynamics. In particular, ribonuclease HI (RNase H), an 18 kD globular protein that hydrolyzes the RNA strand of RNA:DNA hybrid substrates, has been extensively studied by NMR spectroscopy to characterize the differences in dynamics between homologs from the mesophilic organism E. coli and the thermophilic organism T. thermophilus . Herein, molecular dynamics simulations are reported for five homologous RNase H proteins of varying thermostabilities and enzymatic activities from organisms of markedly different preferred growth temperatures. For the E. coli and T. thermophilus proteins, strong agreement is obtained between simulated and experimental values for NMR order parameters and for dynamically averaged chemical shifts, suggesting that these simulations can be a productive platform for predicting the effects of individual amino acid residues on dynamic behavior. Analyses of the simulations reveal that a single residue differentiates between two different and otherwise conserved dynamic processes in a region of the protein known to form part of the substrate-binding interface. Additional key residues within these two categories are identified through the temperature-dependence of these conformational processes.

Description: by Rene Rex, Nelli Bill, Kerstin Schmidt-Hohagen, Dietmar Schomburg The Roseobacter clade is a ubiquitous group of marine α - proteobacteria . To gain insight into the versatile metabolism of this clade, we took a constraint-based approach and created a genome-scale metabolic model ( i Dsh827) of Dinoroseobacter shibae DFL12T. Our model is the first accounting for the energy demand of motility, the light-driven ATP generation and experimentally determined specific biomass composition. To cover a large variety of environmental conditions, as well as plasmid and single gene knock-out mutants, we simulated 391,560 different physiological states using flux balance analysis. We analyzed our results with regard to energy metabolism, validated them experimentally, and revealed a pronounced metabolic response to the availability of light. Furthermore, we introduced the energy demand of motility as an important parameter in genome-scale metabolic models. The results of our simulations also gave insight into the changing usage of the two degradation routes for dimethylsulfoniopropionate, an abundant compound in the ocean. A side product of dimethylsulfoniopropionate degradation is dimethyl sulfide, which seeds cloud formation and thus enhances the reflection of sunlight. By our exhaustive simulations, we were able to identify single-gene knock-out mutants, which show an increased production of dimethyl sulfide. In addition to the single-gene knock-out simulations we studied the effect of plasmid loss on the metabolism. Moreover, we explored the possible use of a functioning phosphofructokinase for D. shibae .

Description: by Ruth Nussinov

Description: by Farzad Farkhooi, Anja Froese, Eilif Muller, Randolf Menzel, Martin P. Nawrot Most neurons in peripheral sensory pathways initially respond vigorously when a preferred stimulus is presented, but adapt as stimulation continues. It is unclear how this phenomenon affects stimulus coding in the later stages of sensory processing. Here, we show that a temporally sparse and reliable stimulus representation develops naturally in sequential stages of a sensory network with adapting neurons. As a modeling framework we employ a mean-field approach together with an adaptive population density treatment, accompanied by numerical simulations of spiking neural networks. We find that cellular adaptation plays a critical role in the dynamic reduction of the trial-by-trial variability of cortical spike responses by transiently suppressing self-generated fast fluctuations in the cortical balanced network. This provides an explanation for a widespread cortical phenomenon by a simple mechanism. We further show that in the insect olfactory system cellular adaptation is sufficient to explain the emergence of the temporally sparse and reliable stimulus representation in the mushroom body. Our results reveal a generic, biophysically plausible mechanism that can explain the emergence of a temporally sparse and reliable stimulus representation within a sequential processing architecture.

Description: by Francisco Martínez-Jiménez, George Papadatos, Lun Yang, Iain M. Wallace, Vinod Kumar, Ursula Pieper, Andrej Sali, James R. Brown, John P. Overington, Marc A. Marti-Renom Mycobacterium tuberculosis , the causative agent of tuberculosis (TB), infects an estimated two billion people worldwide and is the leading cause of mortality due to infectious disease. The development of new anti-TB therapeutics is required, because of the emergence of multi-drug resistance strains as well as co-infection with other pathogens, especially HIV. Recently, the pharmaceutical company GlaxoSmithKline published the results of a high-throughput screen (HTS) of their two million compound library for anti-mycobacterial phenotypes. The screen revealed 776 compounds with significant activity against the M. tuberculosis H37Rv strain, including a subset of 177 prioritized compounds with high potency and low in vitro cytotoxicity. The next major challenge is the identification of the target proteins. Here, we use a computational approach that integrates historical bioassay data, chemical properties and structural comparisons of selected compounds to propose their potential targets in M. tuberculosis . We predicted 139 target - compound links, providing a necessary basis for further studies to characterize the mode of action of these compounds. The results from our analysis, including the predicted structural models, are available to the wider scientific community in the open source mode, to encourage further development of novel TB therapeutics.

Description: by A. J. K. Phillips, B. D. Fulcher, P. A. Robinson, E. B. Klerman Circadian rhythms are fundamental to life. In mammals, these rhythms are generated by pacemaker neurons in the suprachiasmatic nucleus (SCN) of the hypothalamus. The SCN is remarkably consistent in structure and function between species, yet mammalian rest/activity patterns are extremely diverse, including diurnal, nocturnal, and crepuscular behaviors. Two mechanisms have been proposed to account for this diversity: (i) modulation of SCN output by downstream nuclei, and (ii) direct effects of light on activity. These two mechanisms are difficult to disentangle experimentally and their respective roles remain unknown. To address this, we developed a computational model to simulate the two mechanisms and their influence on temporal niche. In our model, SCN output is relayed via the subparaventricular zone (SPZ) to the dorsomedial hypothalamus (DMH), and thence to ventrolateral preoptic nuclei (VLPO) and lateral hypothalamus (LHA). Using this model, we generated rich phenotypes that closely resemble experimental data. Modulation of SCN output at the SPZ was found to generate a full spectrum of diurnal-to-nocturnal phenotypes. Intriguingly, we also uncovered a novel mechanism for crepuscular behavior: if DMH/VLPO and DMH/LHA projections act cooperatively, daily activity is unimodal, but if they act competitively, activity can become bimodal. In addition, we successfully reproduced diurnal/nocturnal switching in the rodent Octodon degu using coordinated inversions in both masking and circadian modulation. Finally, the model correctly predicted the SCN lesion phenotype in squirrel monkeys: loss of circadian rhythmicity and emergence of ∼4-h sleep/wake cycles. In capturing these diverse phenotypes, the model provides a powerful new framework for understanding rest/activity patterns and relating them to underlying physiology. Given the ubiquitous effects of temporal organization on all aspects of animal behavior and physiology, this study sheds light on the physiological changes required to orchestrate adaptation to various temporal niches.

Description: by Ioannis N. Melas, Regina Samaga, Leonidas G. Alexopoulos, Steffen Klamt Cross-referencing experimental data with our current knowledge of signaling network topologies is one central goal of mathematical modeling of cellular signal transduction networks. We present a new methodology for data-driven interrogation and training of signaling networks. While most published methods for signaling network inference operate on Bayesian, Boolean, or ODE models, our approach uses integer linear programming (ILP) on interaction graphs to encode constraints on the qualitative behavior of the nodes. These constraints are posed by the network topology and their formulation as ILP allows us to predict the possible qualitative changes (up, down, no effect) of the activation levels of the nodes for a given stimulus. We provide four basic operations to detect and remove inconsistencies between measurements and predicted behavior: (i) find a topology-consistent explanation for responses of signaling nodes measured in a stimulus-response experiment (if none exists, find the closest explanation); (ii) determine a minimal set of nodes that need to be corrected to make an inconsistent scenario consistent; (iii) determine the optimal subgraph of the given network topology which can best reflect measurements from a set of experimental scenarios; (iv) find possibly missing edges that would improve the consistency of the graph with respect to a set of experimental scenarios the most. We demonstrate the applicability of the proposed approach by interrogating a manually curated interaction graph model of EGFR/ErbB signaling against a library of high-throughput phosphoproteomic data measured in primary hepatocytes. Our methods detect interactions that are likely to be inactive in hepatocytes and provide suggestions for new interactions that, if included, would significantly improve the goodness of fit. Our framework is highly flexible and the underlying model requires only easily accessible biological knowledge. All related algorithms were implemented in a freely available toolbox SigNetTrainer making it an appealing approach for various applications.

Description: by Tristan Ursell, Rosanna Man Wah Chau, Susanne Wisen, Devaki Bhaya, Kerwyn Casey Huang The emergent behaviors of communities of genotypically identical cells cannot be easily predicted from the behaviors of individual cells. In many cases, it is thought that direct cell-cell communication plays a critical role in the transition from individual to community behaviors. In the unicellular photosynthetic cyanobacterium Synechocystis sp. PCC 6803, individual cells exhibit light-directed motility (“phototaxis”) over surfaces, resulting in the emergence of dynamic spatial organization of multicellular communities. To probe this striking community behavior, we carried out time-lapse video microscopy coupled with quantitative analysis of single-cell dynamics under varying light conditions. These analyses suggest that cells secrete an extracellular substance that modifies the physical properties of the substrate, leading to enhanced motility and the ability for groups of cells to passively guide one another. We developed a biophysical model that demonstrates that this form of indirect, surface-based communication is sufficient to create distinct motile groups whose shape, velocity, and dynamics qualitatively match our experimental observations, even in the absence of direct cellular interactions or changes in single-cell behavior. Our computational analysis of the predicted community behavior, across a matrix of cellular concentrations and light biases, demonstrates that spatial patterning follows robust scaling laws and provides a useful resource for the generation of testable hypotheses regarding phototactic behavior. In addition, we predict that degradation of the surface modification may account for the secondary patterns occasionally observed after the initial formation of a community structure. Taken together, our modeling and experiments provide a framework to show that the emergent spatial organization of phototactic communities requires modification of the substrate, and this form of surface-based communication could provide insight into the behavior of a wide array of biological communities.

Description: by Jorge F. Mejias, Gary Marsat, Kieran Bol, Leonard Maler, André Longtin Cancellation of redundant information is a highly desirable feature of sensory systems, since it would potentially lead to a more efficient detection of novel information. However, biologically plausible mechanisms responsible for such selective cancellation, and especially those robust to realistic variations in the intensity of the redundant signals, are mostly unknown. In this work, we study, via in vivo experimental recordings and computational models, the behavior of a cerebellar-like circuit in the weakly electric fish which is known to perform cancellation of redundant stimuli. We experimentally observe contrast invariance in the cancellation of spatially and temporally redundant stimuli in such a system. Our model, which incorporates heterogeneously-delayed feedback, bursting dynamics and burst-induced STDP, is in agreement with our in vivo observations. In addition, the model gives insight on the activity of granule cells and parallel fibers involved in the feedback pathway, and provides a strong prediction on the parallel fiber potentiation time scale. Finally, our model predicts the existence of an optimal learning contrast around 15% contrast levels, which are commonly experienced by interacting fish.

Description: by Jing Tang, Leena Karhinen, Tao Xu, Agnieszka Szwajda, Bhagwan Yadav, Krister Wennerberg, Tero Aittokallio A recent trend in drug development is to identify drug combinations or multi-target agents that effectively modify multiple nodes of disease-associated networks. Such polypharmacological effects may reduce the risk of emerging drug resistance by means of attacking the disease networks through synergistic and synthetic lethal interactions. However, due to the exponentially increasing number of potential drug and target combinations, systematic approaches are needed for prioritizing the most potent multi-target alternatives on a global network level. We took a functional systems pharmacology approach toward the identification of selective target combinations for specific cancer cells by combining large-scale screening data on drug treatment efficacies and drug-target binding affinities. Our model-based prediction approach, named TIMMA, takes advantage of the polypharmacological effects of drugs and infers combinatorial drug efficacies through system-level target inhibition networks. Case studies in MCF-7 and MDA-MB-231 breast cancer and BxPC-3 pancreatic cancer cells demonstrated how the target inhibition modeling allows systematic exploration of functional interactions between drugs and their targets to maximally inhibit multiple survival pathways in a given cancer type. The TIMMA prediction results were experimentally validated by means of systematic siRNA-mediated silencing of the selected targets and their pairwise combinations, showing increased ability to identify not only such druggable kinase targets that are essential for cancer survival either individually or in combination, but also synergistic interactions indicative of non-additive drug efficacies. These system-level analyses were enabled by a novel model construction method utilizing maximization and minimization rules, as well as a model selection algorithm based on sequential forward floating search. Compared with an existing computational solution, TIMMA showed both enhanced prediction accuracies in cross validation as well as significant reduction in computation times. Such cost-effective computational-experimental design strategies have the potential to greatly speed-up the drug testing efforts by prioritizing those interventions and interactions warranting further study in individual cancer cases.

Description: by Kartik Subramanian, Mark R. Paul, John J. Tyson The free-living aquatic bacterium, Caulobacter crescentus , exhibits two different morphologies during its life cycle. The morphological change from swarmer cell to stalked cell is a result of changes of function of two bi-functional histidine kinases, PleC and CckA. Here, we describe a detailed molecular mechanism by which the function of PleC changes between phosphatase and kinase state. By mathematical modeling of our proposed molecular interactions, we derive conditions under which PleC, CckA and its response regulators exhibit bistable behavior, thus providing a scenario for robust switching between swarmer and stalked states. Our simulations are in reasonable agreement with in vitro and in vivo experimental observations of wild type and mutant phenotypes. According to our model, the kinase form of PleC is essential for the swarmer-to-stalked transition and to prevent premature development of the swarmer pole. Based on our results, we reconcile some published experimental observations and suggest novel mutants to test our predictions.

Description: by Yun M. Zhao, Anthony R. French Natural killer (NK) cells are innate lymphocytes that provide early host defense against intracellular pathogens, such as viruses. Although NK cell development, homeostasis, and proliferation are regulated by IL-15, the influence of IL-15 receptor (IL-15R)-mediated signaling at the cellular level has not been quantitatively characterized. We developed a mathematical model to analyze the kinetic interactions that control the formation and localization of IL-15/IL-15R complexes. Our computational results demonstrated that IL-15/IL-15R complexes on the cell surface were a key determinant of the magnitude of the IL-15 proliferative signal and that IL-15R occupancy functioned as an effective surrogate measure of receptor signaling. Ligand binding and receptor internalization modulated IL-15R occupancy. Our work supports the hypothesis that the total number and duration of IL-15/IL-15R complexes on the cell surface crosses a quantitative threshold prior to the initiation of NK cell division. Furthermore, our model predicted that the upregulation of IL-15Rα on NK cells substantially increased IL-15R complex formation and accelerated the expansion of dividing NK cells with the greatest impact at low IL-15 concentrations. Model predictions of the threshold requirement for NK cell recruitment to the cell cycle and the subsequent exponential proliferation correlated well with experimental data. In summary, our modeling analysis provides quantitative insight into the regulation of NK cell proliferation at the receptor level and provides a framework for the development of IL-15 based immunotherapies to modulate NK cell proliferation.

Description: by William Pontius, Michael W. Sneddon, Thierry Emonet In many sensory systems, transmembrane receptors are spatially organized in large clusters. Such arrangement may facilitate signal amplification and the integration of multiple stimuli. However, this organization likely also affects the kinetics of signaling since the cytoplasmic enzymes that modulate the activity of the receptors must localize to the cluster prior to receptor modification. Here we examine how these spatial considerations shape signaling dynamics at rest and in response to stimuli. As a model system, we use the chemotaxis pathway of Escherichia coli , a canonical system for the study of how organisms sense, respond, and adapt to environmental stimuli. In bacterial chemotaxis, adaptation is mediated by two enzymes that localize to the clustered receptors and modulate their activity through methylation-demethylation. Using a novel stochastic simulation, we show that distributive receptor methylation is necessary for successful adaptation to stimulus and also leads to large fluctuations in receptor activity in the steady state. These fluctuations arise from noise in the number of localized enzymes combined with saturated modification kinetics between the localized enzymes and the receptor substrate. An analytical model explains how saturated enzyme kinetics and large fluctuations can coexist with an adapted state robust to variation in the expression levels of the pathway constituents, a key requirement to ensure the functionality of individual cells within a population. This contrasts with the well-mixed covalent modification system studied by Goldbeter and Koshland in which mean activity becomes ultrasensitive to protein abundances when the enzymes operate at saturation. Large fluctuations in receptor activity have been quantified experimentally and may benefit the cell by enhancing its ability to explore empty environments and track shallow nutrient gradients. Here we clarify the mechanistic relationship of these large fluctuations to well-studied aspects of the chemotaxis system, precise adaptation and functional robustness.

Description: by Vishal N. Patel, Giridharan Gokulrangan, Salim A. Chowdhury, Yanwen Chen, Andrew E. Sloan, Mehmet Koyutürk, Jill Barnholtz-Sloan, Mark R. Chance To determine a molecular basis for prognostic differences in glioblastoma multiforme (GBM), we employed a combinatorial network analysis framework to exhaustively search for molecular patterns in protein-protein interaction (PPI) networks. We identified a dysregulated molecular signature distinguishing short-term (survival 635 days) survivors of GBM using whole genome expression data from The Cancer Genome Atlas (TCGA). A 50-gene subnetwork signature achieved 80% prediction accuracy when tested against an independent gene expression dataset. Functional annotations for the subnetwork signature included “protein kinase cascade,” “IκB kinase/NFκB cascade,” and “regulation of programmed cell death” – all of which were not significant in signatures of existing subtypes. Finally, we used label-free proteomics to examine how our subnetwork signature predicted protein level expression differences in an independent GBM cohort of 16 patients. We found that the genes discovered using network biology had a higher probability of dysregulated protein expression than either genes exhibiting individual differential expression or genes derived from known GBM subtypes. In particular, the long-term survivor subtype was characterized by increased protein expression of DNM1 and MAPK1 and decreased expression of HSPA9, PSMD3, and CANX. Overall, we demonstrate that the combinatorial analysis of gene expression data constrained by PPIs outlines an approach for the discovery of robust and translatable molecular signatures in GBM.

Description: by Martin D. Vesper, Bert L. de Groot Hemoglobin is the prototypic allosteric protein. Still, its molecular allosteric mechanism is not fully understood. To elucidate the mechanism of cooperativity on an atomistic level, we developed a novel computational technique to analyse the coupling of tertiary and quaternary motions. From Molecular Dynamics simulations showing spontaneous quaternary transitions, we separated the transition trajectories into two orthogonal sets of motions: one consisting of intra-chain motions only (referred to as tertiary-only ) and one consisting of global inter-chain motions only (referred to as quaternary-only ). The two underlying subspaces are orthogonal by construction and their direct sum is the space of full motions. Using Functional Mode Analysis, we were able to identify a collective coordinate within the tertiary-only subspace that is correlated to the most dominant motion within the quaternary-only motions, hence providing direct insight into the allosteric coupling mechanism between tertiary and quaternary conformation changes. This coupling-motion is substantially different from tertiary structure changes between the crystallographic structures of the T- and R-state. We found that hemoglobin's allosteric mechanism of communication between subunits is equally based on hydrogen bonds and steric interactions. In addition, we were able to affect the T-to-R transition rates by choosing different histidine protonation states, thereby providing a possible atomistic explanation for the Bohr effect.

Description: by James J. Lee, Justin Huang, Christopher G. England, Lacey R. McNally, Hermann B. Frieboes A clear contradiction exists between cytotoxic in-vitro studies demonstrating effectiveness of Gemcitabine to curtail pancreatic cancer and in-vivo studies failing to show Gemcitabine as an effective treatment. The outcome of chemotherapy in metastatic stages, where surgery is no longer viable, shows a 5-year survival

Description: by Sven-Eric Schelhorn, Matthias Fischer, Laura Tolosi, Janine Altmüller, Peter Nürnberg, Herbert Pfister, Thomas Lengauer, Frank Berthold In excess of % of human cancer incidents have a viral cofactor. Epidemiological studies of idiopathic human cancers indicate that additional tumor viruses remain to be discovered. Recent advances in sequencing technology have enabled systematic screenings of human tumor transcriptomes for viral transcripts. However, technical problems such as low abundances of viral transcripts in large volumes of sequencing data, viral sequence divergence, and homology between viral and human factors significantly confound identification of tumor viruses. We have developed a novel computational approach for detecting viral transcripts in human cancers that takes the aforementioned confounding factors into account and is applicable to a wide variety of viruses and tumors. We apply the approach to conducting the first systematic search for viruses in neuroblastoma, the most common cancer in infancy. The diverse clinical progression of this disease as well as related epidemiological and virological findings are highly suggestive of a pathogenic cofactor. However, a viral etiology of neuroblastoma is currently contested. We mapped transcriptomes of neuroblastoma as well as positive and negative controls to the human and all known viral genomes in order to detect both known and unknown viruses. Analysis of controls, comparisons with related methods, and statistical estimates demonstrate the high sensitivity of our approach. Detailed investigation of putative viral transcripts within neuroblastoma samples did not provide evidence for the existence of any known human viruses. Likewise, de-novo assembly and analysis of chimeric transcripts did not result in expression signatures associated with novel human pathogens. While confounding factors such as sample dilution or viral clearance in progressed tumors may mask viral cofactors in the data, in principle, this is rendered less likely by the high sensitivity of our approach and the number of biological replicates analyzed. Therefore, our results suggest that frequent viral cofactors of metastatic neuroblastoma are unlikely.

Description: by Fan Jin, Chen Yu, Luhua Lai, Zhirong Liu Intrinsically disordered proteins (IDPs) were found to be widely associated with human diseases and may serve as potential drug design targets. However, drug design targeting IDPs is still in the very early stages. Progress in drug design is usually achieved using experimental screening; however, the structural disorder of IDPs makes it difficult to characterize their interaction with ligands using experiments alone. To better understand the structure of IDPs and their interactions with small molecule ligands, we performed extensive simulations on the c-Myc 370–409 peptide and its binding to a reported small molecule inhibitor, ligand 10074-A4. We found that the conformational space of the apo c-Myc 370–409 peptide was rather dispersed and that the conformations of the peptide were stabilized mainly by charge interactions and hydrogen bonds. Under the binding of the ligand, c-Myc 370–409 remained disordered. The ligand was found to bind to c-Myc 370–409 at different sites along the chain and behaved like a ‘ligand cloud’. In contrast to ligand binding to more rigid target proteins that usually results in a dominant bound structure, ligand binding to IDPs may better be described as ligand clouds around protein clouds. Nevertheless, the binding of the ligand and a non-ligand to the c-Myc 370–409 target could be clearly distinguished. The present study provides insights that will help improve rational drug design that targets IDPs.

Description: by Shuguang Yuan, Rongliang Wu, Dorota Latek, Bartosz Trzaskowski, Slawomir Filipek Sphingosine 1-phosphate (S1P) is a lysophospholipid mediator which activates G protein–coupled sphingosine 1-phosphate receptors and thus evokes a variety of cell and tissue responses including lymphocyte trafficking, endothelial development, integrity, and maturation. We performed five all-atom 700 ns molecular dynamics simulations of the sphingosine 1-phosphate receptor 1 (S1P 1 ) based on recently released crystal structure of that receptor with an antagonist. We found that the initial movements of amino acid residues occurred in the area of highly conserved W269 6.48 in TM6 which is close to the ligand binding location. Those residues located in the central part of the receptor and adjacent to kinks of TM helices comprise of a transmission switch. Side chains movements of those residues were coupled to the movements of water molecules inside the receptor which helped in the gradual opening of intracellular part of the receptor. The most stable parts of the protein were helices TM1 and TM2, while the largest movement was observed for TM7, possibly due to the short intracellular part starting with a helix kink at P 7.50 , which might be the first helix to move at the intracellular side. We show for the first time the detailed view of the concerted action of the transmission switch and Trp (W 6.48 ) rotamer toggle switch leading to redirection of water molecules flow in the central part of the receptor. That event is a prerequisite for subsequent changes in intracellular part of the receptor involving water influx and opening of the receptor structure.

Description: by Thomas Hannagan, Jonathan Grainger In a recent study, Rauschecker et al. convincingly demonstrate that visual words evoke neural activation signals in the Visual Word Form Area that can be classified based on where they were presented in the visual fields. This result goes against the prevailing consensus, and begs an explanation. We show that one of the simplest possible models for word recognition, a multilayer feedforward network, will exhibit precisely the same behavior when trained to recognize words at different locations. The model suggests that the VWFA initially starts with information about location, which is not being suppressed during reading acquisition more than is needed to meet the requirements of location-invariant word recognition. Some new interpretations of Rauschecker et al.'s results are proposed, and three specific predictions are derived to be tested in further studies.

Description: by André Yoshiaki Kashiwabara, Ígor Bonadio, Vitor Onuchic, Felipe Amado, Rafael Mathias, Alan Mitchell Durham Discrete Markovian models can be used to characterize patterns in sequences of values and have many applications in biological sequence analysis, including gene prediction, CpG island detection, alignment, and protein profiling. We present ToPS, a computational framework that can be used to implement different applications in bioinformatics analysis by combining eight kinds of models: (i) independent and identically distributed process; (ii) variable-length Markov chain; (iii) inhomogeneous Markov chain; (iv) hidden Markov model; (v) profile hidden Markov model; (vi) pair hidden Markov model; (vii) generalized hidden Markov model; and (viii) similarity based sequence weighting. The framework includes functionality for training, simulation and decoding of the models. Additionally, it provides two methods to help parameter setting: Akaike and Bayesian information criteria (AIC and BIC). The models can be used stand-alone, combined in Bayesian classifiers, or included in more complex, multi-model, probabilistic architectures using GHMMs. In particular the framework provides a novel, flexible, implementation of decoding in GHMMs that detects when the architecture can be traversed efficiently.

Description: by Biswa Sengupta, Simon B. Laughlin, Jeremy E. Niven A balance between excitatory and inhibitory synaptic currents is thought to be important for several aspects of information processing in cortical neurons in vivo , including gain control, bandwidth and receptive field structure. These factors will affect the firing rate of cortical neurons and their reliability, with consequences for their information coding and energy consumption. Yet how balanced synaptic currents contribute to the coding efficiency and energy efficiency of cortical neurons remains unclear. We used single compartment computational models with stochastic voltage-gated ion channels to determine whether synaptic regimes that produce balanced excitatory and inhibitory currents have specific advantages over other input regimes. Specifically, we compared models with only excitatory synaptic inputs to those with equal excitatory and inhibitory conductances, and stronger inhibitory than excitatory conductances (i.e. approximately balanced synaptic currents). Using these models, we show that balanced synaptic currents evoke fewer spikes per second than excitatory inputs alone or equal excitatory and inhibitory conductances. However, spikes evoked by balanced synaptic inputs are more informative (bits/spike), so that spike trains evoked by all three regimes have similar information rates (bits/s). Consequently, because spikes dominate the energy consumption of our computational models, approximately balanced synaptic currents are also more energy efficient than other synaptic regimes. Thus, by producing fewer, more informative spikes approximately balanced synaptic currents in cortical neurons can promote both coding efficiency and energy efficiency.

Description: by Joshua L. Payne, Andreas Wagner Gene regulatory circuits drive the development, physiology, and behavior of organisms from bacteria to humans. The phenotypes or functions of such circuits are embodied in the gene expression patterns they form. Regulatory circuits are typically multifunctional, forming distinct gene expression patterns in different embryonic stages, tissues, or physiological states. Any one circuit with a single function can be realized by many different regulatory genotypes. Multifunctionality presumably constrains this number, but we do not know to what extent. We here exhaustively characterize a genotype space harboring millions of model regulatory circuits and all their possible functions. As a circuit's number of functions increases, the number of genotypes with a given number of functions decreases exponentially but can remain very large for a modest number of functions. However, the sets of circuits that can form any one set of functions becomes increasingly fragmented. As a result, historical contingency becomes widespread in circuits with many functions. Whether a circuit can acquire an additional function in the course of its evolution becomes increasingly dependent on the function it already has. Circuits with many functions also become increasingly brittle and sensitive to mutation. These observations are generic properties of a broad class of circuits and independent of any one circuit genotype or phenotype.

Description: by Nanye Long, Samuel P. Dickson, Jessica M. Maia, Hee Shin Kim, Qianqian Zhu, Andrew S. Allen Although many methods are available to test sequence variants for association with complex diseases and traits, methods that specifically seek to identify causal variants are less developed. Here we develop and evaluate a Bayesian hierarchical regression method that incorporates prior information on the likelihood of variant causality through weighting of variant effects. By simulation studies using both simulated and real sequence variants, we compared a standard single variant test for analyzing variant-disease association with the proposed method using different weighting schemes. We found that by leveraging linkage disequilibrium of variants with known GWAS signals and sequence conservation (phastCons), the proposed method provides a powerful approach for detecting causal variants while controlling false positives.

Description: by Drew H. Bryant, Mark Moll, Paul W. Finn, Lydia E. Kavraki The protein kinases are a large family of enzymes that play fundamental roles in propagating signals within the cell. Because of the high degree of binding site similarity shared among protein kinases, designing drug compounds with high specificity among the kinases has proven difficult. However, computational approaches to comparing the 3-dimensional geometry and physicochemical properties of key binding site residue positions have been shown to be informative of inhibitor selectivity. The Combinatorial Clustering Of Residue Position Subsets (ccorps) method, introduced here, provides a semi-supervised learning approach for identifying structural features that are correlated with a given set of annotation labels. Here, ccorps is applied to the problem of identifying structural features of the kinase atp binding site that are informative of inhibitor binding. ccorps is demonstrated to make perfect or near-perfect predictions for the binding affinity profile of 8 of the 38 kinase inhibitors studied, while only having overall poor predictive ability for 1 of the 38 compounds. Additionally, ccorps is shown to identify shared structural features across phylogenetically diverse groups of kinases that are correlated with binding affinity for particular inhibitors; such instances of structural similarity among phylogenetically diverse kinases are also shown to not be rare among kinases. Finally, these function-specific structural features may serve as potential starting points for the development of highly specific kinase inhibitors.

Description: by Christian P. Koch, Anna M. Perna, Max Pillong, Nickolay K. Todoroff, Paul Wrede, Gerd Folkers, Jan A. Hiss, Gisbert Schneider Designed peptides that bind to major histocompatibility protein I (MHC-I) allomorphs bear the promise of representing epitopes that stimulate a desired immune response. A rigorous bioinformatical exploration of sequence patterns hidden in peptides that bind to the mouse MHC-I allomorph H-2K b is presented. We exemplify and validate these motif findings by systematically dissecting the epitope SIINFEKL and analyzing the resulting fragments for their binding potential to H-2K b in a thermal denaturation assay. The results demonstrate that only fragments exclusively retaining the carboxy- or amino-terminus of the reference peptide exhibit significant binding potential, with the N-terminal pentapeptide SIINF as shortest ligand. This study demonstrates that sophisticated machine-learning algorithms excel at extracting fine-grained patterns from peptide sequence data and predicting MHC-I binding peptides, thereby considerably extending existing linear prediction models and providing a fresh view on the computer-based molecular design of future synthetic vaccines. The server for prediction is available at http://modlab-cadd.ethz.ch (SLiDER tool, MHC-I version 2012).

Description: by Nicholas Ketz, Srinimisha G. Morkonda, Randall C. O'Reilly The learning mechanism in the hippocampus has almost universally been assumed to be Hebbian in nature, where individual neurons in an engram join together with synaptic weight increases to support facilitated recall of memories later. However, it is also widely known that Hebbian learning mechanisms impose significant capacity constraints, and are generally less computationally powerful than learning mechanisms that take advantage of error signals. We show that the differential phase relationships of hippocampal subfields within the overall theta rhythm enable a powerful form of error-driven learning, which results in significantly greater capacity, as shown in computer simulations. In one phase of the theta cycle, the bidirectional connectivity between CA1 and entorhinal cortex can be trained in an error-driven fashion to learn to effectively encode the cortical inputs in a compact and sparse form over CA1. In a subsequent portion of the theta cycle, the system attempts to recall an existing memory, via the pathway from entorhinal cortex to CA3 and CA1. Finally the full theta cycle completes when a strong target encoding representation of the current input is imposed onto the CA1 via direct projections from entorhinal cortex. The difference between this target encoding and the attempted recall of the same representation on CA1 constitutes an error signal that can drive the learning of CA3 to CA1 synapses. This CA3 to CA1 pathway is critical for enabling full reinstatement of recalled hippocampal memories out in cortex. Taken together, these new learning dynamics enable a much more robust, high-capacity model of hippocampal learning than was available previously under the classical Hebbian model.

Description: by Jörg Bornschein, Marc Henniges, Jörg Lücke Simple cells in primary visual cortex were famously found to respond to low-level image components such as edges. Sparse coding and independent component analysis (ICA) emerged as the standard computational models for simple cell coding because they linked their receptive fields to the statistics of visual stimuli. However, a salient feature of image statistics, occlusions of image components, is not considered by these models. Here we ask if occlusions have an effect on the predicted shapes of simple cell receptive fields. We use a comparative approach to answer this question and investigate two models for simple cells: a standard linear model and an occlusive model. For both models we simultaneously estimate optimal receptive fields, sparsity and stimulus noise. The two models are identical except for their component superposition assumption. We find the image encoding and receptive fields predicted by the models to differ significantly. While both models predict many Gabor-like fields, the occlusive model predicts a much sparser encoding and high percentages of ‘globular’ receptive fields. This relatively new center-surround type of simple cell response is observed since reverse correlation is used in experimental studies. While high percentages of ‘globular’ fields can be obtained using specific choices of sparsity and overcompleteness in linear sparse coding, no or only low proportions are reported in the vast majority of studies on linear models (including all ICA models). Likewise, for the here investigated linear model and optimal sparsity, only low proportions of ‘globular’ fields are observed. In comparison, the occlusive model robustly infers high proportions and can match the experimentally observed high proportions of ‘globular’ fields well. Our computational study, therefore, suggests that ‘globular’ fields may be evidence for an optimal encoding of visual occlusions in primary visual cortex.

Description: by Elzbieta Petelenz-Kurdziel, Clemens Kuehn, Bodil Nordlander, Dagmara Klein, Kuk-Ki Hong, Therese Jacobson, Peter Dahl, Jörg Schaber, Jens Nielsen, Stefan Hohmann, Edda Klipp We provide an integrated dynamic view on a eukaryotic osmolyte system, linking signaling with regulation of gene expression, metabolic control and growth. Adaptation to osmotic changes enables cells to adjust cellular activity and turgor pressure to an altered environment. The yeast Saccharomyces cerevisiae adapts to hyperosmotic stress by activating the HOG signaling cascade, which controls glycerol accumulation. The Hog1 kinase stimulates transcription of genes encoding enzymes required for glycerol production (Gpd1, Gpp2) and glycerol import (Stl1) and activates a regulatory enzyme in glycolysis (Pfk26/27). In addition, glycerol outflow is prevented by closure of the Fps1 glycerol facilitator. In order to better understand the contributions to glycerol accumulation of these different mechanisms and how redox and energy metabolism as well as biomass production are maintained under such conditions we collected an extensive dataset. Over a period of 180 min after hyperosmotic shock we monitored in wild type and different mutant cells the concentrations of key metabolites and proteins relevant for osmoadaptation. The dataset was used to parameterize an ODE model that reproduces the generated data very well. A detailed computational analysis using time-dependent response coefficients showed that Pfk26/27 contributes to rerouting glycolytic flux towards lower glycolysis. The transient growth arrest following hyperosmotic shock further adds to redirecting almost all glycolytic flux from biomass towards glycerol production. Osmoadaptation is robust to loss of individual adaptation pathways because of the existence and upregulation of alternative routes of glycerol accumulation. For instance, the Stl1 glycerol importer contributes to glycerol accumulation in a mutant with diminished glycerol production capacity. In addition, our observations suggest a role for trehalose accumulation in osmoadaptation and that Hog1 probably directly contributes to the regulation of the Fps1 glycerol facilitator. Taken together, we elucidated how different metabolic adaptation mechanisms cooperate and provide hypotheses for further experimental studies.

Description: by Balázs Könnyű, S. Kashif Sadiq, Tamás Turányi, Rita Hírmondó, Barbara Müller, Hans-Georg Kräusslich, Peter V. Coveney, Viktor Müller Proteolytic processing of Gag and Gag-Pol polyproteins by the viral protease (PR) is crucial for the production of infectious HIV-1, and inhibitors of the viral PR are an integral part of current antiretroviral therapy. The process has several layers of complexity (multiple cleavage sites and substrates; multiple enzyme forms; PR auto-processing), which calls for a systems level approach to identify key vulnerabilities and optimal treatment strategies. Here we present the first full reaction kinetics model of proteolytic processing by HIV-1 PR, taking into account all canonical cleavage sites within Gag and Gag-Pol, intermediate products and enzyme forms, enzyme dimerization, the initial auto-cleavage of full-length Gag-Pol as well as self-cleavage of PR. The model allows us to identify the rate limiting step of virion maturation and the parameters with the strongest effect on maturation kinetics. Using the modelling framework, we predict interactions and compensatory potential between individual cleavage rates and drugs, characterize the time course of the process, explain the steep dose response curves associated with PR inhibitors and gain new insights into drug action. While the results of the model are subject to limitations arising from the simplifying assumptions used and from the uncertainties in the parameter estimates, the developed framework provides an extendable open-access platform to incorporate new data and hypotheses in the future.

Description: by Simon Garnier, Tucker Murphy, Matthew Lutz, Edward Hurme, Simon Leblanc, Iain D. Couzin Robustness and adaptability are central to the functioning of biological systems, from gene networks to animal societies. Yet the mechanisms by which living organisms achieve both stability to perturbations and sensitivity to input are poorly understood. Here, we present an integrated study of a living architecture in which army ants interconnect their bodies to span gaps. We demonstrate that these self-assembled bridges are a highly effective means of maintaining traffic flow over unpredictable terrain. The individual-level rules responsible depend only on locally-estimated traffic intensity and the number of neighbours to which ants are attached within the structure. We employ a parameterized computational model to reveal that bridges are tuned to be maximally stable in the face of regular, periodic fluctuations in traffic. However analysis of the model also suggests that interactions among ants give rise to feedback processes that result in bridges being highly responsive to sudden interruptions in traffic. Subsequent field experiments confirm this prediction and thus the dual nature of stability and flexibility in living bridges. Our study demonstrates the importance of robust and adaptive modular architecture to efficient traffic organisation and reveals general principles regarding the regulation of form in biological self-assemblies.

Description: by Brad Busse, Stephen Smith Synapses of the mammalian central nervous system are highly diverse in function and molecular composition. Synapse diversity per se may be critical to brain function, since memory and homeostatic mechanisms are thought to be rooted primarily in activity-dependent plastic changes in specific subsets of individual synapses. Unfortunately, the measurement of synapse diversity has been restricted by the limitations of methods capable of measuring synapse properties at the level of individual synapses. Array tomography is a new high-resolution, high-throughput proteomic imaging method that has the potential to advance the measurement of unit-level synapse diversity across large and diverse synapse populations. Here we present an automated feature extraction and classification algorithm designed to quantify synapses from high-dimensional array tomographic data too voluminous for manual analysis. We demonstrate the use of this method to quantify laminar distributions of synapses in mouse somatosensory cortex and validate the classification process by detecting the presence of known but uncommon proteomic profiles. Such classification and quantification will be highly useful in identifying specific subpopulations of synapses exhibiting plasticity in response to perturbations from the environment or the sensory periphery.

Description: by Himadri Mukhopadhyay, Shaun-Paul Cordoba, Philip K. Maini, P. Anton van der Merwe, Omer Dushek Receptor phosphorylation is thought to be tightly regulated because phosphorylated receptors initiate signaling cascades leading to cellular activation. The T cell antigen receptor (TCR) on the surface of T cells is phosphorylated by the kinase Lck and dephosphorylated by the phosphatase CD45 on multiple immunoreceptor tyrosine-based activation motifs (ITAMs). Intriguingly, Lck sequentially phosphorylates ITAMs and ZAP-70, a cytosolic kinase, binds to phosphorylated ITAMs with differential affinities. The purpose of multiple ITAMs, their sequential phosphorylation, and the differential ZAP-70 affinities are unknown. Here, we use a systems model to show that this signaling architecture produces emergent ultrasensitivity resulting in switch-like responses at the scale of individual TCRs. Importantly, this switch-like response is an emergent property, so that removal of multiple ITAMs, sequential phosphorylation, or differential affinities abolishes the switch. We propose that highly regulated TCR phosphorylation is achieved by an emergent switch-like response and use the systems model to design novel chimeric antigen receptors for therapy.

Description: by Clive G. Bowsher, Margaritis Voliotis, Peter S. Swain Cells live in changing, dynamic environments. To understand cellular decision-making, we must therefore understand how fluctuating inputs are processed by noisy biomolecular networks. Here we present a general methodology for analyzing the fidelity with which different statistics of a fluctuating input are represented, or encoded, in the output of a signaling system over time. We identify two orthogonal sources of error that corrupt perfect representation of the signal: dynamical error, which occurs when the network responds on average to other features of the input trajectory as well as to the signal of interest, and mechanistic error, which occurs because biochemical reactions comprising the signaling mechanism are stochastic. Trade-offs between these two errors can determine the system's fidelity. By developing mathematical approaches to derive dynamics conditional on input trajectories we can show, for example, that increased biochemical noise (mechanistic error) can improve fidelity and that both negative and positive feedback degrade fidelity, for standard models of genetic autoregulation. For a group of cells, the fidelity of the collective output exceeds that of an individual cell and negative feedback then typically becomes beneficial. We can also predict the dynamic signal for which a given system has highest fidelity and, conversely, how to modify the network design to maximize fidelity for a given dynamic signal. Our approach is general, has applications to both systems and synthetic biology, and will help underpin studies of cellular behavior in natural, dynamic environments.

Description: by Kirill S. Korolev Species expand their geographical ranges following an environmental change, long range dispersal, or a new adaptation. Range expansions not only bring an ecological change, but also affect the evolution of the expanding species. Although the dynamics of deleterious, neutral, and beneficial mutations have been extensively studied in expanding populations, the fate of alleles under frequency-dependent selection remains largely unexplored. The dynamics of cooperative alleles are particularly interesting because selection can be both frequency and density dependent, resulting in a coupling between population and evolutionary dynamics. This coupling leads to an increase in the frequency of cooperators at the expansion front, and, under certain conditions, the entire front can be taken over by cooperators. Thus, a mixed population wave can split into an expansion wave of only cooperators followed by an invasion wave of defectors. After the splitting, cooperators increase in abundance by expanding into new territories faster than they are invaded by defectors. Our results not only provide an explanation for the maintenance of cooperation but also elucidate the effect of eco-evolutionary feedback on the maintenance of genetic diversity during range expansions. When cooperators do not split away, we find that defectors can spread much faster with cooperators than they would be able to on their own or by invading cooperators. This enhanced rate of expansion in mixed waves could counterbalance the loss of genetic diversity due to the founder effect for mutations under frequency-dependent selection. Although we focus on cooperator-defector interactions, our analysis could also be relevant for other systems described by reaction-diffusion equations.

Description: by Christoph Schmal, Peter Reimann, Dorothee Staiger The circadian clock controls many physiological processes in higher plants and causes a large fraction of the genome to be expressed with a 24h rhythm. The transcripts encoding the RNA-binding proteins At GRP7 ( Arabidopsis thaliana Glycine Rich Protein 7 ) and At GRP8 oscillate with evening peaks. The circadian clock components CCA1 and LHY negatively affect AtGRP7 expression at the level of transcription. At GRP7 and At GRP8, in turn, negatively auto-regulate and reciprocally cross-regulate post-transcriptionally: high protein levels promote the generation of an alternative splice form that is rapidly degraded. This clock-regulated feedback loop has been proposed to act as a molecular slave oscillator in clock output. While mathematical models describing the circadian core oscillator in Arabidopsis thaliana were introduced recently, we propose here the first model of a circadian slave oscillator. We define the slave oscillator in terms of ordinary differential equations and identify the model's parameters by an optimization procedure based on experimental results. The model successfully reproduces the pertinent experimental findings such as waveforms, phases, and half-lives of the time-dependent concentrations. Furthermore, we obtain insights into possible mechanisms underlying the observed experimental dynamics: the negative auto-regulation and reciprocal cross-regulation via alternative splicing could be responsible for the sharply peaking waveforms of the AtGRP7 and AtGRP8 mRNA. Moreover, our results suggest that the AtGRP8 transcript oscillations are subordinated to those of AtGRP7 due to a higher impact of At GRP7 protein on alternative splicing of its own and of the AtGRP8 pre-mRNA compared to the impact of At GRP8 protein. Importantly, a bifurcation analysis provides theoretical evidence that the slave oscillator could be a toggle switch, arising from the reciprocal cross-regulation at the post-transcriptional level. In view of this, transcriptional repression of AtGRP7 and AtGRP8 by LHY and CCA1 induces oscillations of the toggle switch, leading to the observed high-amplitude oscillations of AtGRP7 mRNA.

Description: by Oscar Franzén, Jon Jerlström-Hultqvist, Elin Einarsson, Johan Ankarklev, Marcela Ferella, Björn Andersson, Staffan G. Svärd Giardia intestinalis is a common cause of diarrheal disease and it consists of eight genetically distinct genotypes or assemblages (A-H). Only assemblages A and B infect humans and are suggested to represent two different Giardia species. Correlations exist between assemblage type and host-specificity and to some extent symptoms. Phenotypical differences have been documented between assemblages and genome sequences are available for A, B and E. We have characterized and compared the polyadenylated transcriptomes of assemblages A, B and E. Four genetically different isolates were studied (WB (AI), AS175 (AII), P15 (E) and GS (B)) using paired-end, strand-specific RNA-seq. Most of the genome was transcribed in trophozoites grown in vitro , but at vastly different levels. RNA-seq confirmed many of the present annotations and refined the current genome annotation. Gene expression divergence was found to recapitulate the known phylogeny, and uncovered lineage-specific differences in expression. Polyadenylation sites were mapped for over 70% of the genes and revealed many examples of conserved and unexpectedly long 3′ UTRs. 28 open reading frames were found in a non-transcribed gene cluster on chromosome 5 of the WB isolate. Analysis of allele-specific expression revealed a correlation between allele-dosage and allele expression in the GS isolate. Previously reported cis -splicing events were confirmed and global mapping of cis -splicing identified only one novel intron. These observations can possibly explain differences in host-preference and symptoms, and it will be the basis for further studies of Giardia pathogenesis and biology.

Description: by Syed Abbas Bukhari, Gustavo Caetano-Anollés The spatial arrangements of secondary structures in proteins, irrespective of their connectivity, depict the overall shape and organization of protein domains. These features have been used in the CATH and SCOP classifications to hierarchically partition fold space and define the architectural make up of proteins. Here we use phylogenomic methods and a census of CATH structures in hundreds of genomes to study the origin and diversification of protein architectures (A) and their associated topologies (T) and superfamilies (H). Phylogenies that describe the evolution of domain structures and proteomes were reconstructed from the structural census and used to generate timelines of domain discovery. Phylogenies of CATH domains at T and H levels of structural abstraction and associated chronologies revealed patterns of reductive evolution, the early rise of Archaea, three epochs in the evolution of the protein world, and patterns of structural sharing between superkingdoms. Phylogenies of proteomes confirmed the early appearance of Archaea. While these findings are in agreement with previous phylogenomic studies based on the SCOP classification, phylogenies unveiled sharing patterns between Archaea and Eukarya that are recent and can explain the canonical bacterial rooting typically recovered from sequence analysis. Phylogenies of CATH domains at A level uncovered general patterns of architectural origin and diversification. The tree of A structures showed that ancient structural designs such as the 3-layer (αβα) sandwich (3.40) or the orthogonal bundle (1.10) are comparatively simpler in their makeup and are involved in basic cellular functions. In contrast, modern structural designs such as prisms , propellers, 2-solenoid, super-roll, clam, trefoil and box are not widely distributed and were probably adopted to perform specialized functions. Our timelines therefore uncover a universal tendency towards protein structural complexity that is remarkable.

Description: by Hermann B. Frieboes, Bryan R. Smith, Yao-Li Chuang, Ken Ito, Allison M. Roettgers, Sanjiv S. Gambhir, Vittorio Cristini Non-Hodgkin's lymphoma is a disseminated, highly malignant cancer, with resistance to drug treatment based on molecular- and tissue-scale characteristics that are intricately linked. A critical element of molecular resistance has been traced to the loss of functionality in proteins such as the tumor suppressor p53 . We investigate the tissue-scale physiologic effects of this loss by integrating in vivo and immunohistological data with computational modeling to study the spatiotemporal physical dynamics of lymphoma growth. We compare between drug-sensitive Eμ-myc Arf-/- and drug-resistant Eμ-myc p53-/- lymphoma cell tumors grown in live mice. Initial values for the model parameters are obtained in part by extracting values from the cellular-scale from whole-tumor histological staining of the tumor-infiltrated inguinal lymph node in vivo . We compare model-predicted tumor growth with that observed from intravital microscopy and macroscopic imaging in vivo , finding that the model is able to accurately predict lymphoma growth. A critical physical mechanism underlying drug-resistant phenotypes may be that the Eμ-myc p53-/- cells seem to pack more closely within the tumor than the Eμ-myc Arf-/- cells, thus possibly exacerbating diffusion gradients of oxygen, leading to cell quiescence and hence resistance to cell-cycle specific drugs. Tighter cell packing could also maintain steeper gradients of drug and lead to insufficient toxicity. The transport phenomena within the lymphoma may thus contribute in nontrivial, complex ways to the difference in drug sensitivity between Eμ-myc Arf-/- and Eμ-myc p53-/- tumors, beyond what might be solely expected from loss of functionality at the molecular scale. We conclude that computational modeling tightly integrated with experimental data gives insight into the dynamics of Non-Hodgkin's lymphoma and provides a platform to generate confirmable predictions of tumor growth.

Description: by Yuguo Yu, Thomas S. McTavish, Michael L. Hines, Gordon M. Shepherd, Cesare Valenti, Michele Migliore In the olfactory bulb, lateral inhibition mediated by granule cells has been suggested to modulate the timing of mitral cell firing, thereby shaping the representation of input odorants. Current experimental techniques, however, do not enable a clear study of how the mitral-granule cell network sculpts odor inputs to represent odor information spatially and temporally. To address this critical step in the neural basis of odor recognition, we built a biophysical network model of mitral and granule cells, corresponding to 1/100th of the real system in the rat, and used direct experimental imaging data of glomeruli activated by various odors. The model allows the systematic investigation and generation of testable hypotheses of the functional mechanisms underlying odor representation in the olfactory bulb circuit. Specifically, we demonstrate that lateral inhibition emerges within the olfactory bulb network through recurrent dendrodendritic synapses when constrained by a range of balanced excitatory and inhibitory conductances. We find that the spatio-temporal dynamics of lateral inhibition plays a critical role in building the glomerular-related cell clusters observed in experiments, through the modulation of synaptic weights during odor training. Lateral inhibition also mediates the development of sparse and synchronized spiking patterns of mitral cells related to odor inputs within the network, with the frequency of these synchronized spiking patterns also modulated by the sniff cycle.

Description: by Marco Di Stefano, Angelo Rosa, Vincenzo Belcastro, Diego di Bernardo, Cristian Micheletti The connection between chromatin nuclear organization and gene activity is vividly illustrated by the observation that transcriptional coregulation of certain genes appears to be directly influenced by their spatial proximity. This fact poses the more general question of whether it is at all feasible that the numerous genes that are coregulated on a given chromosome, especially those at large genomic distances, might become proximate inside the nucleus. This problem is studied here using steered molecular dynamics simulations in order to enforce the colocalization of thousands of knowledge-based gene sequences on a model for the gene-rich human chromosome 19. Remarkably, it is found that most () gene pairs can be brought simultaneously into contact. This is made possible by the low degree of intra-chromosome entanglement and the large number of cliques in the gene coregulatory network. A clique is a set of genes coregulated all together as a group. The constrained conformations for the model chromosome 19 are further shown to be organized in spatial macrodomains that are similar to those inferred from recent HiC measurements. The findings indicate that gene coregulation and colocalization are largely compatible and that this relationship can be exploited to draft the overall spatial organization of the chromosome in vivo . The more general validity and implications of these findings could be investigated by applying to other eukaryotic chromosomes the general and transferable computational strategy introduced here.

Description: by Marta I. Garrido, Maneesh Sahani, Raymond J. Dolan We constantly look for patterns in the environment that allow us to learn its key regularities. These regularities are fundamental in enabling us to make predictions about what is likely to happen next. The physiological study of regularity extraction has focused primarily on repetitive sequence-based rules within the sensory environment, or on stimulus-outcome associations in the context of reward-based decision-making. Here we ask whether we implicitly encode non-sequential stochastic regularities, and detect violations therein. We addressed this question using a novel experimental design and both behavioural and magnetoencephalographic (MEG) metrics associated with responses to pure-tone sounds with frequencies sampled from a Gaussian distribution. We observed that sounds in the tail of the distribution evoked a larger response than those that fell at the centre. This response resembled the mismatch negativity (MMN) evoked by surprising or unlikely events in traditional oddball paradigms. Crucially, responses to physically identical outliers were greater when the distribution was narrower. These results show that humans implicitly keep track of the uncertainty induced by apparently random distributions of sensory events. Source reconstruction suggested that the statistical-context-sensitive responses arose in a temporo-parietal network, areas that have been associated with attention orientation to unexpected events. Our results demonstrate a very early neurophysiological marker of the brain's ability to implicitly encode complex statistical structure in the environment. We suggest that this sensitivity provides a computational basis for our ability to make perceptual inferences in noisy environments and to make decisions in an uncertain world.

Description: by Emily V. Chambers, Wendy A. Bickmore, Colin A. Semple Several recent studies have examined different aspects of mammalian higher order chromatin structure – replication timing, lamina association and Hi-C inter-locus interactions — and have suggested that most of these features of genome organisation are conserved over evolution. However, the extent of evolutionary divergence in higher order structure has not been rigorously measured across the mammalian genome, and until now little has been known about the characteristics of any divergent loci present. Here, we generate a dataset combining multiple measurements of chromatin structure and organisation over many embryonic cell types for both human and mouse that, for the first time, allows a comprehensive assessment of the extent of structural divergence between mammalian genomes. Comparison of orthologous regions confirms that all measurable facets of higher order structure are conserved between human and mouse, across the vast majority of the detectably orthologous genome. This broad similarity is observed in spite of many loci possessing cell type specific structures. However, we also identify hundreds of regions (from 100 Kb to 2.7 Mb in size) showing consistent evidence of divergence between these species, constituting at least 10% of the orthologous mammalian genome and encompassing many hundreds of human and mouse genes. These regions show unusual shifts in human GC content, are unevenly distributed across both genomes, and are enriched in human subtelomeric regions. Divergent regions are also relatively enriched for genes showing divergent expression patterns between human and mouse ES cells, implying these regions cause divergent regulation. Particular divergent loci are strikingly enriched in genes implicated in vertebrate development, suggesting important roles for structural divergence in the evolution of mammalian developmental programmes. These data suggest that, though relatively rare in the mammalian genome, divergence in higher order chromatin structure has played important roles during evolution.

Description: by Matthew R. Nassar, Joshua I. Gold Fitting models to behavior is commonly used to infer the latent computational factors responsible for generating behavior. However, the complexity of many behaviors can handicap the interpretation of such models. Here we provide perspectives on problems that can arise when interpreting parameter fits from models that provide incomplete descriptions of behavior. We illustrate these problems by fitting commonly used and neurophysiologically motivated reinforcement-learning models to simulated behavioral data sets from learning tasks. These model fits can pass a host of standard goodness-of-fit tests and other model-selection diagnostics even when the models do not provide a complete description of the behavioral data. We show that such incomplete models can be misleading by yielding biased estimates of the parameters explicitly included in the models. This problem is particularly pernicious when the neglected factors are unknown and therefore not easily identified by model comparisons and similar methods. An obvious conclusion is that a parsimonious description of behavioral data does not necessarily imply an accurate description of the underlying computations. Moreover, general goodness-of-fit measures are not a strong basis to support claims that a particular model can provide a generalized understanding of the computations that govern behavior. To help overcome these challenges, we advocate the design of tasks that provide direct reports of the computational variables of interest. Such direct reports complement model-fitting approaches by providing a more complete, albeit possibly more task-specific, representation of the factors that drive behavior. Computational models then provide a means to connect such task-specific results to a more general algorithmic understanding of the brain.

Description: by Purushottam D. Dixit, Sergei Maslov In addition to their biological function, protein complexes reduce the exposure of the constituent proteins to the risk of undesired oligomerization by reducing the concentration of the free monomeric state. We interpret this reduced risk as a stabilization of the functional state of the protein. We estimate that protein-protein interactions can account for of additional stabilization; a substantial contribution to intrinsic stability. We hypothesize that proteins in the interaction network act as evolutionary capacitors which allows their binding partners to explore regions of the sequence space which correspond to less stable proteins. In the interaction network of baker's yeast, we find that statistically proteins that receive higher energetic benefits from the interaction network are more likely to misfold. A simplified fitness landscape wherein the fitness of an organism is inversely proportional to the total concentration of unfolded proteins provides an evolutionary justification for the proposed trends. We conclude by outlining clear biophysical experiments to test our predictions.

Description: by Andrei Zinovyev, Inna Kuperstein, Emmanuel Barillot, Wolf-Dietrich Heyer Systematic analysis of synthetic lethality (SL) constitutes a critical tool for systems biology to decipher molecular pathways. The most accepted mechanistic explanation of SL is that the two genes function in parallel, mutually compensatory pathways, known as between-pathway SL. However, recent genome-wide analyses in yeast identified a significant number of within-pathway negative genetic interactions. The molecular mechanisms leading to within-pathway SL are not fully understood. Here, we propose a novel mechanism leading to within-pathway SL involving two genes functioning in a single non-essential pathway. This type of SL termed within-reversible-pathway SL involves reversible pathway steps, catalyzed by different enzymes in the forward and backward directions, and kinetic trapping of a potentially toxic intermediate. Experimental data with recombinational DNA repair genes validate the concept. Mathematical modeling recapitulates the possibility of kinetic trapping and revealed the potential contributions of synthetic, dosage-lethal interactions in such a genetic system as well as the possibility of within-pathway positive masking interactions. Analysis of yeast gene interaction and pathway data suggests broad applicability of this novel concept. These observations extend the canonical interpretation of synthetic-lethal or synthetic-sick interactions with direct implications to reconstruct molecular pathways and improve therapeutic approaches to diseases such as cancer.

Description: by Tina Toni, Bruce Tidor Biological systems are inherently variable, with their dynamics influenced by intrinsic and extrinsic sources. These systems are often only partially characterized, with large uncertainties about specific sources of extrinsic variability and biochemical properties. Moreover, it is not yet well understood how different sources of variability combine and affect biological systems in concert. To successfully design biomedical therapies or synthetic circuits with robust performance, it is crucial to account for uncertainty and effects of variability. Here we introduce an efficient modeling and simulation framework to study systems that are simultaneously subject to multiple sources of variability, and apply it to make design decisions on small genetic networks that play a role of basic design elements of synthetic circuits. Specifically, the framework was used to explore the effect of transcriptional and post-transcriptional autoregulation on fluctuations in protein expression in simple genetic networks. We found that autoregulation could either suppress or increase the output variability, depending on specific noise sources and network parameters. We showed that transcriptional autoregulation was more successful than post-transcriptional in suppressing variability across a wide range of intrinsic and extrinsic magnitudes and sources. We derived the following design principles to guide the design of circuits that best suppress variability: (i) high protein cooperativity and low miRNA cooperativity, (ii) imperfect complementarity between miRNA and mRNA was preferred to perfect complementarity, and (iii) correlated expression of mRNA and miRNA – for example, on the same transcript – was best for suppression of protein variability. Results further showed that correlations in kinetic parameters between cells affected the ability to suppress variability, and that variability in transient states did not necessarily follow the same principles as variability in the steady state. Our model and findings provide a general framework to guide design principles in synthetic biology.

Description: by Simon Garnier, Maud Combe, Christian Jost, Guy Theraulaz Interactions between individuals and the structure of their environment play a crucial role in shaping self-organized collective behaviors. Recent studies have shown that ants crossing asymmetrical bifurcations in a network of galleries tend to follow the branch that deviates the least from their incoming direction. At the collective level, the combination of this tendency and the pheromone-based recruitment results in a greater likelihood of selecting the shortest path between the colony's nest and a food source in a network containing asymmetrical bifurcations. It was not clear however what the origin of this behavioral bias is. Here we propose that it results from a simple interaction between the behavior of the ants and the geometry of the network, and that it does not require the ability to measure the angle of the bifurcation. We tested this hypothesis using groups of ant-like robots whose perceptual and cognitive abilities can be fully specified. We programmed them only to lay down and follow light trails, avoid obstacles and move according to a correlated random walk, but not to use more sophisticated orientation methods. We recorded the behavior of the robots in networks of galleries presenting either only symmetrical bifurcations or a combination of symmetrical and asymmetrical bifurcations. Individual robots displayed the same pattern of branch choice as individual ants when crossing a bifurcation, suggesting that ants do not actually measure the geometry of the bifurcations when travelling along a pheromone trail. Finally at the collective level, the group of robots was more likely to select one of the possible shorter paths between two designated areas when moving in an asymmetrical network, as observed in ants. This study reveals the importance of the shape of trail networks for foraging in ants and emphasizes the underestimated role of the geometrical properties of transportation networks in general.

Description: by Michael A. Carlin, Mounya Elhilali The processing characteristics of neurons in the central auditory system are directly shaped by and reflect the statistics of natural acoustic environments, but the principles that govern the relationship between natural sound ensembles and observed responses in neurophysiological studies remain unclear. In particular, accumulating evidence suggests the presence of a code based on sustained neural firing rates, where central auditory neurons exhibit strong, persistent responses to their preferred stimuli. Such a strategy can indicate the presence of ongoing sounds, is involved in parsing complex auditory scenes, and may play a role in matching neural dynamics to varying time scales in acoustic signals. In this paper, we describe a computational framework for exploring the influence of a code based on sustained firing rates on the shape of the spectro-temporal receptive field (STRF), a linear kernel that maps a spectro-temporal acoustic stimulus to the instantaneous firing rate of a central auditory neuron. We demonstrate the emergence of richly structured STRFs that capture the structure of natural sounds over a wide range of timescales, and show how the emergent ensembles resemble those commonly reported in physiological studies. Furthermore, we compare ensembles that optimize a sustained firing code with one that optimizes a sparse code, another widely considered coding strategy, and suggest how the resulting population responses are not mutually exclusive. Finally, we demonstrate how the emergent ensembles contour the high-energy spectro-temporal modulations of natural sounds, forming a discriminative representation that captures the full range of modulation statistics that characterize natural sound ensembles. These findings have direct implications for our understanding of how sensory systems encode the informative components of natural stimuli and potentially facilitate multi-sensory integration.

Description: by James B. Gilmore, Anthony D. Kelleher, David A. Cooper, John M. Murray A recent investigation of the effect of different antiretroviral drug classes on first phase dynamics of HIV RNA plasma virus levels has indicated that drugs acting at stages closer to viral production, such as the integrase inhibitor raltegravir, can produce a steeper first phase decay slope that may not be due to drug efficacy. Moreover it was found that for most drug classes the first phase transitions from a faster (phase IA) to a slightly slower decay region (phase IB) before the start of the usual second phase. Neither of these effects has been explained to date. We use a mathematical model that incorporates the different stages of the HIV viral life cycle in CD4+ T cells: viral entry, reverse transcription, integration, and viral production, to investigate the intracellular HIV mechanisms responsible for these complex plasma virus decay dynamics. We find differences in the phase IA slope across drug classes arise from a higher death rate of cells when they enter the productively infected stage post-integration, with a half-life of approximately 8 hours in this stage, whereas cells in earlier stages of the infection cycle have half-lives similar to uninfected cells. This implies any immune clearance is predominantly limited to the productive infection stage. We also show that the slowing of phase IA to phase IB at day 2 to 4 of monotherapy, depending on drug class, is a result of new rounds of infection. The level at which this slowing occurs is a better indicator of drug efficacy than the slope of the initial decay.

Description: by Christopher P. Said, David J. Heeger Binocular rivalry and cross-orientation suppression are well-studied forms of competition in visual cortex, but models of these two types of competition are in tension with one another. Binocular rivalry occurs during the presentation of dichoptic grating stimuli, where two orthogonal gratings presented separately to the two eyes evoke strong alternations in perceptual dominance. Cross-orientation suppression occurs during the presentation of plaid stimuli, where the responses to a component grating presented to both eyes is weakened by the presence of a superimposed orthogonal grating. Conventional models of rivalry that rely on strong competition between orientation-selective neurons incorrectly predict rivalry between the components of plaids. Lowering the inhibitory weights in such models reduces rivalry for plaids, but also reduces it for dichoptic gratings. Using an exhaustive grid search, we show that this problem cannot be solved simply by adjusting the parameters of the model. Instead, we propose a robust class of models that rely on ocular opponency neurons, previously proposed as a mechanism for efficient stereo coding, to yield rivalry only for dichoptic gratings, not for plaids. This class of models reconciles models of binocular rivalry with the divisive normalization framework that has been used to explain cross-orientation. Our model makes novel predictions that we confirmed with psychophysical tests.

Description: by Annabelle Ballesta, Jonathan Lopez, Nikolay Popgeorgiev, Philippe Gonzalo, Marie Doumic, Germain Gillet Src tyrosine kinases are deregulated in numerous cancers and may favor tumorigenesis and tumor progression. We previously described that Src activation in NIH-3T3 mouse fibroblasts promoted cell resistance to apoptosis. Indeed, Src was found to accelerate the degradation of the pro-apoptotic BH3-only protein Bik and compromised Bax activation as well as subsequent mitochondrial outer membrane permeabilization. The present study undertook a systems biomedicine approach to design optimal anticancer therapeutic strategies using Src-transformed and parental fibroblasts as a biological model. First, a mathematical model of Bik kinetics was designed and fitted to biological data. It guided further experimental investigation that showed that Bik total amount remained constant during staurosporine exposure, and suggested that Bik protein might undergo activation to induce apoptosis. Then, a mathematical model of the mitochondrial pathway of apoptosis was designed and fitted to experimental results. It showed that Src inhibitors could circumvent resistance to apoptosis in Src-transformed cells but gave no specific advantage to parental cells. In addition, it predicted that inhibitors of Bcl-2 antiapoptotic proteins such as ABT-737 should not be used in this biological system in which apoptosis resistance relied on the deficiency of an apoptosis accelerator but not on the overexpression of an apoptosis inhibitor, which was experimentally verified. Finally, we designed theoretically optimal therapeutic strategies using the data-calibrated model. All of them relied on the observed Bax overexpression in Src-transformed cells compared to parental fibroblasts. Indeed, they all involved Bax downregulation such that Bax levels would still be high enough to induce apoptosis in Src-transformed cells but not in parental ones. Efficacy of this counterintuitive therapeutic strategy was further experimentally validated. Thus, the use of Bax inhibitors might be an unexpected way to specifically target cancer cells with deregulated Src tyrosine kinase activity.

Description: by Masayo Inoue, Kunihiko Kaneko Cells generally adapt to environmental changes by first exhibiting an immediate response and then gradually returning to their original state to achieve homeostasis. Although simple network motifs consisting of a few genes have been shown to exhibit such adaptive dynamics, they do not reflect the complexity of real cells, where the expression of a large number of genes activates or represses other genes, permitting adaptive behaviors. Here, we investigated the responses of gene regulatory networks containing many genes that have undergone numerical evolution to achieve high fitness due to the adaptive response of only a single target gene; this single target gene responds to changes in external inputs and later returns to basal levels. Despite setting a single target, most genes showed adaptive responses after evolution. Such adaptive dynamics were not due to common motifs within a few genes; even without such motifs, almost all genes showed adaptation, albeit sometimes partial adaptation, in the sense that expression levels did not always return to original levels. The genes split into two groups: genes in the first group exhibited an initial increase in expression and then returned to basal levels, while genes in the second group exhibited the opposite changes in expression. From this model, genes in the first group received positive input from other genes within the first group, but negative input from genes in the second group, and vice versa. Thus, the adaptation dynamics of genes from both groups were consolidated. This cooperative adaptive behavior was commonly observed if the number of genes involved was larger than the order of ten. These results have implications in the collective responses of gene expression networks in microarray measurements of yeast Saccharomyces cerevisiae and the significance to the biological homeostasis of systems with many components.

Description: by Adria Carbo, Raquel Hontecillas, Barbara Kronsteiner, Monica Viladomiu, Mireia Pedragosa, Pinyi Lu, Casandra W. Philipson, Stefan Hoops, Madhav Marathe, Stephen Eubank, Keith Bisset, Katherine Wendelsdorf, Abdul Jarrah, Yongguo Mei, Josep Bassaganya-Riera Differentiation of CD4+ T cells into effector or regulatory phenotypes is tightly controlled by the cytokine milieu, complex intracellular signaling networks and numerous transcriptional regulators. We combined experimental approaches and computational modeling to investigate the mechanisms controlling differentiation and plasticity of CD4+ T cells in the gut of mice. Our computational model encompasses the major intracellular pathways involved in CD4+ T cell differentiation into T helper 1 (Th1), Th2, Th17 and induced regulatory T cells (iTreg). Our modeling efforts predicted a critical role for peroxisome proliferator-activated receptor gamma (PPARγ) in modulating plasticity between Th17 and iTreg cells. PPARγ regulates differentiation, activation and cytokine production, thereby controlling the induction of effector and regulatory responses, and is a promising therapeutic target for dysregulated immune responses and inflammation. Our modeling efforts predict that following PPARγ activation, Th17 cells undergo phenotype switch and become iTreg cells. This prediction was validated by results of adoptive transfer studies showing an increase of colonic iTreg and a decrease of Th17 cells in the gut mucosa of mice with colitis following pharmacological activation of PPARγ. Deletion of PPARγ in CD4+ T cells impaired mucosal iTreg and enhanced colitogenic Th17 responses in mice with CD4+ T cell-induced colitis. Thus, for the first time we provide novel molecular evidence in vivo demonstrating that PPARγ in addition to regulating CD4+ T cell differentiation also plays a major role controlling Th17 and iTreg plasticity in the gut mucosa.

Description: by Domenico Campolo, Ferdinan Widjaja, Hong Xu, Wei Tech Ang, Etienne Burdet This work introduces a coordinate-independent method to analyse movement variability of tasks performed with hand-held tools, such as a pen or a surgical scalpel. We extend the classical uncontrolled manifold (UCM) approach by exploiting the geometry of rigid body motions, used to describe tool configurations. In particular, we analyse variability during a static pointing task with a hand-held tool, where subjects are asked to keep the tool tip in steady contact with another object. In this case the tool is redundant with respect to the task, as subjects control position/orientation of the tool, i.e. 6 degrees-of-freedom (dof), to maintain the tool tip position (3dof) steady. To test the new method, subjects performed a pointing task with and without arm support. The additional dof introduced in the unsupported condition, injecting more variability into the system, represented a resource to minimise variability in the task space via coordinated motion. The results show that all of the seven subjects channeled more variability along directions not directly affecting the task (UCM), consistent with previous literature but now shown in a coordinate-independent way. Variability in the unsupported condition was only slightly larger at the endpoint but much larger in the UCM.

Description: by Lin-Tai Da, Fátima Pardo Avila, Dong Wang, Xuhui Huang The dynamics of the PP i release during the transcription elongation of bacterial RNA polymerase and its effects on the Trigger Loop (TL) opening motion are still elusive. Here, we built a Markov State Model (MSM) from extensive all-atom molecular dynamics (MD) simulations to investigate the mechanism of the PP i release. Our MSM has identified a simple two-state mechanism for the PP i release instead of a more complex four-state mechanism observed in RNA polymerase II (Pol II). We observed that the PP i release in bacterial RNA polymerase occurs at sub-microsecond timescale, which is ∼3-fold faster than that in Pol II. After escaping from the active site, the (Mg-PP i ) 2− group passes through a single elongated metastable region where several positively charged residues on the secondary channel provide favorable interactions. Surprisingly, we found that the PP i release is not coupled with the TL unfolding but correlates tightly with the side-chain rotation of the TL residue R1239. Our work sheds light on the dynamics underlying the transcription elongation of the bacterial RNA polymerase.

Description: by Izzet B. Yildiz, Katharina von Kriegstein, Stefan J. Kiebel Our knowledge about the computational mechanisms underlying human learning and recognition of sound sequences, especially speech, is still very limited. One difficulty in deciphering the exact means by which humans recognize speech is that there are scarce experimental findings at a neuronal, microscopic level. Here, we show that our neuronal-computational understanding of speech learning and recognition may be vastly improved by looking at an animal model, i.e., the songbird, which faces the same challenge as humans: to learn and decode complex auditory input, in an online fashion. Motivated by striking similarities between the human and songbird neural recognition systems at the macroscopic level, we assumed that the human brain uses the same computational principles at a microscopic level and translated a birdsong model into a novel human sound learning and recognition model with an emphasis on speech. We show that the resulting Bayesian model with a hierarchy of nonlinear dynamical systems can learn speech samples such as words rapidly and recognize them robustly, even in adverse conditions. In addition, we show that recognition can be performed even when words are spoken by different speakers and with different accents—an everyday situation in which current state-of-the-art speech recognition models often fail. The model can also be used to qualitatively explain behavioral data on human speech learning and derive predictions for future experiments.

Description: by John C. Williams, Jianjin Xu, Zhongju Lu, Aleksandra Klimas, Xuxin Chen, Christina M. Ambrosi, Ira S. Cohen, Emilia Entcheva Channelrhodospin-2 (ChR2), a light-sensitive ion channel, and its variants have emerged as new excitatory optogenetic tools not only in neuroscience, but also in other areas, including cardiac electrophysiology. An accurate quantitative model of ChR2 is necessary for in silico prediction of the response to optical stimulation in realistic tissue/organ settings. Such a model can guide the rational design of new ion channel functionality tailored to different cell types/tissues. Focusing on one of the most widely used ChR2 mutants (H134R) with enhanced current, we collected a comprehensive experimental data set of the response of this ion channel to different irradiances and voltages, and used these data to develop a model of ChR2 with empirically-derived voltage- and irradiance- dependence, where parameters were fine-tuned via simulated annealing optimization. This ChR2 model offers: 1) accurate inward rectification in the current-voltage response across irradiances; 2) empirically-derived voltage- and light-dependent kinetics (activation, deactivation and recovery from inactivation); and 3) accurate amplitude and morphology of the response across voltage and irradiance settings. Temperature-scaling factors (Q 10 ) were derived and model kinetics was adjusted to physiological temperatures. Using optical action potential clamp, we experimentally validated model-predicted ChR2 behavior in guinea pig ventricular myocytes. The model was then incorporated in a variety of cardiac myocytes, including human ventricular, atrial and Purkinje cell models. We demonstrate the ability of ChR2 to trigger action potentials in human cardiomyocytes at relatively low light levels, as well as the differential response of these cells to light, with the Purkinje cells being most easily excitable and ventricular cells requiring the highest irradiance at all pulse durations. This new experimentally-validated ChR2 model will facilitate virtual experimentation in neural and cardiac optogenetics at the cell and organ level and provide guidance for the development of in vivo tools.

Description: by Kobi Snitz, Adi Yablonka, Tali Weiss, Idan Frumin, Rehan M. Khan, Noam Sobel To understand the brain mechanisms of olfaction we must understand the rules that govern the link between odorant structure and odorant perception. Natural odors are in fact mixtures made of many molecules, and there is currently no method to look at the molecular structure of such odorant-mixtures and predict their smell. In three separate experiments, we asked 139 subjects to rate the pairwise perceptual similarity of 64 odorant-mixtures ranging in size from 4 to 43 mono-molecular components. We then tested alternative models to link odorant-mixture structure to odorant-mixture perceptual similarity. Whereas a model that considered each mono-molecular component of a mixture separately provided a poor prediction of mixture similarity, a model that represented the mixture as a single structural vector provided consistent correlations between predicted and actual perceptual similarity (r≥0.49, p

Description: by Laurence T. Hunt, Mark W. Woolrich, Matthew F. S. Rushworth, Timothy E. J. Behrens A central question in cognitive neuroscience regards the means by which options are compared and decisions are resolved during value-guided choice. It is clear that several component processes are needed; these include identifying options, a value-based comparison, and implementation of actions to execute the decision. What is less clear is the temporal precedence and functional organisation of these component processes in the brain. Competing models of decision making have proposed that value comparison may occur in the space of alternative actions, or in the space of abstract goods. We hypothesized that the signals observed might in fact depend upon the framing of the decision. We recorded magnetoencephalographic data from humans performing value-guided choices in which two closely related trial types were interleaved. In the first trial type, each option was revealed separately, potentially causing subjects to estimate each action's value as it was revealed and perform comparison in action-space. In the second trial type, both options were presented simultaneously, potentially leading to comparison in abstract goods-space prior to commitment to a specific action. Distinct activity patterns (in distinct brain regions) on the two trial types demonstrated that the observed frame of reference used for decision making indeed differed, despite the information presented being formally identical, between the two trial types. This provides a potential reconciliation of conflicting accounts of value-guided choice.

Description: by Ryan Tasseff, Anjali Bheda-Malge, Teresa DiColandrea, Charles C. Bascom, Robert J. Isfort, Richard Gelinas The hair cycle is a dynamic process where follicles repeatedly move through phases of growth, retraction, and relative quiescence. This process is an example of temporal and spatial biological complexity. Understanding of the hair cycle and its regulation would shed light on many other complex systems relevant to biological and medical research. Currently, a systematic characterization of gene expression and summarization within the context of a mathematical model is not yet available. Given the cyclic nature of the hair cycle, we felt it was important to consider a subset of genes with periodic expression. To this end, we combined several mathematical approaches with high-throughput, whole mouse skin, mRNA expression data to characterize aspects of the dynamics and the possible cell populations corresponding to potentially periodic patterns. In particular two gene clusters, demonstrating properties of out-of-phase synchronized expression, were identified. A mean field, phase coupled oscillator model was shown to quantitatively recapitulate the synchronization observed in the data. Furthermore, we found only one configuration of positive-negative coupling to be dynamically stable, which provided insight on general features of the regulation. Subsequent bifurcation analysis was able to identify and describe alternate states based on perturbation of system parameters. A 2-population mixture model and cell type enrichment was used to associate the two gene clusters to features of background mesenchymal populations and rapidly expanding follicular epithelial cells. Distinct timing and localization of expression was also shown by RNA and protein imaging for representative genes. Taken together, the evidence suggests that synchronization between expanding epithelial and background mesenchymal cells may be maintained, in part, by inhibitory regulation, and potential mediators of this regulation were identified. Furthermore, the model suggests that impairing this negative regulation will drive a bifurcation which may represent transition into a pathological state such as hair miniaturization.

Description: by Varsha Dhankani, J. Nathan Kutz, Joshua T. Schiffer Herpes simplex virus-2 (HSV-2) is a chronic reactivating infection that leads to recurrent shedding episodes in the genital tract. A minority of episodes are prolonged, and associated with development of painful ulcers. However, currently, available tools poorly predict viral trajectories and timing of reactivations in infected individuals. We employed principal components analysis (PCA) and singular value decomposition (SVD) to interpret HSV-2 genital tract shedding time series data, as well as simulation output from a stochastic spatial mathematical model. Empirical and model-derived, time-series data gathered over 〉30 days consists of multiple complex episodes that could not be reduced to a manageable number of descriptive features with PCA and SVD. However, single HSV-2 shedding episodes, even those with prolonged duration and complex morphologies consisting of multiple erratic peaks, were consistently described using a maximum of four dominant features. Modeled and clinical episodes had equivalent distributions of dominant features, implying similar dynamics in real and simulated episodes. We applied linear discriminant analysis (LDA) to simulation output and identified that local immune cell density at the viral reactivation site had a predictive effect on episode duration, though longer term shedding suggested chaotic dynamics and could not be predicted based on spatial patterns of immune cell density. These findings suggest that HSV-2 shedding patterns within an individual are impossible to predict over weeks or months, and that even highly complex single HSV-2 episodes can only be partially predicted based on spatial distribution of immune cell density.

Description: by Marcin J. Skwark, Daniele Raimondi, Mirco Michel, Arne Elofsson Given sufficient large protein families, and using a global statistical inference approach, it is possible to obtain sufficient accuracy in protein residue contact predictions to predict the structure of many proteins. However, these approaches do not consider the fact that the contacts in a protein are neither randomly, nor independently distributed, but actually follow precise rules governed by the structure of the protein and thus are interdependent. Here, we present PconsC2, a novel method that uses a deep learning approach to identify protein-like contact patterns to improve contact predictions. A substantial enhancement can be seen for all contacts independently on the number of aligned sequences, residue separation or secondary structure type, but is largest for β -sheet containing proteins. In addition to being superior to earlier methods based on statistical inferences, in comparison to state of the art methods using machine learning, PconsC2 is superior for families with more than 100 effective sequence homologs. The improved contact prediction enables improved structure prediction.

Description: by Veronika Boskova, Sebastian Bonhoeffer, Tanja Stadler Quantifying epidemiological dynamics is crucial for understanding and forecasting the spread of an epidemic. The coalescent and the birth-death model are used interchangeably to infer epidemiological parameters from the genealogical relationships of the pathogen population under study, which in turn are inferred from the pathogen genetic sequencing data. To compare the performance of these widely applied models, we performed a simulation study. We simulated phylogenetic trees under the constant rate birth-death model and the coalescent model with a deterministic exponentially growing infected population. For each tree, we re-estimated the epidemiological parameters using both a birth-death and a coalescent based method, implemented as an MCMC procedure in BEAST v2.0. In our analyses that estimate the growth rate of an epidemic based on simulated birth-death trees, the point estimates such as the maximum a posteriori/maximum likelihood estimates are not very different. However, the estimates of uncertainty are very different. The birth-death model had a higher coverage than the coalescent model, i.e. contained the true value in the highest posterior density (HPD) interval more often (2–13% vs. 31–75% error). The coverage of the coalescent decreases with decreasing basic reproductive ratio and increasing sampling probability of infecteds. We hypothesize that the biases in the coalescent are due to the assumption of deterministic rather than stochastic population size changes. Both methods performed reasonably well when analyzing trees simulated under the coalescent. The methods can also identify other key epidemiological parameters as long as one of the parameters is fixed to its true value. In summary, when using genetic data to estimate epidemic dynamics, our results suggest that the birth-death method will be less sensitive to population fluctuations of early outbreaks than the coalescent method that assumes a deterministic exponentially growing infected population.

Description: by Junaid Hassan, Linda L. Bergaust, I. David Wheat, Lars R. Bakken In response to impending anoxic conditions, denitrifying bacteria sustain respiratory metabolism by producing enzymes for reducing nitrogen oxyanions/-oxides (NO x ) to N 2 (denitrification). Since denitrifying bacteria are non-fermentative, the initial production of denitrification proteome depends on energy from aerobic respiration. Thus, if a cell fails to synthesise a minimum of denitrification proteome before O 2 is completely exhausted, it will be unable to produce it later due to energy-limitation. Such entrapment in anoxia is recently claimed to be a major phenomenon in batch cultures of the model organism Paracoccus denitrificans on the basis of measured e − -flow rates to O 2 and NO x . Here we constructed a dynamic model and explicitly simulated actual kinetics of recruitment of the cells to denitrification to directly and more accurately estimate the recruited fraction (). Transcription of nirS is pivotal for denitrification, for it triggers a cascade of events leading to the synthesis of a full-fledged denitrification proteome. The model is based on the hypothesis that nirS has a low probability (, h −1 ) of initial transcription, but once initiated, the transcription is greatly enhanced through positive feedback by NO, resulting in the recruitment of the transcribing cell to denitrification. We assume that the recruitment is initiated as [O 2 ] falls below a critical threshold and terminates (assuming energy-limitation) as [O 2 ] exhausts. With  = 0.005 h −1 , the model robustly simulates observed denitrification kinetics for a range of culture conditions. The resulting (fraction of the cells recruited to denitrification) falls within 0.038–0.161. In contrast, if the recruitment of the entire population is assumed, the simulated denitrification kinetics deviate grossly from those observed. The phenomenon can be understood as a ‘bet-hedging strategy’: switching to denitrification is a gain if anoxic spell lasts long but is a waste of energy if anoxia turns out to be a ‘false alarm’.

Description: by Thomas R. Caulfield, Fabienne C. Fiesel, Elisabeth L. Moussaud-Lamodière, Daniel F. A. R. Dourado, Samuel C. Flores, Wolfdieter Springer Loss-of-function mutations in PINK1 or PARKIN are the most common causes of autosomal recessive Parkinson's disease. Both gene products, the Ser/Thr kinase PINK1 and the E3 Ubiquitin ligase Parkin, functionally cooperate in a mitochondrial quality control pathway. Upon stress, PINK1 activates Parkin and enables its translocation to and ubiquitination of damaged mitochondria to facilitate their clearance from the cell. Though PINK1-dependent phosphorylation of Ser65 is an important initial step, the molecular mechanisms underlying the activation of Parkin's enzymatic functions remain unclear. Using molecular modeling, we generated a complete structural model of human Parkin at all atom resolution. At steady state, the Ub ligase is maintained inactive in a closed, auto-inhibited conformation that results from intra-molecular interactions. Evidently, Parkin has to undergo major structural rearrangements in order to unleash its catalytic activity. As a spark, we have modeled PINK1-dependent Ser65 phosphorylation in silico and provide the first molecular dynamics simulation of Parkin conformations along a sequential unfolding pathway that could release its intertwined domains and enable its catalytic activity. We combined free (unbiased) molecular dynamics simulation, Monte Carlo algorithms, and minimal-biasing methods with cell-based high content imaging and biochemical assays. Phosphorylation of Ser65 results in widening of a newly defined cleft and dissociation of the regulatory N-terminal UBL domain. This motion propagates through further opening conformations that allow binding of an Ub-loaded E2 co-enzyme. Subsequent spatial reorientation of the catalytic centers of both enzymes might facilitate the transfer of the Ub moiety to charge Parkin. Our structure-function study provides the basis to elucidate regulatory mechanisms and activity of the neuroprotective Parkin. This may open up new avenues for the development of small molecule Parkin activators through targeted drug design.

Description: by The PLOS Computational Biology Staff

Description: by The PLOS Computational Biology Staff

Description: by Adviti Naik, Damjana Rozman, Aleš Belič Current state-of-the-art mathematical models to investigate complex biological processes, in particular liver-associated pathologies, have limited expansiveness, flexibility, representation of integrated regulation and rely on the availability of detailed kinetic data. We generated the SteatoNet, a multi-pathway, multi-tissue model and in silico platform to investigate hepatic metabolism and its associated deregulations. SteatoNet is based on object-oriented modelling, an approach most commonly applied in automotive and process industries, whereby individual objects correspond to functional entities. Objects were compiled to feature two novel hepatic modelling aspects: the interaction of hepatic metabolic pathways with extra-hepatic tissues and the inclusion of transcriptional and post-transcriptional regulation. SteatoNet identification at normalised steady state circumvents the need for constraining kinetic parameters. Validation and identification of flux disturbances that have been proven experimentally in liver patients and animal models highlights the ability of SteatoNet to effectively describe biological behaviour. SteatoNet identifies crucial pathway branches (transport of glucose, lipids and ketone bodies) where changes in flux distribution drive the healthy liver towards hepatic steatosis, the primary stage of non-alcoholic fatty liver disease. Cholesterol metabolism and its transcription regulators are highlighted as novel steatosis factors. SteatoNet thus serves as an intuitive in silico platform to identify systemic changes associated with complex hepatic metabolic disorders.

Description: by Xiaolin Tang, Mourad Bendjennat, Saveez Saffarian Vesicular stomatitis virus (VSV) is the prototype for negative sense non segmented (NNS) RNA viruses which include potent human and animal pathogens such as Rabies, Ebola and measles. The polymerases of NNS RNA viruses only initiate transcription at or near the 3′ end of their genome template. We measured the dissociation constant of VSV polymerases from their whole genome template to be 20 pM. Given this low dissociation constant, initiation and sustainability of transcription becomes nontrivial. To explore possible mechanisms, we simulated the first hour of transcription using Monte Carlo methods and show that a one-time initial dissociation of all polymerases during entry is not sufficient to sustain transcription. We further show that efficient transcription requires a sliding mechanism for non-transcribing polymerases and can be realized with different polymerase-polymerase interactions and distinct template topologies. In conclusion, we highlight a model in which collisions between transcribing and sliding non-transcribing polymerases result in release of the non-transcribing polymerases allowing for redistribution of polymerases between separate templates during transcription and suggest specific experiments to further test these mechanisms.

Description: by The PLOS Computational Biology Staff

Description: by Ganesh Shahane, Chirag Parsania, Durba Sengupta, Manali Joshi The human β 2 -adrenergic receptor (β 2 AR), a member of the G-protein coupled receptor (GPCR) family, is expressed in bronchial smooth muscle cells. Upon activation by agonists, β 2 AR causes bronchodilation and relief in asthma patients. The N-terminal polymorphism of β 2 AR at the 16 th position, Arg16Gly, has warranted a lot of attention since it is linked to variations in response to albuterol (agonist) treatment. Although the β 2 AR is one of the well-studied GPCRs, the N-terminus which harbors this mutation, is absent in all available experimental structures. The goal of this work was to study the molecular level differences between the N-terminal variants using structural modeling and atomistic molecular dynamics simulations. Our simulations reveal that the N-terminal region of the Arg variant shows greater dynamics than the Gly variant, leading to differential placement. Further, the position and dynamics of the N-terminal region, further, affects the ligand binding-site accessibility. Interestingly, long-range effects are also seen at the ligand binding site, which is marginally larger in the Gly as compared to the Arg variant resulting in the preferential docking of albuterol to the Gly variant. This study thus reveals key differences between the variants providing a molecular framework towards understanding the variable drug response in asthma patients.

Description: by Krishna Praneeth Kilambi, Kavan Reddy, Jeffrey J. Gray Protein-protein interactions depend on a host of environmental factors. Local pH conditions influence the interactions through the protonation states of the ionizable residues that can change upon binding. In this work, we present a pH-sensitive docking approach, pHDock, that can sample side-chain protonation states of five ionizable residues (Asp, Glu, His, Tyr, Lys) on-the-fly during the docking simulation. pHDock produces successful local docking funnels in approximately half (79/161) the protein complexes, including 19 cases where standard RosettaDock fails. pHDock also performs better than the two control cases comprising docking at pH 7.0 or using fixed, predetermined protonation states. On average, the top-ranked pHDock structures have lower interface RMSDs and recover more native interface residue-residue contacts and hydrogen bonds compared to RosettaDock. Addition of backbone flexibility using a computationally-generated conformational ensemble further improves native contact and hydrogen bond recovery in the top-ranked structures. Although pHDock is designed to improve docking, it also successfully predicts a large pH-dependent binding affinity change in the Fc–FcRn complex, suggesting that it can be exploited to improve affinity predictions. The approaches in the study contribute to the goal of structural simulations of whole-cell protein-protein interactions including all the environmental factors, and they can be further expanded for pH-sensitive protein design.

Description: by Tadashi Ando, Jeffrey Skolnick DNA binding proteins efficiently search for their cognitive sites on long genomic DNA by combining 3D diffusion and 1D diffusion (sliding) along the DNA. Recent experimental results and theoretical analyses revealed that the proteins show a rotation-coupled sliding along DNA helical pitch. Here, we performed Brownian dynamics simulations using newly developed coarse-grained protein and DNA models for evaluating how hydrodynamic interactions between the protein and DNA molecules, binding affinity of the protein to DNA, and DNA fluctuations affect the one dimensional diffusion of the protein on the DNA. Our results indicate that intermolecular hydrodynamic interactions reduce 1D diffusivity by 30%. On the other hand, structural fluctuations of DNA give rise to steric collisions between the CG-proteins and DNA, resulting in faster 1D sliding of the protein. Proteins with low binding affinities consistent with experimental estimates of non-specific DNA binding show hopping along the CG-DNA. This hopping significantly increases sliding speed. These simulation studies provide additional insights into the mechanism of how DNA binding proteins find their target sites on the genome.

Description: by Krzysztof Puszynski, Alberto Gandolfi, Alberto d'Onofrio In this work we investigate, by means of a computational stochastic model, how tumor cells with wild-type p53 gene respond to the drug Nutlin, an agent that interferes with the Mdm2-mediated p53 regulation. In particular, we show how the stochastic gene-switching controlled by p53 can explain experimental dose-response curves, i.e., the observed inter-cell variability of the cell viability under Nutlin action. The proposed model describes in some detail the regulation network of p53, including the negative feedback loop mediated by Mdm2 and the positive loop mediated by PTEN, as well as the reversible inhibition of Mdm2 caused by Nutlin binding. The fate of the individual cell is assumed to be decided by the rising of nuclear-phosphorylated p53 over a certain threshold. We also performed in silico experiments to evaluate the dose-response curve after a single drug dose delivered in mice, or after its fractionated administration. Our results suggest that dose-splitting may be ineffective at low doses and effective at high doses. This complex behavior can be due to the interplay among the existence of a threshold on the p53 level for its cell activity, the nonlinearity of the relationship between the bolus dose and the peak of active p53, and the relatively fast elimination of the drug.

Description: by Sven Jahnke, Raoul-Martin Memmesheimer, Marc Timme Reliable signal transmission constitutes a key requirement for neural circuit function. The propagation of synchronous pulse packets through recurrent circuits is hypothesized to be one robust form of signal transmission and has been extensively studied in computational and theoretical works. Yet, although external or internally generated oscillations are ubiquitous across neural systems, their influence on such signal propagation is unclear. Here we systematically investigate the impact of oscillations on propagating synchrony. We find that for standard, additive couplings and a net excitatory effect of oscillations, robust propagation of synchrony is enabled in less prominent feed-forward structures than in systems without oscillations. In the presence of non-additive coupling (as mediated by fast dendritic spikes), even balanced oscillatory inputs may enable robust propagation. Here, emerging resonances create complex locking patterns between oscillations and spike synchrony. Interestingly, these resonances make the circuits capable of selecting specific pathways for signal transmission. Oscillations may thus promote reliable transmission and, in co-action with dendritic nonlinearities, provide a mechanism for information processing by selectively gating and routing of signals. Our results are of particular interest for the interpretation of sharp wave/ripple complexes in the hippocampus, where previously learned spike patterns are replayed in conjunction with global high-frequency oscillations. We suggest that the oscillations may serve to stabilize the replay.

Description: by Michael M. H. Graf, Lin Zhixiong, Urban Bren, Dietmar Haltrich, Wilfred F. van Gunsteren, Chris Oostenbrink The flavoenzyme pyranose dehydrogenase (PDH) from the litter decomposing fungus Agaricus meleagris oxidizes many different carbohydrates occurring during lignin degradation. This promiscuous substrate specificity makes PDH a promising catalyst for bioelectrochemical applications. A generalized approach to simulate all 32 possible aldohexopyranoses in the course of one or a few molecular dynamics (MD) simulations is reported. Free energy calculations according to the one-step perturbation (OSP) method revealed the solvation free energies (ΔG solv ) of all 32 aldohexopyranoses in water, which have not yet been reported in the literature. The free energy difference between β- and α-anomers (ΔG β-α ) of all d-stereoisomers in water were compared to experimental values with a good agreement. Moreover, the free-energy differences (ΔG) of the 32 stereoisomers bound to PDH in two different poses were calculated from MD simulations. The relative binding free energies (ΔΔG bind ) were calculated and, where available, compared to experimental values, approximated from K m values. The agreement was very good for one of the poses, in which the sugars are positioned in the active site for oxidation at C1 or C2. Distance analysis between hydrogens of the monosaccharide and the reactive N5-atom of the flavin adenine dinucleotide (FAD) revealed that oxidation is possible at HC1 or HC2 for pose A, and at HC3 or HC4 for pose B. Experimentally detected oxidation products could be rationalized for the majority of monosaccharides by combining ΔΔG bind and a reweighted distance analysis. Furthermore, several oxidation products were predicted for sugars that have not yet been tested experimentally, directing further analyses. This study rationalizes the relationship between binding free energies and substrate promiscuity in PDH, providing novel insights for its applicability in bioelectrochemistry. The results suggest that a similar approach could be applied to study promiscuity of other enzymes.