Brain Microtubule Electrical Oscillations-Empirical Mode Decomposition Analysis.

Microtubules (MTs) are essential cytoskeletal polymers of eukaryote cells implicated in various cell functions, including cell division, cargo transfer, and cell signaling. MTs also are highly charged polymers that generate electrical oscillations that may underlie their ability to act as nonlinear transmission lines. However, the oscillatory composition and time-frequency differences of the MT electrical oscillations have not been identified. Here, we applied the Empirical Mode Decomposition (EMD) to bovine brain MT sheet recordings to determine the number and fundamental frequencies of the Intrinsic Modes Functions (IMF) and evaluate their energetic contribution to the electrical signal. As previously reported, raw signals were obtained from cow brain MTs (Cantero et al. Sci Rep 6:27143, 2016), sampled, filtered, and subjected to signal decomposition from representative experiments. Filtered signals (200 Hz) allowed us to identify either six or seven IMFs. The reconstructed tracings faithfully resembled the original signals, with identifiable frequency peaks. To extend the analysis to obtain time-frequency information and the energy implicated in each IMF, we applied the Hilbert-Huang Transform (HHT) and the Continuous Wavelet Transform (CWT) to the same samples. The analyses disclosed the presence of more fundamental frequency peaks than initially reported and evidenced the advantages and disadvantages of each transform. The study indicates that the EMD is a robust approach to quantifying signal decomposition of brain MT oscillations and suggests novel similarities with human brain wave electroencephalogram (EEG) recordings. The evidence points to the potentially fundamental role of MT oscillations in brain electrical activity.

Electrical recordings from dendritic spines of adult mouse hippocampus and effect of the actin cytoskeleton.

Dendritic spines (DS) are tiny protrusions implicated in excitatory postsynaptic responses in the CNS. To achieve their function, DS concentrate a high density of ion channels and dynamic actin networks in a tiny specialized compartment. However, to date there is no direct information on DS ionic conductances. Here, we used several experimental techniques to obtain direct electrical information from DS of the adult mouse hippocampus. First, we optimized a method to isolate DS from the dissected hippocampus. Second, we used the lipid bilayer membrane (BLM) reconstitution and patch clamping techniques and obtained heretofore unavailable electrical phenotypes on ion channels present in the DS membrane. Third, we also patch clamped DS directly in cultured adult mouse hippocampal neurons, to validate the electrical information observed with the isolated preparation. Electron microscopy and immunochemistry of PDS-95 and NMDA receptors and intrinsic actin networks confirmed the enrichment of the isolated DS preparation, showing open and closed DS, and multi-headed DS. The preparation was used to identify single channel activities and "whole-DS" electrical conductance. We identified NMDA and Ca2+-dependent intrinsic electrical activity in isolated DS and in situ DS of cultured adult mouse hippocampal neurons. In situ recordings in the presence of local NMDA, showed that individual DS intrinsic electrical activity often back-propagated to the dendrite from which it sprouted. The DS electrical oscillations were modulated by changes in actin cytoskeleton dynamics by addition of the F-actin disrupter agent, cytochalasin D, and exogenous actin-binding proteins. The data indicate that DS are elaborate excitable electrical devices, whose activity is a functional interplay between ion channels and the underlying actin networks. The data argue in favor of the active contribution of individual DS to the electrical activity of neurons at the level of both the membrane conductance and cytoskeletal signaling.

A vectorial tree distance measure.

A vectorial distance measure for trees is presented. Given two trees, we define a Tree-Alignment (T-Alignment). We T-align the trees from their centers outwards, starting from the root-branches, to make the next level as similar as possible. The algorithm is recursive; condition on the T-alignment of the root-branches we T-align the sub-branches, thereafter each T-alignment is conditioned on the previous one. We define a minimal T-alignment under a lexicographic order which follows the intuition that the differences between the two trees constitutes a vector. Given such a minimal T-alignment, the difference in the number of branches calculated at any level defines the entry of the distance vector at that level. We compare our algorithm to other well-known tree distance measures in the task of clustering sets of phylogenetic trees. We use the TreeSimGM simulator for generating stochastic phylogenetic trees. The vectorial tree distance (VTD) can successfully separate symmetric from asymmetric trees, and hierarchical from non-hierarchical trees. We also test the algorithm as a classifier of phylogenetic trees extracted from two members of the fungi kingdom, mushrooms and mildews, thus showimg that the algorithm can separate real world phylogenetic trees. The Matlab code can be accessed via: https://gitlab.com/avner.priel/vectorial-tree-distance .

