On propagation in networks, promising models beyond network diffusion to describe degenerative brain diseases and traumatic brain injuries.

Introduction: Connections among neurons form one of the most amazing and effective network in nature. At higher level, also the functional structures of the brain is organized as a network. It is therefore natural to use modern techniques of network analysis to describe the structures of networks in the brain. Many studies have been conducted in this area, showing that the structure of the neuronal network is complex, with a small-world topology, modularity and the presence of hubs. Other studies have been conducted to investigate the dynamical processes occurring in brain networks, analyzing local and large-scale network dynamics. Recently, network diffusion dynamics have been proposed as a model for the progression of brain degenerative diseases and for traumatic brain injuries. Methods: In this paper, the dynamics of network diffusion is re-examined and reaction-diffusion models on networks is introduced in order to better describe the degenerative dynamics in the brain. Results: Numerical simulations of the dynamics of injuries in the brain connectome are presented. Different choices of reaction term and initial condition provide very different phenomenologies, showing how network propagation models are highly flexible. Discussion: The uniqueness of this research lies in the fact that it is the first time that reaction-diffusion dynamics have been applied to the connectome to model the evolution of neurodegenerative diseases or traumatic brain injury. In addition, the generality of these models allows the introduction of non-constant diffusion and different reaction terms with non-constant parameters, allowing a more precise definition of the pathology to be studied.

NIAPU: network-informed adaptive positive-unlabeled learning for disease gene identification.

MOTIVATION: Gene-disease associations are fundamental for understanding disease etiology and developing effective interventions and treatments. Identifying genes not yet associated with a disease due to a lack of studies is a challenging task in which prioritization based on prior knowledge is an important element. The computational search for new candidate disease genes may be eased by positive-unlabeled learning, the machine learning (ML) setting in which only a subset of instances are labeled as positive while the rest of the dataset is unlabeled. In this work, we propose a set of effective network-based features to be used in a novel Markov diffusion-based multi-class labeling strategy for putative disease gene discovery. RESULTS: The performances of the new labeling algorithm and the effectiveness of the proposed features have been tested on 10 different disease datasets using three ML algorithms. The new features have been compared against classical topological and functional/ontological features and a set of network- and biological-derived features already used in gene discovery tasks. The predictive power of the integrated methodology in searching for new disease genes has been found to be competitive against state-of-the-art algorithms. AVAILABILITY AND IMPLEMENTATION: The source code of NIAPU can be accessed at https://github.com/AndMastro/NIAPU. The source data used in this study are available online on the respective websites. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Targeting SARS-CoV-2 nsp13 Helicase and Assessment of Druggability Pockets: Identification of Two Potent Inhibitors by a Multi-Site In Silico Drug Repurposing Approach.

The SARS-CoV-2 non-structural protein 13 (nsp13) helicase is an essential enzyme for viral replication and has been identified as an attractive target for the development of new antiviral drugs. In detail, the helicase catalyzes the unwinding of double-stranded DNA or RNA in a 5' to 3' direction and acts in concert with the replication-transcription complex (nsp7/nsp8/nsp12). In this work, bioinformatics and computational tools allowed us to perform a detailed conservation analysis of the SARS-CoV-2 helicase genome and to further predict the druggable enzyme's binding pockets. Thus, a structure-based virtual screening was used to identify valuable compounds that are capable of recognizing multiple nsp13 pockets. Starting from a database of around 4000 drugs already approved by the Food and Drug Administration (FDA), we chose 14 shared compounds capable of recognizing three out of four sites. Finally, by means of visual inspection analysis and based on their commercial availability, five promising compounds were submitted to in vitro assays. Among them, PF-03715455 was able to block both the unwinding and NTPase activities of nsp13 in a micromolar range.

Network Proximity-Based Drug Repurposing Strategy for Early and Late Stages of Primary Biliary Cholangitis.

