Search results for "Computer simulation"
showing 10 items of 1054 documents
Wavelength selection of rippling patterns in myxobacteria
2016
Rippling patterns of myxobacteria appear in starving colonies before they aggregate to form fruiting bodies. These periodic traveling cell density waves arise from the coordination of individual cell reversals, resulting from an internal clock regulating them, and from contact signaling during bacterial collisions. Here we revisit a mathematical model of rippling in myxobacteria due to Igoshin et al.\ [Proc. Natl. Acad. Sci. USA {\bf 98}, 14913 (2001) and Phys. Rev. E {\bf 70}, 041911 (2004)]. Bacteria in this model are phase oscillators with an extra internal phase through which they are coupled to a mean-field of oppositely moving bacteria. Previously, patterns for this model were obtaine…
The role of spatial structure in the evolution of viral innate immunity evasion: A diffusion-reaction cellular automaton model
2020
Most viruses have evolved strategies for preventing interferon (IFN) secretion and evading innate immunity. Recent work has shown that viral shutdown of IFN secretion can be viewed as a social trait, since the ability of a given virus to evade IFN-mediated immunity depends on the phenotype of neighbor viruses. Following this idea, we investigate the role of spatial structure in the evolution of innate immunity evasion. For this, we model IFN signaling and viral spread using a spatially explicit approximation that combines a diffusion-reaction model and cellular automaton. Our results indicate that the benefits of preventing IFN secretion for a virus are strongly determined by spatial struct…
Epigenetic mutations can both help and hinder adaptive evolution.
2015
Epigenetic variation is being integrated into our understanding of adaptation, yet we lack models on how epigenetic mutations affect evolution that includes de novo genetic change. We model the effects of epigenetic mutations on the dynamics and endpoints of adaptive walks—a process where a series of beneficial mutations move a population towards a fitness optimum. We use an individual-based model of an asexual population, where mutational effects are drawn from Fisher's geometric model. We find cases where epigenetic mutations speed adaptation or result in populations with higher fitness. However, we also find cases where they slow adaptation or result in populations with lower fitness. Th…
A Simple Method to Predict Blood-Brain Barrier Permeability of Drug- Like Compounds Using Classification Trees
2017
Background: To know the ability of a compound to penetrate the blood-brain barrier (BBB) is a challenging task; despite the numerous efforts realized to predict/measure BBB passage, they still have several drawbacks. Methods: The prediction of the permeability through the BBB is carried out using classification trees. A large data set of 497 compounds (recently published) is selected to develop the tree model. Results: The best model shows an accuracy higher than 87.6% for training set; the model was also validated using 10-fold cross-validation procedure and through a test set achieving accuracy values of 86.1% and 87.9%, correspondingly. We give a brief explanation, in structural terms, o…
Assessment of in vivo organ-uptake and in silico prediction of CYP mediated metabolism of DA-Phen, a new dopaminergic agent
2017
Abstract The drug development process strives to predict metabolic fate of a drug candidate, together with its uptake in major organs, whether they act as target, deposit or metabolism sites, to the aim of establish a relationship between the pharmacodynamics and the pharmacokinetics and highlight the potential toxicity of the drug candidate. The present study was aimed at evaluating the in vivo uptake of 2-Amino-N-[2-(3,4-dihydroxy-phenyl)-ethyl]-3-phenyl-propionamide (DA-Phen) − a new dopaminergic neurotransmission modulator, in target and non-target organs of animal subjects and integrating these data with SMARTCyp results, an in silico method that predicts the sites of cytochrome P450-m…
Assessing statistical significance in multivariable genome wide association analysis
2016
Motivation: Although Genome Wide Association Studies (GWAS) genotype a very large number of single nucleotide polymorphisms (SNPs), the data are often analyzed one SNP at a time. The low predictive power of single SNPs, coupled with the high significance threshold needed to correct for multiple testing, greatly decreases the power of GWAS. Results: We propose a procedure in which all the SNPs are analyzed in a multiple generalized linear model, and we show its use for extremely high-dimensional datasets. Our method yields P-values for assessing significance of single SNPs or groups of SNPs while controlling for all other SNPs and the family wise error rate (FWER). Thus, our method tests whe…
Stagewise pseudo-value regression for time-varying effects on the cumulative incidence
2015
In a competing risks setting, the cumulative incidence of an event of interest describes the absolute risk for this event as a function of time. For regression analysis, one can either choose to model all competing events by separate cause-specific hazard models or directly model the association between covariates and the cumulative incidence of one of the events. With a suitable link function, direct regression models allow for a straightforward interpretation of covariate effects on the cumulative incidence. In practice, where data can be right-censored, these regression models are implemented using a pseudo-value approach. For a grid of time points, the possibly unobserved binary event s…
L1-Penalized Censored Gaussian Graphical Model
2018
Graphical lasso is one of the most used estimators for inferring genetic networks. Despite its diffusion, there are several fields in applied research where the limits of detection of modern measurement technologies make the use of this estimator theoretically unfounded, even when the assumption of a multivariate Gaussian distribution is satisfied. Typical examples are data generated by polymerase chain reactions and flow cytometer. The combination of censoring and high-dimensionality make inference of the underlying genetic networks from these data very challenging. In this article, we propose an $\ell_1$-penalized Gaussian graphical model for censored data and derive two EM-like algorithm…
Model selection for factorial Gaussian graphical models with an application to dynamic regulatory networks.
2016
Abstract Factorial Gaussian graphical Models (fGGMs) have recently been proposed for inferring dynamic gene regulatory networks from genomic high-throughput data. In the search for true regulatory relationships amongst the vast space of possible networks, these models allow the imposition of certain restrictions on the dynamic nature of these relationships, such as Markov dependencies of low order – some entries of the precision matrix are a priori zeros – or equal dependency strengths across time lags – some entries of the precision matrix are assumed to be equal. The precision matrix is then estimated by l 1-penalized maximum likelihood, imposing a further constraint on the absolute value…
Simulation-based estimation of branching models for LTR retrotransposons
2017
Abstract Motivation LTR retrotransposons are mobile elements that are able, like retroviruses, to copy and move inside eukaryotic genomes. In the present work, we propose a branching model for studying the propagation of LTR retrotransposons in these genomes. This model allows us to take into account both the positions and the degradation level of LTR retrotransposons copies. In our model, the duplication rate is also allowed to vary with the degradation level. Results Various functions have been implemented in order to simulate their spread and visualization tools are proposed. Based on these simulation tools, we have developed a first method to evaluate the parameters of this propagation …