A single nucleotide variant of human PARP1 determines response to PARP inhibitors.

The introduction of novel cancer drugs and innovative treatments brings great hope for cancer patients, but also an urgent need to match drugs to suitable patients, since certain drugs that benefit one patient may actually harm others. The newly developed poly-ADP ribose polymerase (PARP) inhibitors (PARPis) are a group of pharmacological enzyme inhibitors used clinically for multiple indications. Several forms of cancer tend to be PARP dependent, making PARP an attractive target for cancer therapy. Specifically, PARPis are commonly used in BRCA-associated breast cancers patients, since unrepaired single-strand breaks are converted into double-strand breaks and BRCA-associated tumors cannot repair them by homologous recombination so that PARPi leads to tumor cell death, by a mechanism called "Synthetic Lethality". Unfortunately, not all patients respond to PARPi, and it is not currently possible to predict who will or will not respond. Here, we present a specific genomic marker, which reflects a single-nucleotide polymorphism of human PARP1 and correlates in vitro with response to PARPi, throughout all indications. In addition, we report that this SNP is associated with re-shaping mRNA, and mRNA levels, and influences the final protein structure to expose new binding sites while hiding others. The status of the SNP is therefore critical to patients' care, as it relates responses to PARPi to the PARP1-SNP carried.

Network Representation of T-Cell Repertoire- A Novel Tool to Analyze Immune Response to Cancer Formation.

The T cell repertoire potentially presents complexity compatible, or greater than, that of the human brain. T cell based immune response is involved with practically every part of human physiology, and high-throughput biology needed to follow the T-cell repertoire has made great leaps with the advent of massive parallel sequencing [1]. Nevertheless, tools to handle and observe the dynamics of this complexity have only recently started to emerge [e.g., 2, 3, 4] in parallel with sequencing technologies. Here, we present a network-based view of the dynamics of the T cell repertoire, during the course of mammary tumors development in a mouse model. The transition from the T cell receptor as a feature, to network-based clustering, followed by network-based temporal analyses, provides novel insights to the workings of the system and provides novel tools to observe cancer progression via the perspective of the immune system. The crux of the approach here is at the network-motivated clustering. The purpose of the clustering step is not merely data reduction and exposing structures, but rather to detect hubs, or attractors, within the T cell receptor repertoire that might shed light on the behavior of the immune system as a dynamic network. The Clone-Attractor is in fact an extension of the clone concept, i.e., instead of looking at particular clones we observe the extended clonal network by assigning clusters to graph nodes and edges to adjacent clusters (editing distance metric). Viewing the system as dynamical brings to the fore the notion of an attractors landscape, hence the possibility to chart this space and map the sample state at a given time to a vector in this large space. Based on this representation we applied two different methods to demonstrate its effectiveness in identifying changes in the repertoire that correlate with changes in the phenotype: (1) network analysis of the TCR repertoire in which two measures were calculated and demonstrated the ability to differentiate control from transgenic samples, and, (2) machine learning classifier capable of both stratifying control and trangenic samples, as well as to stratify pre-cancer and cancer samples.

Metabolic control of type 1 regulatory T cell differentiation by AHR and HIF1-α.

Our understanding of the pathways that regulate lymphocyte metabolism, as well as the effects of metabolism and its products on the immune response, is still limited. We report that a metabolic program controlled by the transcription factors hypoxia inducible factor-1α (HIF1-α) and aryl hydrocarbon receptor (AHR) supports the differentiation of type 1 regulatory T cell (Tr1) cells. HIF1-α controls the early metabolic reprograming of Tr1 cells. At later time points, AHR promotes HIF1-α degradation and takes control of Tr1 cell metabolism. Extracellular ATP (eATP) and hypoxia, linked to inflammation, trigger AHR inactivation by HIF1-α and inhibit Tr1 cell differentiation. Conversely, CD39 promotes Tr1 cell differentiation by depleting eATP. CD39 also contributes to Tr1 suppressive activity by generating adenosine in cooperation with CD73 expressed by responder T cells and antigen-presenting cells. These results suggest that HIF1-α and AHR integrate immunological, metabolic and environmental signals to regulate the immune response.