Primary biliary cholangitis (PBC) is a chronic, cholestatic, immune-mediated, and progressive liver disorder. Treatment to preventing the disease from advancing into later and irreversible stages is still an unmet clinical need. Accordingly, we set up a drug repurposing framework to find potential therapeutic agents targeting relevant pathways derived from an expanded pool of genes involved in different stages of PBC. Starting with updated human protein-protein interaction data and genes specifically involved in the early and late stages of PBC, a network medicine approach was used to provide a PBC "proximity" or "involvement" gene ranking using network diffusion algorithms and machine learning models. The top genes in the proximity ranking, when combined with the original PBC-related genes, resulted in a final dataset of the genes most involved in PBC disease. Finally, a drug repurposing strategy was implemented by mining and utilizing dedicated drug-gene interaction and druggable genome information knowledge bases (e.g., the DrugBank repository). We identified several potential drug candidates interacting with PBC pathways after performing an over-representation analysis on our initial 1121-seed gene list and the resulting disease-associated (algorithm-obtained) genes. The mechanism and potential therapeutic applications of such drugs were then thoroughly discussed, with a particular emphasis on different stages of PBC disease. We found that interleukin/EGFR/TNF-alpha inhibitors, branched-chain amino acids, geldanamycin, tauroursodeoxycholic acid, genistein, antioestrogens, curcumin, antineovascularisation agents, enzyme/protease inhibitors, and antirheumatic agents are promising drugs targeting distinct stages of PBC. We developed robust and transparent selection mechanisms for prioritizing already approved medicinal products or investigational products for repurposing based on recognized unmet medical needs in PBC, as well as solid preliminary data to achieve this goal.

Robust estimation of time-dependent precision matrix with application to the cryptocurrency market.

Most financial signals show time dependency that, combined with noisy and extreme events, poses serious problems in the parameter estimations of statistical models. Moreover, when addressing asset pricing, portfolio selection, and investment strategies, accurate estimates of the relationship among assets are as necessary as are delicate in a time-dependent context. In this regard, fundamental tools that increasingly attract research interests are precision matrix and graphical models, which are able to obtain insights into the joint evolution of financial quantities. In this paper, we present a robust divergence estimator for a time-varying precision matrix that can manage both the extreme events and time-dependency that affect financial time series. Furthermore, we provide an algorithm to handle parameter estimations that uses the "maximization-minimization" approach. We apply the methodology to synthetic data to test its performances. Then, we consider the cryptocurrency market as a real data application, given its remarkable suitability for the proposed method because of its volatile and unregulated nature.

Evaluation of HIV-1 integrase variability by combining computational and probabilistic approaches.

This study aimed at updating previous data on HIV-1 integrase variability, by using effective bioinformatics methods combining different statistical instruments from simple entropy and mutation rate to more specific approaches such as Hellinger distance. A total of 2133 HIV-1 integrase sequences were analyzed in: i) 1460 samples from drug-naïve [DN] individuals; ii) 386 samples from drug-experienced but INI-naïve [IN] individuals; iii) 287 samples from INI-experienced [IE] individuals. Within the three groups, 76 amino acid positions were highly conserved (≤0.2% variation, Hellinger distance: <0.25%), with 35 fully invariant positions; while, 80 positions were conserved (>0.2% to <1% variation, Hellinger distance: <1%). The H12-H16-C40-C43 and D64-D116-E152 motifs were all well conserved. Some residues were affected by dramatic changes in their mutation distributions, especially between DN and IE samples (Hellinger distance ≥1%). In particular, 15 positions (D6, S24, V31, S39, L74, A91, S119, T122, T124, T125, V126, K160, N222, S230, C280) showed a significant decrease of mutation rate in IN and/or IE samples compared to DN samples. Conversely, 8 positions showed significantly higher mutation rate in samples from treated individuals (IN and/or IE) compared to DN. Some of these positions, such as E92, T97, G140, Y143, Q148 and N155, were already known to be associated with resistance to integrase inhibitors; other positions including S24, M154, V165 and D270 are not yet documented to be associated with resistance. Our study confirms the high conservation of HIV-1 integrase and identified highly invariant positions using robust and innovative methods. The role of novel mutations located in the critical region of HIV-1 integrase deserves further investigation.

The farther the better: Investigating how distance from human self affects the propensity of a peptide to be presented on cell surface by MHC class I molecules, the case of Trypanosoma cruzi.

