Power Law Tails In Phylogenetic Systems

Covariance analysis of protein sequence alignments uses coevolving pairs of sequence positions to predict features of protein structure and function. However, current methods ignore the phylogenetic relationships between sequences, potentially corrupting the identification of covarying positions. Here, we use random matrix theory to demonstrate the existence of a power law tail that distinguishes the spectrum of covariance caused by phylogeny from that caused by structural interactions. The power law is essentially independent of the phylogenetic tree topology, depending on just two parameters - the sequence length, and the average branch length. We demonstrate that these power law tails are ubiquitous in the large protein sequence alignments used to predict contacts in 3D structure, as predicted by our theory. This suggests that to decouple phylogenetic effects from the interactions between sequence distal sites that control biological function, it is necessary to remove or downweight the eigenvectors of the covariance matrix with largest eigenvalues. We confirm that truncating these eigenvectors improves contact prediction.This work was supported by a Next Generation fellowship (to L.J.C.), a Marie Curie Career Integration Grant [Evo-Couplings, Grant 631609], and an Engineering and Physical Sciences Research Council PhD studentship (to C.Q.)

Predicting protein-ligand affinity with a random matrix framework

Rapid determination of whether a candidate compound will bind to a particular target receptor remains a stumbling block in drug discovery. We use an approach inspired by random matrix theory to decompose the known ligand set of a target in terms of orthogonal "signals" of salient chemical features, and distinguish these from the much larger set of ligand chemical features that are not relevant for binding to that particular target receptor. After removing the noise caused by finite sampling, we show that the similarity of an unknown ligand to the remaining, cleaned chemical features is a robust predictor of ligand-target affinity, performing as well or better than any algorithm in the published literature. We interpret our algorithm as deriving a model for the binding energy between a target receptor and the set of known ligands, where the underlying binding energy model is related to the classic Ising model in statistical physics.This research was funded by a grant from Roche Pharmaceuticals. A.A.L. acknowledges the support of a Fulbright Fellowship. L.J.C. was supported by a Next Generation Fellowship, and a Marie Curie Career Integration Grant (Evo-Couplings, Grant 631609). M.P.B. is an investigator of the Simons Foundation, and also acknowledges support from the National Science Foundation through DMS-1411694

Intracellular Electric Field and pH Optimize Protein Localization and Movement

Mammalian cell function requires timely and accurate transmission of information from the cell membrane (CM) to the nucleus (N). These pathways have been intensively investigated and many critical components and interactions have been identified. However, the physical forces that control movement of these proteins have received scant attention. Thus, transduction pathways are typically presented schematically with little regard to spatial constraints that might affect the underlying dynamics necessary for protein-protein interactions and molecular movement from the CM to the N. We propose messenger protein localization and movements are highly regulated and governed by Coulomb interactions between: 1. A recently discovered, radially directed E-field from the NM into the CM and 2. Net protein charge determined by its isoelectric point, phosphorylation state, and the cytosolic pH. These interactions, which are widely applied in elecrophoresis, provide a previously unknown mechanism for localization of messenger proteins within the cytoplasm as well as rapid shuttling between the CM and N. Here we show these dynamics optimize the speed, accuracy and efficiency of transduction pathways even allowing measurement of the location and timing of ligand binding at the CM –previously unknown components of intracellular information flow that are, nevertheless, likely necessary for detecting spatial gradients and temporal fluctuations in ligand concentrations within the environment. The model has been applied to the RAF-MEK-ERK pathway and scaffolding protein KSR1 using computer simulations and in-vitro experiments. The computer simulations predicted distinct distributions of phosphorylated and unphosphorylated components of this transduction pathway which were experimentally confirmed in normal breast epithelial cells (HMEC)

Plants used during pregnancy, childbirth and postpartum healthcare in Lao PDR: A comparative study of the Brou, Saek and Kry ethnic groups

<p>Abstract</p> <p>Background</p> <p>In many Southeast Asian cultures the activities and diet during the postpartum period are culturally dictated and a period of confinement is observed. Plants play an important role in recovery during the postpartum period in diet, traditional medicine, steam bath and mother roasting (where mother and child placed on a bed above a brazier with charcoal embers on which aromatic plants are laid). This research focuses on the use of plants during pregnancy, parturition, postpartum recovery and infant healthcare among three ethnic groups, the Brou, Saek and Kry. It aims to identify culturally important traditions that may facilitate implementation of culturally appropriate healthcare.</p> <p>Methods</p> <p>Data were collected in 10 different villages in Khammouane province, Lao PDR, through group and individual interviews with women by female interviewers.</p> <p>Results</p> <p>A total of 55 different plant species are used in women's healthcare, of which over 90% are used in postpartum recovery. Consensus Analysis rejects the hypothesis that the three ethnic groups belong to a single culture for postpartum plant use, and multidimensional scaling reveals non-overlapping clusters per ethnic group.</p> <p>Conclusion</p> <p>Medicinal plant use is common among the Brou, Saek and Kry to facilitate childbirth, alleviate menstruation problems, assist recovery after miscarriage, mitigate postpartum haemorrhage, aid postpartum recovery, and for use in infant care. The wealth of novel insights into plant use and preparation will help to understand culturally important practices such as confinement, dietary restrictions, mother roasting and herbal steam baths and their incorporation into modern healthcare.</p

Close-range remote sensing of Saturn's rings during Cassini's ring-grazing orbits and Grand Finale