Keeping time: could quantum beating in microtubules be the basis for the neural synchrony related to consciousness?

This paper discusses the possibility of quantum coherent oscillations playing a role in neuronal signaling. Consciousness correlates strongly with coherent neural oscillations, however the mechanisms by which neurons synchronize are not fully elucidated. Recent experimental evidence of quantum beats in light-harvesting complexes of plants (LHCII) and bacteria provided a stimulus for seeking similar effects in important structures found in animal cells, especially in neurons. We argue that microtubules (MTs), which play critical roles in all eukaryotic cells, possess structural and functional characteristics that are consistent with quantum coherent excitations in the aromatic groups of their tryptophan residues. Furthermore we outline the consequences of these findings on neuronal processes including the emergence of consciousness.

Model of ionic currents through microtubule nanopores and the lumen.

It has been suggested that microtubules and other cytoskeletal filaments may act as electrical transmission lines. An electrical circuit model of the microtubule is constructed incorporating features of its cylindrical structure with nanopores in its walls. This model is used to study how ionic conductance along the lumen is affected by flux through the nanopores, both with and without an external potential applied across its two ends. Based on the results of Brownian dynamics simulations, the nanopores were found to have asymmetric inner and outer conductances, manifested as nonlinear IV curves. Our simulations indicate that a combination of this asymmetry and an internal voltage source arising from the motion of the C-terminal tails causes cations to be pumped across the microtubule wall and propagate in both directions down the microtubule through the lumen, returning to the bulk solution through its open ends. This effect is demonstrated to add directly to the longitudinal current through the lumen resulting from an external voltage source applied across the two ends of the microtubule. The predicted persistent currents directed through the microtubule wall and along the lumen could be significant in directing the dissipation of weak, endogenous potential gradients toward one end of the microtubule within the cellular environment.

Reversibility and efficiency in coding protein information.

Why the genetic code has a fixed length? Protein information is transferred by coding each amino acid using codons whose length equals 3 for all amino acids. Hence the most probable and the least probable amino acid get a codeword with an equal length. Moreover, the distributions of amino acids found in nature are not uniform and therefore the efficiency of such codes is sub-optimal. The origins of these apparently non-efficient codes are yet unclear. In this paper we propose an a priori argument for the energy efficiency of such codes resulting from their reversibility, in contrast to their time inefficiency. Such codes are reversible in the sense that a primitive processor, reading three letters in each step, can always reverse its operation, undoing its process. We examine the codes for the distributions of amino acids that exist in nature and show that they could not be both time efficient and reversible. We investigate a family of Zipf-type distributions and present their efficient (non-fixed length) prefix code, their graphs, and the condition for their reversibility. We prove that for a large family of such distributions, if the code is time efficient, it could not be reversible. In other words, if pre-biotic processes demand reversibility, the protein code could not be time efficient. The benefits of reversibility are clear: reversible processes are adiabatic, namely, they dissipate a very small amount of energy. Such processes must be done slowly enough; therefore time efficiency is non-important. It is reasonable to assume that early biochemical complexes were more prone towards energy efficiency, where forward and backward processes were almost symmetrical.

Microtubule ionic conduction and its implications for higher cognitive functions.

The neuronal cytoskeleton has been hypothesized to play a role in higher cognitive functions including learning, memory and consciousness. Experimental evidence suggests that both microtubules and actin filaments act as biological electrical wires that can transmit and amplify electric signals via the flow of condensed ion clouds. The potential transmission of electrical signals via the cytoskeleton is of extreme importance to the electrical activity of neurons in general. In this regard, the unique structure, geometry and electrostatics of microtubules are discussed with the expected impact on their specific functions within the neuron. Electric circuit models of ionic flow along microtubules are discussed in the context of experimental data, and the specific importance of both the tubulin C-terminal tail regions, and the nano-pore openings lining the microtubule wall is elucidated. Overall, these recent results suggest that ions, condensed around the surface of the major filaments of the cytoskeleton, flow along and through microtubules in the presence of potential differences, thus acting as transmission lines propagating intracellular signals in a given cell. The significance of this conductance to the functioning of the electrically active neuron, and to higher cognitive function is also discussed.

Neural cytoskeleton capabilities for learning and memory.