More than twenty years ago the reverse vaccinology paradigm came to light trying to design new vaccines based on the analysis of genomic information in order to select those pathogen peptides able to trigger an immune response. In this context, focusing on the proteome of Trypanosoma cruzi, we investigated the link between the probabilities for pathogen peptides to be presented on a cell surface and their distance from human self. We found a reasonable but, as far as we know, undiscovered property: the farther the distance between a peptide and the human-self the higher the probability for that peptide to be presented on a cell surface. We also found that the most distant peptides from human self bind, on average, a broader collection of HLAs than expected, implying a potential immunological role in a large portion of individuals. Finally, introducing a novel quantitative indicator for a peptide to measure its potential immunological role, we proposed a pool of peptides that could be potential epitopes and that can be suitable for experimental testing. The software to compute peptide classes according to the distance from human self is free available at http://www.iasi.cnr.it/~dsantoni/nullomers.

Designing a Network Proximity-Based Drug Repurposing Strategy for COVID-19.

The ongoing COVID-19 pandemic still requires fast and effective efforts from all fronts, including epidemiology, clinical practice, molecular medicine, and pharmacology. A comprehensive molecular framework of the disease is needed to better understand its pathological mechanisms, and to design successful treatments able to slow down and stop the impressive pace of the outbreak and harsh clinical symptomatology, possibly via the use of readily available, off-the-shelf drugs. This work engages in providing a wider picture of the human molecular landscape of the SARS-CoV-2 infection via a network medicine approach as the ground for a drug repurposing strategy. Grounding on prior knowledge such as experimentally validated host proteins known to be viral interactors, tissue-specific gene expression data, and using network analysis techniques such as network propagation and connectivity significance, the host molecular reaction network to the viral invasion is explored and exploited to infer and prioritize candidate target genes, and finally to propose drugs to be repurposed for the treatment of COVID-19. Ranks of potential target genes have been obtained for coherent groups of tissues/organs, potential and distinct sites of interaction between the virus and the organism. The normalization and the aggregation of the different scores allowed to define a preliminary, restricted list of genes candidates as pharmacological targets for drug repurposing, with the aim of contrasting different phases of the virus infection and viral replication cycle.

In the search of potential epitopes for Wuhan seafood market pneumonia virus using high order nullomers.

Alarms periodically emerge for viral pneumonia infections due to coronavirus. In all cases, these are zoonoses passing the barrier between species and infect humans. The legitimate concern of the international community is due to the fact that the new identified coronavirus, named SARS-CoV-2 (previously called 2019-nCoV), has a quite high mortality rate, around 2%, and a strong ability to spread, with an estimated reproduction number higher than 2. Even though all countries are doing their utmost to stop the pandemic, the only reliable solution to tackle the infection is the rapid development of a vaccine. For this purpose, the means of bioinformatics, applied in the context of reverse-vaccinology paradigm, can be of fundamental help to select the most promising peptides able to trigger an effective immune response. In this short report, using the concept of nullomer and introducing a distance from human self, we provide a list of peptides that could deserve experimental investigation in the view of a potential vaccine for SARS-CoV-2.

Reaction fronts in persistent random walks with demographic stochasticity.

Standard reaction-diffusion systems are characterized by infinite velocities and no persistence in the movement of individuals, two conditions that are violated when considering living organisms. Here we consider a discrete particle model in which individuals move following a persistent random walk with finite speed and grow with logistic dynamics. We show that, when the number of individuals is very large, the individual-based model is well described by the continuous reactive Cattaneo equation (RCE), but for smaller values of the carrying capacity important finite-population effects arise. The effects of fluctuations on the propagation speed are investigated both considering the RCE with a cutoff in the reaction term and by means of numerical simulations of the individual-based model. Finally, a more general Lévy walk process for the transport of individuals is examined and an expression for the front speed of the resulting traveling wave is proposed.

Investigating transcription factor synergism in humans.

Proteins are the core and the engine of every process in cells thus the study of mechanisms that drive the regulation of protein expression, is essential. Transcription factors play a central role in this extremely complex task and they synergically co-operate in order to provide a fine tuning of protein expressions. In the present study, we designed a mathematically well-founded procedure to investigate the mutual positioning of transcription factors binding sites related to a given couple of transcription factors in order to evaluate the possible association between them. We obtained a list of highly related transcription factors couples, whose binding site occurrences significantly group together for a given set of gene promoters, identifying the biological contexts in which the couples are involved in and the processes they should contribute to regulate.