Saturn’s rings are an accessible exemplar of an astrophysical disk, tracing the Saturn system’s dynamical processes and history. We present close-range remote-sensing observations of the main rings from the Cassini spacecraft. We find detailed sculpting of the rings by embedded masses, and banded texture belts throughout the rings. Saturn-orbiting streams of material impact the F ring. There are fine-scaled correlations among optical depth, spectral properties, and temperature in the B ring, but anticorrelations within strong density waves in the A ring. There is no spectral distinction between plateaux and the rest of the C ring, whereas the region outward of the Keeler gap is spectrally distinct from nearby regions. These results likely indicate that radial stratification of particle physical properties, rather than compositional differences, is responsible for producing these ring structures. © 2019 American Association for the Advancement of Science. All rights reserved

Neuropeptide Signaling Differentially Affects Phase Maintenance and Rhythm Generation in SCN and Extra-SCN Circadian Oscillators

Circadian rhythms in physiology and behavior are coordinated by the brain's dominant circadian pacemaker located in the suprachiasmatic nuclei (SCN) of the hypothalamus. Vasoactive intestinal polypeptide (VIP) and its receptor, VPAC2, play important roles in the functioning of the SCN pacemaker. Mice lacking VPAC2 receptors (Vipr2−/−) express disrupted behavioral and metabolic rhythms and show altered SCN neuronal activity and clock gene expression. Within the brain, the SCN is not the only site containing endogenous circadian oscillators, nor is it the only site of VPAC2 receptor expression; both VPAC2 receptors and rhythmic clock gene/protein expression have been noted in the arcuate (Arc) and dorsomedial (DMH) nuclei of the mediobasal hypothalamus, and in the pituitary gland. The functional role of VPAC2 receptors in rhythm generation and maintenance in these tissues is, however, unknown. We used wild type (WT) and Vipr2−/− mice expressing a luciferase reporter (PER2::LUC) to investigate whether circadian rhythms in the clock gene protein PER2 in these extra-SCN tissues were compromised by the absence of the VPAC2 receptor. Vipr2−/− SCN cultures expressed significantly lower amplitude PER2::LUC oscillations than WT SCN. Surprisingly, in Vipr2−/− Arc/ME/PT complex (Arc, median eminence and pars tuberalis), DMH and pituitary, the period, amplitude and rate of damping of rhythms were not significantly different to WT. Intriguingly, while we found WT SCN and Arc/ME/PT tissues to maintain a consistent circadian phase when cultured, the phase of corresponding Vipr2−/− cultures was reset by cull/culture procedure. These data demonstrate that while the main rhythm parameters of extra-SCN circadian oscillations are maintained in Vipr2−/− mice, the ability of these oscillators to resist phase shifts is compromised. These deficiencies may contribute towards the aberrant behavior and metabolism associated with Vipr2−/− animals. Further, our data indicate a link between circadian rhythm strength and the ability of tissues to resist circadian phase resetting

The Passive Yet Successful Way of Planktonic Life: Genomic and Experimental Analysis of the Ecology of a Free-Living Polynucleobacter Population

Background: The bacterial taxon Polynucleobacter necessarius subspecies asymbioticus represents a group of planktonic freshwater bacteria with cosmopolitan and ubiquitous distribution in standing freshwater habitats. These bacteria comprise,1 % to 70 % (on average about 20%) of total bacterioplankton cells in various freshwater habitats. The ubiquity of this taxon was recently explained by intra-taxon ecological diversification, i.e. specialization of lineages to specific environmental conditions; however, details on specific adaptations are not known. Here we investigated by means of genomic and experimental analyses the ecological adaptation of a persistent population dwelling in a small acidic pond. Findings: The investigated population (F10 lineage) contributed on average 11 % to total bacterioplankton in the pond during the vegetation periods (ice-free period, usually May to November). Only a low degree of genetic diversification of the population could be revealed. These bacteria are characterized by a small genome size (2.1 Mb), a relatively small number of genes involved in transduction of environmental signals, and the lack of motility and quorum sensing. Experiments indicated that these bacteria live as chemoorganotrophs by mainly utilizing low-molecular-weight substrates derived from photooxidation of humic substances. Conclusions: Evolutionary genome streamlining resulted in a highly passive lifestyle so far only known among free-living bacteria from pelagic marine taxa dwelling in environmentally stable nutrient-poor off-shore systems. Surprisingly, such a lifestyle is also successful in a highly dynamic and nutrient-richer environment such as the water column of the investigate

A Threshold Equation for Action Potential Initiation

In central neurons, the threshold for spike initiation can depend on the stimulus and varies between cells and between recording sites in a given cell, but it is unclear what mechanisms underlie this variability. Properties of ionic channels are likely to play a role in threshold modulation. We examined in models the influence of Na channel activation, inactivation, slow voltage-gated channels and synaptic conductances on spike threshold. We propose a threshold equation which quantifies the contribution of all these mechanisms. It provides an instantaneous time-varying value of the threshold, which applies to neurons with fluctuating inputs. We deduce a differential equation for the threshold, similar to the equations of gating variables in the Hodgkin-Huxley formalism, which describes how the spike threshold varies with the membrane potential, depending on channel properties. We find that spike threshold depends logarithmically on Na channel density, and that Na channel inactivation and K channels can dynamically modulate it in an adaptive way: the threshold increases with membrane potential and after every action potential. Our equation was validated with simulations of a previously published multicompartemental model of spike initiation. Finally, we observed that threshold variability in models depends crucially on the shape of the Na activation function near spike initiation (about −55 mV), while its parameters are adjusted near half-activation voltage (about −30 mV), which might explain why many models exhibit little threshold variability, contrary to experimental observations. We conclude that ionic channels can account for large variations in spike threshold