This paper proposes a physical model involving the key structures within the neural cytoskeleton as major players in molecular-level processing of information required for learning and memory storage. In particular, actin filaments and microtubules are macromolecules having highly charged surfaces that enable them to conduct electric signals. The biophysical properties of these filaments relevant to the conduction of ionic current include a condensation of counterions on the filament surface and a nonlinear complex physical structure conducive to the generation of modulated waves. Cytoskeletal filaments are often directly connected with both ionotropic and metabotropic types of membrane-embedded receptors, thereby linking synaptic inputs to intracellular functions. Possible roles for cable-like, conductive filaments in neurons include intracellular information processing, regulating developmental plasticity, and mediating transport. The cytoskeletal proteins form a complex network capable of emergent information processing, and they stand to intervene between inputs to and outputs from neurons. In this manner, the cytoskeletal matrix is proposed to work with neuronal membrane and its intrinsic components (e.g., ion channels, scaffolding proteins, and adaptor proteins), especially at sites of synaptic contacts and spines. An information processing model based on cytoskeletal networks is proposed that may underlie certain types of learning and memory.

Effect of calcium on electrical energy transfer by microtubules.

Microtubules (MTs) are important cytoskeletal superstructures implicated in neuronal morphology and function, which are involved in vesicle trafficking, neurite formation and differentiation and other morphological changes. The structural and functional properties of MTs depend on their high intrinsic charge density and functional regulation by the MT depolymerising properties of changes in Ca(2 + ) concentration. Recently, we reported on remarkable properties of isolated MTs, which behave as biomolecular transistors capable of amplifying electrical signals (Priel et al., Biophys J 90:4639-4643, 2006). Here, we demonstrate that MT-bathing (cytoplasmic) Ca(2 + ) concentrations modulate the electrodynamic properties of MTs. Electrical amplification by MTs was exponentially dependent on the Ca(2 + ) concentration between 10( - 7) and 10( - 2) M. However, the electrical connectivity (coupling) of MTs was optimal at a narrower window of Ca(2 + ) concentrations. We observed that while raising bathing Ca(2 + ) concentration increased electrical amplification by MTs, energy transfer was highest in the presence of ethylene glycol tetraacetic acid (lowest Ca(2 + ) concentration). Our data indicate that Ca(2 + ) is an important modulator of electrical amplification by MTs, supporting the hypothesis that this divalent cation, which adsorbs onto the polymer's surface, plays an important role as a regulator of the electrical properties of MTs. The Ca(2 + )-dependent ability of MTs to modulate and amplify electrical signals may provide a novel means of cell signaling, likely contributing to neuronal function.

A biopolymer transistor: electrical amplification by microtubules.

Microtubules (MTs) are important cytoskeletal structures engaged in a number of specific cellular activities, including vesicular traffic, cell cyto-architecture and motility, cell division, and information processing within neuronal processes. MTs have also been implicated in higher neuronal functions, including memory and the emergence of "consciousness". How MTs handle and process electrical information, however, is heretofore unknown. Here we show new electrodynamic properties of MTs. Isolated, taxol-stabilized MTs behave as biomolecular transistors capable of amplifying electrical information. Electrical amplification by MTs can lead to the enhancement of dynamic information, and processivity in neurons can be conceptualized as an "ionic-based" transistor, which may affect, among other known functions, neuronal computational capabilities.

Transitions in microtubule C-termini conformations as a possible dendritic signaling phenomenon.

We model the dynamical states of the C-termini of tubulin dimers that comprise neuronal microtubules. We use molecular dynamics and other computational tools to explore the time-dependent behavior of conformational states of a C-terminus of tubulin within a microtubule and assume that each C-terminus interacts via screened Coulomb forces with the surface of a tubulin dimer, with neighboring C-termini and also with any adjacent microtubule-associated protein 2 (MAP2). Each C-terminus can either bind to the tubulin surface via one of the several positively charged regions or can be allowed to explore the space available in the solution surrounding the dimer. We find that the preferential orientation of each C-terminus is away from the tubulin surface but binding to the surface may also take place, albeit at a lower probability. The results of our model suggest that perturbations generated by the C-termini interactions with counterions surrounding a MAP2 may propagate over distances greater than those between adjacent microtubules. Thus, the MAP2 structure is able to act as a kind of biological wire (or a cable) transmitting local electrostatic perturbations resulting in ionic concentration gradients from one microtubule to another. We briefly discuss the implications the current dynamic modeling may have on synaptic activation and potentiation.