Nullomers and High Order Nullomers in Genomic Sequences.

A nullomer is an oligomer that does not occur as a subsequence in a given DNA sequence, i.e. it is an absent word of that sequence. The importance of nullomers in several applications, from drug discovery to forensic practice, is now debated in the literature. Here, we investigated the nature of nullomers, whether their absence in genomes has just a statistical explanation or it is a peculiar feature of genomic sequences. We introduced an extension of the notion of nullomer, namely high order nullomers, which are nullomers whose mutated sequences are still nullomers. We studied different aspects of them: comparison with nullomers of random sequences, CpG distribution and mean helical rise. In agreement with previous results we found that the number of nullomers in the human genome is much larger than expected by chance. Nevertheless antithetical results were found when considering a random DNA sequence preserving dinucleotide frequencies. The analysis of CpG frequencies in nullomers and high order nullomers revealed, as expected, a high CpG content but it also highlighted a strong dependence of CpG frequencies on the dinucleotide position, suggesting that nullomers have their own peculiar structure and are not simply sequences whose CpG frequency is biased. Furthermore, phylogenetic trees were built on eleven species based on both the similarities between the dinucleotide frequencies and the number of nullomers two species share, showing that nullomers are fairly conserved among close species. Finally the study of mean helical rise of nullomers sequences revealed significantly high mean rise values, reinforcing the hypothesis that those sequences have some peculiar structural features. The obtained results show that nullomers are the consequence of the peculiar structure of DNA (also including biased CpG frequency and CpGs islands), so that the hypermutability model, also taking into account CpG islands, seems to be not sufficient to explain nullomer phenomenon. Finally, high order nullomers could emphasize those features that already make simple nullomers useful in several applications.

Linear and anomalous front propagation in systems with non-Gaussian diffusion: The importance of tails.

We investigate front propagation in systems with diffusive and subdiffusive behavior. The scaling behavior of moments of the diffusive problem, both in the standard and in the anomalous cases, is not enough to determine the features of the reactive front. In fact, the shape of the bulk of the probability distribution of the transport process, which determines the diffusive properties, is important just for preasymptotic behavior of front propagation, while the precise shape of the tails of the probability distribution determines asymptotic behavior of front propagation.

Natural vs. random protein sequences: Discovering combinatorics properties on amino acid words.

Casual mutations and natural selection have driven the evolution of protein amino acid sequences that we observe at present in nature. The question about which is the dominant force of proteins evolution is still lacking of an unambiguous answer. Casual mutations tend to randomize protein sequences while, in order to have the correct functionality, one expects that selection mechanisms impose rigid constraints on amino acid sequences. Moreover, one also has to consider that the space of all possible amino acid sequences is so astonishingly large that it could be reasonable to have a well tuned amino acid sequence indistinguishable from a random one. In order to study the possibility to discriminate between random and natural amino acid sequences, we introduce different measures of association between pairs of amino acids in a sequence, and apply them to a dataset of 1047 natural protein sequences and 10,470 random sequences, carefully generated in order to preserve the relative length and amino acid distribution of the natural proteins. We analyze the multidimensional measures with machine learning techniques and show that, to a reasonable extent, natural protein sequences can be differentiated from random ones.

Extinction dynamics of a discrete population in an oasis.

Understanding the conditions ensuring the persistence of a population is an issue of primary importance in population biology. The first theoretical approach to the problem dates back to the 1950s with the Kierstead, Slobodkin, and Skellam (KiSS) model, namely a continuous reaction-diffusion equation for a population growing on a patch of finite size L surrounded by a deadly environment with infinite mortality, i.e., an oasis in a desert. The main outcome of the model is that only patches above a critical size allow for population persistence. Here we introduce an individual-based analog of the KiSS model to investigate the effects of discreteness and demographic stochasticity. In particular, we study the average time to extinction both above and below the critical patch size of the continuous model and investigate the quasistationary distribution of the number of individuals for patch sizes above the critical threshold.

Reaction spreading on percolating clusters.

Reaction-diffusion processes in two-dimensional percolating structures are investigated. Two different problems are addressed: reaction spreading on a percolating cluster and front propagation through a percolating channel. For reaction spreading, numerical data and analytical estimates show a power-law behavior of the reaction product as M(t)~t(d(l)), where d(l) is the connectivity dimension. In a percolating channel, a statistically stationary traveling wave develops. The speed and the width of the traveling wave are numerically computed. While the front speed is a low-fluctuating quantity and its behavior can be understood using a simple theoretical argument, the front width is a high-fluctuating quantity showing a power-law behavior as a function of the size of the channel.

Front speed in reactive compressible stirred media.

We investigated a nonlinear advection-diffusion-reaction equation for a passive scalar field. The purpose is to understand how the compressibility can affect the front dynamics and the bulk burning rate. We study two classes of flows: periodic shear flow and cellular flow, analyzing the system by varying the extent of compressibility and the reaction rate. We find that the bulk burning rate v(f) in a shear flow increases with compressibility intensity ε, following the relation Δv(f)~ε(2). Furthermore, the faster the reaction is, the more important the difference is with respect to the laminar case. The effect has been quantitatively measured, and it turns out to be generally small. For the cellular flow, two extreme cases have been investigated, with the whole perturbation situated either in the center of the vortex or in the periphery. The dependence in this case does not show a monotonic scaling with different behavior in the two cases. The enhancing remains modest and is always less than 20%.

Reaction spreading on graphs.

We study reaction-diffusion processes on graphs through an extension of the standard reaction-diffusion equation starting from first principles. We focus on reaction spreading, i.e., on the time evolution of the reaction product M(t). At variance with pure diffusive processes, characterized by the spectral dimension d{s}, the important quantity for reaction spreading is found to be the connectivity dimension d{l}. Numerical data, in agreement with analytical estimates based on the features of n independent random walkers on the graph, show that M(t)∼t{d{l}}. In the case of Erdös-Renyi random graphs, the reaction product is characterized by an exponential growth M(t)e{αt} with α proportional to ln(k), where (k) is the average degree of the graph.

Invasions in heterogeneous habitats in the presence of advection.

We investigate invasions from a biological reservoir to an initially empty, heterogeneous habitat in the presence of advection. The habitat consists of a periodic alternation of favorable and unfavorable patches. In the latter the population dies at fixed rate. In the former it grows either with the logistic or with an Allee effect type dynamics, where the population has to overcome a threshold to grow. We study the conditions for successful invasions and the speed of the invasion process, which is numerically and analytically investigated in several limits. Generically advection enhances the downstream invasion speed but decreases the population size of the invading species, and can even inhibit the invasion process. Remarkably, however, the rate of population increase, which quantifies the invasion efficiency, is maximized by an optimal advection velocity. In models with Allee effect, differently from the logistic case, above a critical unfavorable patch size the population localizes in a favorable patch, being unable to invade the habitat. However, we show that advection, when intense enough, may activate the invasion process.

A model of ischemia-induced neuroblast activation in the adult subventricular zone.

We have developed a rat brain organotypic culture model, in which tissue slices contain cortex-subventricular zone-striatum regions, to model neuroblast activity in response to in vitro ischemia. Neuroblast activation has been described in terms of two main parameters, proliferation and migration from the subventricular zone into the injured cortex. We observed distinct phases of neuroblast activation as is known to occur after in vivo ischemia. Thus, immediately after oxygen/glucose deprivation (6-24 hours), neuroblasts reduce their proliferative and migratory activity, whereas, at longer time points after the insult (2 to 5 days), they start to proliferate and migrate into the damaged cortex. Antagonism of ionotropic receptors for extracellular ATP during and after the insult unmasks an early activation of neuroblasts in the subventricular zone, which responded with a rapid and intense migration of neuroblasts into the damaged cortex (within 24 hours). The process is further enhanced by elevating the production of the chemoattractant SDf-1alpha and may also be boosted by blocking the activation of microglia. This organotypic model which we have developed is an excellent in vitro system to study neurogenesis after ischemia and other neurodegenerative diseases. Its application has revealed a SOS response to oxygen/glucose deprivation, which is inhibited by unfavorable conditions due to the ischemic environment. Finally, experimental quantifications have allowed us to elaborate a mathematical model to describe neuroblast activation and to develop a computer simulation which should have promising applications for the screening of drug candidates for novel therapies of ischemia-related pathologies.