Genetic interactions contribute less than additive effects to quantitative trait variation in yeast

PubMed Central

Bloom, Joshua S.; Kotenko, Iulia; Sadhu, Meru J.; Treusch, Sebastian; Albert, Frank W.; Kruglyak, Leonid

2015-01-01

Genetic mapping studies of quantitative traits typically focus on detecting loci that contribute additively to trait variation. Genetic interactions are often proposed as a contributing factor to trait variation, but the relative contribution of interactions to trait variation is a subject of debate. Here we use a very large cross between two yeast strains to accurately estimate the fraction of phenotypic variance due to pairwise QTL–QTL interactions for 20 quantitative traits. We find that this fraction is 9% on average, substantially less than the contribution of additive QTL (43%). Statistically significant QTL–QTL pairs typically have small individual effect sizes, but collectively explain 40% of the pairwise interaction variance. We show that pairwise interaction variance is largely explained by pairs of loci at least one of which has a significant additive effect. These results refine our understanding of the genetic architecture of quantitative traits and help guide future mapping studies. PMID:26537231

Frequency-Dependent Selection: The High Potential for Permanent Genetic Variation in the Diallelic, Pairwise Interaction Model

PubMed Central

Asmussen, M. A.; Basnayake, E.

1990-01-01

A detailed analytic and numerical study is made of the potential for permanent genetic variation in frequency-dependent models based on pairwise interactions among genotypes at a single diallelic locus. The full equilibrium structure and qualitative gene-frequency dynamics are derived analytically for a symmetric model, in which pairwise fitnesses are chiefly determined by the genetic similarity of the individuals involved. This is supplemented by an extensive numerical investigation of the general model, the symmetric model, and nine other special cases. Together the results show that there is a high potential for permanent genetic diversity in the pairwise interaction model, and provide insight into the extent to which various forms of genotypic interactions enhance or reduce this potential. Technically, although two stable polymorphic equilibria are possible, the increased likelihood of maintaining both alleles, and the poor performance of protected polymorphism conditions as a measure of this likelihood, are primarily due to a greater variety and frequency of equilibrium patterns with one stable polymorphic equilibrium, in conjunction with a disproportionately large domain of attraction for stable internal equilibria. PMID:2341034

Bispectral pairwise interacting source analysis for identifying systems of cross-frequency interacting brain sources from electroencephalographic or magnetoencephalographic signals

NASA Astrophysics Data System (ADS)

Chella, Federico; Pizzella, Vittorio; Zappasodi, Filippo; Nolte, Guido; Marzetti, Laura

2016-05-01

Brain cognitive functions arise through the coordinated activity of several brain regions, which actually form complex dynamical systems operating at multiple frequencies. These systems often consist of interacting subsystems, whose characterization is of importance for a complete understanding of the brain interaction processes. To address this issue, we present a technique, namely the bispectral pairwise interacting source analysis (biPISA), for analyzing systems of cross-frequency interacting brain sources when multichannel electroencephalographic (EEG) or magnetoencephalographic (MEG) data are available. Specifically, the biPISA makes it possible to identify one or many subsystems of cross-frequency interacting sources by decomposing the antisymmetric components of the cross-bispectra between EEG or MEG signals, based on the assumption that interactions are pairwise. Thanks to the properties of the antisymmetric components of the cross-bispectra, biPISA is also robust to spurious interactions arising from mixing artifacts, i.e., volume conduction or field spread, which always affect EEG or MEG functional connectivity estimates. This method is an extension of the pairwise interacting source analysis (PISA), which was originally introduced for investigating interactions at the same frequency, to the study of cross-frequency interactions. The effectiveness of this approach is demonstrated in simulations for up to three interacting source pairs and for real MEG recordings of spontaneous brain activity. Simulations show that the performances of biPISA in estimating the phase difference between the interacting sources are affected by the increasing level of noise rather than by the number of the interacting subsystems. The analysis of real MEG data reveals an interaction between two pairs of sources of central mu and beta rhythms, localizing in the proximity of the left and right central sulci.

Searching for Collective Behavior in a Large Network of Sensory Neurons

PubMed Central

Tkačik, Gašper; Marre, Olivier; Amodei, Dario; Schneidman, Elad; Bialek, William; Berry, Michael J.

2014-01-01

Maximum entropy models are the least structured probability distributions that exactly reproduce a chosen set of statistics measured in an interacting network. Here we use this principle to construct probabilistic models which describe the correlated spiking activity of populations of up to 120 neurons in the salamander retina as it responds to natural movies. Already in groups as small as 10 neurons, interactions between spikes can no longer be regarded as small perturbations in an otherwise independent system; for 40 or more neurons pairwise interactions need to be supplemented by a global interaction that controls the distribution of synchrony in the population. Here we show that such “K-pairwise” models—being systematic extensions of the previously used pairwise Ising models—provide an excellent account of the data. We explore the properties of the neural vocabulary by: 1) estimating its entropy, which constrains the population's capacity to represent visual information; 2) classifying activity patterns into a small set of metastable collective modes; 3) showing that the neural codeword ensembles are extremely inhomogenous; 4) demonstrating that the state of individual neurons is highly predictable from the rest of the population, allowing the capacity for error correction. PMID:24391485

Quantum Spin Dynamics with Pairwise-Tunable, Long-Range Interactions

DTIC Science & Technology

2016-08-05

rection of the arrows. Dashed (dotted) lines mark the NNN hopping terms (coefficients ±t2). NNNN long -range hopping along curved lines are included to...Quantum spin dynamics with pairwise-tunable, long -range interactions C.-L. Hunga,b,1,2, Alejandro González-Tudelac,1,2, J. Ignacio Ciracc, and H. J...atoms) that interact by way of a variety of processes, such as atomic collisions. Such pro- cesses typically lead to short -range, nearest-neighbor

Pairwise Interaction Extended Point-Particle (PIEP) model for multiphase jets and sedimenting particles

NASA Astrophysics Data System (ADS)

Liu, Kai; Balachandar, S.

2017-11-01

We perform a series of Euler-Lagrange direct numerical simulations (DNS) for multiphase jets and sedimenting particles. The forces the flow exerts on the particles in these two-way coupled simulations are computed using the Basset-Bousinesq-Oseen (BBO) equations. These forces do not explicitly account for particle-particle interactions, even though such pairwise interactions induced by the perturbations from neighboring particles may be important especially when the particle volume fraction is high. Such effects have been largely unaddressed in the literature. Here, we implement the Pairwise Interaction Extended Point-Particle (PIEP) model to simulate the effect of neighboring particle pairs. A simple collision model is also applied to avoid unphysical overlapping of solid spherical particles. The simulation results indicate that the PIEP model provides a more elaborative and complicated movement of the dispersed phase (droplets and particles). Office of Naval Research (ONR) Multidisciplinary University Research Initiative (MURI) project N00014-16-1-2617.

Statistical mechanics of letters in words

PubMed Central

Stephens, Greg J.; Bialek, William

2013-01-01

We consider words as a network of interacting letters, and approximate the probability distribution of states taken on by this network. Despite the intuition that the rules of English spelling are highly combinatorial and arbitrary, we find that maximum entropy models consistent with pairwise correlations among letters provide a surprisingly good approximation to the full statistics of words, capturing ~92% of the multi-information in four-letter words and even “discovering” words that were not represented in the data. These maximum entropy models incorporate letter interactions through a set of pairwise potentials and thus define an energy landscape on the space of possible words. Guided by the large letter redundancy we seek a lower-dimensional encoding of the letter distribution and show that distinctions between local minima in the landscape account for ~68% of the four-letter entropy. We suggest that these states provide an effective vocabulary which is matched to the frequency of word use and much smaller than the full lexicon. PMID:20866490

Bioinformatic prediction and in vivo validation of residue-residue interactions in human proteins

NASA Astrophysics Data System (ADS)

Jordan, Daniel; Davis, Erica; Katsanis, Nicholas; Sunyaev, Shamil

2014-03-01

Identifying residue-residue interactions in protein molecules is important for understanding both protein structure and function in the context of evolutionary dynamics and medical genetics. Such interactions can be difficult to predict using existing empirical or physical potentials, especially when residues are far from each other in sequence space. Using a multiple sequence alignment of 46 diverse vertebrate species we explore the space of allowed sequences for orthologous protein families. Amino acid changes that are known to damage protein function allow us to identify specific changes that are likely to have interacting partners. We fit the parameters of the continuous-time Markov process used in the alignment to conclude that these interactions are primarily pairwise, rather than higher order. Candidates for sites under pairwise epistasis are predicted, which can then be tested by experiment. We report the results of an initial round of in vivo experiments in a zebrafish model that verify the presence of multiple pairwise interactions predicted by our model. These experimentally validated interactions are novel, distant in sequence, and are not readily explained by known biochemical or biophysical features.

Self-Diffusion of Drops in a Dilute Sheared Emulsion

NASA Technical Reports Server (NTRS)

Loewenberg, Michael; Hinch, E. J.

1996-01-01

Self-diffusion coefficients that describe cross-flow migration of non-Brownian drops in a dilute sheared emulsion were obtained by trajectory calculations. A boundary integral formulation was used to describe pairwise interactions between deformable drops; interactions between undeformed drops were described with mobility functions for spherical drops. The results indicate that drops have large anisotropic self-diffusivities which depend strongly on the drop viscosity and modestly on the shear-rate. Pairwise interactions between drops in shear-flow do not appreciably promote drop breakup.

Weighted projected networks: mapping hypergraphs to networks.

PubMed

López, Eduardo

2013-05-01

Many natural, technological, and social systems incorporate multiway interactions, yet are characterized and measured on the basis of weighted pairwise interactions. In this article, I propose a family of models in which pairwise interactions originate from multiway interactions, by starting from ensembles of hypergraphs and applying projections that generate ensembles of weighted projected networks. I calculate analytically the statistical properties of weighted projected networks, and suggest ways these could be used beyond theoretical studies. Weighted projected networks typically exhibit weight disorder along links even for very simple generating hypergraph ensembles. Also, as the size of a hypergraph changes, a signature of multiway interaction emerges on the link weights of weighted projected networks that distinguishes them from fundamentally weighted pairwise networks. This signature could be used to search for hidden multiway interactions in weighted network data. I find the percolation threshold and size of the largest component for hypergraphs of arbitrary uniform rank, translate the results into projected networks, and show that the transition is second order. This general approach to network formation has the potential to shed new light on our understanding of weighted networks.

Pairwise Maximum Entropy Models for Studying Large Biological Systems: When They Can Work and When They Can't

PubMed Central

Roudi, Yasser; Nirenberg, Sheila; Latham, Peter E.

2009-01-01

One of the most critical problems we face in the study of biological systems is building accurate statistical descriptions of them. This problem has been particularly challenging because biological systems typically contain large numbers of interacting elements, which precludes the use of standard brute force approaches. Recently, though, several groups have reported that there may be an alternate strategy. The reports show that reliable statistical models can be built without knowledge of all the interactions in a system; instead, pairwise interactions can suffice. These findings, however, are based on the analysis of small subsystems. Here, we ask whether the observations will generalize to systems of realistic size, that is, whether pairwise models will provide reliable descriptions of true biological systems. Our results show that, in most cases, they will not. The reason is that there is a crossover in the predictive power of pairwise models: If the size of the subsystem is below the crossover point, then the results have no predictive power for large systems. If the size is above the crossover point, then the results may have predictive power. This work thus provides a general framework for determining the extent to which pairwise models can be used to predict the behavior of large biological systems. Applied to neural data, the size of most systems studied so far is below the crossover point. PMID:19424487

Template-based protein-protein docking exploiting pairwise interfacial residue restraints.

PubMed

Xue, Li C; Rodrigues, João P G L M; Dobbs, Drena; Honavar, Vasant; Bonvin, Alexandre M J J

2017-05-01

Although many advanced and sophisticated ab initio approaches for modeling protein-protein complexes have been proposed in past decades, template-based modeling (TBM) remains the most accurate and widely used approach, given a reliable template is available. However, there are many different ways to exploit template information in the modeling process. Here, we systematically evaluate and benchmark a TBM method that uses conserved interfacial residue pairs as docking distance restraints [referred to as alpha carbon-alpha carbon (CA-CA)-guided docking]. We compare it with two other template-based protein-protein modeling approaches, including a conserved non-pairwise interfacial residue restrained docking approach [referred to as the ambiguous interaction restraint (AIR)-guided docking] and a simple superposition-based modeling approach. Our results show that, for most cases, the CA-CA-guided docking method outperforms both superposition with refinement and the AIR-guided docking method. We emphasize the superiority of the CA-CA-guided docking on cases with medium to large conformational changes, and interactions mediated through loops, tails or disordered regions. Our results also underscore the importance of a proper refinement of superimposition models to reduce steric clashes. In summary, we provide a benchmarked TBM protocol that uses conserved pairwise interface distance as restraints in generating realistic 3D protein-protein interaction models, when reliable templates are available. The described CA-CA-guided docking protocol is based on the HADDOCK platform, which allows users to incorporate additional prior knowledge of the target system to further improve the quality of the resulting models. © The Author 2016. Published by Oxford University Press.

Enhanced identification of synergistic and antagonistic emergent interactions among three or more drugs

PubMed Central

White, Cynthia; Mao, Zhiyuan; Savage, Van M.

2016-01-01

Interactions among drugs play a critical role in the killing efficacy of multi-drug treatments. Recent advances in theory and experiment for three-drug interactions enable the search for emergent interactions—ones not predictable from pairwise interactions. Previous work has shown it is easier to detect synergies and antagonisms among pairwise interactions when a rescaling method is applied to the interaction metric. However, no study has carefully examined whether new types of normalization might be needed for emergence. Here, we propose several rescaling methods for enhancing the classification of the higher order drug interactions based on our conceptual framework. To choose the rescaling that best separates synergism, antagonism and additivity, we conducted bacterial growth experiments in the presence of single, pairwise and triple-drug combinations among 14 antibiotics. We found one of our rescaling methods is far better at distinguishing synergistic and antagonistic emergent interactions than any of the other methods. Using our new method, we find around 50% of emergent interactions are additive, much less than previous reports of greater than 90% additivity. We conclude that higher order emergent interactions are much more common than previously believed, and we argue these findings for drugs suggest that appropriate rescaling is crucial to infer higher order interactions. PMID:27278366

Intransitivity is infrequent and fails to promote annual plant coexistence without pairwise niche differences.

PubMed

Godoy, Oscar; Stouffer, Daniel B; Kraft, Nathan J B; Levine, Jonathan M

2017-05-01

Intransitive competition is often projected to be a widespread mechanism of species coexistence in ecological communities. However, it is unknown how much of the coexistence we observe in nature results from this mechanism when species interactions are also stabilized by pairwise niche differences. We combined field-parameterized models of competition among 18 annual plant species with tools from network theory to quantify the prevalence of intransitive competitive relationships. We then analyzed the predicted outcome of competitive interactions with and without pairwise niche differences. Intransitive competition was found for just 15-19% of the 816 possible triplets, and this mechanism was never sufficient to stabilize the coexistence of the triplet when the pair-wise niche differences between competitors were removed. Of the transitive and intransitive triplets, only four were predicted to coexist and these were more similar in multidimensional trait space defined by 11 functional traits than non-coexisting triplets. Our results argue that intransitive competition may be less frequent than recently posed, and that even when it does operate, pairwise niche differences may be key to possible coexistence. © 2017 by the Ecological Society of America.

The effective colloid interaction in the Asakura-Oosawa model. Assessment of non-pairwise terms from the virial expansion.

PubMed

Santos, Andrés; López de Haro, Mariano; Fiumara, Giacomo; Saija, Franz

2015-06-14

The relevance of neglecting three- and four-body interactions in the coarse-grained version of the Asakura-Oosawa model is examined. A mapping between the first few virial coefficients of the binary nonadditive hard-sphere mixture representative of this model and those arising from the coarse-grained (pairwise) depletion potential approximation allows for a quantitative evaluation of the effect of such interactions. This turns out to be especially important for large size ratios and large reservoir polymer packing fractions.

Evolutionary relationships can be more important than abiotic conditions in predicting the outcome of plant-plant interactions

PubMed Central

Soliveres, Santiago; Torices, Rubén; Maestre, Fernando T.

2015-01-01

Positive and negative plant-plant interactions are major processes shaping plant communities. They are affected by environmental conditions and evolutionary relationships among the interacting plants. However, the generality of these factors as drivers of pairwise plant interactions and their combined effects remain virtually unknown. We conducted an observational study to assess how environmental conditions (altitude, temperature, irradiance and rainfall), the dispersal mechanism of beneficiary species and evolutionary relationships affected the co-occurrence of pairwise interactions in 11 Stipa tenacissima steppes located along an environmental gradient in Spain. We studied 197 pairwise plant-plant interactions involving the two major nurse plants (the resprouting shrub Quercus coccifera and the tussock grass S. tenacissima) found in these communities. The relative importance of the studied factors varied with the nurse species considered. None of the factors studied were good predictors of the co-ocurrence between S. tenacissima and its neighbours. However, both the dispersal mechanism of the beneficiary species and the phylogenetic distance between interacting species were crucial factors affecting the co-occurrence between Q. coccifera and its neighbours, while climatic conditions (irradiance) played a secondary role. Values of phylogenetic distance between 207-272.8 Myr led to competition, while values outside this range or fleshy-fruitness in the beneficiary species led to positive interactions. The low importance of environmental conditions as a general driver of pairwise interactions was caused by the species-specific response to changes in either rainfall or radiation. This result suggests that factors other than climatic conditions must be included in theoretical models aimed to generally predict the outcome of plant-plant interactions. Our study helps to improve current theory on plant-plant interactions and to understand how these interactions can respond to expected modifications in species composition and climate associated to ongoing global environmental change. PMID:25914426

COMPUTING THERAPY FOR PRECISION MEDICINE: COLLABORATIVE FILTERING INTEGRATES AND PREDICTS MULTI-ENTITY INTERACTIONS.

PubMed

Regenbogen, Sam; Wilkins, Angela D; Lichtarge, Olivier

2016-01-01

Biomedicine produces copious information it cannot fully exploit. Specifically, there is considerable need to integrate knowledge from disparate studies to discover connections across domains. Here, we used a Collaborative Filtering approach, inspired by online recommendation algorithms, in which non-negative matrix factorization (NMF) predicts interactions among chemicals, genes, and diseases only from pairwise information about their interactions. Our approach, applied to matrices derived from the Comparative Toxicogenomics Database, successfully recovered Chemical-Disease, Chemical-Gene, and Disease-Gene networks in 10-fold cross-validation experiments. Additionally, we could predict each of these interaction matrices from the other two. Integrating all three CTD interaction matrices with NMF led to good predictions of STRING, an independent, external network of protein-protein interactions. Finally, this approach could integrate the CTD and STRING interaction data to improve Chemical-Gene cross-validation performance significantly, and, in a time-stamped study, it predicted information added to CTD after a given date, using only data prior to that date. We conclude that collaborative filtering can integrate information across multiple types of biological entities, and that as a first step towards precision medicine it can compute drug repurposing hypotheses.

COMPUTING THERAPY FOR PRECISION MEDICINE: COLLABORATIVE FILTERING INTEGRATES AND PREDICTS MULTI-ENTITY INTERACTIONS

PubMed Central

REGENBOGEN, SAM; WILKINS, ANGELA D.; LICHTARGE, OLIVIER

2015-01-01

Biomedicine produces copious information it cannot fully exploit. Specifically, there is considerable need to integrate knowledge from disparate studies to discover connections across domains. Here, we used a Collaborative Filtering approach, inspired by online recommendation algorithms, in which non-negative matrix factorization (NMF) predicts interactions among chemicals, genes, and diseases only from pairwise information about their interactions. Our approach, applied to matrices derived from the Comparative Toxicogenomics Database, successfully recovered Chemical-Disease, Chemical-Gene, and Disease-Gene networks in 10-fold cross-validation experiments. Additionally, we could predict each of these interaction matrices from the other two. Integrating all three CTD interaction matrices with NMF led to good predictions of STRING, an independent, external network of protein-protein interactions. Finally, this approach could integrate the CTD and STRING interaction data to improve Chemical-Gene cross-validation performance significantly, and, in a time-stamped study, it predicted information added to CTD after a given date, using only data prior to that date. We conclude that collaborative filtering can integrate information across multiple types of biological entities, and that as a first step towards precision medicine it can compute drug repurposing hypotheses. PMID:26776170

Whole Protein Native Fitness Potentials

NASA Astrophysics Data System (ADS)

Faraggi, Eshel; Kloczkowski, Andrzej

2013-03-01

Protein structure prediction can be separated into two tasks: sample the configuration space of the protein chain, and assign a fitness between these hypothetical models and the native structure of the protein. One of the more promising developments in this area is that of knowledge based energy functions. However, standard approaches using pair-wise interactions have shown shortcomings demonstrated by the superiority of multi-body-potentials. These shortcomings are due to residue pair-wise interaction being dependent on other residues along the chain. We developed a method that uses whole protein information filtered through machine learners to score protein models based on their likeness to native structures. For all models we calculated parameters associated with the distance to the solvent and with distances between residues. These parameters, in addition to energy estimates obtained by using a four-body-potential, DFIRE, and RWPlus were used as training for machine learners to predict the fitness of the models. Testing on CASP 9 targets showed that our method is superior to DFIRE, RWPlus, and the four-body potential, which are considered standards in the field.

Diversity Increases Indirect Interactions, Attenuates the Intensity of Competition, and Promotes Coexistence.

PubMed

Aschehoug, Erik T; Callaway, Ragan M

2015-10-01

A fundamental assumption of coexistence theory is that competition inevitably decreases species diversity. Consequently, in the quest to understand the ecological regulators of diversity, there has been a great deal of focus on processes with the potential to reduce competitive exclusion. However, the notion that competition must decrease diversity is largely based on the outcome of two-species interaction experiments and models, despite the fact that species rarely interact only in pairs in natural systems. In a field experiment, we found that competition among native perennial plants in multispecies assemblages was far weaker than competition between those same species in pairwise arrangements and that indirect interactions appeared to weaken direct competitive effects. These results suggest that community assembly theory based on pairwise approaches may overestimate the strength of competition and likelihood of competitive exclusion in species-rich communities. We also found that Centaurea stoebe, a North American invader, retained strong competitive effects when competing against North American natives in both pairwise and multispecies assemblages. Our experimental results support an emerging body of theory suggesting that complex networks of competing species may generate strong indirect interactions that can maintain diversity and that ecological differentiation may not be necessary to attenuate competition.

OxfordGrid: a web interface for pairwise comparative map views.

PubMed

Yang, Hongyu; Gingle, Alan R

2005-12-01

OxfordGrid is a web application and database schema for storing and interactively displaying genetic map data in a comparative, dot-plot, fashion. Its display is composed of a matrix of cells, each representing a pairwise comparison of mapped probe data for two linkage groups or chromosomes. These are arranged along the axes with one forming grid columns and the other grid rows with the degree and pattern of synteny/colinearity between the two linkage groups manifested in the cell's dot density and structure. A mouse click over the selected grid cell launches an image map-based display for the selected cell. Both individual and linear groups of mapped probes can be selected and displayed. Also, configurable links can be used to access other web resources for mapped probe information. OxfordGrid is implemented in C#/ASP.NET and the package, including MySQL schema creation scripts, is available at ftp://cggc.agtec.uga.edu/OxfordGrid/.

Everyday bat vocalizations contain information about emitter, addressee, context, and behavior

PubMed Central

Prat, Yosef; Taub, Mor; Yovel, Yossi

2016-01-01

Animal vocal communication is often diverse and structured. Yet, the information concealed in animal vocalizations remains elusive. Several studies have shown that animal calls convey information about their emitter and the context. Often, these studies focus on specific types of calls, as it is rarely possible to probe an entire vocal repertoire at once. In this study, we continuously monitored Egyptian fruit bats for months, recording audio and video around-the-clock. We analyzed almost 15,000 vocalizations, which accompanied the everyday interactions of the bats, and were all directed toward specific individuals, rather than broadcast. We found that bat vocalizations carry ample information about the identity of the emitter, the context of the call, the behavioral response to the call, and even the call’s addressee. Our results underline the importance of studying the mundane, pairwise, directed, vocal interactions of animals. PMID:28005079

Metabolic network prediction through pairwise rational kernels.

PubMed

Roche-Lima, Abiel; Domaratzki, Michael; Fristensky, Brian

2014-09-26

Metabolic networks are represented by the set of metabolic pathways. Metabolic pathways are a series of biochemical reactions, in which the product (output) from one reaction serves as the substrate (input) to another reaction. Many pathways remain incompletely characterized. One of the major challenges of computational biology is to obtain better models of metabolic pathways. Existing models are dependent on the annotation of the genes. This propagates error accumulation when the pathways are predicted by incorrectly annotated genes. Pairwise classification methods are supervised learning methods used to classify new pair of entities. Some of these classification methods, e.g., Pairwise Support Vector Machines (SVMs), use pairwise kernels. Pairwise kernels describe similarity measures between two pairs of entities. Using pairwise kernels to handle sequence data requires long processing times and large storage. Rational kernels are kernels based on weighted finite-state transducers that represent similarity measures between sequences or automata. They have been effectively used in problems that handle large amount of sequence information such as protein essentiality, natural language processing and machine translations. We create a new family of pairwise kernels using weighted finite-state transducers (called Pairwise Rational Kernel (PRK)) to predict metabolic pathways from a variety of biological data. PRKs take advantage of the simpler representations and faster algorithms of transducers. Because raw sequence data can be used, the predictor model avoids the errors introduced by incorrect gene annotations. We then developed several experiments with PRKs and Pairwise SVM to validate our methods using the metabolic network of Saccharomyces cerevisiae. As a result, when PRKs are used, our method executes faster in comparison with other pairwise kernels. Also, when we use PRKs combined with other simple kernels that include evolutionary information, the accuracy values have been improved, while maintaining lower construction and execution times. The power of using kernels is that almost any sort of data can be represented using kernels. Therefore, completely disparate types of data can be combined to add power to kernel-based machine learning methods. When we compared our proposal using PRKs with other similar kernel, the execution times were decreased, with no compromise of accuracy. We also proved that by combining PRKs with other kernels that include evolutionary information, the accuracy can also also be improved. As our proposal can use any type of sequence data, genes do not need to be properly annotated, avoiding accumulation errors because of incorrect previous annotations.

Acoustically mediated long-range interaction among multiple spherical particles exposed to a plane standing wave

NASA Astrophysics Data System (ADS)

Zhang, Shenwei; Qiu, Chunyin; Wang, Mudi; Ke, Manzhu; Liu, Zhengyou

2016-11-01

In this work, we study the acoustically mediated interaction forces among multiple well-separated spherical particles trapped in the same node or antinode plane of a standing wave. An analytical expression of the acoustic interaction force is derived, which is accurate even for the particles beyond the Rayleigh limit. Interestingly, the multi-particle system can be decomposed into a series of independent two-particle systems described by pairwise interactions. Each pairwise interaction is a long-range interaction, as characterized by a soft oscillatory attenuation (at the power exponent of n = -1 or -2). The vector additivity of the acoustic interaction force, which is not well expected considering the nonlinear nature of the acoustic radiation force, is greatly useful for exploring a system consisting of a large number of particles. The capability of self-organizing a big particle cluster can be anticipated through such acoustically controllable long-range interaction.

Time-Frequency Analysis Reveals Pairwise Interactions in Insect Swarms

NASA Astrophysics Data System (ADS)

Puckett, James G.; Ni, Rui; Ouellette, Nicholas T.

2015-06-01

The macroscopic emergent behavior of social animal groups is a classic example of dynamical self-organization, and is thought to arise from the local interactions between individuals. Determining these interactions from empirical data sets of real animal groups, however, is challenging. Using multicamera imaging and tracking, we studied the motion of individual flying midges in laboratory mating swarms. By performing a time-frequency analysis of the midge trajectories, we show that the midge behavior can be segmented into two distinct modes: one that is independent and composed of low-frequency maneuvers, and one that consists of higher-frequency nearly harmonic oscillations conducted in synchrony with another midge. We characterize these pairwise interactions, and make a hypothesis as to their biological function.

Hi-C 2.0: An optimized Hi-C procedure for high-resolution genome-wide mapping of chromosome conformation.

PubMed

Belaghzal, Houda; Dekker, Job; Gibcus, Johan H

2017-07-01

Chromosome conformation capture-based methods such as Hi-C have become mainstream techniques for the study of the 3D organization of genomes. These methods convert chromatin interactions reflecting topological chromatin structures into digital information (counts of pair-wise interactions). Here, we describe an updated protocol for Hi-C (Hi-C 2.0) that integrates recent improvements into a single protocol for efficient and high-resolution capture of chromatin interactions. This protocol combines chromatin digestion and frequently cutting enzymes to obtain kilobase (kb) resolution. It also includes steps to reduce random ligation and the generation of uninformative molecules, such as unligated ends, to improve the amount of valid intra-chromosomal read pairs. This protocol allows for obtaining information on conformational structures such as compartment and topologically associating domains, as well as high-resolution conformational features such as DNA loops. Copyright © 2017 Elsevier Inc. All rights reserved.

Pairwise additivity of energy components in protein-ligand binding: The HIV II protease-Indinavir case

NASA Astrophysics Data System (ADS)

Ucisik, Melek N.; Dashti, Danial S.; Faver, John C.; Merz, Kenneth M.

2011-08-01

An energy expansion (binding energy decomposition into n-body interaction terms for n ≥ 2) to express the receptor-ligand binding energy for the fragmented HIV II protease-Indinavir system is described to address the role of cooperativity in ligand binding. The outcome of this energy expansion is compared to the total receptor-ligand binding energy at the Hartree-Fock, density functional theory, and semiempirical levels of theory. We find that the sum of the pairwise interaction energies approximates the total binding energy to ˜82% for HF and to >95% for both the M06-L density functional and PM6-DH2 semiempirical method. The contribution of the three-body interactions amounts to 18.7%, 3.8%, and 1.4% for HF, M06-L, and PM6-DH2, respectively. We find that the expansion can be safely truncated after n = 3. That is, the contribution of the interactions involving more than three parties to the total binding energy of Indinavir to the HIV II protease receptor is negligible. Overall, we find that the two-body terms represent a good approximation to the total binding energy of the system, which points to pairwise additivity in the present case. This basic principle of pairwise additivity is utilized in fragment-based drug design approaches and our results support its continued use. The present results can also aid in the validation of non-bonded terms contained within common force fields and in the correction of systematic errors in physics-based score functions.

Rigid-Docking Approaches to Explore Protein-Protein Interaction Space.

PubMed

Matsuzaki, Yuri; Uchikoga, Nobuyuki; Ohue, Masahito; Akiyama, Yutaka

Protein-protein interactions play core roles in living cells, especially in the regulatory systems. As information on proteins has rapidly accumulated on publicly available databases, much effort has been made to obtain a better picture of protein-protein interaction networks using protein tertiary structure data. Predicting relevant interacting partners from their tertiary structure is a challenging task and computer science methods have the potential to assist with this. Protein-protein rigid docking has been utilized by several projects, docking-based approaches having the advantages that they can suggest binding poses of predicted binding partners which would help in understanding the interaction mechanisms and that comparing docking results of both non-binders and binders can lead to understanding the specificity of protein-protein interactions from structural viewpoints. In this review we focus on explaining current computational prediction methods to predict pairwise direct protein-protein interactions that form protein complexes.

cDF Theory Software for mesoscopic modeling of equilibrium and transport phenomena

DOE Office of Scientific and Technical Information (OSTI.GOV)

2015-12-01

The approach is based on classical Density Functional Theory ((cDFT) coupled with the Poisson-Nernst-Planck (PNP) transport kinetics model and quantum mechanical description of short-range interaction and elementary transport processes. The model we proposed and implemented is fully atomistic, taking into account pairwise short-range and manybody long-range interactions. But in contrast to standard molecular dynamics (MD) simulations, where long-range manybody interactions are evaluated as a sum of pair-wise atom-atom contributions, we include them analytically based on wellestablished theories of electrostatic and excluded volume interactions in multicomponent systems. This feature of the PNP/cDFT approach allows us to reach well beyond the length-scalesmore » accessible to MD simulations, while retaining the essential physics of interatomic interactions from first principles and in a parameter-free fashion.« less

Interspecific Tests of Allelism Reveal the Evolutionary Timing and Pattern of Accumulation of Reproductive Isolation Mutations

PubMed Central

Sherman, Natasha A.; Victorine, Anna; Wang, Richard J.; Moyle, Leonie C.

2014-01-01

Despite extensive theory, little is known about the empirical accumulation and evolutionary timing of mutations that contribute to speciation. Here we combined QTL (Quantitative Trait Loci) analyses of reproductive isolation, with information on species evolutionary relationships, to reconstruct the order and timing of mutations contributing to reproductive isolation between three plant (Solanum) species. To evaluate whether reproductive isolation QTL that appear to coincide in more than one species pair are homologous, we used cross-specific tests of allelism and found evidence for both homologous and lineage-specific (non-homologous) alleles at these co-localized loci. These data, along with isolation QTL unique to single species pairs, indicate that >85% of isolation-causing mutations arose later in the history of divergence between species. Phylogenetically explicit analyses of these data support non-linear models of accumulation of hybrid incompatibility, although the specific best-fit model differs between seed (pairwise interactions) and pollen (multi-locus interactions) sterility traits. Our findings corroborate theory that predicts an acceleration (‘snowballing’) in the accumulation of isolation loci as lineages progressively diverge, and suggest different underlying genetic bases for pollen versus seed sterility. Pollen sterility in particular appears to be due to complex genetic interactions, and we show this is consistent with a snowball model where later arising mutations are more likely to be involved in pairwise or multi-locus interactions that specifically involve ancestral alleles, compared to earlier arising mutations. PMID:25211473

Comorbidities in the diseasome are more apparent than real: What Bayesian filtering reveals about the comorbidities of depression

PubMed Central

Bolgar, Bence; Deakin, Bill

2017-01-01

Comorbidity patterns have become a major source of information to explore shared mechanisms of pathogenesis between disorders. In hypothesis-free exploration of comorbid conditions, disease-disease networks are usually identified by pairwise methods. However, interpretation of the results is hindered by several confounders. In particular a very large number of pairwise associations can arise indirectly through other comorbidity associations and they increase exponentially with the increasing breadth of the investigated diseases. To investigate and filter this effect, we computed and compared pairwise approaches with a systems-based method, which constructs a sparse Bayesian direct multimorbidity map (BDMM) by systematically eliminating disease-mediated comorbidity relations. Additionally, focusing on depression-related parts of the BDMM, we evaluated correspondence with results from logistic regression, text-mining and molecular-level measures for comorbidities such as genetic overlap and the interactome-based association score. We used a subset of the UK Biobank Resource, a cross-sectional dataset including 247 diseases and 117,392 participants who filled out a detailed questionnaire about mental health. The sparse comorbidity map confirmed that depressed patients frequently suffer from both psychiatric and somatic comorbid disorders. Notably, anxiety and obesity show strong and direct relationships with depression. The BDMM identified further directly co-morbid somatic disorders, e.g. irritable bowel syndrome, fibromyalgia, or migraine. Using the subnetwork of depression and metabolic disorders for functional analysis, the interactome-based system-level score showed the best agreement with the sparse disease network. This indicates that these epidemiologically strong disease-disease relations have improved correspondence with expected molecular-level mechanisms. The substantially fewer number of comorbidity relations in the BDMM compared to pairwise methods implies that biologically meaningful comorbid relations may be less frequent than earlier pairwise methods suggested. The computed interactive comprehensive multimorbidity views over the diseasome are available on the web at Co=MorNet: bioinformatics.mit.bme.hu/UKBNetworks. PMID:28644851

Hierarchical group dynamics in pigeon flocks.

PubMed

Nagy, Máté; Akos, Zsuzsa; Biro, Dora; Vicsek, Tamás

2010-04-08

Animals that travel together in groups display a variety of fascinating motion patterns thought to be the result of delicate local interactions among group members. Although the most informative way of investigating and interpreting collective movement phenomena would be afforded by the collection of high-resolution spatiotemporal data from moving individuals, such data are scarce and are virtually non-existent for long-distance group motion within a natural setting because of the associated technological difficulties. Here we present results of experiments in which track logs of homing pigeons flying in flocks of up to 10 individuals have been obtained by high-resolution lightweight GPS devices and analysed using a variety of correlation functions inspired by approaches common in statistical physics. We find a well-defined hierarchy among flock members from data concerning leading roles in pairwise interactions, defined on the basis of characteristic delay times between birds' directional choices. The average spatial position of a pigeon within the flock strongly correlates with its place in the hierarchy, and birds respond more quickly to conspecifics perceived primarily through the left eye-both results revealing differential roles for birds that assume different positions with respect to flock-mates. From an evolutionary perspective, our results suggest that hierarchical organization of group flight may be more efficient than an egalitarian one, at least for those flock sizes that permit regular pairwise interactions among group members, during which leader-follower relationships are consistently manifested.

A configuration space of homologous proteins conserving mutual information and allowing a phylogeny inference based on pair-wise Z-score probabilities.

PubMed

Bastien, Olivier; Ortet, Philippe; Roy, Sylvaine; Maréchal, Eric

2005-03-10

Popular methods to reconstruct molecular phylogenies are based on multiple sequence alignments, in which addition or removal of data may change the resulting tree topology. We have sought a representation of homologous proteins that would conserve the information of pair-wise sequence alignments, respect probabilistic properties of Z-scores (Monte Carlo methods applied to pair-wise comparisons) and be the basis for a novel method of consistent and stable phylogenetic reconstruction. We have built up a spatial representation of protein sequences using concepts from particle physics (configuration space) and respecting a frame of constraints deduced from pair-wise alignment score properties in information theory. The obtained configuration space of homologous proteins (CSHP) allows the representation of real and shuffled sequences, and thereupon an expression of the TULIP theorem for Z-score probabilities. Based on the CSHP, we propose a phylogeny reconstruction using Z-scores. Deduced trees, called TULIP trees, are consistent with multiple-alignment based trees. Furthermore, the TULIP tree reconstruction method provides a solution for some previously reported incongruent results, such as the apicomplexan enolase phylogeny. The CSHP is a unified model that conserves mutual information between proteins in the way physical models conserve energy. Applications include the reconstruction of evolutionary consistent and robust trees, the topology of which is based on a spatial representation that is not reordered after addition or removal of sequences. The CSHP and its assigned phylogenetic topology, provide a powerful and easily updated representation for massive pair-wise genome comparisons based on Z-score computations.

Collective translational and rotational Monte Carlo cluster move for general pairwise interaction

NASA Astrophysics Data System (ADS)

Růžička, Štěpán; Allen, Michael P.

2014-09-01

Virtual move Monte Carlo is a cluster algorithm which was originally developed for strongly attractive colloidal, molecular, or atomistic systems in order to both approximate the collective dynamics and avoid sampling of unphysical kinetic traps. In this paper, we present the algorithm in the form, which selects the moving cluster through a wider class of virtual states and which is applicable to general pairwise interactions, including hard-core repulsion. The newly proposed way of selecting the cluster increases the acceptance probability by up to several orders of magnitude, especially for rotational moves. The results have their applications in simulations of systems interacting via anisotropic potentials both to enhance the sampling of the phase space and to approximate the dynamics.

Particle-based simulations of self-motile suspensions

NASA Astrophysics Data System (ADS)

Hinz, Denis F.; Panchenko, Alexander; Kim, Tae-Yeon; Fried, Eliot

2015-11-01

A simple model for simulating flows of active suspensions is investigated. The approach is based on dissipative particle dynamics. While the model is potentially applicable to a wide range of self-propelled particle systems, the specific class of self-motile bacterial suspensions is considered as a modeling scenario. To mimic the rod-like geometry of a bacterium, two dissipative particle dynamics particles are connected by a stiff harmonic spring to form an aggregate dissipative particle dynamics molecule. Bacterial motility is modeled through a constant self-propulsion force applied along the axis of each such aggregate molecule. The model accounts for hydrodynamic interactions between self-propelled agents through the pairwise dissipative interactions conventional to dissipative particle dynamics. Numerical simulations are performed using a customized version of the open-source software package LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) software package. Detailed studies of the influence of agent concentration, pairwise dissipative interactions, and Stokes friction on the statistics of the system are provided. The simulations are used to explore the influence of hydrodynamic interactions in active suspensions. For high agent concentrations in combination with dominating pairwise dissipative forces, strongly correlated motion patterns and a fluid-like spectral distributions of kinetic energy are found. In contrast, systems dominated by Stokes friction exhibit weaker spatial correlations of the velocity field. These results indicate that hydrodynamic interactions may play an important role in the formation of spatially extended structures in active suspensions.

Relating Diseases by Integrating Gene Associations and Information Flow through Protein Interaction Network

PubMed Central

Hamaneh, Mehdi Bagheri; Yu, Yi-Kuo

2014-01-01

Identifying similar diseases could potentially provide deeper understanding of their underlying causes, and may even hint at possible treatments. For this purpose, it is necessary to have a similarity measure that reflects the underpinning molecular interactions and biological pathways. We have thus devised a network-based measure that can partially fulfill this goal. Our method assigns weights to all proteins (and consequently their encoding genes) by using information flow from a disease to the protein interaction network and back. Similarity between two diseases is then defined as the cosine of the angle between their corresponding weight vectors. The proposed method also provides a way to suggest disease-pathway associations by using the weights assigned to the genes to perform enrichment analysis for each disease. By calculating pairwise similarities between 2534 diseases, we show that our disease similarity measure is strongly correlated with the probability of finding the diseases in the same disease family and, more importantly, sharing biological pathways. We have also compared our results to those of MimMiner, a text-mining method that assigns pairwise similarity scores to diseases. We find the results of the two methods to be complementary. It is also shown that clustering diseases based on their similarities and performing enrichment analysis for the cluster centers significantly increases the term association rate, suggesting that the cluster centers are better representatives for biological pathways than the diseases themselves. This lends support to the view that our similarity measure is a good indicator of relatedness of biological processes involved in causing the diseases. Although not needed for understanding this paper, the raw results are available for download for further study at ftp://ftp.ncbi.nlm.nih.gov/pub/qmbpmn/DiseaseRelations/. PMID:25360770

Relating diseases by integrating gene associations and information flow through protein interaction network.

PubMed

Hamaneh, Mehdi Bagheri; Yu, Yi-Kuo

2014-01-01

Identifying similar diseases could potentially provide deeper understanding of their underlying causes, and may even hint at possible treatments. For this purpose, it is necessary to have a similarity measure that reflects the underpinning molecular interactions and biological pathways. We have thus devised a network-based measure that can partially fulfill this goal. Our method assigns weights to all proteins (and consequently their encoding genes) by using information flow from a disease to the protein interaction network and back. Similarity between two diseases is then defined as the cosine of the angle between their corresponding weight vectors. The proposed method also provides a way to suggest disease-pathway associations by using the weights assigned to the genes to perform enrichment analysis for each disease. By calculating pairwise similarities between 2534 diseases, we show that our disease similarity measure is strongly correlated with the probability of finding the diseases in the same disease family and, more importantly, sharing biological pathways. We have also compared our results to those of MimMiner, a text-mining method that assigns pairwise similarity scores to diseases. We find the results of the two methods to be complementary. It is also shown that clustering diseases based on their similarities and performing enrichment analysis for the cluster centers significantly increases the term association rate, suggesting that the cluster centers are better representatives for biological pathways than the diseases themselves. This lends support to the view that our similarity measure is a good indicator of relatedness of biological processes involved in causing the diseases. Although not needed for understanding this paper, the raw results are available for download for further study at ftp://ftp.ncbi.nlm.nih.gov/pub/qmbpmn/DiseaseRelations/.

Pairwise Force SPH Model for Real-Time Multi-Interaction Applications.

PubMed

Yang, Tao; Martin, Ralph R; Lin, Ming C; Chang, Jian; Hu, Shi-Min

2017-10-01

In this paper, we present a novel pairwise-force smoothed particle hydrodynamics (PF-SPH) model to enable simulation of various interactions at interfaces in real time. Realistic capture of interactions at interfaces is a challenging problem for SPH-based simulations, especially for scenarios involving multiple interactions at different interfaces. Our PF-SPH model can readily handle multiple types of interactions simultaneously in a single simulation; its basis is to use a larger support radius than that used in standard SPH. We adopt a novel anisotropic filtering term to further improve the performance of interaction forces. The proposed model is stable; furthermore, it avoids the particle clustering problem which commonly occurs at the free surface. We show how our model can be used to capture various interactions. We also consider the close connection between droplets and bubbles, and show how to animate bubbles rising in liquid as well as bubbles in air. Our method is versatile, physically plausible and easy-to-implement. Examples are provided to demonstrate the capabilities and effectiveness of our approach.

Measuring pair-wise molecular interactions in a complex mixture

NASA Astrophysics Data System (ADS)

Chakraborty, Krishnendu; Varma, Manoj M.; Venkatapathi, Murugesan

2016-03-01

Complex biological samples such as serum contain thousands of proteins and other molecules spanning up to 13 orders of magnitude in concentration. Present measurement techniques do not permit the analysis of all pair-wise interactions between the components of such a complex mixture to a given target molecule. In this work we explore the use of nanoparticle tags which encode the identity of the molecule to obtain the statistical distribution of pair-wise interactions using their Localized Surface Plasmon Resonance (LSPR) signals. The nanoparticle tags are chosen such that the binding between two molecules conjugated to the respective nanoparticle tags can be recognized by the coupling of their LSPR signals. This numerical simulation is done by DDA to investigate this approach using a reduced system consisting of three nanoparticles (a gold ellipsoid with aspect ratio 2.5 and short axis 16 nm, and two silver ellipsoids with aspect ratios 3 and 2 and short axes 8 nm and 10 nm respectively) and the set of all possible dimers formed between them. Incident light was circularly polarized and all possible particle and dimer orientations were considered. We observed that minimum peak separation between two spectra is 5 nm while maximum is 184nm.

Design, Implementation and Deployment of PAIRwise

ERIC Educational Resources Information Center

Knight, Allan; Almeroth, Kevin; Bimber, Bruce

2008-01-01

Increased access to the Internet has dramatically increased the sources from which students can deliberately or accidentally copy information. This article discusses our motivation to design, implement, and deploy an Internet based plagiarism detection system, called PAIRwise, to address this growing problem. We give details as to how we detect…

Detecting epistasis with the marginal epistasis test in genetic mapping studies of quantitative traits

PubMed Central

Zeng, Ping; Mukherjee, Sayan; Zhou, Xiang

2017-01-01

Epistasis, commonly defined as the interaction between multiple genes, is an important genetic component underlying phenotypic variation. Many statistical methods have been developed to model and identify epistatic interactions between genetic variants. However, because of the large combinatorial search space of interactions, most epistasis mapping methods face enormous computational challenges and often suffer from low statistical power due to multiple test correction. Here, we present a novel, alternative strategy for mapping epistasis: instead of directly identifying individual pairwise or higher-order interactions, we focus on mapping variants that have non-zero marginal epistatic effects—the combined pairwise interaction effects between a given variant and all other variants. By testing marginal epistatic effects, we can identify candidate variants that are involved in epistasis without the need to identify the exact partners with which the variants interact, thus potentially alleviating much of the statistical and computational burden associated with standard epistatic mapping procedures. Our method is based on a variance component model, and relies on a recently developed variance component estimation method for efficient parameter inference and p-value computation. We refer to our method as the “MArginal ePIstasis Test”, or MAPIT. With simulations, we show how MAPIT can be used to estimate and test marginal epistatic effects, produce calibrated test statistics under the null, and facilitate the detection of pairwise epistatic interactions. We further illustrate the benefits of MAPIT in a QTL mapping study by analyzing the gene expression data of over 400 individuals from the GEUVADIS consortium. PMID:28746338

Validating two-dimensional leadership models on three-dimensionally structured fish schools

PubMed Central

Nagy, Máté; Holbrook, Robert I.; Biro, Dora; Burt de Perera, Theresa

2017-01-01

Identifying leader–follower interactions is crucial for understanding how a group decides where or when to move, and how this information is transferred between members. Although many animal groups have a three-dimensional structure, previous studies investigating leader–follower interactions have often ignored vertical information. This raises the question of whether commonly used two-dimensional leader–follower analyses can be used justifiably on groups that interact in three dimensions. To address this, we quantified the individual movements of banded tetra fish (Astyanax mexicanus) within shoals by computing the three-dimensional trajectories of all individuals using a stereo-camera technique. We used these data firstly to identify and compare leader–follower interactions in two and three dimensions, and secondly to analyse leadership with respect to an individual's spatial position in three dimensions. We show that for 95% of all pairwise interactions leadership identified through two-dimensional analysis matches that identified through three-dimensional analysis, and we reveal that fish attend to the same shoalmates for vertical information as they do for horizontal information. Our results therefore highlight that three-dimensional analyses are not always required to identify leader–follower relationships in species that move freely in three dimensions. We discuss our results in terms of the importance of taking species' sensory capacities into account when studying interaction networks within groups. PMID:28280582

Unjamming in models with analytic pairwise potentials

NASA Astrophysics Data System (ADS)

Kooij, Stefan; Lerner, Edan

2017-06-01

Canonical models for studying the unjamming scenario in systems of soft repulsive particles assume pairwise potentials with a sharp cutoff in the interaction range. The sharp cutoff renders the potential nonanalytic but makes it possible to describe many properties of the solid in terms of the coordination number z , which has an unambiguous definition in these cases. Pairwise potentials without a sharp cutoff in the interaction range have not been studied in this context, but should in fact be considered to understand the relevance of the unjamming phenomenology in systems where such a cutoff is not present. In this work we explore two systems with such interactions: an inverse power law and an exponentially decaying pairwise potential, with the control parameters being the exponent (of the inverse power law) for the former and the number density for the latter. Both systems are shown to exhibit the characteristic features of the unjamming transition, among which are the vanishing of the shear-to-bulk modulus ratio and the emergence of an excess of low-frequency vibrational modes. We establish a relation between the pressure-to-bulk modulus ratio and the distance to unjamming in each of our model systems. This allows us to predict the dependence of other key observables on the distance to unjamming. Our results provide the means for a quantitative estimation of the proximity of generic glass-forming models to the unjamming transition in the absence of a clear-cut definition of the coordination number and highlight the general irrelevance of nonaffine contributions to the bulk modulus.

Unjamming in models with analytic pairwise potentials.

PubMed

Kooij, Stefan; Lerner, Edan

2017-06-01

Canonical models for studying the unjamming scenario in systems of soft repulsive particles assume pairwise potentials with a sharp cutoff in the interaction range. The sharp cutoff renders the potential nonanalytic but makes it possible to describe many properties of the solid in terms of the coordination number z, which has an unambiguous definition in these cases. Pairwise potentials without a sharp cutoff in the interaction range have not been studied in this context, but should in fact be considered to understand the relevance of the unjamming phenomenology in systems where such a cutoff is not present. In this work we explore two systems with such interactions: an inverse power law and an exponentially decaying pairwise potential, with the control parameters being the exponent (of the inverse power law) for the former and the number density for the latter. Both systems are shown to exhibit the characteristic features of the unjamming transition, among which are the vanishing of the shear-to-bulk modulus ratio and the emergence of an excess of low-frequency vibrational modes. We establish a relation between the pressure-to-bulk modulus ratio and the distance to unjamming in each of our model systems. This allows us to predict the dependence of other key observables on the distance to unjamming. Our results provide the means for a quantitative estimation of the proximity of generic glass-forming models to the unjamming transition in the absence of a clear-cut definition of the coordination number and highlight the general irrelevance of nonaffine contributions to the bulk modulus.

ICAP - An Interactive Cluster Analysis Procedure for analyzing remotely sensed data

NASA Technical Reports Server (NTRS)

Wharton, S. W.; Turner, B. J.

1981-01-01

An Interactive Cluster Analysis Procedure (ICAP) was developed to derive classifier training statistics from remotely sensed data. ICAP differs from conventional clustering algorithms by allowing the analyst to optimize the cluster configuration by inspection, rather than by manipulating process parameters. Control of the clustering process alternates between the algorithm, which creates new centroids and forms clusters, and the analyst, who can evaluate and elect to modify the cluster structure. Clusters can be deleted, or lumped together pairwise, or new centroids can be added. A summary of the cluster statistics can be requested to facilitate cluster manipulation. The principal advantage of this approach is that it allows prior information (when available) to be used directly in the analysis, since the analyst interacts with ICAP in a straightforward manner, using basic terms with which he is more likely to be familiar. Results from testing ICAP showed that an informed use of ICAP can improve classification, as compared to an existing cluster analysis procedure.

Pair-Wise and Many-Body Dispersive Interactions Coupled to an Optimally Tuned Range-Separated Hybrid Functional.

PubMed

Agrawal, Piyush; Tkatchenko, Alexandre; Kronik, Leeor

2013-08-13

We propose a nonempirical, pair-wise or many-body dispersion-corrected, optimally tuned range-separated hybrid functional. This functional retains the advantages of the optimal-tuning approach in the prediction of the electronic structure. At the same time, it gains accuracy in the prediction of binding energies for dispersively bound systems, as demonstrated on the S22 and S66 benchmark sets of weakly bound dimers.

Generalized priority-queue network dynamics: Impact of team and hierarchy

NASA Astrophysics Data System (ADS)

Cho, Won-Kuk; Min, Byungjoon; Goh, K.-I.; Kim, I.-M.

2010-06-01

We study the effect of team and hierarchy on the waiting-time dynamics of priority-queue networks. To this end, we introduce generalized priority-queue network models incorporating interaction rules based on team-execution and hierarchy in decision making, respectively. It is numerically found that the waiting-time distribution exhibits a power law for long waiting times in both cases, yet with different exponents depending on the team size and the position of queue nodes in the hierarchy, respectively. The observed power-law behaviors have in many cases a corresponding single or pairwise-interacting queue dynamics, suggesting that the pairwise interaction may constitute a major dynamic consequence in the priority-queue networks. It is also found that the reciprocity of influence is a relevant factor for the priority-queue network dynamics.

Coevolution study of mitochondria respiratory chain proteins: toward the understanding of protein--protein interaction.

PubMed

Yang, Ming; Ge, Yan; Wu, Jiayan; Xiao, Jingfa; Yu, Jun

2011-05-20

Coevolution can be seen as the interdependency between evolutionary histories. In the context of protein evolution, functional correlation proteins are ever-present coordinated evolutionary characters without disruption of organismal integrity. As to complex system, there are two forms of protein--protein interactions in vivo, which refer to inter-complex interaction and intra-complex interaction. In this paper, we studied the difference of coevolution characters between inter-complex interaction and intra-complex interaction using "Mirror tree" method on the respiratory chain (RC) proteins. We divided the correlation coefficients of every pairwise RC proteins into two groups corresponding to the binary protein--protein interaction in intra-complex and the binary protein--protein interaction in inter-complex, respectively. A dramatical discrepancy is detected between the coevolution characters of the two sets of protein interactions (Wilcoxon test, p-value = 4.4 × 10(-6)). Our finding reveals some critical information on coevolutionary study and assists the mechanical investigation of protein--protein interaction. Furthermore, the results also provide some unique clue for supramolecular organization of protein complexes in the mitochondrial inner membrane. More detailed binding sites map and genome information of nuclear encoded RC proteins will be extraordinary valuable for the further mitochondria dynamics study. Copyright © 2011. Published by Elsevier Ltd.

Distantly related lipocalins share two conserved clusters of hydrophobic residues: use in homology modeling

PubMed Central

Adam, Benoit; Charloteaux, Benoit; Beaufays, Jerome; Vanhamme, Luc; Godfroid, Edmond; Brasseur, Robert; Lins, Laurence

2008-01-01

Background Lipocalins are widely distributed in nature and are found in bacteria, plants, arthropoda and vertebra. In hematophagous arthropods, they are implicated in the successful accomplishment of the blood meal, interfering with platelet aggregation, blood coagulation and inflammation and in the transmission of disease parasites such as Trypanosoma cruzi and Borrelia burgdorferi. The pairwise sequence identity is low among this family, often below 30%, despite a well conserved tertiary structure. Under the 30% identity threshold, alignment methods do not correctly assign and align proteins. The only safe way to assign a sequence to that family is by experimental determination. However, these procedures are long and costly and cannot always be applied. A way to circumvent the experimental approach is sequence and structure analyze. To further help in that task, the residues implicated in the stabilisation of the lipocalin fold were determined. This was done by analyzing the conserved interactions for ten lipocalins having a maximum pairwise identity of 28% and various functions. Results It was determined that two hydrophobic clusters of residues are conserved by analysing the ten lipocalin structures and sequences. One cluster is internal to the barrel, involving all strands and the 310 helix. The other is external, involving four strands and the helix lying parallel to the barrel surface. These clusters are also present in RaHBP2, a unusual "outlier" lipocalin from tick Rhipicephalus appendiculatus. This information was used to assess assignment of LIR2 a protein from Ixodes ricinus and to build a 3D model that helps to predict function. FTIR data support the lipocalin fold for this protein. Conclusion By sequence and structural analyzes, two conserved clusters of hydrophobic residues in interactions have been identified in lipocalins. Since the residues implicated are not conserved for function, they should provide the minimal subset necessary to confer the lipocalin fold. This information has been used to assign LIR2 to lipocalins and to investigate its structure/function relationship. This study could be applied to other protein families with low pairwise similarity, such as the structurally related fatty acid binding proteins or avidins. PMID:18190694

Using Replicates in Information Retrieval Evaluation.

PubMed

Voorhees, Ellen M; Samarov, Daniel; Soboroff, Ian

2017-09-01

This article explores a method for more accurately estimating the main effect of the system in a typical test-collection-based evaluation of information retrieval systems, thus increasing the sensitivity of system comparisons. Randomly partitioning the test document collection allows for multiple tests of a given system and topic (replicates). Bootstrap ANOVA can use these replicates to extract system-topic interactions-something not possible without replicates-yielding a more precise value for the system effect and a narrower confidence interval around that value. Experiments using multiple TREC collections demonstrate that removing the topic-system interactions substantially reduces the confidence intervals around the system effect as well as increases the number of significant pairwise differences found. Further, the method is robust against small changes in the number of partitions used, against variability in the documents that constitute the partitions, and the measure of effectiveness used to quantify system effectiveness.

SVM-dependent pairwise HMM: an application to protein pairwise alignments.

PubMed

Orlando, Gabriele; Raimondi, Daniele; Khan, Taushif; Lenaerts, Tom; Vranken, Wim F

2017-12-15

Methods able to provide reliable protein alignments are crucial for many bioinformatics applications. In the last years many different algorithms have been developed and various kinds of information, from sequence conservation to secondary structure, have been used to improve the alignment performances. This is especially relevant for proteins with highly divergent sequences. However, recent works suggest that different features may have different importance in diverse protein classes and it would be an advantage to have more customizable approaches, capable to deal with different alignment definitions. Here we present Rigapollo, a highly flexible pairwise alignment method based on a pairwise HMM-SVM that can use any type of information to build alignments. Rigapollo lets the user decide the optimal features to align their protein class of interest. It outperforms current state of the art methods on two well-known benchmark datasets when aligning highly divergent sequences. A Python implementation of the algorithm is available at http://ibsquare.be/rigapollo. wim.vranken@vub.be. Supplementary data are available at Bioinformatics online. © The Author (2017). Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com

Multilevel summation with B-spline interpolation for pairwise interactions in molecular dynamics simulations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hardy, David J., E-mail: dhardy@illinois.edu; Schulten, Klaus; Wolff, Matthew A.

2016-03-21

The multilevel summation method for calculating electrostatic interactions in molecular dynamics simulations constructs an approximation to a pairwise interaction kernel and its gradient, which can be evaluated at a cost that scales linearly with the number of atoms. The method smoothly splits the kernel into a sum of partial kernels of increasing range and decreasing variability with the longer-range parts interpolated from grids of increasing coarseness. Multilevel summation is especially appropriate in the context of dynamics and minimization, because it can produce continuous gradients. This article explores the use of B-splines to increase the accuracy of the multilevel summation methodmore » (for nonperiodic boundaries) without incurring additional computation other than a preprocessing step (whose cost also scales linearly). To obtain accurate results efficiently involves technical difficulties, which are overcome by a novel preprocessing algorithm. Numerical experiments demonstrate that the resulting method offers substantial improvements in accuracy and that its performance is competitive with an implementation of the fast multipole method in general and markedly better for Hamiltonian formulations of molecular dynamics. The improvement is great enough to establish multilevel summation as a serious contender for calculating pairwise interactions in molecular dynamics simulations. In particular, the method appears to be uniquely capable for molecular dynamics in two situations, nonperiodic boundary conditions and massively parallel computation, where the fast Fourier transform employed in the particle–mesh Ewald method falls short.« less

Multilevel summation with B-spline interpolation for pairwise interactions in molecular dynamics simulations.

PubMed

Hardy, David J; Wolff, Matthew A; Xia, Jianlin; Schulten, Klaus; Skeel, Robert D

2016-03-21

The multilevel summation method for calculating electrostatic interactions in molecular dynamics simulations constructs an approximation to a pairwise interaction kernel and its gradient, which can be evaluated at a cost that scales linearly with the number of atoms. The method smoothly splits the kernel into a sum of partial kernels of increasing range and decreasing variability with the longer-range parts interpolated from grids of increasing coarseness. Multilevel summation is especially appropriate in the context of dynamics and minimization, because it can produce continuous gradients. This article explores the use of B-splines to increase the accuracy of the multilevel summation method (for nonperiodic boundaries) without incurring additional computation other than a preprocessing step (whose cost also scales linearly). To obtain accurate results efficiently involves technical difficulties, which are overcome by a novel preprocessing algorithm. Numerical experiments demonstrate that the resulting method offers substantial improvements in accuracy and that its performance is competitive with an implementation of the fast multipole method in general and markedly better for Hamiltonian formulations of molecular dynamics. The improvement is great enough to establish multilevel summation as a serious contender for calculating pairwise interactions in molecular dynamics simulations. In particular, the method appears to be uniquely capable for molecular dynamics in two situations, nonperiodic boundary conditions and massively parallel computation, where the fast Fourier transform employed in the particle-mesh Ewald method falls short.

Multilevel summation with B-spline interpolation for pairwise interactions in molecular dynamics simulations

NASA Astrophysics Data System (ADS)

Hardy, David J.; Wolff, Matthew A.; Xia, Jianlin; Schulten, Klaus; Skeel, Robert D.

2016-03-01

The multilevel summation method for calculating electrostatic interactions in molecular dynamics simulations constructs an approximation to a pairwise interaction kernel and its gradient, which can be evaluated at a cost that scales linearly with the number of atoms. The method smoothly splits the kernel into a sum of partial kernels of increasing range and decreasing variability with the longer-range parts interpolated from grids of increasing coarseness. Multilevel summation is especially appropriate in the context of dynamics and minimization, because it can produce continuous gradients. This article explores the use of B-splines to increase the accuracy of the multilevel summation method (for nonperiodic boundaries) without incurring additional computation other than a preprocessing step (whose cost also scales linearly). To obtain accurate results efficiently involves technical difficulties, which are overcome by a novel preprocessing algorithm. Numerical experiments demonstrate that the resulting method offers substantial improvements in accuracy and that its performance is competitive with an implementation of the fast multipole method in general and markedly better for Hamiltonian formulations of molecular dynamics. The improvement is great enough to establish multilevel summation as a serious contender for calculating pairwise interactions in molecular dynamics simulations. In particular, the method appears to be uniquely capable for molecular dynamics in two situations, nonperiodic boundary conditions and massively parallel computation, where the fast Fourier transform employed in the particle-mesh Ewald method falls short.

Collective decision dynamics in the presence of external drivers

NASA Astrophysics Data System (ADS)

Bassett, Danielle S.; Alderson, David L.; Carlson, Jean M.

2012-09-01

We develop a sequence of models describing information transmission and decision dynamics for a network of individual agents subject to multiple sources of influence. Our general framework is set in the context of an impending natural disaster, where individuals, represented by nodes on the network, must decide whether or not to evacuate. Sources of influence include a one-to-many externally driven global broadcast as well as pairwise interactions, across links in the network, in which agents transmit either continuous opinions or binary actions. We consider both uniform and variable threshold rules on the individual opinion as baseline models for decision making. Our results indicate that (1) social networks lead to clustering and cohesive action among individuals, (2) binary information introduces high temporal variability and stagnation, and (3) information transmission over the network can either facilitate or hinder action adoption, depending on the influence of the global broadcast relative to the social network. Our framework highlights the essential role of local interactions between agents in predicting collective behavior of the population as a whole.

Pairwise-additive hydrophobic effect for alkanes in water

PubMed Central

Wu, Jianzhong; Prausnitz, John M.

2008-01-01

Pairwise additivity of the hydrophobic effect is indicated by reliable experimental Henry's constants for a large number of linear and branched low-molecular-weight alkanes in water. Pairwise additivity suggests that the hydrophobic effect is primarily a local phenomenon and that the hydrophobic interaction may be represented by a semiempirical force field. By representing the hydrophobic potential between two methane molecules as a linear function of the overlap volume of the hydration layers, we find that the contact value of the hydrophobic potential (−0.72 kcal/mol) is smaller than that from quantum mechanics simulations (−2.8 kcal/mol) but is close to that from classical molecular dynamics (−0.5∼−0.9 kcal/mol). PMID:18599448

Method for stationarity-segmentation of spike train data with application to the Pearson cross-correlation.

PubMed

Quiroga-Lombard, Claudio S; Hass, Joachim; Durstewitz, Daniel

2013-07-01

Correlations among neurons are supposed to play an important role in computation and information coding in the nervous system. Empirically, functional interactions between neurons are most commonly assessed by cross-correlation functions. Recent studies have suggested that pairwise correlations may indeed be sufficient to capture most of the information present in neural interactions. Many applications of correlation functions, however, implicitly tend to assume that the underlying processes are stationary. This assumption will usually fail for real neurons recorded in vivo since their activity during behavioral tasks is heavily influenced by stimulus-, movement-, or cognition-related processes as well as by more general processes like slow oscillations or changes in state of alertness. To address the problem of nonstationarity, we introduce a method for assessing stationarity empirically and then "slicing" spike trains into stationary segments according to the statistical definition of weak-sense stationarity. We examine pairwise Pearson cross-correlations (PCCs) under both stationary and nonstationary conditions and identify another source of covariance that can be differentiated from the covariance of the spike times and emerges as a consequence of residual nonstationarities after the slicing process: the covariance of the firing rates defined on each segment. Based on this, a correction of the PCC is introduced that accounts for the effect of segmentation. We probe these methods both on simulated data sets and on in vivo recordings from the prefrontal cortex of behaving rats. Rather than for removing nonstationarities, the present method may also be used for detecting significant events in spike trains.

Efficient conformational space exploration in ab initio protein folding simulation.

PubMed

Ullah, Ahammed; Ahmed, Nasif; Pappu, Subrata Dey; Shatabda, Swakkhar; Ullah, A Z M Dayem; Rahman, M Sohel

2015-08-01

Ab initio protein folding simulation largely depends on knowledge-based energy functions that are derived from known protein structures using statistical methods. These knowledge-based energy functions provide us with a good approximation of real protein energetics. However, these energy functions are not very informative for search algorithms and fail to distinguish the types of amino acid interactions that contribute largely to the energy function from those that do not. As a result, search algorithms frequently get trapped into the local minima. On the other hand, the hydrophobic-polar (HP) model considers hydrophobic interactions only. The simplified nature of HP energy function makes it limited only to a low-resolution model. In this paper, we present a strategy to derive a non-uniform scaled version of the real 20×20 pairwise energy function. The non-uniform scaling helps tackle the difficulty faced by a real energy function, whereas the integration of 20×20 pairwise information overcomes the limitations faced by the HP energy function. Here, we have applied a derived energy function with a genetic algorithm on discrete lattices. On a standard set of benchmark protein sequences, our approach significantly outperforms the state-of-the-art methods for similar models. Our approach has been able to explore regions of the conformational space which all the previous methods have failed to explore. Effectiveness of the derived energy function is presented by showing qualitative differences and similarities of the sampled structures to the native structures. Number of objective function evaluation in a single run of the algorithm is used as a comparison metric to demonstrate efficiency.

gRINN: a tool for calculation of residue interaction energies and protein energy network analysis of molecular dynamics simulations.

PubMed

Serçinoglu, Onur; Ozbek, Pemra

2018-05-25

Atomistic molecular dynamics (MD) simulations generate a wealth of information related to the dynamics of proteins. If properly analyzed, this information can lead to new insights regarding protein function and assist wet-lab experiments. Aiming to identify interactions between individual amino acid residues and the role played by each in the context of MD simulations, we present a stand-alone software called gRINN (get Residue Interaction eNergies and Networks). gRINN features graphical user interfaces (GUIs) and a command-line interface for generating and analyzing pairwise residue interaction energies and energy correlations from protein MD simulation trajectories. gRINN utilizes the features of NAMD or GROMACS MD simulation packages and automatizes the steps necessary to extract residue-residue interaction energies from user-supplied simulation trajectories, greatly simplifying the analysis for the end-user. A GUI, including an embedded molecular viewer, is provided for visualization of interaction energy time-series, distributions, an interaction energy matrix, interaction energy correlations and a residue correlation matrix. gRINN additionally offers construction and analysis of Protein Energy Networks, providing residue-based metrics such as degrees, betweenness-centralities, closeness centralities as well as shortest path analysis. gRINN is free and open to all users without login requirement at http://grinn.readthedocs.io.

Post-weaning social and cognitive performance of piglets raised pre-weaning either in a complex multi-suckling group housing system or in a conventional system with a crated sow.

PubMed

van Nieuwamerongen, S E; Mendl, M; Held, S; Soede, N M; Bolhuis, J E

2017-09-01

We studied the social and cognitive performance of piglets raised pre-weaning either in a conventional system with a sow in a farrowing crate (FC) or in a multi-suckling (MS) system in which 5 sows and their piglets could interact in a more physically enriched and spacious environment. After weaning at 4 weeks of age, 8 groups of 4 litter-mates per pre-weaning housing treatment were studied under equal and enriched post-weaning housing conditions. From each pen, one pair consisting of a dominant and a submissive pig was selected, based on a feed competition test (FCT) 2 weeks post-weaning. This pair was used in an informed forager test (IFT) which measured aspects of spatial learning and foraging strategies in a competitive context. During individual training, submissive (informed) pigs learned to remember a bait location in a testing arena with 8 buckets (the same bucket was baited in a search visit and a subsequent relocation visit), whereas dominant (non-informed) pigs always found the bait in a random bucket (search visits only). After learning their task, the informed pigs' individual search visit was followed by a pairwise relocation visit in which they were accompanied by the non-informed pig. Effects of pre-weaning housing treatment were not distinctly present regarding the occurrence of aggression in the FCT and the learning performance during individual training in the IFT. During paired visits, informed and non-informed pigs changed their behaviour in response to being tested pairwise instead of individually, but MS and FC pigs showed few distinct behavioural differences.

Modularity, pollination systems, and interaction turnover in plant-pollinator networks across space.

PubMed

Carstensen, Daniel W; Sabatino, Malena; Morellato, Leonor Patricia C

2016-05-01

Mutualistic interaction networks have been shown to be structurally conserved over space and time while pairwise interactions show high variability. In such networks, modularity is the division of species into compartments, or modules, where species within modules share more interactions with each other than they do with species from other modules. Such a modular structure is common in mutualistic networks and several evolutionary and ecological mechanisms have been proposed as underlying drivers. One prominent explanation is the existence of pollination syndromes where flowers tend to attract certain pollinators as determined by a set of traits. We investigate the modularity of seven community level plant-pollinator networks sampled in rupestrian grasslands, or campos rupestres, in SE Brazil. Defining pollination systems as corresponding groups of flower syndromes and pollinator functional groups, we test the two hypotheses that (1) interacting species from the same pollination system are more often assigned to the same module than interacting species from different pollination systems and; that (2) interactions between species from the same pollination system are more consistent across space than interactions between species from different pollination systems. Specifically we ask (1) whether networks are consistently modular across space; (2) whether interactions among species of the same pollination system occur more often inside modules, compared to interactions among species of different pollination systems, and finally; (3) whether the spatial variation in interaction identity, i.e., spatial interaction rewiring, is affected by trait complementarity among species as indicated by pollination systems. We confirm that networks are consistently modular across space and that interactions within pollination systems principally occur inside modules. Despite a strong tendency, we did not find a significant effect of pollination systems on the spatial consistency of pairwise interactions. These results indicate that the spatial rewiring of interactions could be constrained by pollination systems, resulting in conserved network structures in spite of high variation in pairwise interactions. Our findings suggest a relevant role of pollination systems in structuring plant-pollinator networks and we argue that structural patterns at the sub-network level can help us to fully understand how and why interactions vary across space and time.

Anisotropic interaction rules in circular motions of pigeon flocks: An empirical study based on sparse Bayesian learning

NASA Astrophysics Data System (ADS)

Chen, Duxin; Xu, Bowen; Zhu, Tao; Zhou, Tao; Zhang, Hai-Tao

2017-08-01

Coordination shall be deemed to the result of interindividual interaction among natural gregarious animal groups. However, revealing the underlying interaction rules and decision-making strategies governing highly coordinated motion in bird flocks is still a long-standing challenge. Based on analysis of high spatial-temporal resolution GPS data of three pigeon flocks, we extract the hidden interaction principle by using a newly emerging machine learning method, namely the sparse Bayesian learning. It is observed that the interaction probability has an inflection point at pairwise distance of 3-4 m closer than the average maximum interindividual distance, after which it decays strictly with rising pairwise metric distances. Significantly, the density of spatial neighbor distribution is strongly anisotropic, with an evident lack of interactions along individual velocity. Thus, it is found that in small-sized bird flocks, individuals reciprocally cooperate with a variational number of neighbors in metric space and tend to interact with closer time-varying neighbors, rather than interacting with a fixed number of topological ones. Finally, extensive numerical investigation is conducted to verify both the revealed interaction and decision-making principle during circular flights of pigeon flocks.

Registration of segmented histological images using thin plate splines and belief propagation

NASA Astrophysics Data System (ADS)

Kybic, Jan

2014-03-01

We register images based on their multiclass segmentations, for cases when correspondence of local features cannot be established. A discrete mutual information is used as a similarity criterion. It is evaluated at a sparse set of location on the interfaces between classes. A thin-plate spline regularization is approximated by pairwise interactions. The problem is cast into a discrete setting and solved efficiently by belief propagation. Further speedup and robustness is provided by a multiresolution framework. Preliminary experiments suggest that our method can provide similar registration quality to standard methods at a fraction of the computational cost.

Numerical approach for finite volume three-body interaction

NASA Astrophysics Data System (ADS)

Guo, Peng; Gasparian, Vladimir

2018-01-01

In the present work, we study a numerical approach to one dimensional finite volume three-body interaction, the method is demonstrated by considering a toy model of three spinless particles interacting with pair-wise δ -function potentials. The numerical results are compared with the exact solutions of three spinless bosons interaction when the strength of short-range interactions are set equal for all pairs.

Classification of forest-based ecotourism areas in Pocahontas County of West Virginia using GIS and pairwise comparison method

Treesearch

Ishwar Dhami; Jinyang. Deng

2012-01-01

Many previous studies have examined ecotourism primarily from the perspective of tourists while largely ignoring ecotourism destinations. This study used geographical information system (GIS) and pairwise comparison to identify forest-based ecotourism areas in Pocahontas County, West Virginia. The study adopted the criteria and scores developed by Boyd and Butler (1994...

Pairwise Classifier Ensemble with Adaptive Sub-Classifiers for fMRI Pattern Analysis.

PubMed

Kim, Eunwoo; Park, HyunWook

2017-02-01

The multi-voxel pattern analysis technique is applied to fMRI data for classification of high-level brain functions using pattern information distributed over multiple voxels. In this paper, we propose a classifier ensemble for multiclass classification in fMRI analysis, exploiting the fact that specific neighboring voxels can contain spatial pattern information. The proposed method converts the multiclass classification to a pairwise classifier ensemble, and each pairwise classifier consists of multiple sub-classifiers using an adaptive feature set for each class-pair. Simulated and real fMRI data were used to verify the proposed method. Intra- and inter-subject analyses were performed to compare the proposed method with several well-known classifiers, including single and ensemble classifiers. The comparison results showed that the proposed method can be generally applied to multiclass classification in both simulations and real fMRI analyses.

GetReal in network meta-analysis: a review of the methodology.

PubMed

Efthimiou, Orestis; Debray, Thomas P A; van Valkenhoef, Gert; Trelle, Sven; Panayidou, Klea; Moons, Karel G M; Reitsma, Johannes B; Shang, Aijing; Salanti, Georgia

2016-09-01

Pairwise meta-analysis is an established statistical tool for synthesizing evidence from multiple trials, but it is informative only about the relative efficacy of two specific interventions. The usefulness of pairwise meta-analysis is thus limited in real-life medical practice, where many competing interventions may be available for a certain condition and studies informing some of the pairwise comparisons may be lacking. This commonly encountered scenario has led to the development of network meta-analysis (NMA). In the last decade, several applications, methodological developments, and empirical studies in NMA have been published, and the area is thriving as its relevance to public health is increasingly recognized. This article presents a review of the relevant literature on NMA methodology aiming to pinpoint the developments that have appeared in the field. Copyright © 2016 John Wiley & Sons, Ltd. Copyright © 2016 John Wiley & Sons, Ltd.

A general transformation to canonical form for potentials in pairwise interatomic interactions.

PubMed

Walton, Jay R; Rivera-Rivera, Luis A; Lucchese, Robert R; Bevan, John W

2015-06-14

A generalized formulation of explicit force-based transformations is introduced to investigate the concept of a canonical potential in both fundamental chemical and intermolecular bonding. Different classes of representative ground electronic state pairwise interatomic interactions are referenced to a chosen canonical potential illustrating application of such transformations. Specifically, accurately determined potentials of the diatomic molecules H2, H2(+), HF, LiH, argon dimer, and one-dimensional dissociative coordinates in Ar-HBr, OC-HF, and OC-Cl2 are investigated throughout their bound potentials. Advantages of the current formulation for accurately evaluating equilibrium dissociation energies and a fundamentally different unified perspective on nature of intermolecular interactions will be emphasized. In particular, this canonical approach has significance to previous assertions that there is no very fundamental distinction between van der Waals bonding and covalent bonding or for that matter hydrogen and halogen bonds.

Solution to urn models of pairwise interaction with application to social, physical, and biological sciences

NASA Astrophysics Data System (ADS)

Pickering, William; Lim, Chjan

2017-07-01

We investigate a family of urn models that correspond to one-dimensional random walks with quadratic transition probabilities that have highly diverse applications. Well-known instances of these two-urn models are the Ehrenfest model of molecular diffusion, the voter model of social influence, and the Moran model of population genetics. We also provide a generating function method for diagonalizing the corresponding transition matrix that is valid if and only if the underlying mean density satisfies a linear differential equation and express the eigenvector components as terms of ordinary hypergeometric functions. The nature of the models lead to a natural extension to interaction between agents in a general network topology. We analyze the dynamics on uncorrelated heterogeneous degree sequence networks and relate the convergence times to the moments of the degree sequences for various pairwise interaction mechanisms.

Network Analysis of Protein Adaptation: Modeling the Functional Impact of Multiple Mutations

PubMed Central

Beleva Guthrie, Violeta; Masica, David L; Fraser, Andrew; Federico, Joseph; Fan, Yunfan; Camps, Manel; Karchin, Rachel

2018-01-01

Abstract The evolution of new biochemical activities frequently involves complex dependencies between mutations and rapid evolutionary radiation. Mutation co-occurrence and covariation have previously been used to identify compensating mutations that are the result of physical contacts and preserve protein function and fold. Here, we model pairwise functional dependencies and higher order interactions that enable evolution of new protein functions. We use a network model to find complex dependencies between mutations resulting from evolutionary trade-offs and pleiotropic effects. We present a method to construct these networks and to identify functionally interacting mutations in both extant and reconstructed ancestral sequences (Network Analysis of Protein Adaptation). The time ordering of mutations can be incorporated into the networks through phylogenetic reconstruction. We apply NAPA to three distantly homologous β-lactamase protein clusters (TEM, CTX-M-3, and OXA-51), each of which has experienced recent evolutionary radiation under substantially different selective pressures. By analyzing the network properties of each protein cluster, we identify key adaptive mutations, positive pairwise interactions, different adaptive solutions to the same selective pressure, and complex evolutionary trajectories likely to increase protein fitness. We also present evidence that incorporating information from phylogenetic reconstruction and ancestral sequence inference can reduce the number of spurious links in the network, whereas preserving overall network community structure. The analysis does not require structural or biochemical data. In contrast to function-preserving mutation dependencies, which are frequently from structural contacts, gain-of-function mutation dependencies are most commonly between residues distal in protein structure. PMID:29522102

Controllability of symmetric spin networks

NASA Astrophysics Data System (ADS)

Albertini, Francesca; D'Alessandro, Domenico

2018-05-01

We consider a network of n spin 1/2 systems which are pairwise interacting via Ising interaction and are controlled by the same electro-magnetic control field. Such a system presents symmetries since the Hamiltonian is unchanged if we permute two spins. This prevents full (operator) controllability, in that not every unitary evolution can be obtained. We prove however that controllability is verified if we restrict ourselves to unitary evolutions which preserve the above permutation invariance. For low dimensional cases, n = 2 and n = 3, we provide an analysis of the Lie group of available evolutions and give explicit control laws to transfer between two arbitrary permutation invariant states. This class of states includes highly entangled states such as Greenberger-Horne-Zeilinger (GHZ) states and W states, which are of interest in quantum information.

SubVis: an interactive R package for exploring the effects of multiple substitution matrices on pairwise sequence alignment

PubMed Central

Coan, Heather B.; Youker, Robert T.

2017-01-01

Understanding how proteins mutate is critical to solving a host of biological problems. Mutations occur when an amino acid is substituted for another in a protein sequence. The set of likelihoods for amino acid substitutions is stored in a matrix and input to alignment algorithms. The quality of the resulting alignment is used to assess the similarity of two or more sequences and can vary according to assumptions modeled by the substitution matrix. Substitution strategies with minor parameter variations are often grouped together in families. For example, the BLOSUM and PAM matrix families are commonly used because they provide a standard, predefined way of modeling substitutions. However, researchers often do not know if a given matrix family or any individual matrix within a family is the most suitable. Furthermore, predefined matrix families may inaccurately reflect a particular hypothesis that a researcher wishes to model or otherwise result in unsatisfactory alignments. In these cases, the ability to compare the effects of one or more custom matrices may be needed. This laborious process is often performed manually because the ability to simultaneously load multiple matrices and then compare their effects on alignments is not readily available in current software tools. This paper presents SubVis, an interactive R package for loading and applying multiple substitution matrices to pairwise alignments. Users can simultaneously explore alignments resulting from multiple predefined and custom substitution matrices. SubVis utilizes several of the alignment functions found in R, a common language among protein scientists. Functions are tied together with the Shiny platform which allows the modification of input parameters. Information regarding alignment quality and individual amino acid substitutions is displayed with the JavaScript language which provides interactive visualizations for revealing both high-level and low-level alignment information. PMID:28674656

InfAcrOnt: calculating cross-ontology term similarities using information flow by a random walk.

PubMed

Cheng, Liang; Jiang, Yue; Ju, Hong; Sun, Jie; Peng, Jiajie; Zhou, Meng; Hu, Yang

2018-01-19

Since the establishment of the first biomedical ontology Gene Ontology (GO), the number of biomedical ontology has increased dramatically. Nowadays over 300 ontologies have been built including extensively used Disease Ontology (DO) and Human Phenotype Ontology (HPO). Because of the advantage of identifying novel relationships between terms, calculating similarity between ontology terms is one of the major tasks in this research area. Though similarities between terms within each ontology have been studied with in silico methods, term similarities across different ontologies were not investigated as deeply. The latest method took advantage of gene functional interaction network (GFIN) to explore such inter-ontology similarities of terms. However, it only used gene interactions and failed to make full use of the connectivity among gene nodes of the network. In addition, all existent methods are particularly designed for GO and their performances on the extended ontology community remain unknown. We proposed a method InfAcrOnt to infer similarities between terms across ontologies utilizing the entire GFIN. InfAcrOnt builds a term-gene-gene network which comprised ontology annotations and GFIN, and acquires similarities between terms across ontologies through modeling the information flow within the network by random walk. In our benchmark experiments on sub-ontologies of GO, InfAcrOnt achieves a high average area under the receiver operating characteristic curve (AUC) (0.9322 and 0.9309) and low standard deviations (1.8746e-6 and 3.0977e-6) in both human and yeast benchmark datasets exhibiting superior performance. Meanwhile, comparisons of InfAcrOnt results and prior knowledge on pair-wise DO-HPO terms and pair-wise DO-GO terms show high correlations. The experiment results show that InfAcrOnt significantly improves the performance of inferring similarities between terms across ontologies in benchmark set.

Comparative analysis of microarray data in Arabidopsis transcriptome during compatible interactions with plant viruses

USDA-ARS?s Scientific Manuscript database

To analyze transcriptome response to virus infection, we have assembled currently available microarray data on changes in gene expression levels in compatible Arabidopsis-virus interactions. We used the mean r (Pearson’s correlation coefficient) for neighboring pairs to estimate pairwise local simil...

Referring to the social performance promotes cooperation in spatial prisoner's dilemma games

NASA Astrophysics Data System (ADS)

Shigaki, Keizo; Tanimoto, Jun; Wang, Zhen; Kokubo, Satoshi; Hagishima, Aya; Ikegaya, Naoki

2012-09-01

We propose a new pairwise Fermi updating rule by considering a social average payoff when an agent copies a neighbor's strategy. In the update rule, a focal agent compares her payoff with the social average payoff of the same strategy that her pairwise opponent has. This concept might be justified by the fact that people reference global and, somehow, statistical information, not local information when imitating social behaviors. We presume several possible ways for the social average. Simulation results prove that the social average of some limited agents realizes more significant cooperation than that of the entire population.

Pairwise adaptive thermostats for improved accuracy and stability in dissipative particle dynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Leimkuhler, Benedict, E-mail: b.leimkuhler@ed.ac.uk; Shang, Xiaocheng, E-mail: x.shang@brown.edu

2016-11-01

We examine the formulation and numerical treatment of dissipative particle dynamics (DPD) and momentum-conserving molecular dynamics. We show that it is possible to improve both the accuracy and the stability of DPD by employing a pairwise adaptive Langevin thermostat that precisely matches the dynamical characteristics of DPD simulations (e.g., autocorrelation functions) while automatically correcting thermodynamic averages using a negative feedback loop. In the low friction regime, it is possible to replace DPD by a simpler momentum-conserving variant of the Nosé–Hoover–Langevin method based on thermostatting only pairwise interactions; we show that this method has an extra order of accuracy for anmore » important class of observables (a superconvergence result), while also allowing larger timesteps than alternatives. All the methods mentioned in the article are easily implemented. Numerical experiments are performed in both equilibrium and nonequilibrium settings; using Lees–Edwards boundary conditions to induce shear flow.« less

Systematic chemical-genetic and chemical-chemical interaction datasets for prediction of compound synergism

PubMed Central

Wildenhain, Jan; Spitzer, Michaela; Dolma, Sonam; Jarvik, Nick; White, Rachel; Roy, Marcia; Griffiths, Emma; Bellows, David S.; Wright, Gerard D.; Tyers, Mike

2016-01-01

The network structure of biological systems suggests that effective therapeutic intervention may require combinations of agents that act synergistically. However, a dearth of systematic chemical combination datasets have limited the development of predictive algorithms for chemical synergism. Here, we report two large datasets of linked chemical-genetic and chemical-chemical interactions in the budding yeast Saccharomyces cerevisiae. We screened 5,518 unique compounds against 242 diverse yeast gene deletion strains to generate an extended chemical-genetic matrix (CGM) of 492,126 chemical-gene interaction measurements. This CGM dataset contained 1,434 genotype-specific inhibitors, termed cryptagens. We selected 128 structurally diverse cryptagens and tested all pairwise combinations to generate a benchmark dataset of 8,128 pairwise chemical-chemical interaction tests for synergy prediction, termed the cryptagen matrix (CM). An accompanying database resource called ChemGRID was developed to enable analysis, visualisation and downloads of all data. The CGM and CM datasets will facilitate the benchmarking of computational approaches for synergy prediction, as well as chemical structure-activity relationship models for anti-fungal drug discovery. PMID:27874849

Identification of multi-loci hubs from 4C-seq demonstrates the functional importance of simultaneous interactions.

PubMed

Jiang, Tingting; Raviram, Ramya; Snetkova, Valentina; Rocha, Pedro P; Proudhon, Charlotte; Badri, Sana; Bonneau, Richard; Skok, Jane A; Kluger, Yuval

2016-10-14

Use of low resolution single cell DNA FISH and population based high resolution chromosome conformation capture techniques have highlighted the importance of pairwise chromatin interactions in gene regulation. However, it is unlikely that associations involving regulatory elements act in isolation of other interacting partners that also influence their impact. Indeed, the influence of multi-loci interactions remains something of an enigma as beyond low-resolution DNA FISH we do not have the appropriate tools to analyze these. Here we present a method that uses standard 4C-seq data to identify multi-loci interactions from the same cell. We demonstrate the feasibility of our method using 4C-seq data sets that identify known pairwise and novel tri-loci interactions involving the Tcrb and Igk antigen receptor enhancers. We further show that the three Igk enhancers, MiEκ, 3'Eκ and Edκ, interact simultaneously in this super-enhancer cluster, which add to our previous findings showing that loss of one element decreases interactions between all three elements as well as reducing their transcriptional output. These findings underscore the functional importance of simultaneous interactions and provide new insight into the relationship between enhancer elements. Our method opens the door for studying multi-loci interactions and their impact on gene regulation in other biological settings. © The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research.

Identification of multi-loci hubs from 4C-seq demonstrates the functional importance of simultaneous interactions

PubMed Central

Jiang, Tingting; Raviram, Ramya; Snetkova, Valentina; Rocha, Pedro P.; Proudhon, Charlotte; Badri, Sana; Bonneau, Richard; Skok, Jane A.; Kluger, Yuval

2016-01-01

Use of low resolution single cell DNA FISH and population based high resolution chromosome conformation capture techniques have highlighted the importance of pairwise chromatin interactions in gene regulation. However, it is unlikely that associations involving regulatory elements act in isolation of other interacting partners that also influence their impact. Indeed, the influence of multi-loci interactions remains something of an enigma as beyond low-resolution DNA FISH we do not have the appropriate tools to analyze these. Here we present a method that uses standard 4C-seq data to identify multi-loci interactions from the same cell. We demonstrate the feasibility of our method using 4C-seq data sets that identify known pairwise and novel tri-loci interactions involving the Tcrb and Igk antigen receptor enhancers. We further show that the three Igk enhancers, MiEκ, 3′Eκ and Edκ, interact simultaneously in this super-enhancer cluster, which add to our previous findings showing that loss of one element decreases interactions between all three elements as well as reducing their transcriptional output. These findings underscore the functional importance of simultaneous interactions and provide new insight into the relationship between enhancer elements. Our method opens the door for studying multi-loci interactions and their impact on gene regulation in other biological settings. PMID:27439714

Collective memory in primate conflict implied by temporal scaling collapse.

PubMed

Lee, Edward D; Daniels, Bryan C; Krakauer, David C; Flack, Jessica C

2017-09-01

In biological systems, prolonged conflict is costly, whereas contained conflict permits strategic innovation and refinement. Causes of variation in conflict size and duration are not well understood. We use a well-studied primate society model system to study how conflicts grow. We find conflict duration is a 'first to fight' growth process that scales superlinearly, with the number of possible pairwise interactions. This is in contrast with a 'first to fail' process that characterizes peaceful durations. Rescaling conflict distributions reveals a universal curve, showing that the typical time scale of correlated interactions exceeds nearly all individual fights. This temporal correlation implies collective memory across pairwise interactions beyond those assumed in standard models of contagion growth or iterated evolutionary games. By accounting for memory, we make quantitative predictions for interventions that mitigate or enhance the spread of conflict. Managing conflict involves balancing the efficient use of limited resources with an intervention strategy that allows for conflict while keeping it contained and controlled. © 2017 The Author(s).

DockTrina: docking triangular protein trimers.

PubMed

Popov, Petr; Ritchie, David W; Grudinin, Sergei

2014-01-01

In spite of the abundance of oligomeric proteins within a cell, the structural characterization of protein-protein interactions is still a challenging task. In particular, many of these interactions involve heteromeric complexes, which are relatively difficult to determine experimentally. Hence there is growing interest in using computational techniques to model such complexes. However, assembling large heteromeric complexes computationally is a highly combinatorial problem. Nonetheless the problem can be simplified greatly by considering interactions between protein trimers. After dimers and monomers, triangular trimers (i.e. trimers with pair-wise contacts between all three pairs of proteins) are the most frequently observed quaternary structural motifs according to the three-dimensional (3D) complex database. This article presents DockTrina, a novel protein docking method for modeling the 3D structures of nonsymmetrical triangular trimers. The method takes as input pair-wise contact predictions from a rigid body docking program. It then scans and scores all possible combinations of pairs of monomers using a very fast root mean square deviation test. Finally, it ranks the predictions using a scoring function which combines triples of pair-wise contact terms and a geometric clash penalty term. The overall approach takes less than 2 min per complex on a modern desktop computer. The method is tested and validated using a benchmark set of 220 bound and seven unbound protein trimer structures. DockTrina will be made available at http://nano-d.inrialpes.fr/software/docktrina. Copyright © 2013 Wiley Periodicals, Inc.

Evaluation of lattice sums by the Poisson sum formula

NASA Technical Reports Server (NTRS)

Ray, R. D.

1975-01-01

The Poisson sum formula was applied to the problem of summing pairwise interactions between an observer molecule and a semi-infinite regular array of solid state molecules. The transformed sum is often much more rapidly convergent than the original sum, and forms a Fourier series in the solid surface coordinates. The method is applicable to a variety of solid state structures and functional forms of the pairwise potential. As an illustration of the method, the electric field above the (100) face of the CsCl structure is calculated and compared to earlier results obtained by direct summation.

Tools for Protecting the Privacy of Specific Individuals in Video

NASA Astrophysics Data System (ADS)

Chen, Datong; Chang, Yi; Yan, Rong; Yang, Jie

2007-12-01

This paper presents a system for protecting the privacy of specific individuals in video recordings. We address the following two problems: automatic people identification with limited labeled data, and human body obscuring with preserved structure and motion information. In order to address the first problem, we propose a new discriminative learning algorithm to improve people identification accuracy using limited training data labeled from the original video and imperfect pairwise constraints labeled from face obscured video data. We employ a robust face detection and tracking algorithm to obscure human faces in the video. Our experiments in a nursing home environment show that the system can obtain a high accuracy of people identification using limited labeled data and noisy pairwise constraints. The study result indicates that human subjects can perform reasonably well in labeling pairwise constraints with the face masked data. For the second problem, we propose a novel method of body obscuring, which removes the appearance information of the people while preserving rich structure and motion information. The proposed approach provides a way to minimize the risk of exposing the identities of the protected people while maximizing the use of the captured data for activity/behavior analysis.

Random Partition Distribution Indexed by Pairwise Information

PubMed Central

Dahl, David B.; Day, Ryan; Tsai, Jerry W.

2017-01-01

We propose a random partition distribution indexed by pairwise similarity information such that partitions compatible with the similarities are given more probability. The use of pairwise similarities, in the form of distances, is common in some clustering algorithms (e.g., hierarchical clustering), but we show how to use this type of information to define a prior partition distribution for flexible Bayesian modeling. A defining feature of the distribution is that it allocates probability among partitions within a given number of subsets, but it does not shift probability among sets of partitions with different numbers of subsets. Our distribution places more probability on partitions that group similar items yet keeps the total probability of partitions with a given number of subsets constant. The distribution of the number of subsets (and its moments) is available in closed-form and is not a function of the similarities. Our formulation has an explicit probability mass function (with a tractable normalizing constant) so the full suite of MCMC methods may be used for posterior inference. We compare our distribution with several existing partition distributions, showing that our formulation has attractive properties. We provide three demonstrations to highlight the features and relative performance of our distribution. PMID:29276318

Active learning for semi-supervised clustering based on locally linear propagation reconstruction.

PubMed

Chang, Chin-Chun; Lin, Po-Yi

2015-03-01

The success of semi-supervised clustering relies on the effectiveness of side information. To get effective side information, a new active learner learning pairwise constraints known as must-link and cannot-link constraints is proposed in this paper. Three novel techniques are developed for learning effective pairwise constraints. The first technique is used to identify samples less important to cluster structures. This technique makes use of a kernel version of locally linear embedding for manifold learning. Samples neither important to locally linear propagation reconstructions of other samples nor on flat patches in the learned manifold are regarded as unimportant samples. The second is a novel criterion for query selection. This criterion considers not only the importance of a sample to expanding the space coverage of the learned samples but also the expected number of queries needed to learn the sample. To facilitate semi-supervised clustering, the third technique yields inferred must-links for passing information about flat patches in the learned manifold to semi-supervised clustering algorithms. Experimental results have shown that the learned pairwise constraints can capture the underlying cluster structures and proven the feasibility of the proposed approach. Copyright © 2014 Elsevier Ltd. All rights reserved.

Multipole Algorithms for Molecular Dynamics Simulation on High Performance Computers.

NASA Astrophysics Data System (ADS)

Elliott, William Dewey

1995-01-01

A fundamental problem in modeling large molecular systems with molecular dynamics (MD) simulations is the underlying N-body problem of computing the interactions between all pairs of N atoms. The simplest algorithm to compute pair-wise atomic interactions scales in runtime {cal O}(N^2), making it impractical for interesting biomolecular systems, which can contain millions of atoms. Recently, several algorithms have become available that solve the N-body problem by computing the effects of all pair-wise interactions while scaling in runtime less than {cal O}(N^2). One algorithm, which scales {cal O}(N) for a uniform distribution of particles, is called the Greengard-Rokhlin Fast Multipole Algorithm (FMA). This work describes an FMA-like algorithm called the Molecular Dynamics Multipole Algorithm (MDMA). The algorithm contains several features that are new to N-body algorithms. MDMA uses new, efficient series expansion equations to compute general 1/r^{n } potentials to arbitrary accuracy. In particular, the 1/r Coulomb potential and the 1/r^6 portion of the Lennard-Jones potential are implemented. The new equations are based on multivariate Taylor series expansions. In addition, MDMA uses a cell-to-cell interaction region of cells that is closely tied to worst case error bounds. The worst case error bounds for MDMA are derived in this work also. These bounds apply to other multipole algorithms as well. Several implementation enhancements are described which apply to MDMA as well as other N-body algorithms such as FMA and tree codes. The mathematics of the cell -to-cell interactions are converted to the Fourier domain for reduced operation count and faster computation. A relative indexing scheme was devised to locate cells in the interaction region which allows efficient pre-computation of redundant information and prestorage of much of the cell-to-cell interaction. Also, MDMA was integrated into the MD program SIgMA to demonstrate the performance of the program over several simulation timesteps. One MD application described here highlights the utility of including long range contributions to Lennard-Jones potential in constant pressure simulations. Another application shows the time dependence of long range forces in a multiple time step MD simulation.

BiGGER: a new (soft) docking algorithm for predicting protein interactions.

PubMed

Palma, P N; Krippahl, L; Wampler, J E; Moura, J J

2000-06-01

A new computationally efficient and automated "soft docking" algorithm is described to assist the prediction of the mode of binding between two proteins, using the three-dimensional structures of the unbound molecules. The method is implemented in a software package called BiGGER (Bimolecular Complex Generation with Global Evaluation and Ranking) and works in two sequential steps: first, the complete 6-dimensional binding spaces of both molecules is systematically searched. A population of candidate protein-protein docked geometries is thus generated and selected on the basis of the geometric complementarity and amino acid pairwise affinities between the two molecular surfaces. Most of the conformational changes observed during protein association are treated in an implicit way and test results are equally satisfactory, regardless of starting from the bound or the unbound forms of known structures of the interacting proteins. In contrast to other methods, the entire molecular surfaces are searched during the simulation, using absolutely no additional information regarding the binding sites. In a second step, an interaction scoring function is used to rank the putative docked structures. The function incorporates interaction terms that are thought to be relevant to the stabilization of protein complexes. These include: geometric complementarity of the surfaces, explicit electrostatic interactions, desolvation energy, and pairwise propensities of the amino acid side chains to contact across the molecular interface. The relative functional contribution of each of these interaction terms to the global scoring function has been empirically adjusted through a neural network optimizer using a learning set of 25 protein-protein complexes of known crystallographic structures. In 22 out of 25 protein-protein complexes tested, near-native docked geometries were found with C(alpha) RMS deviations < or =4.0 A from the experimental structures, of which 14 were found within the 20 top ranking solutions. The program works on widely available personal computers and takes 2 to 8 hours of CPU time to run any of the docking tests herein presented. Finally, the value and limitations of the method for the study of macromolecular interactions, not yet revealed by experimental techniques, are discussed.

Highly viscous antibody solutions are a consequence of network formation caused by domain-domain electrostatic complementarities: insights from coarse-grained simulations.

PubMed

Buck, Patrick M; Chaudhri, Anuj; Kumar, Sandeep; Singh, Satish K

2015-01-05

Therapeutic monoclonal antibody (mAb) candidates that form highly viscous solutions at concentrations above 100 mg/mL can lead to challenges in bioprocessing, formulation development, and subcutaneous drug delivery. Earlier studies of mAbs with concentration-dependent high viscosity have indicated that mAbs with negatively charged Fv regions have a dipole-like quality that increases the likelihood of reversible self-association. This suggests that weak electrostatic intermolecular interactions can form transient antibody networks that participate in resistance to solution deformation under shear stress. Here this hypothesis is explored by parametrizing a coarse-grained (CG) model of an antibody using the domain charges from four different mAbs that have had their concentration-dependent viscosity behaviors previously determined. Multicopy molecular dynamics simulations were performed for these four CG mAbs at several concentrations to understand the effect of surface charge on mass diffusivity, pairwise interactions, and electrostatic network formation. Diffusion coefficients computed from simulations were in qualitative agreement with experimentally determined viscosities for all four mAbs. Contact analysis revealed an overall greater number of pairwise interactions for the two mAbs in this study with high concentration viscosity issues. Further, using equilibrated solution trajectories, the two mAbs with high concentration viscosity issues quantitatively formed more features of an electrostatic network than the other mAbs. The change in the number of these network features as a function of concentration is related to the number of pairwise interactions formed by electrostatic complementarities between antibody domains. Thus, transient antibody network formation caused by domain-domain electrostatic complementarities is the most probable origin of high concentration viscosity for mAbs in this study.

ICAP: An Interactive Cluster Analysis Procedure for analyzing remotely sensed data. [to classify the radiance data to produce a thematic map

NASA Technical Reports Server (NTRS)

Wharton, S. W.

1980-01-01

An Interactive Cluster Analysis Procedure (ICAP) was developed to derive classifier training statistics from remotely sensed data. The algorithm interfaces the rapid numerical processing capacity of a computer with the human ability to integrate qualitative information. Control of the clustering process alternates between the algorithm, which creates new centroids and forms clusters and the analyst, who evaluate and elect to modify the cluster structure. Clusters can be deleted or lumped pairwise, or new centroids can be added. A summary of the cluster statistics can be requested to facilitate cluster manipulation. The ICAP was implemented in APL (A Programming Language), an interactive computer language. The flexibility of the algorithm was evaluated using data from different LANDSAT scenes to simulate two situations: one in which the analyst is assumed to have no prior knowledge about the data and wishes to have the clusters formed more or less automatically; and the other in which the analyst is assumed to have some knowledge about the data structure and wishes to use that information to closely supervise the clustering process. For comparison, an existing clustering method was also applied to the two data sets.

Pairwise contact energy statistical potentials can help to find probability of point mutations.

PubMed

Saravanan, K M; Suvaithenamudhan, S; Parthasarathy, S; Selvaraj, S

2017-01-01

To adopt a particular fold, a protein requires several interactions between its amino acid residues. The energetic contribution of these residue-residue interactions can be approximated by extracting statistical potentials from known high resolution structures. Several methods based on statistical potentials extracted from unrelated proteins are found to make a better prediction of probability of point mutations. We postulate that the statistical potentials extracted from known structures of similar folds with varying sequence identity can be a powerful tool to examine probability of point mutation. By keeping this in mind, we have derived pairwise residue and atomic contact energy potentials for the different functional families that adopt the (α/β) 8 TIM-Barrel fold. We carried out computational point mutations at various conserved residue positions in yeast Triose phosphate isomerase enzyme for which experimental results are already reported. We have also performed molecular dynamics simulations on a subset of point mutants to make a comparative study. The difference in pairwise residue and atomic contact energy of wildtype and various point mutations reveals probability of mutations at a particular position. Interestingly, we found that our computational prediction agrees with the experimental studies of Silverman et al. (Proc Natl Acad Sci 2001;98:3092-3097) and perform better prediction than i Mutant and Cologne University Protein Stability Analysis Tool. The present work thus suggests deriving pairwise contact energy potentials and molecular dynamics simulations of functionally important folds could help us to predict probability of point mutations which may ultimately reduce the time and cost of mutation experiments. Proteins 2016; 85:54-64. © 2016 Wiley Periodicals, Inc. © 2016 Wiley Periodicals, Inc.

The Influence of Higher-Order Epistasis on Biological Fitness Landscape Topography

NASA Astrophysics Data System (ADS)

Weinreich, Daniel M.; Lan, Yinghong; Jaffe, Jacob; Heckendorn, Robert B.

2018-07-01

The effect of a mutation on the organism often depends on what other mutations are already present in its genome. Geneticists refer to such mutational interactions as epistasis. Pairwise epistatic effects have been recognized for over a century, and their evolutionary implications have received theoretical attention for nearly as long. However, pairwise epistatic interactions themselves can vary with genomic background. This is called higher-order epistasis, and its consequences for evolution are much less well understood. Here, we assess the influence that higher-order epistasis has on the topography of 16 published, biological fitness landscapes. We find that on average, their effects on fitness landscape declines with order, and suggest that notable exceptions to this trend may deserve experimental scrutiny. We conclude by highlighting opportunities for further theoretical and experimental work dissecting the influence that epistasis of all orders has on fitness landscape topography and on the efficiency of evolution by natural selection.

The Influence of Higher-Order Epistasis on Biological Fitness Landscape Topography

NASA Astrophysics Data System (ADS)

Weinreich, Daniel M.; Lan, Yinghong; Jaffe, Jacob; Heckendorn, Robert B.

2018-02-01

The effect of a mutation on the organism often depends on what other mutations are already present in its genome. Geneticists refer to such mutational interactions as epistasis. Pairwise epistatic effects have been recognized for over a century, and their evolutionary implications have received theoretical attention for nearly as long. However, pairwise epistatic interactions themselves can vary with genomic background. This is called higher-order epistasis, and its consequences for evolution are much less well understood. Here, we assess the influence that higher-order epistasis has on the topography of 16 published, biological fitness landscapes. We find that on average, their effects on fitness landscape declines with order, and suggest that notable exceptions to this trend may deserve experimental scrutiny. We conclude by highlighting opportunities for further theoretical and experimental work dissecting the influence that epistasis of all orders has on fitness landscape topography and on the efficiency of evolution by natural selection.

Hierarchical structure of biological systems

PubMed Central

Alcocer-Cuarón, Carlos; Rivera, Ana L; Castaño, Victor M

2014-01-01

A general theory of biological systems, based on few fundamental propositions, allows a generalization of both Wierner and Berthalanffy approaches to theoretical biology. Here, a biological system is defined as a set of self-organized, differentiated elements that interact pair-wise through various networks and media, isolated from other sets by boundaries. Their relation to other systems can be described as a closed loop in a steady-state, which leads to a hierarchical structure and functioning of the biological system. Our thermodynamical approach of hierarchical character can be applied to biological systems of varying sizes through some general principles, based on the exchange of energy information and/or mass from and within the systems. PMID:24145961

Hierarchical structure of biological systems: a bioengineering approach.

PubMed

Alcocer-Cuarón, Carlos; Rivera, Ana L; Castaño, Victor M

2014-01-01

A general theory of biological systems, based on few fundamental propositions, allows a generalization of both Wierner and Berthalanffy approaches to theoretical biology. Here, a biological system is defined as a set of self-organized, differentiated elements that interact pair-wise through various networks and media, isolated from other sets by boundaries. Their relation to other systems can be described as a closed loop in a steady-state, which leads to a hierarchical structure and functioning of the biological system. Our thermodynamical approach of hierarchical character can be applied to biological systems of varying sizes through some general principles, based on the exchange of energy information and/or mass from and within the systems.

PTM Along Track Algorithm to Maintain Spacing During Same Direction Pair-Wise Trajectory Management Operations

NASA Technical Reports Server (NTRS)

Carreno, Victor A.

2015-01-01

Pair-wise Trajectory Management (PTM) is a cockpit based delegated responsibility separation standard. When an air traffic service provider gives a PTM clearance to an aircraft and the flight crew accepts the clearance, the flight crew will maintain spacing and separation from a designated aircraft. A PTM along track algorithm will receive state information from the designated aircraft and from the own ship to produce speed guidance for the flight crew to maintain spacing and separation

Adaptive multi-view clustering based on nonnegative matrix factorization and pairwise co-regularization

NASA Astrophysics Data System (ADS)

Zhang, Tianzhen; Wang, Xiumei; Gao, Xinbo

2018-04-01

Nowadays, several datasets are demonstrated by multi-view, which usually include shared and complementary information. Multi-view clustering methods integrate the information of multi-view to obtain better clustering results. Nonnegative matrix factorization has become an essential and popular tool in clustering methods because of its interpretation. However, existing nonnegative matrix factorization based multi-view clustering algorithms do not consider the disagreement between views and neglects the fact that different views will have different contributions to the data distribution. In this paper, we propose a new multi-view clustering method, named adaptive multi-view clustering based on nonnegative matrix factorization and pairwise co-regularization. The proposed algorithm can obtain the parts-based representation of multi-view data by nonnegative matrix factorization. Then, pairwise co-regularization is used to measure the disagreement between views. There is only one parameter to auto learning the weight values according to the contribution of each view to data distribution. Experimental results show that the proposed algorithm outperforms several state-of-the-arts algorithms for multi-view clustering.

Thermal Entanglement Between Atoms in the Four-Cavity Linear Chain Coupled by Single-Mode Fibers

NASA Astrophysics Data System (ADS)

Wang, Jun-Biao; Zhang, Guo-Feng

2018-05-01

Natural thermal entanglement between atoms of a linear arranged four coupled cavities system is studied. The results show that there is no thermal pairwise entanglement between atoms if atom-field interaction strength f or fiber-cavity coupling constant J equals to zero, both f and J can induce thermal pairwise entanglement in a certain range. Numerical simulations show that the nearest neighbor concurrence C A B is always greater than alternate concurrence C A C in the same condition. In addition, the effect of temperature T on the entanglement of alternate qubits is much stronger than the nearest neighbor qubits.

Emergent behaviors of the Schrödinger-Lohe model on cooperative-competitive networks

NASA Astrophysics Data System (ADS)

Huh, Hyungjin; Ha, Seung-Yeal; Kim, Dohyun

2017-12-01

We present several sufficient frameworks leading to the emergent behaviors of the coupled Schrödinger-Lohe (S-L) model under the same one-body external potential on cooperative-competitive networks. The S-L model was first introduced as a possible phenomenological model exhibiting quantum synchronization and its emergent dynamics on all-to-all cooperative networks has been treated via two distinct approaches, Lyapunov functional approach and the finite-dimensional reduction based on pairwise correlations. In this paper, we further generalize the finite-dimensional dynamical systems approach for pairwise correlation functions on cooperative-competitive networks and provide several sufficient frameworks leading to the collective exponential synchronization. For small systems consisting of three and four quantum subsystem, we also show that the system for pairwise correlations can be reduced to the Lotka-Volterra model with cooperative and competitive interactions, in which lots of interesting dynamical patterns appear, e.g., existence of closed orbits and limit-cycles.

Trophic interaction modifications: an empirical and theoretical framework.

PubMed

Terry, J Christopher D; Morris, Rebecca J; Bonsall, Michael B

2017-10-01

Consumer-resource interactions are often influenced by other species in the community. At present these 'trophic interaction modifications' are rarely included in ecological models despite demonstrations that they can drive system dynamics. Here, we advocate and extend an approach that has the potential to unite and represent this key group of non-trophic interactions by emphasising the change to trophic interactions induced by modifying species. We highlight the opportunities this approach brings in comparison to frameworks that coerce trophic interaction modifications into pairwise relationships. To establish common frames of reference and explore the value of the approach, we set out a range of metrics for the 'strength' of an interaction modification which incorporate increasing levels of contextual information about the system. Through demonstrations in three-species model systems, we establish that these metrics capture complimentary aspects of interaction modifications. We show how the approach can be used in a range of empirical contexts; we identify as specific gaps in current understanding experiments with multiple levels of modifier species and the distributions of modifications in networks. The trophic interaction modification approach we propose can motivate and unite empirical and theoretical studies of system dynamics, providing a route to confront ecological complexity. © 2017 The Authors. Ecology Letters published by CNRS and John Wiley & Sons Ltd.

A Small World of Neuronal Synchrony

PubMed Central

Yu, Shan; Huang, Debin; Singer, Wolf

2008-01-01

A small-world network has been suggested to be an efficient solution for achieving both modular and global processing—a property highly desirable for brain computations. Here, we investigated functional networks of cortical neurons using correlation analysis to identify functional connectivity. To reconstruct the interaction network, we applied the Ising model based on the principle of maximum entropy. This allowed us to assess the interactions by measuring pairwise correlations and to assess the strength of coupling from the degree of synchrony. Visual responses were recorded in visual cortex of anesthetized cats, simultaneously from up to 24 neurons. First, pairwise correlations captured most of the patterns in the population's activity and, therefore, provided a reliable basis for the reconstruction of the interaction networks. Second, and most importantly, the resulting networks had small-world properties; the average path lengths were as short as in simulated random networks, but the clustering coefficients were larger. Neurons differed considerably with respect to the number and strength of interactions, suggesting the existence of “hubs” in the network. Notably, there was no evidence for scale-free properties. These results suggest that cortical networks are optimized for the coexistence of local and global computations: feature detection and feature integration or binding. PMID:18400792

Σ 3 (111 ) grain boundary of body-centered cubic Ti-Mo and Ti-V alloys: First-principles and model calculations

NASA Astrophysics Data System (ADS)

Yan, Jia-Yi; Ehteshami, Hossein; Korzhavyi, Pavel A.; Borgenstam, Annika

2017-07-01

The energetics and atomic structures of Σ 3 [1 1 ¯0 ] (111 ) grain boundary (GB) of body-centered cubic (bcc) Ti-Mo and Ti-V alloys are investigated using density-functional-theory calculations and virtual crystal approximation. The electron density in bcc structure and the atomic displacements and excess energy of the GB are correlated to bcc-ω phase stability. Model calculations based on pairwise interplanar interactions successfully reproduce the chemical part of GB energy. The chemical GB energy can be expressed as a sum of excess pairwise interactions between bcc (111) layers, which are obtained from Gaussian elimination of the total energies of a number of periodic structures. The energy associated with the relaxation near the GB is solved by numerical minimization using the derivatives of the excess interactions. Anharmonic interlayer interactions are necessary for obtaining accurate relaxation energy and excess GB volume from model calculations. The effect of GB on vibrational spectrum is also investigated. Segregation energies of B and Y to a substitutional site on the GB plane are calculated. Preliminary results suggest that Y tends to segregate, while B tends to antisegregate.

Effective empirical corrections for basis set superposition error in the def2-SVPD basis: gCP and DFT-C

NASA Astrophysics Data System (ADS)

Witte, Jonathon; Neaton, Jeffrey B.; Head-Gordon, Martin

2017-06-01

With the aim of mitigating the basis set error in density functional theory (DFT) calculations employing local basis sets, we herein develop two empirical corrections for basis set superposition error (BSSE) in the def2-SVPD basis, a basis which—when stripped of BSSE—is capable of providing near-complete-basis DFT results for non-covalent interactions. Specifically, we adapt the existing pairwise geometrical counterpoise (gCP) approach to the def2-SVPD basis, and we develop a beyond-pairwise approach, DFT-C, which we parameterize across a small set of intermolecular interactions. Both gCP and DFT-C are evaluated against the traditional Boys-Bernardi counterpoise correction across a set of 3402 non-covalent binding energies and isomerization energies. We find that the DFT-C method represents a significant improvement over gCP, particularly for non-covalently-interacting molecular clusters. Moreover, DFT-C is transferable among density functionals and can be combined with existing functionals—such as B97M-V—to recover large-basis results at a fraction of the cost.

An efficient algorithm for pairwise local alignment of protein interaction networks

DOE PAGES

Chen, Wenbin; Schmidt, Matthew; Tian, Wenhong; ...

2015-04-01

Recently, researchers seeking to understand, modify, and create beneficial traits in organisms have looked for evolutionarily conserved patterns of protein interactions. Their conservation likely means that the proteins of these conserved functional modules are important to the trait's expression. In this paper, we formulate the problem of identifying these conserved patterns as a graph optimization problem, and develop a fast heuristic algorithm for this problem. We compare the performance of our network alignment algorithm to that of the MaWISh algorithm [Koyuturk M, Kim Y, Topkara U, Subramaniam S, Szpankowski W, Grama A, Pairwise alignment of protein interaction networks, J Computmore » Biol 13(2): 182-199, 2006.], which bases its search algorithm on a related decision problem formulation. We find that our algorithm discovers conserved modules with a larger number of proteins in an order of magnitude less time. In conclusion, the protein sets found by our algorithm correspond to known conserved functional modules at comparable precision and recall rates as those produced by the MaWISh algorithm.« less

The structure of pairwise correlation in mouse primary visual cortex reveals functional organization in the absence of an orientation map.

PubMed

Denman, Daniel J; Contreras, Diego

2014-10-01

Neural responses to sensory stimuli are not independent. Pairwise correlation can reduce coding efficiency, occur independent of stimulus representation, or serve as an additional channel of information, depending on the timescale of correlation and the method of decoding. Any role for correlation depends on its magnitude and structure. In sensory areas with maps, like the orientation map in primary visual cortex (V1), correlation is strongly related to the underlying functional architecture, but it is unclear whether this correlation structure is an essential feature of the system or arises from the arrangement of cells in the map. We assessed the relationship between functional architecture and pairwise correlation by measuring both synchrony and correlated spike count variability in mouse V1, which lacks an orientation map. We observed significant pairwise synchrony, which was organized by distance and relative orientation preference between cells. We also observed nonzero correlated variability in both the anesthetized (0.16) and awake states (0.18). Our results indicate that the structure of pairwise correlation is maintained in the absence of an underlying anatomical organization and may be an organizing principle of the mammalian visual system preserved by nonrandom connectivity within local networks. © The Author 2013. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.

Kernel machine methods for integrative analysis of genome-wide methylation and genotyping studies.

PubMed

Zhao, Ni; Zhan, Xiang; Huang, Yen-Tsung; Almli, Lynn M; Smith, Alicia; Epstein, Michael P; Conneely, Karen; Wu, Michael C

2018-03-01

Many large GWAS consortia are expanding to simultaneously examine the joint role of DNA methylation in addition to genotype in the same subjects. However, integrating information from both data types is challenging. In this paper, we propose a composite kernel machine regression model to test the joint epigenetic and genetic effect. Our approach works at the gene level, which allows for a common unit of analysis across different data types. The model compares the pairwise similarities in the phenotype to the pairwise similarities in the genotype and methylation values; and high correspondence is suggestive of association. A composite kernel is constructed to measure the similarities in the genotype and methylation values between pairs of samples. We demonstrate through simulations and real data applications that the proposed approach can correctly control type I error, and is more robust and powerful than using only the genotype or methylation data in detecting trait-associated genes. We applied our method to investigate the genetic and epigenetic regulation of gene expression in response to stressful life events using data that are collected from the Grady Trauma Project. Within the kernel machine testing framework, our methods allow for heterogeneity in effect sizes, nonlinear, and interactive effects, as well as rapid P-value computation. © 2017 WILEY PERIODICALS, INC.

Using Replicates in Information Retrieval Evaluation

PubMed Central

VOORHEES, ELLEN M.; SAMAROV, DANIEL; SOBOROFF, IAN

2018-01-01

This article explores a method for more accurately estimating the main effect of the system in a typical test-collection-based evaluation of information retrieval systems, thus increasing the sensitivity of system comparisons. Randomly partitioning the test document collection allows for multiple tests of a given system and topic (replicates). Bootstrap ANOVA can use these replicates to extract system-topic interactions—something not possible without replicates—yielding a more precise value for the system effect and a narrower confidence interval around that value. Experiments using multiple TREC collections demonstrate that removing the topic-system interactions substantially reduces the confidence intervals around the system effect as well as increases the number of significant pairwise differences found. Further, the method is robust against small changes in the number of partitions used, against variability in the documents that constitute the partitions, and the measure of effectiveness used to quantify system effectiveness. PMID:29905334

Modelling particles moving in a potential field with pairwise interactions and an application

Treesearch

D. R. Brillinger; Haiganoush Preisler; M. J. Wisdom

2011-01-01

Motions of particles in fields characterized by real-valued potential functions, are considered. Three particular expressions for potential functions are studied. One, U, depends on the ith particle’s location, ri(t) at times t

Learning a constrained conditional random field for enhanced segmentation of fallen trees in ALS point clouds

NASA Astrophysics Data System (ADS)

Polewski, Przemyslaw; Yao, Wei; Heurich, Marco; Krzystek, Peter; Stilla, Uwe

2018-06-01

In this study, we present a method for improving the quality of automatic single fallen tree stem segmentation in ALS data by applying a specialized constrained conditional random field (CRF). The entire processing pipeline is composed of two steps. First, short stem segments of equal length are detected and a subset of them is selected for further processing, while in the second step the chosen segments are merged to form entire trees. The first step is accomplished using the specialized CRF defined on the space of segment labelings, capable of finding segment candidates which are easier to merge subsequently. To achieve this, the CRF considers not only the features of every candidate individually, but incorporates pairwise spatial interactions between adjacent segments into the model. In particular, pairwise interactions include a collinearity/angular deviation probability which is learned from training data as well as the ratio of spatial overlap, whereas unary potentials encode a learned probabilistic model of the laser point distribution around each segment. Each of these components enters the CRF energy with its own balance factor. To process previously unseen data, we first calculate the subset of segments for merging on a grid of balance factors by minimizing the CRF energy. Then, we perform the merging and rank the balance configurations according to the quality of their resulting merged trees, obtained from a learned tree appearance model. The final result is derived from the top-ranked configuration. We tested our approach on 5 plots from the Bavarian Forest National Park using reference data acquired in a field inventory. Compared to our previous segment selection method without pairwise interactions, an increase in detection correctness and completeness of up to 7 and 9 percentage points, respectively, was observed.

Pairwise Trajectory Management (PTM): Concept Overview

NASA Technical Reports Server (NTRS)

Jones, Kenneth M.; Graff, Thomas J.; Chartrand, Ryan C.; Carreno, Victor; Kibler, Jennifer L.

2017-01-01

Pairwise Trajectory Management (PTM) is an Interval Management (IM) concept that utilizes airborne and ground-based capabilities to enable the implementation of airborne pairwise spacing capabilities in oceanic regions. The goal of PTM is to use airborne surveillance and tools to manage an "at or greater than" inter-aircraft spacing. Due to the precision of Automatic Dependent Surveillance-Broadcast (ADS-B) information and the use of airborne spacing guidance, the PTM minimum spacing distance will be less than distances a controller can support with current automation systems that support oceanic operations. Ground tools assist the controller in evaluating the traffic picture and determining appropriate PTM clearances to be issued. Avionics systems provide guidance information that allows the flight crew to conform to the PTM clearance issued by the controller. The combination of a reduced minimum distance and airborne spacing management will increase the capacity and efficiency of aircraft operations at a given altitude or volume of airspace. This paper provides an overview of the proposed application, description of a few key scenarios, high level discussion of expected air and ground equipment and procedure changes, overview of a potential flight crew human-machine interface that would support PTM operations and some initial PTM benefits results.

POEM: Identifying Joint Additive Effects on Regulatory Circuits.

PubMed

Botzman, Maya; Nachshon, Aharon; Brodt, Avital; Gat-Viks, Irit

2016-01-01

Expression Quantitative Trait Locus (eQTL) mapping tackles the problem of identifying variation in DNA sequence that have an effect on the transcriptional regulatory network. Major computational efforts are aimed at characterizing the joint effects of several eQTLs acting in concert to govern the expression of the same genes. Yet, progress toward a comprehensive prediction of such joint effects is limited. For example, existing eQTL methods commonly discover interacting loci affecting the expression levels of a module of co-regulated genes. Such "modularization" approaches, however, are focused on epistatic relations and thus have limited utility for the case of additive (non-epistatic) effects. Here we present POEM (Pairwise effect On Expression Modules), a methodology for identifying pairwise eQTL effects on gene modules. POEM is specifically designed to achieve high performance in the case of additive joint effects. We applied POEM to transcription profiles measured in bone marrow-derived dendritic cells across a population of genotyped mice. Our study reveals widespread additive, trans-acting pairwise effects on gene modules, characterizes their organizational principles, and highlights high-order interconnections between modules within the immune signaling network. These analyses elucidate the central role of additive pairwise effect in regulatory circuits, and provide computational tools for future investigations into the interplay between eQTLs. The software described in this article is available at csgi.tau.ac.il/POEM/.

POEM: Identifying Joint Additive Effects on Regulatory Circuits

PubMed Central

Botzman, Maya; Nachshon, Aharon; Brodt, Avital; Gat-Viks, Irit

2016-01-01

Motivation: Expression Quantitative Trait Locus (eQTL) mapping tackles the problem of identifying variation in DNA sequence that have an effect on the transcriptional regulatory network. Major computational efforts are aimed at characterizing the joint effects of several eQTLs acting in concert to govern the expression of the same genes. Yet, progress toward a comprehensive prediction of such joint effects is limited. For example, existing eQTL methods commonly discover interacting loci affecting the expression levels of a module of co-regulated genes. Such “modularization” approaches, however, are focused on epistatic relations and thus have limited utility for the case of additive (non-epistatic) effects. Results: Here we present POEM (Pairwise effect On Expression Modules), a methodology for identifying pairwise eQTL effects on gene modules. POEM is specifically designed to achieve high performance in the case of additive joint effects. We applied POEM to transcription profiles measured in bone marrow-derived dendritic cells across a population of genotyped mice. Our study reveals widespread additive, trans-acting pairwise effects on gene modules, characterizes their organizational principles, and highlights high-order interconnections between modules within the immune signaling network. These analyses elucidate the central role of additive pairwise effect in regulatory circuits, and provide computational tools for future investigations into the interplay between eQTLs. Availability: The software described in this article is available at csgi.tau.ac.il/POEM/. PMID:27148351

Dynamic density functional theory with hydrodynamic interactions and fluctuations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Donev, Aleksandar, E-mail: donev@courant.nyu.edu; Vanden-Eijnden, Eric, E-mail: eve2@courant.nyu.edu

2014-06-21

We derive a closed equation for the empirical concentration of colloidal particles in the presence of both hydrodynamic and direct interactions. The ensemble average of our functional Langevin equation reproduces known deterministic Dynamic Density Functional Theory (DDFT) [M. Rex and H. Löwen, “Dynamical density functional theory with hydrodynamic interactions and colloids in unstable traps,” Phys. Rev. Lett. 101(14), 148302 (2008)], and, at the same time, it also describes the microscopic fluctuations around the mean behavior. We suggest separating the ideal (non-interacting) contribution from additional corrections due to pairwise interactions. We find that, for an incompressible fluid and in the absencemore » of direct interactions, the mean concentration follows Fick's law just as for uncorrelated walkers. At the same time, the nature of the stochastic terms in fluctuating DDFT is shown to be distinctly different for hydrodynamically-correlated and uncorrelated walkers. This leads to striking differences in the behavior of the fluctuations around Fick's law, even in the absence of pairwise interactions. We connect our own prior work [A. Donev, T. G. Fai, and E. Vanden-Eijnden, “A reversible mesoscopic model of diffusion in liquids: from giant fluctuations to Fick's law,” J. Stat. Mech.: Theory Exp. (2014) P04004] on fluctuating hydrodynamics of diffusion in liquids to the DDFT literature, and demonstrate that the fluid cannot easily be eliminated from consideration if one wants to describe the collective diffusion in colloidal suspensions.« less

SANA NetGO: a combinatorial approach to using Gene Ontology (GO) terms to score network alignments.

PubMed

Hayes, Wayne B; Mamano, Nil

2018-04-15

Gene Ontology (GO) terms are frequently used to score alignments between protein-protein interaction (PPI) networks. Methods exist to measure GO similarity between proteins in isolation, but proteins in a network alignment are not isolated: each pairing is dependent on every other via the alignment itself. Existing measures fail to take into account the frequency of GO terms across networks, instead imposing arbitrary rules on when to allow GO terms. Here we develop NetGO, a new measure that naturally weighs infrequent, informative GO terms more heavily than frequent, less informative GO terms, without arbitrary cutoffs, instead downweighting GO terms according to their frequency in the networks being aligned. This is a global measure applicable only to alignments, independent of pairwise GO measures, in the same sense that the edge-based EC or S3 scores are global measures of topological similarity independent of pairwise topological similarities. We demonstrate the superiority of NetGO in alignments of predetermined quality and show that NetGO correlates with alignment quality better than any existing GO-based alignment measures. We also demonstrate that NetGO provides a measure of taxonomic similarity between species, consistent with existing taxonomic measuresa feature not shared with existing GObased network alignment measures. Finally, we re-score alignments produced by almost a dozen aligners from a previous study and show that NetGO does a better job at separating good alignments from bad ones. Available as part of SANA. whayes@uci.edu. Supplementary data are available at Bioinformatics online.

Depth Reconstruction from Single Images Using a Convolutional Neural Network and a Condition Random Field Model.

PubMed

Liu, Dan; Liu, Xuejun; Wu, Yiguang

2018-04-24

This paper presents an effective approach for depth reconstruction from a single image through the incorporation of semantic information and local details from the image. A unified framework for depth acquisition is constructed by joining a deep Convolutional Neural Network (CNN) and a continuous pairwise Conditional Random Field (CRF) model. Semantic information and relative depth trends of local regions inside the image are integrated into the framework. A deep CNN network is firstly used to automatically learn a hierarchical feature representation of the image. To get more local details in the image, the relative depth trends of local regions are incorporated into the network. Combined with semantic information of the image, a continuous pairwise CRF is then established and is used as the loss function of the unified model. Experiments on real scenes demonstrate that the proposed approach is effective and that the approach obtains satisfactory results.

Evolutionary games in the multiverse.

PubMed

Gokhale, Chaitanya S; Traulsen, Arne

2010-03-23

Evolutionary game dynamics of two players with two strategies has been studied in great detail. These games have been used to model many biologically relevant scenarios, ranging from social dilemmas in mammals to microbial diversity. Some of these games may, in fact, take place between a number of individuals and not just between two. Here we address one-shot games with multiple players. As long as we have only two strategies, many results from two-player games can be generalized to multiple players. For games with multiple players and more than two strategies, we show that statements derived for pairwise interactions no longer hold. For two-player games with any number of strategies there can be at most one isolated internal equilibrium. For any number of players with any number of strategies , there can be at most isolated internal equilibria. Multiplayer games show a great dynamical complexity that cannot be captured based on pairwise interactions. Our results hold for any game and can easily be applied to specific cases, such as public goods games or multiplayer stag hunts.

Graph Curvature for Differentiating Cancer Networks

PubMed Central

Sandhu, Romeil; Georgiou, Tryphon; Reznik, Ed; Zhu, Liangjia; Kolesov, Ivan; Senbabaoglu, Yasin; Tannenbaum, Allen

2015-01-01

Cellular interactions can be modeled as complex dynamical systems represented by weighted graphs. The functionality of such networks, including measures of robustness, reliability, performance, and efficiency, are intrinsically tied to the topology and geometry of the underlying graph. Utilizing recently proposed geometric notions of curvature on weighted graphs, we investigate the features of gene co-expression networks derived from large-scale genomic studies of cancer. We find that the curvature of these networks reliably distinguishes between cancer and normal samples, with cancer networks exhibiting higher curvature than their normal counterparts. We establish a quantitative relationship between our findings and prior investigations of network entropy. Furthermore, we demonstrate how our approach yields additional, non-trivial pair-wise (i.e. gene-gene) interactions which may be disrupted in cancer samples. The mathematical formulation of our approach yields an exact solution to calculating pair-wise changes in curvature which was computationally infeasible using prior methods. As such, our findings lay the foundation for an analytical approach to studying complex biological networks. PMID:26169480

Statistical Mechanics of the US Supreme Court

NASA Astrophysics Data System (ADS)

Lee, Edward D.; Broedersz, Chase P.; Bialek, William

2015-07-01

We build simple models for the distribution of voting patterns in a group, using the Supreme Court of the United States as an example. The maximum entropy model consistent with the observed pairwise correlations among justices' votes, an Ising spin glass, agrees quantitatively with the data. While all correlations (perhaps surprisingly) are positive, the effective pairwise interactions in the spin glass model have both signs, recovering the intuition that ideologically opposite justices negatively influence each another. Despite the competing interactions, a strong tendency toward unanimity emerges from the model, organizing the voting patterns in a relatively simple "energy landscape." Besides unanimity, other energy minima in this landscape, or maxima in probability, correspond to prototypical voting states, such as the ideological split or a tightly correlated, conservative core. The model correctly predicts the correlation of justices with the majority and gives us a measure of their influence on the majority decision. These results suggest that simple models, grounded in statistical physics, can capture essential features of collective decision making quantitatively, even in a complex political context.

Microstructural Dynamics and Rheology of Suspensions of Rigid Fibers

NASA Astrophysics Data System (ADS)

Butler, Jason E.; Snook, Braden

2018-01-01

The dynamics and rheology of suspensions of rigid, non-Brownian fibers in Newtonian fluids are reviewed. Experiments, theories, and computer simulations are considered, with an emphasis on suspensions at semidilute and concentrated conditions. In these suspensions, interactions between the particles strongly influence the microstructure and rheological properties of the suspension. The interactions can arise from hydrodynamic disturbances, giving multibody interactions at long ranges and pairwise lubrication forces over short distances. For concentrated suspensions, additional interactions due to excluded volume (contacts) and adhesive forces are addressed. The relative importance of the various interactions as a function of fiber concentration is assessed.

Relative importance of evolutionary dynamics depends on the composition of microbial predator-prey community.

PubMed

Friman, Ville-Petri; Dupont, Alessandra; Bass, David; Murrell, David J; Bell, Thomas

2016-06-01

Community dynamics are often studied in subsets of pairwise interactions. Scaling pairwise interactions back to the community level is, however, problematic because one given interaction might not reflect ecological and evolutionary outcomes of other functionally similar species interactions or capture the emergent eco-evolutionary dynamics arising only in more complex communities. Here we studied this experimentally by exposing Pseudomonas fluorescens SBW25 prey bacterium to four different protist predators (Tetrahymena pyriformis, Tetrahymena vorax, Chilomonas paramecium and Acanthamoeba polyphaga) in all possible single-predator, two-predator and four-predator communities for hundreds of prey generations covering both ecological and evolutionary timescales. We found that only T. pyriformis selected for prey defence in single-predator communities. Although T. pyriformis selection was constrained in the presence of the intraguild predator, T. vorax, T. pyriformis selection led to evolution of specialised prey defence strategies in the presence of C. paramecium or A. polyphaga. At the ecological level, adapted prey populations were phenotypically more diverse, less stable and less productive compared with non-adapted prey populations. These results suggest that predator community composition affects the relative importance of ecological and evolutionary processes and can crucially determine when rapid evolution has the potential to change ecological properties of microbial communities.

Relative importance of evolutionary dynamics depends on the composition of microbial predator–prey community

PubMed Central

Friman, Ville-Petri; Dupont, Alessandra; Bass, David; Murrell, David J; Bell, Thomas

2016-01-01

Community dynamics are often studied in subsets of pairwise interactions. Scaling pairwise interactions back to the community level is, however, problematic because one given interaction might not reflect ecological and evolutionary outcomes of other functionally similar species interactions or capture the emergent eco-evolutionary dynamics arising only in more complex communities. Here we studied this experimentally by exposing Pseudomonas fluorescens SBW25 prey bacterium to four different protist predators (Tetrahymena pyriformis, Tetrahymena vorax, Chilomonas paramecium and Acanthamoeba polyphaga) in all possible single-predator, two-predator and four-predator communities for hundreds of prey generations covering both ecological and evolutionary timescales. We found that only T. pyriformis selected for prey defence in single-predator communities. Although T. pyriformis selection was constrained in the presence of the intraguild predator, T. vorax, T. pyriformis selection led to evolution of specialised prey defence strategies in the presence of C. paramecium or A. polyphaga. At the ecological level, adapted prey populations were phenotypically more diverse, less stable and less productive compared with non-adapted prey populations. These results suggest that predator community composition affects the relative importance of ecological and evolutionary processes and can crucially determine when rapid evolution has the potential to change ecological properties of microbial communities. PMID:26684728

Spatial heterogeneity of plant-soil feedback affects root interactions and interspecific competition.

PubMed

Hendriks, Marloes; Ravenek, Janneke M; Smit-Tiekstra, Annemiek E; van der Paauw, Jan Willem; de Caluwe, Hannie; van der Putten, Wim H; de Kroon, Hans; Mommer, Liesje

2015-08-01

Plant-soil feedback is receiving increasing interest as a factor influencing plant competition and species coexistence in grasslands. However, we do not know how spatial distribution of plant-soil feedback affects plant below-ground interactions. We investigated the way in which spatial heterogeneity of soil biota affects competitive interactions in grassland plant species. We performed a pairwise competition experiment combined with heterogeneous distribution of soil biota using four grassland plant species and their soil biota. Patches were applied as quadrants of 'own' and 'foreign' soils from all plant species in all pairwise combinations. To evaluate interspecific root responses, species-specific root biomass was quantified using real-time PCR. All plant species suffered negative soil feedback, but strength was species-specific, reflected by a decrease in root growth in own compared with foreign soil. Reduction in root growth in own patches by the superior plant competitor provided opportunities for inferior competitors to increase root biomass in these patches. These patterns did not cascade into above-ground effects during our experiment. We show that root distributions can be determined by spatial heterogeneity of soil biota, affecting plant below-ground competitive interactions. Thus, spatial heterogeneity of soil biota may contribute to plant species coexistence in species-rich grasslands. © 2015 The Authors. New Phytologist © 2015 New Phytologist Trust.

Epiregulin (EREG) and human V-ATPase (TCIRG1): genetic variation, ethnicity and pulmonary tuberculosis susceptibility in Guinea-Bissau and The Gambia

PubMed Central

White, Marquitta J.; Tacconelli, Alessandra; Chen, Jane S.; Wejse, Christian; Hill, Philip C.; Gomez, Victor F; Velez-Edwards, Digna R.; Østergaard, Lars J.; Hu, Ting; Moore, Jason H.; Novelli, Giuseppe; Scott, William K.; Williams, Scott M.; Sirugo, Giorgio

2017-01-01

We analyzed two West African samples (Guinea-Bissau: n = 289 cases, 322 controls; The Gambia: n = 240 cases, 248 controls) to evaluate single nucleotide polymorphisms (SNPs) in Epiregulin (EREG) and V-ATPase (T cell immune regulator 1, TCIRG1) using single and multi-locus analyses to determine whether previously described associations with pulmonary tuberculosis (PTB) in Vietnamese and Italians would replicate in African populations. We did not detect any significant single locus or haplotype associations in either sample. We also performed exploratory pairwise interaction analyses using Visualization of Statistical Epistasis Networks (ViSEN), a novel method to detect only interactions among multiple variables, to elucidate possible interaction effects between SNPs and demographic factors. Although we found no strong evidence of marginal effects, there were several significant pairwise interactions that were identified in either the Guinea-Bissau or The Gambia samples, two of which replicated across populations. Our results indicate that the effects of EREG and TCIRG1 variants on PTB susceptibility, to the extent that they exist, are dependent on gene-gene interactions in West African populations as detected with ViSEN. In addition, epistatic effects are likely to be influenced by inter- and intra-population differences in genetic or environmental context and/or the mycobacterial lineages causing disease. PMID:24898387

Polymer chain collapse induced by many-body dipole correlations.

PubMed

Budkov, Yu A; Kalikin, N N; Kolesnikov, A L

2017-04-01

We present a simple analytical theory of a flexible polymer chain dissolved in a good solvent, carrying permanent freely oriented dipoles on the monomers. We take into account the dipole correlations within the random phase approximation (RPA), as well as a dielectric heterogeneity in the internal polymer volume relative to the bulk solution. We demonstrate that the dipole correlations of monomers can be taken into account as pairwise ones only when the polymer chain is in a coil conformation. In this case the dipole correlations manifest themselves through the Keesom interactions of the permanent dipoles. On the other hand, the dielectric heterogeneity effect (dielectric mismatch effect) leads to the effective interaction between the monomers of the polymeric coil. Both of these effects can be taken into account by renormalizing the second virial coefficient of the monomer-monomer volume interactions. We establish that in the case when the solvent dielectric permittivity exceeds the dielectric permittivity of the polymeric material, the dielectric mismatch effect competes with the dipole attractive interactions, leading to polymer coil expansion. In the opposite case, both the dielectric mismatch effect and the dipole attractive interaction lead to the polymer coil collapse. We analyse the coil-globule transition caused by the dipole correlations of monomers within the many-body theory. We demonstrate that accounting for the dipole correlations higher than the pairwise ones smooths this pure electrostatics driven coil-globule transition of the polymer chain.

Interactions among filamentous fungi Aspergillus niger, Fusarium verticillioides and Clonostachys rosea: fungal biomass, diversity of secreted metabolites and fumonisin production.

PubMed

Chatterjee, Subhankar; Kuang, Yi; Splivallo, Richard; Chatterjee, Paramita; Karlovsky, Petr

2016-05-10

Interactions among fungi colonizing dead organic matter involve exploitation competition and interference competition. Major mechanism of interference competition is antibiosis caused by secreted secondary metabolites. The effect of competition on secondary metabolite production by fungi is however poorly understood. Fungal biomass was rarely monitored in interaction studies; it is not known whether dominance in pairwise interactions follows congruent patterns. Pairwise interactions of three fungal species with different life styles were studied. The saprophyte Aspergillus niger (A.n.), the plant pathogen Fusarium verticillioides (F.v.), and the mycoparasite Clonostachys rosea (C.r.) were grown in single and dual cultures in minimal medium with asparagine as nitrogen source. Competitive fitness shifted with time: in dual C.r./F.v. cultures after 10 d F.v. grew well while C.r. was suppressed; after 20 d C.r. recovered while F.v. became suppressed; and after 30 d most F.v. was destroyed. At certain time points fungal competitive fitness exhibited a rock-paper-scissors pattern: F.v. > A.n., A.n. > C.r., and C.r. > F.v. Most metabolites secreted to the medium at early stages in single and dual cultures were not found at later times. Many metabolites occurring in supernatants of single cultures were suppressed in dual cultures and many new metabolites not occurring in single cultures were found in dual cultures. A. niger showed the greatest ability to suppress the accumulation of metabolites produced by the other fungi. A. niger was also the species with the largest capacity of transforming metabolites produced by other fungi. Fumonisin production by F. verticillioides was suppressed in co-cultures with C. rosea but fumonisin B1 was not degraded by C. rosea nor did it affect the growth of C. rosea up to a concentration of 160 μg/ml. Competitive fitness in pairwise interactions among fungi is incongruent, indicating that species-specific factors and/or effects are involved. Many metabolites secreted by fungi are catabolized by their producers at later growth stages. Diversity of metabolites accumulating in the medium is stimulated by fungus/fungus interactions. C. rosea suppresses the synthesis of fumonisins by F. verticillioides but does not degrade fumonisins.

Potential impact of soil microbiomes on the leaf metabolome and on herbivore feeding behavior

USDA-ARS?s Scientific Manuscript database

: It is known that environmental factors can affect the biosynthesis of leaf metabolites. Similarly, specific pairwise plant-microbe interactions modulate specifically the plant’s metabolome by stimulating production of phytoalexins and other defense-related compounds. However, there is no informati...

Specific Non-Local Interactions Are Not Necessary for Recovering Native Protein Dynamics

PubMed Central

Dasgupta, Bhaskar; Kasahara, Kota; Kamiya, Narutoshi; Nakamura, Haruki; Kinjo, Akira R.

2014-01-01

The elastic network model (ENM) is a widely used method to study native protein dynamics by normal mode analysis (NMA). In ENM we need information about all pairwise distances, and the distance between contacting atoms is restrained to the native value. Therefore ENM requires O(N2) information to realize its dynamics for a protein consisting of N amino acid residues. To see if (or to what extent) such a large amount of specific structural information is required to realize native protein dynamics, here we introduce a novel model based on only O(N) restraints. This model, named the ‘contact number diffusion’ model (CND), includes specific distance restraints for only local (along the amino acid sequence) atom pairs, and semi-specific non-local restraints imposed on each atom, rather than atom pairs. The semi-specific non-local restraints are defined in terms of the non-local contact numbers of atoms. The CND model exhibits the dynamic characteristics comparable to ENM and more correlated with the explicit-solvent molecular dynamics simulation than ENM. Moreover, unrealistic surface fluctuations often observed in ENM were suppressed in CND. On the other hand, in some ligand-bound structures CND showed larger fluctuations of buried protein atoms interacting with the ligand compared to ENM. In addition, fluctuations from CND and ENM show comparable correlations with the experimental B-factor. Although there are some indications of the importance of some specific non-local interactions, the semi-specific non-local interactions are mostly sufficient for reproducing the native protein dynamics. PMID:24625758

Reexamination of the interaction of atoms with a LiF(001) surface

NASA Astrophysics Data System (ADS)

Miraglia, J. E.; Gravielle, M. S.

2017-02-01

Pairwise additive potentials for multielectronic atoms interacting with a LiF(001) surface are revisited by including an improved description of the electron density associated with the different lattice sites, as well as nonlocal electron density contributions. Within this model, the electron distribution around each ionic site of the crystal is described by means of a so-called "onion" approach that accounts for the influence of the Madelung potential. From such densities, binary interatomic potentials are then derived by using well-known nonlocal functionals. Rumpling and long-range contributions due to projectile polarization and van der Waals forces are also included. We apply this pairwise additive approximation to evaluate the interaction potential between closed-shell (He, Ne, Ar, Kr, and Xe) and open-shell (N, S, and Cl) atoms and the LiF surface, analyzing the relative importance of the different contributions. The performance of the proposed potentials is assessed by contrasting angular positions of rainbow and supernumerary rainbow maxima produced by fast grazing incidence with available experimental data. One important result of our model is that both van der Waals contributions and thermal lattice vibrations play a negligible role for normal energies in the eV range.

A density functional approach to ferrogels

NASA Astrophysics Data System (ADS)

Cremer, P.; Heinen, M.; Menzel, A. M.; Löwen, H.

2017-07-01

Ferrogels consist of magnetic colloidal particles embedded in an elastic polymer matrix. As a consequence, their structural and rheological properties are governed by a competition between magnetic particle-particle interactions and mechanical matrix elasticity. Typically, the particles are permanently fixed within the matrix, which makes them distinguishable by their positions. Over time, particle neighbors do not change due to the fixation by the matrix. Here we present a classical density functional approach for such ferrogels. We map the elastic matrix-induced interactions between neighboring colloidal particles distinguishable by their positions onto effective pairwise interactions between indistinguishable particles similar to a ‘pairwise pseudopotential’. Using Monte-Carlo computer simulations, we demonstrate for one-dimensional dipole-spring models of ferrogels that this mapping is justified. We then use the pseudopotential as an input into classical density functional theory of inhomogeneous fluids and predict the bulk elastic modulus of the ferrogel under various conditions. In addition, we propose the use of an ‘external pseudopotential’ when one switches from the viewpoint of a one-dimensional dipole-spring object to a one-dimensional chain embedded in an infinitely extended bulk matrix. Our mapping approach paves the way to describe various inhomogeneous situations of ferrogels using classical density functional concepts of inhomogeneous fluids.

Information theoretic measures of network coordination in high-frequency scalp EEG reveal dynamic patterns associated with seizure termination

PubMed Central

Stamoulis, Catherine; Schomer, Donald L.; Chang, Bernard S.

2013-01-01

How a seizure terminates is still under-studied and, despite its clinical importance, remains an obscure phase of seizure evolution. Recent studies of seizure-related scalp EEGs at frequencies >100 Hz suggest that neural activity, in the form of oscillations and/or neuronal network interactions, may play an important role in preictal/ictal seizure evolution [2, 31]. However, the role of high-frequency activity in seizure termination, is unknown, if it exists at all. Using information theoretic measures of network coordination, this study investigated ictal and immediate postictal neurodynamic interactions encoded in scalp EEGs from a relatively small sample of 8 patients with focal epilepsy and multiple seizures originating in temporal and/or frontal brain regions, at frequencies ≤100 Hz and >100 Hz, respectively. Despite some heterogeneity in the dynamics of these interactions, consistent patterns were also estimated. Specifically, in several seizures, linear or non-linear increase in high-frequency neuronal coordination during ictal intervals, coincided with a corresponding decrease in coordination at frequencies <100 Hz, suggesting a potential interference role of high-frequency activity, to disrupt abnormal ictal synchrony at lower frequencies. These changes in network synchrony started at least 20–30 sec prior to seizure offset, depending on the seizure duration. Opposite patterns were estimated at frequencies ≤100 Hz in several seizures. These results raise the possibility that high-frequency interference may occur in the form of progressive network coordination during the ictal interval, which continues during the postictal interval. This may be one of several possible mechanisms that facilitate seizure termination. In fact, inhibition of pairwise interactions between EEGs by other signals in their spatial neighborhood, quantified by negative interaction information, was estimated at frequencies ≤100 Hz, at least in some seizures. PMID:23608198

NetDecoder: a network biology platform that decodes context-specific biological networks and gene activities.

PubMed

da Rocha, Edroaldo Lummertz; Ung, Choong Yong; McGehee, Cordelia D; Correia, Cristina; Li, Hu

2016-06-02

The sequential chain of interactions altering the binary state of a biomolecule represents the 'information flow' within a cellular network that determines phenotypic properties. Given the lack of computational tools to dissect context-dependent networks and gene activities, we developed NetDecoder, a network biology platform that models context-dependent information flows using pairwise phenotypic comparative analyses of protein-protein interactions. Using breast cancer, dyslipidemia and Alzheimer's disease as case studies, we demonstrate NetDecoder dissects subnetworks to identify key players significantly impacting cell behaviour specific to a given disease context. We further show genes residing in disease-specific subnetworks are enriched in disease-related signalling pathways and information flow profiles, which drive the resulting disease phenotypes. We also devise a novel scoring scheme to quantify key genes-network routers, which influence many genes, key targets, which are influenced by many genes, and high impact genes, which experience a significant change in regulation. We show the robustness of our results against parameter changes. Our network biology platform includes freely available source code (http://www.NetDecoder.org) for researchers to explore genome-wide context-dependent information flow profiles and key genes, given a set of genes of particular interest and transcriptome data. More importantly, NetDecoder will enable researchers to uncover context-dependent drug targets. © The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research.

Natural history matters: how biological constraints shape diversified interactions in pollination networks.

PubMed

Jordano, Pedro

2016-11-01

Species-specific traits constrain the ways organisms interact in nature. Some pairwise interactions among coexisting species simply do not occur; they are impossible to observe despite the fact that partners coexist in the same place. The author discusses these 'forbidden links' of species interaction networks. Photo: a sphingid moth, Manduca sexta visiting a flower of Tocoyena formosa (Rubiaceae) in the Brazilian Cerrado; tongue and corolla tube lengths approximately 100 mm. Courtesy of Felipe Amorim. Sazatornil, F.D., Moré, M., Benitez-Vieyra, S., Cocucci, A.A., Kitching, I.J., Schlumpberger, B.O., Oliveira, P.E., Sazima, M. & Amorim, F.W. (2016) Beyond neutral and forbidden links: morphological matches and the assembly of mutualistic hawkmoth-plant networks. Journal of Animal Ecology, 85, 1586-1594. Species-specific traits and life-history characteristics constrain the ways organisms interact in nature. For example, gape-limited predators are constrained in the sizes of prey they can handle and efficiently consume. When we consider the ubiquity of such constrains, it is evident how hard it can be to be a generalist partner in ecological interactions: a free-living animal or plant cannot simply interact with every available partner it encounters. Some pairwise interactions among coexisting species simply do not occur; they are impossible to observe despite the fact that partners coexist in the same place. Sazatornil et al. () explore the nature of such constraints in the mutualisms among hawkmoths and the plants they pollinate. In this iconic interaction, used by Darwin and Wallace to vividly illustrate the power of natural selection in shaping evolutionary change, both pollinators and plants are sharply constrained in their interaction modes and outcomes. © 2016 The Author. Journal of Animal Ecology © 2016 British Ecological Society.

ScaffoldSeq: Software for characterization of directed evolution populations.

PubMed

Woldring, Daniel R; Holec, Patrick V; Hackel, Benjamin J

2016-07-01

ScaffoldSeq is software designed for the numerous applications-including directed evolution analysis-in which a user generates a population of DNA sequences encoding for partially diverse proteins with related functions and would like to characterize the single site and pairwise amino acid frequencies across the population. A common scenario for enzyme maturation, antibody screening, and alternative scaffold engineering involves naïve and evolved populations that contain diversified regions, varying in both sequence and length, within a conserved framework. Analyzing the diversified regions of such populations is facilitated by high-throughput sequencing platforms; however, length variability within these regions (e.g., antibody CDRs) encumbers the alignment process. To overcome this challenge, the ScaffoldSeq algorithm takes advantage of conserved framework sequences to quickly identify diverse regions. Beyond this, unintended biases in sequence frequency are generated throughout the experimental workflow required to evolve and isolate clones of interest prior to DNA sequencing. ScaffoldSeq software uniquely handles this issue by providing tools to quantify and remove background sequences, cluster similar protein families, and dampen the impact of dominant clones. The software produces graphical and tabular summaries for each region of interest, allowing users to evaluate diversity in a site-specific manner as well as identify epistatic pairwise interactions. The code and detailed information are freely available at http://research.cems.umn.edu/hackel. Proteins 2016; 84:869-874. © 2016 Wiley Periodicals, Inc. © 2016 Wiley Periodicals, Inc.

Minimum and Maximum Entropy Distributions for Binary Systems with Known Means and Pairwise Correlations

DTIC Science & Technology

2017-08-21

distributions, and we discuss some applications for engineered and biological information transmission systems. Keywords: information theory; minimum...of its interpretation as a measure of the amount of information communicable by a neural system to groups of downstream neurons. Previous authors...of the maximum entropy approach. Our results also have relevance for engineered information transmission systems. We show that empirically measured

When do correlations increase with firing rates in recurrent networks?

PubMed Central

2017-01-01

A central question in neuroscience is to understand how noisy firing patterns are used to transmit information. Because neural spiking is noisy, spiking patterns are often quantified via pairwise correlations, or the probability that two cells will spike coincidentally, above and beyond their baseline firing rate. One observation frequently made in experiments, is that correlations can increase systematically with firing rate. Theoretical studies have determined that stimulus-dependent correlations that increase with firing rate can have beneficial effects on information coding; however, we still have an incomplete understanding of what circuit mechanisms do, or do not, produce this correlation-firing rate relationship. Here, we studied the relationship between pairwise correlations and firing rates in recurrently coupled excitatory-inhibitory spiking networks with conductance-based synapses. We found that with stronger excitatory coupling, a positive relationship emerged between pairwise correlations and firing rates. To explain these findings, we used linear response theory to predict the full correlation matrix and to decompose correlations in terms of graph motifs. We then used this decomposition to explain why covariation of correlations with firing rate—a relationship previously explained in feedforward networks driven by correlated input—emerges in some recurrent networks but not in others. Furthermore, when correlations covary with firing rate, this relationship is reflected in low-rank structure in the correlation matrix. PMID:28448499

Simulations of the pairwise kinematic Sunyaev-Zel'dovich signal

DOE PAGES

Flender, Samuel; Bleem, Lindsey; Finkel, Hal; ...

2016-05-26

The pairwise kinematic Sunyaev–Zel'dovich (kSZ) signal from galaxy clusters is a probe of their line of sight momenta, and thus a potentially valuable source of cosmological information. In addition to the momenta, the amplitude of the measured signal depends on the properties of the intracluster gas and observational limitations such as errors in determining cluster centers and redshifts. In this work, we simulate the pairwise kSZ signal of clusters atmore » $$z\\lt 1$$, using the output from a cosmological N-body simulation and including the properties of the intracluster gas via a model that can be varied in post-processing. We find that modifications to the gas profile due to star formation and feedback reduce the pairwise kSZ amplitude of clusters by $$\\sim 50\\%$$, relative to the naive "gas traces mass" assumption. We demonstrate that miscentering can reduce the overall amplitude of the pairwise kSZ signal by up to 10%, while redshift errors can lead to an almost complete suppression of the signal at small separations. We confirm that a high-significance detection is expected from the combination of data from current generation, high-resolution cosmic microwave background experiments, such as the South Pole Telescope, and cluster samples from optical photometric surveys, such as the Dark Energy Survey. As a result, we forecast that future experiments such as Advanced ACTPol in conjunction with data from the Dark Energy Spectroscopic Instrument will yield detection significances of at least $$20\\sigma $$, and up to $$57\\sigma $$ in an optimistic scenario.« less

A kernel regression approach to gene-gene interaction detection for case-control studies.

PubMed

Larson, Nicholas B; Schaid, Daniel J

2013-11-01

Gene-gene interactions are increasingly being addressed as a potentially important contributor to the variability of complex traits. Consequently, attentions have moved beyond single locus analysis of association to more complex genetic models. Although several single-marker approaches toward interaction analysis have been developed, such methods suffer from very high testing dimensionality and do not take advantage of existing information, notably the definition of genes as functional units. Here, we propose a comprehensive family of gene-level score tests for identifying genetic elements of disease risk, in particular pairwise gene-gene interactions. Using kernel machine methods, we devise score-based variance component tests under a generalized linear mixed model framework. We conducted simulations based upon coalescent genetic models to evaluate the performance of our approach under a variety of disease models. These simulations indicate that our methods are generally higher powered than alternative gene-level approaches and at worst competitive with exhaustive SNP-level (where SNP is single-nucleotide polymorphism) analyses. Furthermore, we observe that simulated epistatic effects resulted in significant marginal testing results for the involved genes regardless of whether or not true main effects were present. We detail the benefits of our methods and discuss potential genome-wide analysis strategies for gene-gene interaction analysis in a case-control study design. © 2013 WILEY PERIODICALS, INC.

Self-propulsion and interactions of catalytic particles in a chemically active medium.

PubMed

Banigan, Edward J; Marko, John F

2016-01-01

Enzymatic "machines," such as catalytic rods or colloids, can self-propel and interact by generating gradients of their substrates. We theoretically investigate the behaviors of such machines in a chemically active environment where their catalytic substrates are continuously synthesized and destroyed, as occurs in living cells. We show how the kinetic properties of the medium modulate self-propulsion and pairwise interactions between machines, with the latter controlled by a tunable characteristic interaction range analogous to the Debye screening length in an electrolytic solution. Finally, we discuss the effective force arising between interacting machines and possible biological applications, such as partitioning of bacterial plasmids.

Precision and recall estimates for two-hybrid screens

PubMed Central

Huang, Hailiang; Bader, Joel S.

2009-01-01

Motivation: Yeast two-hybrid screens are an important method to map pairwise protein interactions. This method can generate spurious interactions (false discoveries), and true interactions can be missed (false negatives). Previously, we reported a capture–recapture estimator for bait-specific precision and recall. Here, we present an improved method that better accounts for heterogeneity in bait-specific error rates. Result: For yeast, worm and fly screens, we estimate the overall false discovery rates (FDRs) to be 9.9%, 13.2% and 17.0% and the false negative rates (FNRs) to be 51%, 42% and 28%. Bait-specific FDRs and the estimated protein degrees are then used to identify protein categories that yield more (or fewer) false positive interactions and more (or fewer) interaction partners. While membrane proteins have been suggested to have elevated FDRs, the current analysis suggests that intrinsic membrane proteins may actually have reduced FDRs. Hydrophobicity is positively correlated with decreased error rates and fewer interaction partners. These methods will be useful for future two-hybrid screens, which could use ultra-high-throughput sequencing for deeper sampling of interacting bait–prey pairs. Availability: All software (C source) and datasets are available as supplemental files and at http://www.baderzone.org under the Lesser GPL v. 3 license. Contact: joel.bader@jhu.edu Supplementary information: Supplementary data are available at Bioinformatics online. PMID:19091773

Unified framework for information integration based on information geometry

PubMed Central

Oizumi, Masafumi; Amari, Shun-ichi

2016-01-01

Assessment of causal influences is a ubiquitous and important subject across diverse research fields. Drawn from consciousness studies, integrated information is a measure that defines integration as the degree of causal influences among elements. Whereas pairwise causal influences between elements can be quantified with existing methods, quantifying multiple influences among many elements poses two major mathematical difficulties. First, overestimation occurs due to interdependence among influences if each influence is separately quantified in a part-based manner and then simply summed over. Second, it is difficult to isolate causal influences while avoiding noncausal confounding influences. To resolve these difficulties, we propose a theoretical framework based on information geometry for the quantification of multiple causal influences with a holistic approach. We derive a measure of integrated information, which is geometrically interpreted as the divergence between the actual probability distribution of a system and an approximated probability distribution where causal influences among elements are statistically disconnected. This framework provides intuitive geometric interpretations harmonizing various information theoretic measures in a unified manner, including mutual information, transfer entropy, stochastic interaction, and integrated information, each of which is characterized by how causal influences are disconnected. In addition to the mathematical assessment of consciousness, our framework should help to analyze causal relationships in complex systems in a complete and hierarchical manner. PMID:27930289

Haplotype Reconstruction in Large Pedigrees with Many Untyped Individuals

NASA Astrophysics Data System (ADS)

Li, Xin; Li, Jing

Haplotypes, as they specify the linkage patterns between dispersed genetic variations, provide important information for understanding the genetics of human traits. However haplotypes are not directly available from current genotyping platforms, and hence there are extensive investigations of computational methods to recover such information. Two major computational challenges arising in current family-based disease studies are large family sizes and many ungenotyped family members. Traditional haplotyping methods can neither handle large families nor families with missing members. In this paper, we propose a method which addresses these issues by integrating multiple novel techniques. The method consists of three major components: pairwise identical-bydescent (IBD) inference, global IBD reconstruction and haplotype restoring. By reconstructing the global IBD of a family from pairwise IBD and then restoring the haplotypes based on the inferred IBD, this method can scale to large pedigrees, and more importantly it can handle families with missing members. Compared with existing methods, this method demonstrates much higher power to recover haplotype information, especially in families with many untyped individuals.

An efficient semi-supervised community detection framework in social networks.

PubMed

Li, Zhen; Gong, Yong; Pan, Zhisong; Hu, Guyu

2017-01-01

Community detection is an important tasks across a number of research fields including social science, biology, and physics. In the real world, topology information alone is often inadequate to accurately find out community structure due to its sparsity and noise. The potential useful prior information such as pairwise constraints which contain must-link and cannot-link constraints can be obtained from domain knowledge in many applications. Thus, combining network topology with prior information to improve the community detection accuracy is promising. Previous methods mainly utilize the must-link constraints while cannot make full use of cannot-link constraints. In this paper, we propose a semi-supervised community detection framework which can effectively incorporate two types of pairwise constraints into the detection process. Particularly, must-link and cannot-link constraints are represented as positive and negative links, and we encode them by adding different graph regularization terms to penalize closeness of the nodes. Experiments on multiple real-world datasets show that the proposed framework significantly improves the accuracy of community detection.

Competition can lead to unexpected patterns in tropical ant communities

NASA Astrophysics Data System (ADS)

Ellwood, M. D. Farnon; Blüthgen, Nico; Fayle, Tom M.; Foster, William A.; Menzel, Florian

2016-08-01

Ecological communities are structured by competitive, predatory, mutualistic and parasitic interactions combined with chance events. Separating deterministic from stochastic processes is possible, but finding statistical evidence for specific biological interactions is challenging. We attempt to solve this problem for ant communities nesting in epiphytic bird's nest ferns (Asplenium nidus) in Borneo's lowland rainforest. By recording the frequencies with which each and every single ant species occurred together, we were able to test statistically for patterns associated with interspecific competition. We found evidence for competition, but the resulting co-occurrence pattern was the opposite of what we expected. Rather than detecting species segregation-the classical hallmark of competition-we found species aggregation. Moreover, our approach of testing individual pairwise interactions mostly revealed spatially positive rather than negative associations. Significant negative interactions were only detected among large ants, and among species of the subfamily Ponerinae. Remarkably, the results from this study, and from a corroborating analysis of ant communities known to be structured by competition, suggest that competition within the ants leads to species aggregation rather than segregation. We believe this unexpected result is linked with the displacement of species following asymmetric competition. We conclude that analysing co-occurrence frequencies across complete species assemblages, separately for each species, and for each unique pairwise combination of species, represents a subtle yet powerful way of detecting structure and compartmentalisation in ecological communities.

Volatiles in Inter-Specific Bacterial Interactions

PubMed Central

Tyc, Olaf; Zweers, Hans; de Boer, Wietse; Garbeva, Paolina

2015-01-01

The importance of volatile organic compounds for functioning of microbes is receiving increased research attention. However, to date very little is known on how inter-specific bacterial interactions effect volatiles production as most studies have been focused on volatiles produced by monocultures of well-described bacterial genera. In this study we aimed to understand how inter-specific bacterial interactions affect the composition, production and activity of volatiles. Four phylogenetically different bacterial species namely: Chryseobacterium, Dyella, Janthinobacterium, and Tsukamurella were selected. Earlier results had shown that pairwise combinations of these bacteria induced antimicrobial activity in agar media whereas this was not the case for monocultures. In the current study, we examined if these observations were also reflected by the production of antimicrobial volatiles. Thus, the identity and antimicrobial activity of volatiles produced by the bacteria were determined in monoculture as well in pairwise combinations. Antimicrobial activity of the volatiles was assessed against fungal, oomycetal, and bacterial model organisms. Our results revealed that inter-specific bacterial interactions affected volatiles blend composition. Fungi and oomycetes showed high sensitivity to bacterial volatiles whereas the effect of volatiles on bacteria varied between no effects, growth inhibition to growth promotion depending on the volatile blend composition. In total 35 volatile compounds were detected most of which were sulfur-containing compounds. Two commonly produced sulfur-containing volatile compounds (dimethyl disulfide and dimethyl trisulfide) were tested for their effect on three target bacteria. Here, we display the importance of inter-specific interactions on bacterial volatiles production and their antimicrobial activities. PMID:26733959

Scalable Creation of Long-Lived Multipartite Entanglement

NASA Astrophysics Data System (ADS)

Kaufmann, H.; Ruster, T.; Schmiegelow, C. T.; Luda, M. A.; Kaushal, V.; Schulz, J.; von Lindenfels, D.; Schmidt-Kaler, F.; Poschinger, U. G.

2017-10-01

We demonstrate the deterministic generation of multipartite entanglement based on scalable methods. Four qubits are encoded in 40Ca+, stored in a microstructured segmented Paul trap. These qubits are sequentially entangled by laser-driven pairwise gate operations. Between these, the qubit register is dynamically reconfigured via ion shuttling operations, where ion crystals are separated and merged, and ions are moved in and out of a fixed laser interaction zone. A sequence consisting of three pairwise entangling gates yields a four-ion Greenberger-Horne-Zeilinger state |ψ ⟩=(1 /√{2 })(|0000 ⟩+|1111 ⟩) , and full quantum state tomography reveals a state fidelity of 94.4(3)%. We analyze the decoherence of this state and employ dynamic decoupling on the spatially distributed constituents to maintain 69(5)% coherence at a storage time of 1.1 sec.

Affective Outcomes of Schooling: Full-Information Item Factor Analysis of a Student Questionnaire.

ERIC Educational Resources Information Center

Muraki, Eiji; Engelhard, George, Jr.

Recent developments in dichotomous factor analysis based on multidimensional item response models (Bock and Aitkin, 1981; Muthen, 1978) provide an effective method for exploring the dimensionality of questionnaire items. Implemented in the TESTFACT program, this "full information" item factor analysis accounts not only for the pairwise joint…

[Analysis of variance of repeated data measured by water maze with SPSS].

PubMed

Qiu, Hong; Jin, Guo-qin; Jin, Ru-feng; Zhao, Wei-kang

2007-01-01

To introduce the method of analyzing repeated data measured by water maze with SPSS 11.0, and offer a reference statistical method to clinical and basic medicine researchers who take the design of repeated measures. Using repeated measures and multivariate analysis of variance (ANOVA) process of the general linear model in SPSS and giving comparison among different groups and different measure time pairwise. Firstly, Mauchly's test of sphericity should be used to judge whether there were relations among the repeatedly measured data. If any (P

Prediction of microsleeps using pairwise joint entropy and mutual information between EEG channels.

PubMed

Baseer, Abdul; Weddell, Stephen J; Jones, Richard D

2017-07-01

Microsleeps are involuntary and brief instances of complete loss of responsiveness, typically of 0.5-15 s duration. They adversely affect performance in extended attention-driven jobs and can be fatal. Our aim was to predict microsleeps from 16 channel EEG signals. Two information theoretic concepts - pairwise joint entropy and mutual information - were independently used to continuously extract features from EEG signals. k-nearest neighbor (kNN) with k = 3 was used to calculate both joint entropy and mutual information. Highly correlated features were discarded and the rest were ranked using Fisher score followed by an average of 3-fold cross-validation area under the curve of the receiver operating characteristic (AUC ROC ). Leave-one-out method (LOOM) was performed to test the performance of microsleep prediction system on independent data. The best prediction for 0.25 s ahead was AUCROC, sensitivity, precision, geometric mean (GM), and φ of 0.93, 0.68, 0.33, 0.75, and 0.38 respectively with joint entropy using single linear discriminant analysis (LDA) classifier.

Evolutionary games in the multiverse

PubMed Central

Gokhale, Chaitanya S.; Traulsen, Arne

2010-01-01

Evolutionary game dynamics of two players with two strategies has been studied in great detail. These games have been used to model many biologically relevant scenarios, ranging from social dilemmas in mammals to microbial diversity. Some of these games may, in fact, take place between a number of individuals and not just between two. Here we address one-shot games with multiple players. As long as we have only two strategies, many results from two-player games can be generalized to multiple players. For games with multiple players and more than two strategies, we show that statements derived for pairwise interactions no longer hold. For two-player games with any number of strategies there can be at most one isolated internal equilibrium. For any number of players with any number of strategies , there can be at most isolated internal equilibria. Multiplayer games show a great dynamical complexity that cannot be captured based on pairwise interactions. Our results hold for any game and can easily be applied to specific cases, such as public goods games or multiplayer stag hunts. PMID:20212124

Pairwise gene GO-based measures for biclustering of high-dimensional expression data.

PubMed

Nepomuceno, Juan A; Troncoso, Alicia; Nepomuceno-Chamorro, Isabel A; Aguilar-Ruiz, Jesús S

2018-01-01

Biclustering algorithms search for groups of genes that share the same behavior under a subset of samples in gene expression data. Nowadays, the biological knowledge available in public repositories can be used to drive these algorithms to find biclusters composed of groups of genes functionally coherent. On the other hand, a distance among genes can be defined according to their information stored in Gene Ontology (GO). Gene pairwise GO semantic similarity measures report a value for each pair of genes which establishes their functional similarity. A scatter search-based algorithm that optimizes a merit function that integrates GO information is studied in this paper. This merit function uses a term that addresses the information through a GO measure. The effect of two possible different gene pairwise GO measures on the performance of the algorithm is analyzed. Firstly, three well known yeast datasets with approximately one thousand of genes are studied. Secondly, a group of human datasets related to clinical data of cancer is also explored by the algorithm. Most of these data are high-dimensional datasets composed of a huge number of genes. The resultant biclusters reveal groups of genes linked by a same functionality when the search procedure is driven by one of the proposed GO measures. Furthermore, a qualitative biological study of a group of biclusters show their relevance from a cancer disease perspective. It can be concluded that the integration of biological information improves the performance of the biclustering process. The two different GO measures studied show an improvement in the results obtained for the yeast dataset. However, if datasets are composed of a huge number of genes, only one of them really improves the algorithm performance. This second case constitutes a clear option to explore interesting datasets from a clinical point of view.

Desert ants achieve reliable recruitment across noisy interactions

PubMed Central

Razin, Nitzan; Eckmann, Jean-Pierre; Feinerman, Ofer

2013-01-01

We study how desert ants, Cataglyphis niger, a species that lacks pheromone-based recruitment mechanisms, inform each other about the presence of food. Our results are based on automated tracking that allows us to collect a large database of ant trajectories and interactions. We find that interactions affect an ant's speed within the nest. Fast ants tend to slow down, whereas slow ones increase their speed when encountering a faster ant. Faster ants tend to exit the nest more frequently than slower ones. So, if an ant gains enough speed through encounters with others, then she tends to leave the nest and look for food. On the other hand, we find that the probability for her to leave the nest depends only on her speed, but not on whether she had recently interacted with a recruiter that has found the food. This suggests a recruitment system in which ants communicate their state by very simple interactions. Based on this assumption, we estimate the information-theoretical channel capacity of the ants’ pairwise interactions. We find that the response to the speed of an interacting nest-mate is very noisy. The question is then how random interactions with ants within the nest can be distinguished from those interactions with a recruiter who has found food. Our measurements and model suggest that this distinction does not depend on reliable communication but on behavioural differences between ants that have found the food and those that have not. Recruiters retain high speeds throughout the experiment, regardless of the ants they interact with; non-recruiters communicate with a limited number of nest-mates and adjust their speed following these interactions. These simple rules lead to the formation of a bistable switch on the level of the group that allows the distinction between recruitment and random noise in the nest. A consequence of the mechanism we propose is a negative effect of ant density on exit rates and recruitment success. This is, indeed, confirmed by our measurements. PMID:23486172

Pick-off annihilation of positronium in matter using full correlation single particle potentials: solid He.

PubMed

Zubiaga, A; Tuomisto, F; Puska, M J

2015-01-29

We investigate the modeling of positronium (Ps) states and their pick-off annihilation trapped at open volumes pockets in condensed molecular matter. Our starting point is the interacting many-body system of Ps and a He atom because it is the smallest entity that can mimic the energy gap between the highest occupied and lowest unoccupied molecular orbitals of molecules, and yet the many-body structure of the HePs system can be calculated accurately enough. The exact-diagonalization solution of the HePs system enables us to construct a pairwise full-correlation single-particle potential for the Ps-He interaction, and the total potential in solids is obtained as a superposition of the pairwise potentials. We study in detail Ps states and their pick-off annihilation rates in voids inside solid He and analyze experimental results for Ps-induced voids in liquid He obtaining the radii of the voids. More importantly, we generalize our conclusions by testing the validity of the Tao-Eldrup model, widely used to analyze ortho-Ps annihilation measurements for voids in molecular matter, against our theoretical results for the solid He. Moreover, we discuss the influence of the partial charges of polar molecules and the strength of the van der Waals interaction on the pick-off annihilation rate.

Working with Missing Values

ERIC Educational Resources Information Center

Acock, Alan C.

2005-01-01

Less than optimum strategies for missing values can produce biased estimates, distorted statistical power, and invalid conclusions. After reviewing traditional approaches (listwise, pairwise, and mean substitution), selected alternatives are covered including single imputation, multiple imputation, and full information maximum likelihood…

A molecular-field-based similarity study of non-nucleoside HIV-1 reverse transcriptase inhibitors

NASA Astrophysics Data System (ADS)

Mestres, Jordi; Rohrer, Douglas C.; Maggiora, Gerald M.

1999-01-01

This article describes a molecular-field-based similarity method for aligning molecules by matching their steric and electrostatic fields and an application of the method to the alignment of three structurally diverse non-nucleoside HIV-1 reverse transcriptase inhibitors. A brief description of the method, as implemented in the program MIMIC, is presented, including a discussion of pairwise and multi-molecule similarity-based matching. The application provides an example that illustrates how relative binding orientations of molecules can be determined in the absence of detailed structural information on their target protein. In the particular system studied here, availability of the X-ray crystal structures of the respective ligand-protein complexes provides a means for constructing an 'experimental model' of the relative binding orientations of the three inhibitors. The experimental model is derived by using MIMIC to align the steric fields of the three protein P66 subunit main chains, producing an overlay with a 1.41 Å average rms distance between the corresponding Cα's in the three chains. The inter-chain residue similarities for the backbone structures show that the main-chain conformations are conserved in the region of the inhibitor-binding site, with the major deviations located primarily in the 'finger' and RNase H regions. The resulting inhibitor structure overlay provides an experimental-based model that can be used to evaluate the quality of the direct a priori inhibitor alignment obtained using MIMIC. It is found that the 'best' pairwise alignments do not always correspond to the experimental model alignments. Therefore, simply combining the best pairwise alignments will not necessarily produce the optimal multi-molecule alignment. However, the best simultaneous three-molecule alignment was found to reproduce the experimental inhibitor alignment model. A pairwise consistency index has been derived which gauges the quality of combining the pairwise alignments and aids in efficiently forming the optimal multi-molecule alignment analysis. Two post-alignment procedures are described that provide information on feature-based and field-based pharmacophoric patterns. The former corresponds to traditional pharmacophore models and is derived from the contribution of individual atoms to the total similarity. The latter is based on molecular regions rather than atoms and is constructed by computing the percent contribution to the similarity of individual points in a regular lattice surrounding the molecules, which when contoured and colored visually depict regions of highly conserved similarity. A discussion of how the information provided by each of the procedures is useful in drug design is also presented.

Beyond pairwise strategy updating in the prisoner's dilemma game

NASA Astrophysics Data System (ADS)

Wang, Xiaofeng; Perc, Matjaž; Liu, Yongkui; Chen, Xiaojie; Wang, Long

2012-10-01

In spatial games players typically alter their strategy by imitating the most successful or one randomly selected neighbor. Since a single neighbor is taken as reference, the information stemming from other neighbors is neglected, which begets the consideration of alternative, possibly more realistic approaches. Here we show that strategy changes inspired not only by the performance of individual neighbors but rather by entire neighborhoods introduce a qualitatively different evolutionary dynamics that is able to support the stable existence of very small cooperative clusters. This leads to phase diagrams that differ significantly from those obtained by means of pairwise strategy updating. In particular, the survivability of cooperators is possible even by high temptations to defect and over a much wider uncertainty range. We support the simulation results by means of pair approximations and analysis of spatial patterns, which jointly highlight the importance of local information for the resolution of social dilemmas.

Adsorption interaction in the molecular hydrogen-aluminophosphate AlPO-5 zeolite system

NASA Astrophysics Data System (ADS)

Grenev, I. V.; Gavrilov, V. Yu.

2015-03-01

The adsorption interaction between molecular hydrogen and atoms forming the lattice of AlPO-5 zeolite is studied. The potential of intramolecular interaction is calculated by summing the potentials of individual pairwise H2-O(Al, P) interactions in a fragment of the zeolite structure with a volume of ˜32 nm3. Isopotential surfaces are constructed that allow determination of the shape of zeolite microchannels and the places of the preferential localization of sorbate molecules in the porous space. The calculated and experimental values of the Henry constant of H2 adsorption on AlPO-5 at 77 K are compared.

Molecular simulations of the pairwise interaction of monoclonal antibodies.

PubMed

Lapelosa, Mauro; Patapoff, Thomas W; Zarraga, Isidro E

2014-11-20

Molecular simulations are employed to compute the free energy of pairwise monoclonal antibodies (mAbs) association using a conformational sampling algorithm with a scoring function. The work reported here is aimed at investigating the mAb-mAb association driven by weak interactions with a computational method capable of predicting experimental observations of low binding affinity. The simulations are able to explore the free energy landscape. A steric interaction component, electrostatic interactions, and a nonpolar component of the free energy form the energy scoring function. Electrostatic interactions are calculated by solving the Poisson-Boltzmann equation. The nonpolar component is derived from the van der Waals interactions upon close contact of the protein surfaces. Two mAbs with similar IgG1 framework but with small sequence differences, mAb1 and mAb2, are considered for their different viscosity and propensity to form a weak interacting dimer. mAb1 presents favorable free energy of association at pH 6 with 15 mM of ion concentration reproducing experimental trends of high viscosity and dimer formation at high concentration. Free energy landscape and minimum free energy configurations of the dimer, as well as the second virial coefficient (B22) values are calculated. The energy distributions for mAb1 are obtained, and the most probable configurations are seen to be consistent with experimental measurements. In contrast, mAb2 shows an unfavorable average free energy at the same buffer conditions due to poor electrostatic complementarity, and reversible dimer configurations with favorable free energy are found to be unlikely. Finally, the simulations of the mAb association dynamics provide insights on the self-association responsible for bulk solution behavior and aggregation, which are important to the processing and the quality of biopharmaceuticals.

Exploring the roles of cannot-link constraint in community detection via Multi-variance Mixed Gaussian Generative Model.

PubMed

Yang, Liang; Ge, Meng; Jin, Di; He, Dongxiao; Fu, Huazhu; Wang, Jing; Cao, Xiaochun

2017-01-01

Due to the demand for performance improvement and the existence of prior information, semi-supervised community detection with pairwise constraints becomes a hot topic. Most existing methods have been successfully encoding the must-link constraints, but neglect the opposite ones, i.e., the cannot-link constraints, which can force the exclusion between nodes. In this paper, we are interested in understanding the role of cannot-link constraints and effectively encoding pairwise constraints. Towards these goals, we define an integral generative process jointly considering the network topology, must-link and cannot-link constraints. We propose to characterize this process as a Multi-variance Mixed Gaussian Generative (MMGG) Model to address diverse degrees of confidences that exist in network topology and pairwise constraints and formulate it as a weighted nonnegative matrix factorization problem. The experiments on artificial and real-world networks not only illustrate the superiority of our proposed MMGG, but also, most importantly, reveal the roles of pairwise constraints. That is, though the must-link is more important than cannot-link when either of them is available, both must-link and cannot-link are equally important when both of them are available. To the best of our knowledge, this is the first work on discovering and exploring the importance of cannot-link constraints in semi-supervised community detection.

Exploring the roles of cannot-link constraint in community detection via Multi-variance Mixed Gaussian Generative Model

PubMed Central

Ge, Meng; Jin, Di; He, Dongxiao; Fu, Huazhu; Wang, Jing; Cao, Xiaochun

2017-01-01

Due to the demand for performance improvement and the existence of prior information, semi-supervised community detection with pairwise constraints becomes a hot topic. Most existing methods have been successfully encoding the must-link constraints, but neglect the opposite ones, i.e., the cannot-link constraints, which can force the exclusion between nodes. In this paper, we are interested in understanding the role of cannot-link constraints and effectively encoding pairwise constraints. Towards these goals, we define an integral generative process jointly considering the network topology, must-link and cannot-link constraints. We propose to characterize this process as a Multi-variance Mixed Gaussian Generative (MMGG) Model to address diverse degrees of confidences that exist in network topology and pairwise constraints and formulate it as a weighted nonnegative matrix factorization problem. The experiments on artificial and real-world networks not only illustrate the superiority of our proposed MMGG, but also, most importantly, reveal the roles of pairwise constraints. That is, though the must-link is more important than cannot-link when either of them is available, both must-link and cannot-link are equally important when both of them are available. To the best of our knowledge, this is the first work on discovering and exploring the importance of cannot-link constraints in semi-supervised community detection. PMID:28678864

A participatory approach to designing and enhancing integrated health information technology systems for veterans: protocol.

PubMed

Haun, Jolie N; Nazi, Kim M; Chavez, Margeaux; Lind, Jason D; Antinori, Nicole; Gosline, Robert M; Martin, Tracey L

2015-02-27

The Department of Veterans Affairs (VA) has developed health information technologies (HIT) and resources to improve veteran access to health care programs and services, and to support a patient-centered approach to health care delivery. To improve VA HIT access and meaningful use by veterans, it is necessary to understand their preferences for interacting with various HIT resources to accomplish health management related tasks and to exchange information. The objective of this paper was to describe a novel protocol for: (1) developing a HIT Digital Health Matrix Model; (2) conducting an Analytic Hierarchy Process called pairwise comparison to understand how and why veterans want to use electronic health resources to complete tasks related to health management; and (3) developing visual modeling simulations that depict veterans' preferences for using VA HIT to manage their health conditions and exchange health information. The study uses participatory research methods to understand how veterans prefer to use VA HIT to accomplish health management tasks within a given context, and how they would like to interact with HIT interfaces (eg, look, feel, and function) in the future. This study includes two rounds of veteran focus groups with self-administered surveys and visual modeling simulation techniques. This study will also convene an expert panel to assist in the development of a VA HIT Digital Health Matrix Model, so that both expert panel members and veteran participants can complete an Analytic Hierarchy Process, pairwise comparisons to evaluate and rank the applicability of electronic health resources for a series of health management tasks. This protocol describes the iterative, participatory, and patient-centered process for: (1) developing a VA HIT Digital Health Matrix Model that outlines current VA patient-facing platforms available to veterans, describing their features and relevant contexts for use; and (2) developing visual model simulations based on direct veteran feedback that depict patient preferences for enhancing the synchronization, integration, and standardization of VA patient-facing platforms. Focus group topics include current uses, preferences, facilitators, and barriers to using electronic health resources; recommendations for synchronizing, integrating, and standardizing VA HIT; and preferences on data sharing and delegation within the VA system. This work highlights the practical, technological, and personal factors that facilitate and inhibit use of current VA HIT, and informs an integrated system redesign. The Digital Health Matrix Model and visual modeling simulations use knowledge of veteran preferences and experiences to directly inform enhancements to VA HIT and provide a more holistic and integrated user experience. These efforts are designed to support the adoption and sustained use of VA HIT to support patient self-management and clinical care coordination in ways that are directly aligned with veteran preferences.

A Participatory Approach to Designing and Enhancing Integrated Health Information Technology Systems for Veterans: Protocol

PubMed Central

Nazi, Kim M; Chavez, Margeaux; Lind, Jason D; Antinori, Nicole; Gosline, Robert M; Martin, Tracey L

2015-01-01

Background The Department of Veterans Affairs (VA) has developed health information technologies (HIT) and resources to improve veteran access to health care programs and services, and to support a patient-centered approach to health care delivery. To improve VA HIT access and meaningful use by veterans, it is necessary to understand their preferences for interacting with various HIT resources to accomplish health management related tasks and to exchange information. Objective The objective of this paper was to describe a novel protocol for: (1) developing a HIT Digital Health Matrix Model; (2) conducting an Analytic Hierarchy Process called pairwise comparison to understand how and why veterans want to use electronic health resources to complete tasks related to health management; and (3) developing visual modeling simulations that depict veterans’ preferences for using VA HIT to manage their health conditions and exchange health information. Methods The study uses participatory research methods to understand how veterans prefer to use VA HIT to accomplish health management tasks within a given context, and how they would like to interact with HIT interfaces (eg, look, feel, and function) in the future. This study includes two rounds of veteran focus groups with self-administered surveys and visual modeling simulation techniques. This study will also convene an expert panel to assist in the development of a VA HIT Digital Health Matrix Model, so that both expert panel members and veteran participants can complete an Analytic Hierarchy Process, pairwise comparisons to evaluate and rank the applicability of electronic health resources for a series of health management tasks. Results This protocol describes the iterative, participatory, and patient-centered process for: (1) developing a VA HIT Digital Health Matrix Model that outlines current VA patient-facing platforms available to veterans, describing their features and relevant contexts for use; and (2) developing visual model simulations based on direct veteran feedback that depict patient preferences for enhancing the synchronization, integration, and standardization of VA patient-facing platforms. Focus group topics include current uses, preferences, facilitators, and barriers to using electronic health resources; recommendations for synchronizing, integrating, and standardizing VA HIT; and preferences on data sharing and delegation within the VA system. Conclusions This work highlights the practical, technological, and personal factors that facilitate and inhibit use of current VA HIT, and informs an integrated system redesign. The Digital Health Matrix Model and visual modeling simulations use knowledge of veteran preferences and experiences to directly inform enhancements to VA HIT and provide a more holistic and integrated user experience. These efforts are designed to support the adoption and sustained use of VA HIT to support patient self-management and clinical care coordination in ways that are directly aligned with veteran preferences. PMID:25803324

Evaluating the Quality of Evidence from a Network Meta-Analysis

PubMed Central

Salanti, Georgia; Del Giovane, Cinzia; Chaimani, Anna; Caldwell, Deborah M.; Higgins, Julian P. T.

2014-01-01

Systematic reviews that collate data about the relative effects of multiple interventions via network meta-analysis are highly informative for decision-making purposes. A network meta-analysis provides two types of findings for a specific outcome: the relative treatment effect for all pairwise comparisons, and a ranking of the treatments. It is important to consider the confidence with which these two types of results can enable clinicians, policy makers and patients to make informed decisions. We propose an approach to determining confidence in the output of a network meta-analysis. Our proposed approach is based on methodology developed by the Grading of Recommendations Assessment, Development and Evaluation (GRADE) Working Group for pairwise meta-analyses. The suggested framework for evaluating a network meta-analysis acknowledges (i) the key role of indirect comparisons (ii) the contributions of each piece of direct evidence to the network meta-analysis estimates of effect size; (iii) the importance of the transitivity assumption to the validity of network meta-analysis; and (iv) the possibility of disagreement between direct evidence and indirect evidence. We apply our proposed strategy to a systematic review comparing topical antibiotics without steroids for chronically discharging ears with underlying eardrum perforations. The proposed framework can be used to determine confidence in the results from a network meta-analysis. Judgements about evidence from a network meta-analysis can be different from those made about evidence from pairwise meta-analyses. PMID:24992266

Solitary water wave interactions

NASA Astrophysics Data System (ADS)

Craig, W.; Guyenne, P.; Hammack, J.; Henderson, D.; Sulem, C.

2006-05-01

This article concerns the pairwise nonlinear interaction of solitary waves in the free surface of a body of water lying over a horizontal bottom. Unlike solitary waves in many completely integrable model systems, solitary waves for the full Euler equations do not collide elastically; after interactions, there is a nonzero residual wave that trails the post-collision solitary waves. In this report on new numerical and experimental studies of such solitary wave interactions, we verify that this is the case, both in head-on collisions (the counterpropagating case) and overtaking collisions (the copropagating case), quantifying the degree to which interactions are inelastic. In the situation in which two identical solitary waves undergo a head-on collision, we compare the asymptotic predictions of Su and Mirie [J. Fluid Mech. 98, 509 (1980)] and Byatt-Smith [J. Fluid Mech. 49, 625 (1971)], the wavetank experiments of Maxworthy [J. Fluid Mech. 76, 177 (1976)], and the numerical results of Cooker, Weidman, and Bale [J. Fluid Mech. 342, 141 (1997)] with independent numerical simulations, in which we quantify the phase change, the run-up, and the form of the residual wave and its Fourier signature in both small- and large-amplitude interactions. This updates the prior numerical observations of inelastic interactions in Fenton and Rienecker [J. Fluid Mech. 118, 411 (1982)]. In the case of two nonidentical solitary waves, our precision wavetank experiments are compared with numerical simulations, again observing the run-up, phase lag, and generation of a residual from the interaction. Considering overtaking solitary wave interactions, we compare our experimental observations, numerical simulations, and the asymptotic predictions of Zou and Su [Phys. Fluids 29, 2113 (1986)], and again we quantify the inelastic residual after collisions in the simulations. Geometrically, our numerical simulations of overtaking interactions fit into the three categories of Korteweg-deVries two-soliton solutions defined in Lax [Commun. Pure Appl. Math. 21, 467 (1968)], with, however, a modification in the parameter regime. In all cases we have considered, collisions are seen to be inelastic, although the degree to which interactions depart from elastic is very small. Finally, we give several theoretical results: (i) a relationship between the change in amplitude of solitary waves due to a pairwise collision and the energy carried away from the interaction by the residual component, and (ii) a rigorous estimate of the size of the residual component of pairwise solitary wave collisions. This estimate is consistent with the analytic results of Schneider and Wayne [Commun. Pure Appl. Math. 53, 1475 (2000)], Wright [SIAM J. Math. Anal. 37, 1161 (2005)], and Bona, Colin, and Lannes [Arch. Rat. Mech. Anal. 178, 373 (2005)]. However, in light of our numerical data, both (i) and (ii) indicate a need to reevaluate the asymptotic results in Su and Mirie [J. Fluid Mech. 98, 509 (1980)] and Zou and Su [Phys. Fluids 29, 2113 (1986)].

Scalable Creation of Long-Lived Multipartite Entanglement.

PubMed

Kaufmann, H; Ruster, T; Schmiegelow, C T; Luda, M A; Kaushal, V; Schulz, J; von Lindenfels, D; Schmidt-Kaler, F; Poschinger, U G

2017-10-13

We demonstrate the deterministic generation of multipartite entanglement based on scalable methods. Four qubits are encoded in ^{40}Ca^{+}, stored in a microstructured segmented Paul trap. These qubits are sequentially entangled by laser-driven pairwise gate operations. Between these, the qubit register is dynamically reconfigured via ion shuttling operations, where ion crystals are separated and merged, and ions are moved in and out of a fixed laser interaction zone. A sequence consisting of three pairwise entangling gates yields a four-ion Greenberger-Horne-Zeilinger state |ψ⟩=(1/sqrt[2])(|0000⟩+|1111⟩), and full quantum state tomography reveals a state fidelity of 94.4(3)%. We analyze the decoherence of this state and employ dynamic decoupling on the spatially distributed constituents to maintain 69(5)% coherence at a storage time of 1.1 sec.

Ar(n)HF van der Waals clusters revisited: II. Energetics and HF vibrational frequency shifts from diffusion Monte Carlo calculations on additive and nonadditive potential-energy surfaces for n=1-12.

PubMed

Jiang, Hao; Xu, Minzhong; Hutson, Jeremy M; Bacić, Zlatko

2005-08-01

The ground-state energies and HF vibrational frequency shifts of Ar(n)HF clusters have been calculated on the nonadditive potential-energy surfaces (PESs) for n=2-7 and on the pairwise-additive PESs for the clusters with n=1-12, using the diffusion Monte Carlo (DMC) method. For n>3, the calculations have been performed for the lowest-energy isomer and several higher-lying isomers which are the closest in energy. They provide information about the isomer dependence of the HF redshift, and enable direct comparison with the experimental data recently obtained in helium nanodroplets. The agreement between theory and experiment is excellent, in particular, for the nonadditive DMC redshifts. The relative, incremental redshifts are reproduced accurately even at the lower level of theory, i.e., the DMC and quantum five-dimensional (rigid Ar(n)) calculations on the pairwise-additive PESs. The nonadditive interactions make a significant contribution to the frequency shift, on the order of 10%-12%, and have to be included in the PESs in order for the theory to yield accurate magnitude of the HF redshift. The energy gaps between the DMC ground states of the cluster isomers are very different from the energy separation of their respective minima on the PES, due to the considerable variations in the intermolecular zero-point energy of different Ar(n)HF isomers.

Investigation of genotype x environment interactions for weaning weight for Herefords in three countries.

PubMed

de Mattos, D; Bertrand, J K; Misztal, I

2000-08-01

The objective of this study was to investigate the possibility of genotype x environment interactions for weaning weight (WWT) between different regions of the United States (US) and between Canada (CA), Uruguay (UY), and US for populations of Hereford cattle. Original data were composed of 487,661, 102,986, and 2,322,722 edited weaning weight records from CA, UY, and US, respectively. A total of 359 sires were identified as having progeny across all three countries; 240 of them had at least one progeny with a record in each environment. The data sets within each country were reduced by retaining records from herds with more than 500 WWT records, with an average contemporary group size of greater than nine animals, and that contained WWT records from progeny or maternal grand-progeny of the across-country sires. Data sets within each country were further reduced by randomly selecting among remaining herds. Four regions within US were defined: Upper Plains (UP), Cornbelt (CB), South (S), and Gulf Coast (GC). Similar sampling criteria and common international sires were used to form the within-US regional data sets. A pairwise analysis was done between countries and regions within US (UP-CB vs S-GC, UP vs CB, and S vs GC) for the estimation of (co)variance components and genetic correlation between environments. An accelerated EM-REML algorithm and a multiple-trait animal model that considered WWT as a different trait in each environment were used to estimate parameters in each pairwise analysis. Direct and maternal (in parentheses) estimated genetic correlations for CA vs UY, CA vs US, US vs UY, UP-CB vs S-GC, UP vs CB, and S vs GC were .88 (.84), .86 (.82), .90 (.85), .88 (.87), .88 (.84), and .87 (.85), respectively. The general absence of genotype x country interactions observed in this study, together with a prior study that showed the similarity of genetic and environmental parameters across the three countries, strongly indicates that a joint WWT genetic evaluation for Hereford cattle could be conducted using a model that treated the information from CA, UY, and US as a single population using single population-wide genetic parameters.

On the streaming model for redshift-space distortions

NASA Astrophysics Data System (ADS)

Kuruvilla, Joseph; Porciani, Cristiano

2018-06-01

The streaming model describes the mapping between real and redshift space for 2-point clustering statistics. Its key element is the probability density function (PDF) of line-of-sight pairwise peculiar velocities. Following a kinetic-theory approach, we derive the fundamental equations of the streaming model for ordered and unordered pairs. In the first case, we recover the classic equation while we demonstrate that modifications are necessary for unordered pairs. We then discuss several statistical properties of the pairwise velocities for DM particles and haloes by using a suite of high-resolution N-body simulations. We test the often used Gaussian ansatz for the PDF of pairwise velocities and discuss its limitations. Finally, we introduce a mixture of Gaussians which is known in statistics as the generalised hyperbolic distribution and show that it provides an accurate fit to the PDF. Once inserted in the streaming equation, the fit yields an excellent description of redshift-space correlations at all scales that vastly outperforms the Gaussian and exponential approximations. Using a principal-component analysis, we reduce the complexity of our model for large redshift-space separations. Our results increase the robustness of studies of anisotropic galaxy clustering and are useful for extending them towards smaller scales in order to test theories of gravity and interacting dark-energy models.

Diversity in virus assembly: biology makes things complicated

NASA Astrophysics Data System (ADS)

Zlotnick, Adam

2008-03-01

Icosahedral viruses have an elegance of geometry that implies a general path of assembly. However, structure alone provides insufficient information. Cowpea Chlorotic Mottle Virus (CCMV), an important system for studying virus assembly, consists of 90 coat protein (CP) homodimers condensed around an RNA genome. The crystal structure (Speir et al, 1995) reveals that assembly causes burial of hydrophobic surface and formation of β hexamers, the intertwining of N-termini of the CPs surrounding a quasi-sixfold. This structural view leads to reasonable and erroneous predictions: (i) CCMV capsids are extremely stable, and (ii) β hexamer formation is critical to assembly. Experimentally, we have found that capsids are based on a network of extremely weak (4-5 kT) pairwise interactions and that pentamer formation is the critical step in assembly kinetics. Because of the fragility of CP-Cp interaction, we can redirect assembly to generate and dissociate tubular nanostructures. The dynamic behavior of CCMV reflects the requirements and peculiarities of an evolved biological system; it does not necessarily reflect the behavior predicted from a more static picture of the virus.

Pairwise Trajectory Management (PTM): Concept Description and Documentation

NASA Technical Reports Server (NTRS)

Jones, Kenneth M.; Graff, Thomas J.; Carreno, Victor; Chartrand, Ryan C.; Kibler, Jennifer L.

2018-01-01

Pairwise Trajectory Management (PTM) is an Interval Management (IM) concept that utilizes airborne and ground-based capabilities to enable the implementation of airborne pairwise spacing capabilities in oceanic regions. The goal of PTM is to use airborne surveillance and tools to manage an "at or greater than" inter-aircraft spacing. Due to the accuracy of Automatic Dependent Surveillance-Broadcast (ADS-B) information and the use of airborne spacing guidance, the minimum PTM spacing distance will be less than distances a controller can support with current automation systems that support oceanic operations. Ground tools assist the controller in evaluating the traffic picture and determining appropriate PTM clearances to be issued. Avionics systems provide guidance information that allows the flight crew to conform to the PTM clearance issued by the controller. The combination of a reduced minimum distance and airborne spacing management will increase the capacity and efficiency of aircraft operations at a given altitude or volume of airspace. This document provides an overview of the proposed application, a description of several key scenarios, a high level discussion of expected air and ground equipment and procedure changes, a description of a NASA human-machine interface (HMI) prototype for the flight crew that would support PTM operations, and initial benefits analysis results. Additionally, included as appendices, are the following documents: the PTM Operational Services and Environment Definition (OSED) document and a companion "Future Considerations for the Pairwise Trajectory Management (PTM) Concept: Potential Future Updates for the PTM OSED" paper, a detailed description of the PTM algorithm and PTM Limit Mach rules, initial PTM safety requirements and safety assessment documents, a detailed description of the design, development, and initial evaluations of the proposed flight crew HMI, an overview of the methodology and results of PTM pilot training requirements focus group and human-in-the-loop testing activities, and the PTM Pilot Guide.

The effect of design modifications to the typographical layout of the New York State elementary science learning standards on user preference and process time

NASA Astrophysics Data System (ADS)

Arnold, Jeffery E.

The purpose of this study was to determine the effect of four different design layouts of the New York State elementary science learning standards on user processing time and preference. Three newly developed layouts contained the same information as the standards core curriculum. In this study, the layout of the core guide is referred to as Book. The layouts of the new documents are referred to as Chart, Map, and Tabloid based on the format used to convey content hierarchy information. Most notably, all the new layouts feature larger page sizes, color, page tabs, and an icon based navigation system (IBNS). A convenience sample of 48 New York State educators representing three educator types (16 pre-service teachers, 16 in-service teachers, and 16 administrators) participated in the study. After completing timed tasks accurately, participants scored each layout based on preference. Educator type and layout were the independent variables, and process time and user preference were the dependent variables. A two-factor experimental design with Educator Type as the between variable and with repeated measures on Layout, the within variable, showed a significant difference in process time for Educator Type and Layout. The main effect for Educator Type (F(2, 45) = 8.03, p <.001) was significant with an observed power of .94, and an effect size of .26. The pair-wise comparisons for process time showed that pre-service teachers (p = .02) and administrators (p =.009) completed the assigned tasks more quickly when compared to in-service teachers. The main effect for Layout (F(3, 135) = 4.47, p =.01) was also significant with an observed power of .80, and an effect size of .09. Pair-wise comparisons showed that the newly developed Chart (p = .019) and Map (p = .032) layouts reduced overall process time when compared to the existing state learning standards (Book). The Layout X Educator type interaction was not significant. The same two-factor experimental design on preference, showed the main effect for Layout (F(3, 135) = 28.43, p =.001) was significant. The observed power was 1.0, with an effect size of .39. Pair-wise comparisons for preference scores showed that the Chart (p = .001), Map (p = .001), and Tabloid (p = .001) were preferred over the Book layout. The Layout Type X Educator Type interaction and the main effect for Educator Type were not significant. This study provides evidence that the newly developed design layouts improve usability (as measured by process time and preference scores) of the New York State elementary science learning standard documents. Features in the new layout design, such as the IBNS, may provide a foundation for a visual language and aid users in navigating standard documents across grade level and subject areas. Implications for the next generation of standard documents are presented.

Spatial assignment of symmetry adapted perturbation theory interaction energy components: The atomic SAPT partition

NASA Astrophysics Data System (ADS)

Parrish, Robert M.; Sherrill, C. David

2014-07-01

We develop a physically-motivated assignment of symmetry adapted perturbation theory for intermolecular interactions (SAPT) into atom-pairwise contributions (the A-SAPT partition). The basic precept of A-SAPT is that the many-body interaction energy components are computed normally under the formalism of SAPT, following which a spatially-localized two-body quasiparticle interaction is extracted from the many-body interaction terms. For electrostatics and induction source terms, the relevant quasiparticles are atoms, which are obtained in this work through the iterative stockholder analysis (ISA) procedure. For the exchange, induction response, and dispersion terms, the relevant quasiparticles are local occupied orbitals, which are obtained in this work through the Pipek-Mezey procedure. The local orbital atomic charges obtained from ISA additionally allow the terms involving local orbitals to be assigned in an atom-pairwise manner. Further summation over the atoms of one or the other monomer allows for a chemically intuitive visualization of the contribution of each atom and interaction component to the overall noncovalent interaction strength. Herein, we present the intuitive development and mathematical form for A-SAPT applied in the SAPT0 approximation (the A-SAPT0 partition). We also provide an efficient series of algorithms for the computation of the A-SAPT0 partition with essentially the same computational cost as the corresponding SAPT0 decomposition. We probe the sensitivity of the A-SAPT0 partition to the ISA grid and convergence parameter, orbital localization metric, and induction coupling treatment, and recommend a set of practical choices which closes the definition of the A-SAPT0 partition. We demonstrate the utility and computational tractability of the A-SAPT0 partition in the context of side-on cation-π interactions and the intercalation of DNA by proflavine. A-SAPT0 clearly shows the key processes in these complicated noncovalent interactions, in systems with up to 220 atoms and 2845 basis functions.

Spatial assignment of symmetry adapted perturbation theory interaction energy components: The atomic SAPT partition

DOE Office of Scientific and Technical Information (OSTI.GOV)

Parrish, Robert M.; Sherrill, C. David, E-mail: sherrill@gatech.edu

2014-07-28

We develop a physically-motivated assignment of symmetry adapted perturbation theory for intermolecular interactions (SAPT) into atom-pairwise contributions (the A-SAPT partition). The basic precept of A-SAPT is that the many-body interaction energy components are computed normally under the formalism of SAPT, following which a spatially-localized two-body quasiparticle interaction is extracted from the many-body interaction terms. For electrostatics and induction source terms, the relevant quasiparticles are atoms, which are obtained in this work through the iterative stockholder analysis (ISA) procedure. For the exchange, induction response, and dispersion terms, the relevant quasiparticles are local occupied orbitals, which are obtained in this work throughmore » the Pipek-Mezey procedure. The local orbital atomic charges obtained from ISA additionally allow the terms involving local orbitals to be assigned in an atom-pairwise manner. Further summation over the atoms of one or the other monomer allows for a chemically intuitive visualization of the contribution of each atom and interaction component to the overall noncovalent interaction strength. Herein, we present the intuitive development and mathematical form for A-SAPT applied in the SAPT0 approximation (the A-SAPT0 partition). We also provide an efficient series of algorithms for the computation of the A-SAPT0 partition with essentially the same computational cost as the corresponding SAPT0 decomposition. We probe the sensitivity of the A-SAPT0 partition to the ISA grid and convergence parameter, orbital localization metric, and induction coupling treatment, and recommend a set of practical choices which closes the definition of the A-SAPT0 partition. We demonstrate the utility and computational tractability of the A-SAPT0 partition in the context of side-on cation-π interactions and the intercalation of DNA by proflavine. A-SAPT0 clearly shows the key processes in these complicated noncovalent interactions, in systems with up to 220 atoms and 2845 basis functions.« less

Length-scale crossover of the hydrophobic interaction in a coarse-grained water model

NASA Astrophysics Data System (ADS)

Chaimovich, Aviel; Shell, M. Scott

2013-11-01

It has been difficult to establish a clear connection between the hydrophobic interaction among small molecules typically studied in molecular simulations (a weak, oscillatory force) and that found between large, macroscopic surfaces in experiments (a strong, monotonic force). Here, we show that both types of interaction can emerge with a simple, core-softened water model that captures water's unique pairwise structure. As in hydrophobic hydration, we find that the hydrophobic interaction manifests a length-scale dependence, exhibiting distinct driving forces in the molecular and macroscopic regimes. Moreover, the ability of this simple model to capture both regimes suggests that several features of the hydrophobic force can be understood merely through water's pair correlations.

Length-scale crossover of the hydrophobic interaction in a coarse-grained water model.

PubMed

Chaimovich, Aviel; Shell, M Scott

2013-11-01

It has been difficult to establish a clear connection between the hydrophobic interaction among small molecules typically studied in molecular simulations (a weak, oscillatory force) and that found between large, macroscopic surfaces in experiments (a strong, monotonic force). Here, we show that both types of interaction can emerge with a simple, core-softened water model that captures water's unique pairwise structure. As in hydrophobic hydration, we find that the hydrophobic interaction manifests a length-scale dependence, exhibiting distinct driving forces in the molecular and macroscopic regimes. Moreover, the ability of this simple model to capture both regimes suggests that several features of the hydrophobic force can be understood merely through water's pair correlations.

Information theoretic measures of network coordination in high-frequency scalp EEG reveal dynamic patterns associated with seizure termination.

PubMed

Stamoulis, Catherine; Schomer, Donald L; Chang, Bernard S

2013-08-01

How a seizure terminates is still under-studied and, despite its clinical importance, remains an obscure phase of seizure evolution. Recent studies of seizure-related scalp EEGs at frequencies >100 Hz suggest that neural activity, in the form of oscillations and/or neuronal network interactions, may play an important role in preictal/ictal seizure evolution (Andrade-Valenca et al., 2011; Stamoulis et al., 2012). However, the role of high-frequency activity in seizure termination, is unknown, if it exists at all. Using information theoretic measures of network coordination, this study investigated ictal and immediate postictal neurodynamic interactions encoded in scalp EEGs from a relatively small sample of 8 patients with focal epilepsy and multiple seizures originating in temporal and/or frontal brain regions, at frequencies ≤ 100 Hz and >100 Hz, respectively. Despite some heterogeneity in the dynamics of these interactions, consistent patterns were also estimated. Specifically, in several seizures, linear or non-linear increase in high-frequency neuronal coordination during ictal intervals, coincided with a corresponding decrease in coordination at frequencies <100 Hz, suggesting a potential interference role of high-frequency activity, to disrupt abnormal ictal synchrony at lower frequencies. These changes in network synchrony started at least 20-30s prior to seizure offset, depending on the seizure duration. Opposite patterns were estimated at frequencies ≤ 100 Hz in several seizures. These results raise the possibility that high-frequency interference may occur in the form of progressive network coordination during the ictal interval, which continues during the postictal interval. This may be one of several possible mechanisms that facilitate seizure termination. In fact, inhibition of pairwise interactions between EEGs by other signals in their spatial neighborhood, quantified by negative interaction information, was estimated at frequencies ≤ 100 Hz, at least in some seizures. Copyright © 2013 Elsevier B.V. All rights reserved.

A Computational Study of Rare Gas Clusters: Stepping Stones to the Solid State

ERIC Educational Resources Information Center

Glendening, Eric D.; Halpern, Arthur M.

2012-01-01

An upper-level undergraduate or beginning graduate project is described in which students obtain the Lennard-Jones 6-12 potential parameters for Ne[subscript 2] and Ar[subscript 2] from ab initio calculations and use the results to express pairwise interactions between the atoms in clusters containing up to N = 60 atoms. The students use simulated…

Inverse MDS: Inferring Dissimilarity Structure from Multiple Item Arrangements

PubMed Central

Kriegeskorte, Nikolaus; Mur, Marieke

2012-01-01

The pairwise dissimilarities of a set of items can be intuitively visualized by a 2D arrangement of the items, in which the distances reflect the dissimilarities. Such an arrangement can be obtained by multidimensional scaling (MDS). We propose a method for the inverse process: inferring the pairwise dissimilarities from multiple 2D arrangements of items. Perceptual dissimilarities are classically measured using pairwise dissimilarity judgments. However, alternative methods including free sorting and 2D arrangements have previously been proposed. The present proposal is novel (a) in that the dissimilarity matrix is estimated by “inverse MDS” based on multiple arrangements of item subsets, and (b) in that the subsets are designed by an adaptive algorithm that aims to provide optimal evidence for the dissimilarity estimates. The subject arranges the items (represented as icons on a computer screen) by means of mouse drag-and-drop operations. The multi-arrangement method can be construed as a generalization of simpler methods: It reduces to pairwise dissimilarity judgments if each arrangement contains only two items, and to free sorting if the items are categorically arranged into discrete piles. Multi-arrangement combines the advantages of these methods. It is efficient (because the subject communicates many dissimilarity judgments with each mouse drag), psychologically attractive (because dissimilarities are judged in context), and can characterize continuous high-dimensional dissimilarity structures. We present two procedures for estimating the dissimilarity matrix: a simple weighted-aligned-average of the partial dissimilarity matrices and a computationally intensive algorithm, which estimates the dissimilarity matrix by iteratively minimizing the error of MDS-predictions of the subject’s arrangements. The Matlab code for interactive arrangement and dissimilarity estimation is available from the authors upon request. PMID:22848204

CUBE: Information-optimized parallel cosmological N-body simulation code

NASA Astrophysics Data System (ADS)

Yu, Hao-Ran; Pen, Ue-Li; Wang, Xin

2018-05-01

CUBE, written in Coarray Fortran, is a particle-mesh based parallel cosmological N-body simulation code. The memory usage of CUBE can approach as low as 6 bytes per particle. Particle pairwise (PP) force, cosmological neutrinos, spherical overdensity (SO) halofinder are included.

Machine learning predictions of molecular properties: Accurate many-body potentials and nonlocality in chemical space

DOE PAGES

Hansen, Katja; Biegler, Franziska; Ramakrishnan, Raghunathan; ...

2015-06-04

Simultaneously accurate and efficient prediction of molecular properties throughout chemical compound space is a critical ingredient toward rational compound design in chemical and pharmaceutical industries. Aiming toward this goal, we develop and apply a systematic hierarchy of efficient empirical methods to estimate atomization and total energies of molecules. These methods range from a simple sum over atoms, to addition of bond energies, to pairwise interatomic force fields, reaching to the more sophisticated machine learning approaches that are capable of describing collective interactions between many atoms or bonds. In the case of equilibrium molecular geometries, even simple pairwise force fields demonstratemore » prediction accuracy comparable to benchmark energies calculated using density functional theory with hybrid exchange-correlation functionals; however, accounting for the collective many-body interactions proves to be essential for approaching the “holy grail” of chemical accuracy of 1 kcal/mol for both equilibrium and out-of-equilibrium geometries. This remarkable accuracy is achieved by a vectorized representation of molecules (so-called Bag of Bonds model) that exhibits strong nonlocality in chemical space. The same representation allows us to predict accurate electronic properties of molecules, such as their polarizability and molecular frontier orbital energies.« less

Machine Learning Predictions of Molecular Properties: Accurate Many-Body Potentials and Nonlocality in Chemical Space

PubMed Central

2015-01-01

Simultaneously accurate and efficient prediction of molecular properties throughout chemical compound space is a critical ingredient toward rational compound design in chemical and pharmaceutical industries. Aiming toward this goal, we develop and apply a systematic hierarchy of efficient empirical methods to estimate atomization and total energies of molecules. These methods range from a simple sum over atoms, to addition of bond energies, to pairwise interatomic force fields, reaching to the more sophisticated machine learning approaches that are capable of describing collective interactions between many atoms or bonds. In the case of equilibrium molecular geometries, even simple pairwise force fields demonstrate prediction accuracy comparable to benchmark energies calculated using density functional theory with hybrid exchange-correlation functionals; however, accounting for the collective many-body interactions proves to be essential for approaching the “holy grail” of chemical accuracy of 1 kcal/mol for both equilibrium and out-of-equilibrium geometries. This remarkable accuracy is achieved by a vectorized representation of molecules (so-called Bag of Bonds model) that exhibits strong nonlocality in chemical space. In addition, the same representation allows us to predict accurate electronic properties of molecules, such as their polarizability and molecular frontier orbital energies. PMID:26113956

Statistical Mechanics of US Supreme Court

NASA Astrophysics Data System (ADS)

Lee, Edward; Broedersz, Chase; Bialek, William; Biophysics Theory Group Team

2014-03-01

We build simple models for the distribution of voting patterns in a group, using the Supreme Court of the United States as an example. The least structured, or maximum entropy, model that is consistent with the observed pairwise correlations among justices' votes is equivalent to an Ising spin glass. While all correlations (perhaps surprisingly) are positive, the effective pairwise interactions in the spin glass model have both signs, recovering some of our intuition that justices on opposite sides of the ideological spectrum should have a negative influence on one another. Despite the competing interactions, a strong tendency toward unanimity emerges from the model, and this agrees quantitatively with the data. The model shows that voting patterns are organized in a relatively simple ``energy landscape,'' correctly predicts the extent to which each justice is correlated with the majority, and gives us a measure of the influence that justices exert on one another. These results suggest that simple models, grounded in statistical physics, can capture essential features of collective decision making quantitatively, even in a complex political context. Funded by National Science Foundation Grants PHY-0957573 and CCF-0939370, WM Keck Foundation, Lewis-Sigler Fellowship, Burroughs Wellcome Fund, and Winston Foundation.

Caveats for the spatial arrangement method: Comment on Hout, Goldinger, and Ferguson (2013).

PubMed

Verheyen, Steven; Voorspoels, Wouter; Vanpaemel, Wolf; Storms, Gert

2016-03-01

The gold standard among proximity data collection methods for multidimensional scaling is the (dis)similarity rating of pairwise presented stimuli. A drawback of the pairwise method is its lengthy duration, which may cause participants to change their strategy over time, become fatigued, or disengage altogether. Hout, Goldinger, and Ferguson (2013) recently made a case for the Spatial Arrangement Method (SpAM) as an alternative to the pairwise method, arguing that it is faster and more engaging. SpAM invites participants to directly arrange stimuli on a computer screen such that the interstimuli distances are proportional to psychological proximity. Based on a reanalysis of the Hout et al. (2013), data we identify three caveats for SpAM. An investigation of the distributional characteristics of the SpAM proximity data reveals that the spatial nature of SpAM imposes structure on the data, invoking a bias against featural representations. Individual-differences scaling of the SpAM proximity data reveals that the two-dimensional nature of SpAM allows individuals to only communicate two dimensions of variation among stimuli properly, invoking a bias against high-dimensional scaling representations. Monte Carlo simulations indicate that in order to obtain reliable estimates of the group average, SpAM requires more individuals to be tested. We conclude with an overview of considerations that can inform the choice between SpAM and the pairwise method and offer suggestions on how to overcome their respective limitations. (c) 2016 APA, all rights reserved).

Epistasis between neurochemical gene polymorphisms and risk for ADHD

PubMed Central

Segurado, Ricardo; Bellgrove, Mark A; Manconi, Francesca; Gill, Michael; Hawi, Ziarah

2011-01-01

A number of genes with function related to synaptic neurochemistry have been genetically associated with attention deficit/hyperactivity disorder. However, susceptibility to the development of common psychiatric disorders by single variants acting alone, can so far only explain a small proportion of the heritability of the phenotype. It has been postulated that the unexplained ‘dark heritability' may at least in part be due to epistatic effects, which may account for the small observed marginal associations, and the difficulties with replication of positive findings. We undertook a comprehensive exploration of pair-wise interactions between genetic variants in 24 candidate genic regions involved in monoaminergic catabolism, anabolism, release, re-uptake and signal transmission in a sample of 177 parent-affected child trios using a case-only design and a case–pseudocontrol design using conditional logistic regression. Marker-pairs thresholded on interaction odds ratio (OR) and P-value are presented. We detected a number of interaction ORs >4.0, including an interesting correlation between markers in the ADRA1B and DBH genes in affected individuals, and several further interesting but smaller effects. These effects are no larger than you would expect by chance under the assumption of independence of all pair-wise relations; however, independence is unlikely. Furthermore, the size of these effects is of interest and attempts to replicate these results in other samples are anticipated. PMID:21368917

Spatial interactions between urban areas and cause-specific mortality differentials in France.

PubMed

Ghosn, Walid; Kassie, Daouda; Jougla, Eric; Rican, Stéphane; Rey, Grégoire

2013-11-01

Spatial interactions constitute a challenging but promising approach for investigation of spatial mortality inequalities. Among spatial interactions measures, between-spatial unit migration differentials are a marker of socioeconomic imbalance, but also reflect discrepancies due to other factors. Specifically, this paper asks whether population exchange intensities measure differentials or similarities that are not captured by usual socioeconomic indicators. Urban areas were grouped pairwise by the intensity of connection estimated from a gravity model. The mortality differences for several causes of death were observed to be significantly smaller for strongly connected pairs than for weakly connected pairs even after adjustment on deprivation. © 2013 Published by Elsevier Ltd.

Dynamical pairwise entanglement and two-point correlations in the three-ligand spin-star structure

NASA Astrophysics Data System (ADS)

Motamedifar, M.

2017-10-01

We consider the three-ligand spin-star structure through homogeneous Heisenberg interactions (XXX-3LSSS) in the framework of dynamical pairwise entanglement. It is shown that the time evolution of the central qubit ;one-particle; state (COPS) brings about the generation of quantum W states at periodical time instants. On the contrary, W states cannot be generated from the time evolution of a ligand ;one-particle; state (LOPS). We also investigate the dynamical behavior of two-point quantum correlations as well as the expectation values of the different spin-components for each element in the XXX-3LSSS. It is found that when a W state is generated, the same value of the concurrence between any two arbitrary qubits arises from the xx and yy two-point quantum correlations. On the opposite, zz quantum correlation between any two qubits vanishes at these time instants.

Analysis of Geographic and Pairwise Distances among Chinese Cashmere Goat Populations

PubMed Central

Liu, Jian-Bin; Wang, Fan; Lang, Xia; Zha, Xi; Sun, Xiao-Ping; Yue, Yao-Jing; Feng, Rui-Lin; Yang, Bo-Hui; Guo, Jian

2013-01-01

This study investigated the geographic and pairwise distances of nine Chinese local Cashmere goat populations through the analysis of 20 microsatellite DNA markers. Fluorescence PCR was used to identify the markers, which were selected based on their significance as identified by the Food and Agriculture Organization of the United Nations (FAO) and the International Society for Animal Genetics (ISAG). In total, 206 alleles were detected; the average allele number was 10.30; the polymorphism information content of loci ranged from 0.5213 to 0.7582; the number of effective alleles ranged from 4.0484 to 4.6178; the observed heterozygosity was from 0.5023 to 0.5602 for the practical sample; the expected heterozygosity ranged from 0.5783 to 0.6464; and Allelic richness ranged from 4.7551 to 8.0693. These results indicated that Chinese Cashmere goat populations exhibited rich genetic diversity. Further, the Wright’s F-statistics of subpopulation within total (FST) was 0.1184; the genetic differentiation coefficient (GST) was 0.0940; and the average gene flow (Nm) was 2.0415. All pairwise FST values among the populations were highly significant (p<0.01 or p<0.001), suggesting that the populations studied should all be considered to be separate breeds. Finally, the clustering analysis divided the Chinese Cashmere goat populations into at least four clusters, with the Hexi and Yashan goat populations alone in one cluster. These results have provided useful, practical, and important information for the future of Chinese Cashmere goat breeding. PMID:25049794

Living network meta-analysis compared with pairwise meta-analysis in comparative effectiveness research: empirical study.

PubMed

Nikolakopoulou, Adriani; Mavridis, Dimitris; Furukawa, Toshi A; Cipriani, Andrea; Tricco, Andrea C; Straus, Sharon E; Siontis, George C M; Egger, Matthias; Salanti, Georgia

2018-02-28

To examine whether the continuous updating of networks of prospectively planned randomised controlled trials (RCTs) ("living" network meta-analysis) provides strong evidence against the null hypothesis in comparative effectiveness of medical interventions earlier than the updating of conventional, pairwise meta-analysis. Empirical study of the accumulating evidence about the comparative effectiveness of clinical interventions. Database of network meta-analyses of RCTs identified through searches of Medline, Embase, and the Cochrane Database of Systematic Reviews until 14 April 2015. Network meta-analyses published after January 2012 that compared at least five treatments and included at least 20 RCTs. Clinical experts were asked to identify in each network the treatment comparison of greatest clinical interest. Comparisons were excluded for which direct and indirect evidence disagreed, based on side, or node, splitting test (P<0.10). Cumulative pairwise and network meta-analyses were performed for each selected comparison. Monitoring boundaries of statistical significance were constructed and the evidence against the null hypothesis was considered to be strong when the monitoring boundaries were crossed. A significance level was defined as α=5%, power of 90% (β=10%), and an anticipated treatment effect to detect equal to the final estimate from the network meta-analysis. The frequency and time to strong evidence was compared against the null hypothesis between pairwise and network meta-analyses. 49 comparisons of interest from 44 networks were included; most (n=39, 80%) were between active drugs, mainly from the specialties of cardiology, endocrinology, psychiatry, and rheumatology. 29 comparisons were informed by both direct and indirect evidence (59%), 13 by indirect evidence (27%), and 7 by direct evidence (14%). Both network and pairwise meta-analysis provided strong evidence against the null hypothesis for seven comparisons, but for an additional 10 comparisons only network meta-analysis provided strong evidence against the null hypothesis (P=0.002). The median time to strong evidence against the null hypothesis was 19 years with living network meta-analysis and 23 years with living pairwise meta-analysis (hazard ratio 2.78, 95% confidence interval 1.00 to 7.72, P=0.05). Studies directly comparing the treatments of interest continued to be published for eight comparisons after strong evidence had become evident in network meta-analysis. In comparative effectiveness research, prospectively planned living network meta-analyses produced strong evidence against the null hypothesis more often and earlier than conventional, pairwise meta-analyses. Published by the BMJ Publishing Group Limited. For permission to use (where not already granted under a licence) please go to http://group.bmj.com/group/rights-licensing/permissions.

Living network meta-analysis compared with pairwise meta-analysis in comparative effectiveness research: empirical study

PubMed Central

Nikolakopoulou, Adriani; Mavridis, Dimitris; Furukawa, Toshi A; Cipriani, Andrea; Tricco, Andrea C; Straus, Sharon E; Siontis, George C M; Egger, Matthias

2018-01-01

Abstract Objective To examine whether the continuous updating of networks of prospectively planned randomised controlled trials (RCTs) (“living” network meta-analysis) provides strong evidence against the null hypothesis in comparative effectiveness of medical interventions earlier than the updating of conventional, pairwise meta-analysis. Design Empirical study of the accumulating evidence about the comparative effectiveness of clinical interventions. Data sources Database of network meta-analyses of RCTs identified through searches of Medline, Embase, and the Cochrane Database of Systematic Reviews until 14 April 2015. Eligibility criteria for study selection Network meta-analyses published after January 2012 that compared at least five treatments and included at least 20 RCTs. Clinical experts were asked to identify in each network the treatment comparison of greatest clinical interest. Comparisons were excluded for which direct and indirect evidence disagreed, based on side, or node, splitting test (P<0.10). Outcomes and analysis Cumulative pairwise and network meta-analyses were performed for each selected comparison. Monitoring boundaries of statistical significance were constructed and the evidence against the null hypothesis was considered to be strong when the monitoring boundaries were crossed. A significance level was defined as α=5%, power of 90% (β=10%), and an anticipated treatment effect to detect equal to the final estimate from the network meta-analysis. The frequency and time to strong evidence was compared against the null hypothesis between pairwise and network meta-analyses. Results 49 comparisons of interest from 44 networks were included; most (n=39, 80%) were between active drugs, mainly from the specialties of cardiology, endocrinology, psychiatry, and rheumatology. 29 comparisons were informed by both direct and indirect evidence (59%), 13 by indirect evidence (27%), and 7 by direct evidence (14%). Both network and pairwise meta-analysis provided strong evidence against the null hypothesis for seven comparisons, but for an additional 10 comparisons only network meta-analysis provided strong evidence against the null hypothesis (P=0.002). The median time to strong evidence against the null hypothesis was 19 years with living network meta-analysis and 23 years with living pairwise meta-analysis (hazard ratio 2.78, 95% confidence interval 1.00 to 7.72, P=0.05). Studies directly comparing the treatments of interest continued to be published for eight comparisons after strong evidence had become evident in network meta-analysis. Conclusions In comparative effectiveness research, prospectively planned living network meta-analyses produced strong evidence against the null hypothesis more often and earlier than conventional, pairwise meta-analyses. PMID:29490922

A conditional Granger causality model approach for group analysis in functional MRI

PubMed Central

Zhou, Zhenyu; Wang, Xunheng; Klahr, Nelson J.; Liu, Wei; Arias, Diana; Liu, Hongzhi; von Deneen, Karen M.; Wen, Ying; Lu, Zuhong; Xu, Dongrong; Liu, Yijun

2011-01-01

Granger causality model (GCM) derived from multivariate vector autoregressive models of data has been employed for identifying effective connectivity in the human brain with functional MR imaging (fMRI) and to reveal complex temporal and spatial dynamics underlying a variety of cognitive processes. In the most recent fMRI effective connectivity measures, pairwise GCM has commonly been applied based on single voxel values or average values from special brain areas at the group level. Although a few novel conditional GCM methods have been proposed to quantify the connections between brain areas, our study is the first to propose a viable standardized approach for group analysis of an fMRI data with GCM. To compare the effectiveness of our approach with traditional pairwise GCM models, we applied a well-established conditional GCM to pre-selected time series of brain regions resulting from general linear model (GLM) and group spatial kernel independent component analysis (ICA) of an fMRI dataset in the temporal domain. Datasets consisting of one task-related and one resting-state fMRI were used to investigate connections among brain areas with the conditional GCM method. With the GLM detected brain activation regions in the emotion related cortex during the block design paradigm, the conditional GCM method was proposed to study the causality of the habituation between the left amygdala and pregenual cingulate cortex during emotion processing. For the resting-state dataset, it is possible to calculate not only the effective connectivity between networks but also the heterogeneity within a single network. Our results have further shown a particular interacting pattern of default mode network (DMN) that can be characterized as both afferent and efferent influences on the medial prefrontal cortex (mPFC) and posterior cingulate cortex (PCC). These results suggest that the conditional GCM approach based on a linear multivariate vector autoregressive (MVAR) model can achieve greater accuracy in detecting network connectivity than the widely used pairwise GCM, and this group analysis methodology can be quite useful to extend the information obtainable in fMRI. PMID:21232892

Wiki surveys: open and quantifiable social data collection.

PubMed

Salganik, Matthew J; Levy, Karen E C

2015-01-01

In the social sciences, there is a longstanding tension between data collection methods that facilitate quantification and those that are open to unanticipated information. Advances in technology now enable new, hybrid methods that combine some of the benefits of both approaches. Drawing inspiration from online information aggregation systems like Wikipedia and from traditional survey research, we propose a new class of research instruments called wiki surveys. Just as Wikipedia evolves over time based on contributions from participants, we envision an evolving survey driven by contributions from respondents. We develop three general principles that underlie wiki surveys: they should be greedy, collaborative, and adaptive. Building on these principles, we develop methods for data collection and data analysis for one type of wiki survey, a pairwise wiki survey. Using two proof-of-concept case studies involving our free and open-source website www.allourideas.org, we show that pairwise wiki surveys can yield insights that would be difficult to obtain with other methods.

Wiki Surveys: Open and Quantifiable Social Data Collection

PubMed Central

Salganik, Matthew J.; Levy, Karen E. C.

2015-01-01

In the social sciences, there is a longstanding tension between data collection methods that facilitate quantification and those that are open to unanticipated information. Advances in technology now enable new, hybrid methods that combine some of the benefits of both approaches. Drawing inspiration from online information aggregation systems like Wikipedia and from traditional survey research, we propose a new class of research instruments called wiki surveys. Just as Wikipedia evolves over time based on contributions from participants, we envision an evolving survey driven by contributions from respondents. We develop three general principles that underlie wiki surveys: they should be greedy, collaborative, and adaptive. Building on these principles, we develop methods for data collection and data analysis for one type of wiki survey, a pairwise wiki survey. Using two proof-of-concept case studies involving our free and open-source website www.allourideas.org, we show that pairwise wiki surveys can yield insights that would be difficult to obtain with other methods. PMID:25992565

Experimental characterization of pairwise correlations from triple quantum correlated beams generated by cascaded four-wave mixing processes

NASA Astrophysics Data System (ADS)

Wang, Wei; Cao, Leiming; Lou, Yanbo; Du, Jinjian; Jing, Jietai

2018-01-01

We theoretically and experimentally characterize the performance of the pairwise correlations from triple quantum correlated beams based on the cascaded four-wave mixing (FWM) processes. The pairwise correlations between any two of the beams are theoretically calculated and experimentally measured. The experimental and theoretical results are in good agreement. We find that two of the three pairwise correlations can be in the quantum regime. The other pairwise correlation is always in the classical regime. In addition, we also measure the triple-beam correlation which is always in the quantum regime. Such unbalanced and controllable pairwise correlation structures may be taken as advantages in practical quantum communications, for example, hierarchical quantum secret sharing. Our results also open the way for the classification and application of quantum states generated from the cascaded FWM processes.

Measurements of the pairwise kinematic Sunyaev-Zel'dovich effect with the Atacama Cosmology Telescope and future surveys

NASA Astrophysics Data System (ADS)

Vavagiakis, Eve Marie; De Bernardis, Francesco; Aiola, Simone; Battaglia, Nicholas; Niemack, Michael D.; ACTPol Collaboration

2017-06-01

We have made improved measurements of the kinematic Sunyaev-Zel’dovich (kSZ) effect using data from the Atacama Cosmology Telescope (ACT) and the Baryon Oscillation Spectroscopic Survey (BOSS). We used a map of the Cosmic Microwave Background (CMB) from two seasons of observations each by ACT and the Atacama Cosmology Telescope Polarimeter (ACTPol) receiver. We evaluated the mean pairwise baryon momentum associated with the positions of 50,000 bright galaxies in the BOSS DR11 Large Scale Structure catalog via 600 square degrees of overlapping sky area. The measurement of the kSZ signal arising from the large-scale motions of clusters was made by fitting data to an analytical model. The free parameter of the fit determined the optical depth to microwave photon scattering for the cluster sample. We estimated the covariance matrix of the mean pairwise momentum as a function of galaxy separation using CMB simulations, jackknife evaluation, and bootstrap estimates. The most conservative simulation-based uncertainties gave signal-to-noise estimates between 3.6 and 4.1 for various luminosity cuts. Additionally, we explored a novel approach to estimating cluster optical depths from the average thermal Sunyaev-Zel’dovich (tSZ) signal at the BOSS DR11 catalog positions. Our results were broadly consistent with those obtained from the kSZ signal. In the future, the tSZ signal may provide a valuable probe of cluster optical depths, enabling the extraction of velocities from the kSZ sourced mean pairwise momenta. New CMB maps from three seasons of ACTPol observations with multi-frequency coverage overlap with nearly four times as many DR11 sources and promise to improve statistics and systematics for SZ measurements. With these and other upcoming data, the pairwise kSZ signal is poised to become a powerful new cosmological tool, able to probe large physical scales to inform neutrino physics and test models of modified gravity and dark energy.

Post-Hartree-Fock studies of the He/Mg(0001) interaction: Anti-corrugation, screening, and pairwise additivity

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lara-Castells, María Pilar de, E-mail: Pilar.deLara.Castells@csic.es; Fernández-Perea, Ricardo; Madzharova, Fani

2016-06-28

The adsorption of noble gases on metallic surfaces represents a paradigmatic case of van-der-Waals (vdW) interaction due to the role of screening effects on the corrugation of the interaction potential [J. L. F. Da Silva et al., Phys. Rev. Lett. 90, 066104 (2003)]. The extremely small adsorption energy of He atoms on the Mg(0001) surface (below 3 meV) and the delocalized nature and mobility of the surface electrons make the He/Mg(0001) system particularly challenging, even for state-of-the-art vdW-corrected density functional-based (vdW-DFT) approaches [M. P. de Lara-Castells et al., J. Chem. Phys. 143, 194701 (2015)]. In this work, we meet thismore » challenge by applying two different procedures. First, the dispersion-corrected second-order Möller-Plesset perturbation theory (MP2C) approach is adopted, using bare metal clusters of increasing size. Second, the method of increments [H. Stoll, J. Chem. Phys. 97, 8449 (1992)] is applied at coupled cluster singles and doubles and perturbative triples level, using embedded cluster models of the metal surface. Both approaches provide clear evidences of the anti-corrugation of the interaction potential: the He atom prefers on-top sites, instead of the expected hollow sites. This is interpreted as a signature of the screening of the He atom by the metal for the on-top configuration. The strong screening in the metal is clearly reflected in the relative contribution of successively deeper surface layers to the main dispersion contribution. Aimed to assist future dynamical simulations, a pairwise potential model for the He/surface interaction as a sum of effective He–Mg pair potentials is also presented, as an improvement of the approximation using isolated He–Mg pairs.« less

Biotic interactions reduce microbial carbon use efficiency

NASA Astrophysics Data System (ADS)

Bradford, M.; Maynard, D. S.

2017-12-01

The efficiency by which microbes decompose organic matter governs the amount of carbon that is retained in microbial biomass versus lost to the atmosphere as respiration. This carbon use efficiency (CUE) is affected by various abiotic conditions, such as temperature and nutrient availability. In biogeochemical model simulations, CUE is a key variable regulating how much soil carbon is stored or lost from ecosystems under simulated global changes, such as climate warming. Theoretically, the physiological costs of biotic interactions such as competition should likewise alter CUE, yet the direction and magnitude of these costs are untested. Here we conduct a microcosm experiment to quantify how competitive interactions among saprotrophic fungi alter growth, respiration, and CUE. Free-living decomposer fungi representing a broad range of traits and phylogenies were grown alone, in pairwise competition, and in multi-species (up to 15) communities. By combing culturing and stable carbon isotope approaches, we could resolve the amount of carbon substrate allocated to fungal biomass versus respiration, and so estimate CUE. By then comparing individual performance to community-level outcomes, we show that species interactions induce consistent declines in CUE, regardless of abiotic conditions. Pairwise competition lowers CUE by as much as 25%, with the magnitude of these costs equal to or greater than the observed variation across abiotic conditions. However, depending on the competitive network structure, increasing species richness led to consistent gains or declines in CUE. Our results suggest that the extent to which microbial-mediated carbon fluxes respond to environmental change may be influenced strongly by competitive interactions. As such, knowledge of abiotic conditions and community composition is necessary to confidently project CUE and hence ecosystem carbon dynamics.

Many-body van der Waals interactions in molecules and condensed matter.

PubMed

DiStasio, Robert A; Gobre, Vivekanand V; Tkatchenko, Alexandre

2014-05-28

This work reviews the increasing evidence that many-body van der Waals (vdW) or dispersion interactions play a crucial role in the structure, stability and function of a wide variety of systems in biology, chemistry and physics. Starting with the exact expression for the electron correlation energy provided by the adiabatic connection fluctuation-dissipation theorem, we derive both pairwise and many-body interatomic methods for computing the long-range dispersion energy by considering a model system of coupled quantum harmonic oscillators within the random-phase approximation. By coupling this approach to density functional theory, the resulting many-body dispersion (MBD) method provides an accurate and efficient scheme for computing the frequency-dependent polarizability and many-body vdW energy in molecules and materials with a finite electronic gap. A select collection of applications are presented that ascertain the fundamental importance of these non-bonded interactions across the spectrum of intermolecular (the S22 and S66 benchmark databases), intramolecular (conformational energies of alanine tetrapeptide) and supramolecular (binding energy of the 'buckyball catcher') complexes, as well as molecular crystals (cohesive energies in oligoacenes). These applications demonstrate that electrodynamic response screening and beyond-pairwise many-body vdW interactions--both captured at the MBD level of theory--play a quantitative, and sometimes even qualitative, role in describing the properties considered herein. This work is then concluded with an in-depth discussion of the challenges that remain in the future development of reliable (accurate and efficient) methods for treating many-body vdW interactions in complex materials and provides a roadmap for navigating many of the research avenues that are yet to be explored.

Post-Hartree-Fock studies of the He/Mg(0001) interaction: Anti-corrugation, screening, and pairwise additivity

NASA Astrophysics Data System (ADS)

de Lara-Castells, María Pilar; Fernández-Perea, Ricardo; Madzharova, Fani; Voloshina, Elena

2016-06-01

The adsorption of noble gases on metallic surfaces represents a paradigmatic case of van-der-Waals (vdW) interaction due to the role of screening effects on the corrugation of the interaction potential [J. L. F. Da Silva et al., Phys. Rev. Lett. 90, 066104 (2003)]. The extremely small adsorption energy of He atoms on the Mg(0001) surface (below 3 meV) and the delocalized nature and mobility of the surface electrons make the He/Mg(0001) system particularly challenging, even for state-of-the-art vdW-corrected density functional-based (vdW-DFT) approaches [M. P. de Lara-Castells et al., J. Chem. Phys. 143, 194701 (2015)]. In this work, we meet this challenge by applying two different procedures. First, the dispersion-corrected second-order Möller-Plesset perturbation theory (MP2C) approach is adopted, using bare metal clusters of increasing size. Second, the method of increments [H. Stoll, J. Chem. Phys. 97, 8449 (1992)] is applied at coupled cluster singles and doubles and perturbative triples level, using embedded cluster models of the metal surface. Both approaches provide clear evidences of the anti-corrugation of the interaction potential: the He atom prefers on-top sites, instead of the expected hollow sites. This is interpreted as a signature of the screening of the He atom by the metal for the on-top configuration. The strong screening in the metal is clearly reflected in the relative contribution of successively deeper surface layers to the main dispersion contribution. Aimed to assist future dynamical simulations, a pairwise potential model for the He/surface interaction as a sum of effective He-Mg pair potentials is also presented, as an improvement of the approximation using isolated He-Mg pairs.

Information processing in micro and meso-scale neural circuits during normal and disease states

NASA Astrophysics Data System (ADS)

Luongo, Francisco

Neural computation can occur at multiple spatial and temporal timescales. The sum total of all of these processes is to guide optimal behaviors within the context of the constraints imposed by the physical world. How the circuits of the brain achieves this goal represents a central question in systems neuroscience. Here I explore the many ways in which the circuits of the brain can process information at both the micro and meso scale. Understanding the way information is represented and processed in the brain could shed light on the neuropathology underlying complex neuropsychiatric diseases such as autism and schizophrenia. Chapter 2 establishes an experimental paradigm for assaying patterns of microcircuit activity and examines the role of dopaminergic modulation on prefrontal microcircuits. We find that dopamine type 2 (D2) receptor activation results in an increase in spontaneous activity while dopamine type 1 (D1) activation does not. Chapter 3 of this dissertation presents a study that illustrates how cholingergic activation normally produces what has been suggested as a neural substrate of attention; pairwise decorrelation in microcircuit activity. This study also shows that in two etiologicall distinct mouse models of autism, FMR1 knockout mice and Valproic Acid exposed mice, this ability to decorrelate in the presence of cholinergic activation is lost. This represents a putative microcircuit level biomarker of autism. Chapter 4 examines the structure/function relationship within the prefrontal microcircuit. Spontaneous activity in prefrontal microcircuits is shown to be organized according to a small world architecture. Interestingly, this architecture is important for one concrete function of neuronal microcircuits; the ability to produce temporally stereotyped patterns of activation. In the final chapter, we identify subnetworks in chronic intracranial electrocorticographic (ECoG) recordings using pairwise electrode coherence and dimensionality reduction techniques. We show that we can further reduce the dimensionality of these networks by identifying 'key-interactions' that are informative of the overall subnetwork state at any given point in time. This study highlights that redundancy in ECoG data can be exploited to identify low-dimensional representation of brain-wide subnetworks. Taken together, these studies represent the development of multiple technological and analytical techniques aimed at understanding how information is processed and modulated at emergent circuit and network levels as well as understanding their dysfunction in a neuropsychiatric disease state.

Testing the stress-gradient hypothesis during the restoration of tropical degraded land using the shrub Rhodomyrtus tomentosa as a nurse plant

Treesearch

Nan Liu; Hai Ren; Sufen Yuan; Qinfeng Guo; Long Yang

2013-01-01

The relative importance of facilitation and competition between pairwise plants across abiotic stress gradients as predicted by the stress-gradient hypothesis has been confirmed in arid and temperate ecosystems, but the hypothesis has rarely been tested in tropical systems, particularly across nutrient gradients. The current research examines the interactions between a...

Interaction of Phase Singularities on Spiral Wave Tail: Reconsideration of Capturing the Excitable Gap.

PubMed

Tomii, Naoki; Yamazaki, Masatoshi; Arafune, Tatsuhiko; Kamiya, Kaichiro; Nakazawa, Kazuo; Honjo, Haruo; Shibata, Nitaro; Sakuma, Ichiro

2018-03-09

The action mechanism of stimulation toward spiral waves (SWs) owing to the complex excitation patterns that occur just after point stimulation has not yet been experimentally clarified. This study sought to test our hypothesis that the effect of capturing excitable gap of SW by stimulation can also be explained as the interaction of original phase singularity (PS) and PSs induced by the stimulation on the wave tail (WT) of the original SW. Phase variance analysis was used to quantitatively analyze the post-shock PS trajectories. In a two-dimensional subepicardial layer of Langendorff-perfused rabbit hearts, optical mapping was utilized to record the excitation pattern during stimulation. After SW was induced by S1-S2 shock, single biphasic point stimulation S3 was applied. In 70 of the S1-S2-S3 stimulation episodes applied on six hearts, the original PS was clearly observed just before the S3 point stimulation in 37 episodes. Pairwise PSs were newly induced by the S3 in 20 episodes. The original PS collided with the newly-induced PSs in 16 episodes; otherwise, they did not interact with the original PS. SW shift occurred most efficiently when the S3 shock was applied at the relative refractory period, and PS shifted in the direction of WT. Quantitative tracking of PS clarified that stimulation in desirable conditions induces pairwise PSs on WT and that the collision of PSs causes SW shift along the WT. Results of this study indicate the importance of the interaction of shock-induced excitation with the WT for effective stimulation.

Numerical simulation of the pairwise interaction of deformable cells during migration in a microchannel

NASA Astrophysics Data System (ADS)

Lan, Hongzhi; Khismatullin, Damir B.

2014-07-01

Leukocytes and other circulating cells deform and move relatively to the channel flow in the lateral and translational directions. Their migratory property is important in immune response, hemostasis, cancer progression, delivery of nutrients, and microfluidic technologies such as cell separation and enrichment, and flow cytometry. Using our three-dimensional computational algorithm for multiphase viscoelastic flow, we have investigated the effect of pairwise interaction on the lateral and translational migration of circulating cells in a microchannel. The numerical simulation data show that when two cells with the same size and small separation distance interact, repulsive interaction take place until they reach the same lateral equilibrium position. During this process, they undergo swapping or passing, depending on the initial separation distance between each other. The threshold value of this distance increases with cell deformation, indicating that the cells experiencing larger deformation are more likely to swap. When a series of closely spaced cells with the same size are considered, they generally undergo damped oscillation in both lateral and translational directions until they reach equilibrium positions where they become evenly distributed in the flow direction (self-assembly phenomenon). A series of cells with a large lateral separation distance could collide repeatedly with each other, eventually crossing the centerline and entering the other side of the channel. For a series of cells with different deformability, more deformable cells, upon impact with less deformable cells, move to an equilibrium position closer to the centerline. The results of our study show that the bulk deformation of circulating cells plays a key role in their migration in a microchannel.

Mechanisms of action of Coxiella burnetii effectors inferred from host-pathogen protein interactions.

PubMed

Wallqvist, Anders; Wang, Hao; Zavaljevski, Nela; Memišević, Vesna; Kwon, Keehwan; Pieper, Rembert; Rajagopala, Seesandra V; Reifman, Jaques

2017-01-01

Coxiella burnetii is an obligate Gram-negative intracellular pathogen and the etiological agent of Q fever. Successful infection requires a functional Type IV secretion system, which translocates more than 100 effector proteins into the host cytosol to establish the infection, restructure the intracellular host environment, and create a parasitophorous vacuole where the replicating bacteria reside. We used yeast two-hybrid (Y2H) screening of 33 selected C. burnetii effectors against whole genome human and murine proteome libraries to generate a map of potential host-pathogen protein-protein interactions (PPIs). We detected 273 unique interactions between 20 pathogen and 247 human proteins, and 157 between 17 pathogen and 137 murine proteins. We used orthology to combine the data and create a single host-pathogen interaction network containing 415 unique interactions between 25 C. burnetii and 363 human proteins. We further performed complementary pairwise Y2H testing of 43 out of 91 C. burnetii-human interactions involving five pathogen proteins. We used the combined data to 1) perform enrichment analyses of target host cellular processes and pathways, 2) examine effectors with known infection phenotypes, and 3) infer potential mechanisms of action for four effectors with uncharacterized functions. The host-pathogen interaction profiles supported known Coxiella phenotypes, such as adapting cell morphology through cytoskeletal re-arrangements, protein processing and trafficking, organelle generation, cholesterol processing, innate immune modulation, and interactions with the ubiquitin and proteasome pathways. The generated dataset of PPIs-the largest collection of unbiased Coxiella host-pathogen interactions to date-represents a rich source of information with respect to secreted pathogen effector proteins and their interactions with human host proteins.

Triadic Non-Gaussian teleconnections in the Sea Surface Temperature Field: a source of interannual predictability coming from triadic wave resonances

NASA Astrophysics Data System (ADS)

Pires, Carlos; Trigo, Ricardo; Perdigão, Rui

2015-04-01

Analysis of centennial (1910-2012) time-series of the monthly Sea Surface Temperature anomalies (SSTAs) around the global ocean (extracted from the NOAA ERSST v3b dataset) shows clear evidence of non-Gaussian multivariate PDFs on certain projections, as an indication of both nonlinear correlations and nonlinear teleconnections. Beyond that, we still get statistical non-Gaussian relationships involving sets of three pair-wise uncorrelated variables through the occurrence of statistically significant and cross-validated triadic correlations (TCs),reaching ~30% in certain cases, i.e. non-null third-order cross cumulants between three standardized principal components (PCs) of the SSTA field, which would vanish under multivariate Gaussianity. Further enhanced TCs are obtained in the space of orthogonally rotated standardized PCs by expressing them as a function of the generalized Euler rotation angles and then maximized by gradient-descent methods. There are multiple triads depending of the embedding space of PCs where triads are sought. Furthermore they have no preferred order due to non-unique solutions of the non-linear matricial equations to be solved in the optimization. Triadic correlation is a particular form of the triadic interaction information, defined as the parcel of the mutual information (an Information-Theoretic measure of statistical dependency) which is atributed to triadic statistical synergies, not explained by pair-wise relationships. Spatial patterns of the triad's components generally exhibit wave-like structures in spatial quadrature and satisfying the triadic wave resonance condition. Examples of triads are given in spaces spanned by the leading EOFs of the SSTA field and projecting mostly in the Pacific Ocean (e.g. El Niño, Pacific Decadal Oscillation, North-Pacific Gyre Oscillation and pattrens of waves crossing the Pacific basin). A triadic correlation means a non-null Pearson correlation between the product of any two variables and the remaining third one. This nonlinear correlation may exhibit memory extending to months or years and may even be responsible for some skill recovery at the decadal scale. The triadic cumulant may de decomposed into Fourier cross bi-spectrum terms relying on components satisfying the triadic wave resonance. This holds when the frequency (in cycles per century) of a Fourier component is the sum of frequencies of the other two Fourier components. Therefore, dominant resonances between components interacting constructively, i.e. satisfying the appropriate phase relationship, can be considered as nonlinear sources of predictability on scales ranging from months to decades. The triads and indices derived from them can be used in schemes of long-range forecasting and downscaling.

Multisite Interactions in Lattice-Gas Models

NASA Astrophysics Data System (ADS)

Einstein, T. L.; Sathiyanarayanan, R.

For detailed applications of lattice-gas models to surface systems, multisite interactions often play at least as significant a role as interactions between pairs of adatoms that are separated by a few lattice spacings. We recall that trio (3-adatom, non-pairwise) interactions do not inevitably create phase boundary asymmetries about half coverage. We discuss a sophisticated application to an experimental system and describe refinements in extracting lattice-gas energies from calculations of total energies of several different ordered overlayers. We describe how lateral relaxations complicate matters when there is direct interaction between the adatoms, an issue that is important when examining the angular dependence of step line tensions. We discuss the connector model as an alternative viewpoint and close with a brief account of recent work on organic molecule overlayers.

Using the principle of entropy maximization to infer genetic interaction networks from gene expression patterns.

PubMed

Lezon, Timothy R; Banavar, Jayanth R; Cieplak, Marek; Maritan, Amos; Fedoroff, Nina V

2006-12-12

We describe a method based on the principle of entropy maximization to identify the gene interaction network with the highest probability of giving rise to experimentally observed transcript profiles. In its simplest form, the method yields the pairwise gene interaction network, but it can also be extended to deduce higher-order interactions. Analysis of microarray data from genes in Saccharomyces cerevisiae chemostat cultures exhibiting energy metabolic oscillations identifies a gene interaction network that reflects the intracellular communication pathways that adjust cellular metabolic activity and cell division to the limiting nutrient conditions that trigger metabolic oscillations. The success of the present approach in extracting meaningful genetic connections suggests that the maximum entropy principle is a useful concept for understanding living systems, as it is for other complex, nonequilibrium systems.

Enthalpic parameters of interaction between diglycylglycine and polyatomic alcohols in aqueous solutions

NASA Astrophysics Data System (ADS)

Mezhevoi, I. N.; Badelin, V. G.

2015-12-01

Integral enthalpies of solution Δsol H m of diglycylglycine in aqueous solutions of glycerol, ethylene glycol, and 1,2-propylene glycol are measured via solution calorimetry. The experimental data are used to calculate the standard enthalpies of solution (Δsol H°) and transfer (Δtr H°) of the tripeptide from water to aqueous solutions of polyatomic alcohols. The enthalpic pairwise coefficients h xy of interactions between the tripeptide and polyatomic alcohol molecules are calculated using the McMillan-Mayer solution theory and are found to have positive values. The findings are discussed using the theory of estimating various types of interactions in ternary systems and the effect the structural features of interacting biomolecules have on the thermochemical parameters of diglycylglycine dissolution.

SANSparallel: interactive homology search against Uniprot

PubMed Central

Somervuo, Panu; Holm, Liisa

2015-01-01

Proteins evolve by mutations and natural selection. The network of sequence similarities is a rich source for mining homologous relationships that inform on protein structure and function. There are many servers available to browse the network of homology relationships but one has to wait up to a minute for results. The SANSparallel webserver provides protein sequence database searches with immediate response and professional alignment visualization by third-party software. The output is a list, pairwise alignment or stacked alignment of sequence-similar proteins from Uniprot, UniRef90/50, Swissprot or Protein Data Bank. The stacked alignments are viewed in Jalview or as sequence logos. The database search uses the suffix array neighborhood search (SANS) method, which has been re-implemented as a client-server, improved and parallelized. The method is extremely fast and as sensitive as BLAST above 50% sequence identity. Benchmarks show that the method is highly competitive compared to previously published fast database search programs: UBLAST, DIAMOND, LAST, LAMBDA, RAPSEARCH2 and BLAT. The web server can be accessed interactively or programmatically at http://ekhidna2.biocenter.helsinki.fi/cgi-bin/sans/sans.cgi. It can be used to make protein functional annotation pipelines more efficient, and it is useful in interactive exploration of the detailed evidence supporting the annotation of particular proteins of interest. PMID:25855811

Human Retroviruses: Methods and Protocols

PubMed Central

Zhao, Gongpu; Zhang, Peijun

2015-01-01

Summary After virus fusion with a target cell, the viral core is released into the host cell cytoplasm and undergoes a controlled disassembly process, termed uncoating, before or as reverse transcription takes place. The cellular protein TRIM5α is a host cell restriction factor that blocks HIV-1 infection in rhesus macaque cells by targeting the viral capsid and inducing premature uncoating. The molecular mechanism of the interaction between capsid and TRIM5α remains unclear. Here, we describe an approach that utilizes cryo-electron microscopy (cryoEM) to examine the structural changes exerted on HIV-1 capsid (CA) assembly by TRIM5α binding. The TRIM5α interaction sites on CA assembly were further dissected by combining cryoEM with pair-wise cysteine mutations that crosslink CA either within a CA hexamer or between CA hexamers. Based on the structural information from cryoEM and crosslinking results from in vitro CA assemblies and purified intact HIV-1 cores, we demonstrate that direct binding of TRIM5α CC-SPRY domains to the viral capsid results in disruption and fragmentation of the surface lattice of HIV-1 capsid, specifically at inter-hexamer interfaces. The method described here can be easily adopted to study other important interactions in multi-protein complexes. PMID:24158810

Impaired inference in a case of developmental amnesia.

PubMed

D'Angelo, Maria C; Rosenbaum, R Shayna; Ryan, Jennifer D

2016-10-01

Amnesia is associated with impairments in relational memory, which is critically supported by the hippocampus. By adapting the transitivity paradigm, we previously showed that age-related impairments in inference were mitigated when judgments could be predicated on known pairwise relations, however, such advantages were not observed in the adult-onset amnesic case D.A. Here, we replicate and extend this finding in a developmental amnesic case (N.C.), who also shows impaired relational learning and transitive expression. Unlike D.A., N.C.'s damage affected the extended hippocampal system and diencephalic structures, and does not extend to neocortical areas that are affected in D.A. Critically, despite their differences in etiology and affected structures, N.C. and D.A. perform similarly on the task. N.C. showed intact pairwise knowledge, suggesting that he is able to use existing semantic information, but this semantic knowledge was insufficient to support transitive expression. The present results suggest a critical role for regions connected to the hippocampus and/or medial prefrontal cortex in inference beyond learning of pairwise relations. © 2016 The Authors Hippocampus Published by Wiley Periodicals, Inc. © 2016 The Authors. Wiley Periodicals, Inc.

The use of many-body expansions and geometry optimizations in fragment-based methods.

PubMed

Fedorov, Dmitri G; Asada, Naoya; Nakanishi, Isao; Kitaura, Kazuo

2014-09-16

Conspectus Chemists routinely work with complex molecular systems: solutions, biochemical molecules, and amorphous and composite materials provide some typical examples. The questions one often asks are what are the driving forces for a chemical phenomenon? How reasonable are our views of chemical systems in terms of subunits, such as functional groups and individual molecules? How can one quantify the difference in physicochemical properties of functional units found in a different chemical environment? Are various effects on functional units in molecular systems additive? Can they be represented by pairwise potentials? Are there effects that cannot be represented in a simple picture of pairwise interactions? How can we obtain quantitative values for these effects? Many of these questions can be formulated in the language of many-body effects. They quantify the properties of subunits (fragments), referred to as one-body properties, pairwise interactions (two-body properties), couplings of two-body interactions described by three-body properties, and so on. By introducing the notion of fragments in the framework of quantum chemistry, one obtains two immense benefits: (a) chemists can finally relate to quantum chemistry, which now speaks their language, by discussing chemically interesting subunits and their interactions and (b) calculations become much faster due to a reduced computational scaling. For instance, the somewhat academic sounding question of the importance of three-body effects in water clusters is actually another way of asking how two hydrogen bonds affect each other, when they involve three water molecules. One aspect of this is the many-body charge transfer (CT), because the charge transfers in the two hydrogen bonds are coupled to each other (not independent). In this work, we provide a generalized view on the use of many-body expansions in fragment-based methods, focusing on the general aspects of the property expansion and a contraction of a many-body expansion in a formally two-body series, as exemplified in the development of the fragment molecular orbital (FMO) method. Fragment-based methods have been very successful in delivering the properties of fragments, as well as the fragment interactions, providing insights into complex chemical processes in large molecular systems. We briefly review geometry optimizations performed with fragment-based methods and present an efficient geometry optimization method based on the combination of FMO with molecular mechanics (MM), applied to the complex of a subunit of protein kinase 2 (CK2) with a ligand. FMO results are discussed in comparison with experimental and MM-optimized structures.

Event-chain Monte Carlo algorithms for three- and many-particle interactions

NASA Astrophysics Data System (ADS)

Harland, J.; Michel, M.; Kampmann, T. A.; Kierfeld, J.

2017-02-01

We generalize the rejection-free event-chain Monte Carlo algorithm from many-particle systems with pairwise interactions to systems with arbitrary three- or many-particle interactions. We introduce generalized lifting probabilities between particles and obtain a general set of equations for lifting probabilities, the solution of which guarantees maximal global balance. We validate the resulting three-particle event-chain Monte Carlo algorithms on three different systems by comparison with conventional local Monte Carlo simulations: i) a test system of three particles with a three-particle interaction that depends on the enclosed triangle area; ii) a hard-needle system in two dimensions, where needle interactions constitute three-particle interactions of the needle end points; iii) a semiflexible polymer chain with a bending energy, which constitutes a three-particle interaction of neighboring chain beads. The examples demonstrate that the generalization to many-particle interactions broadens the applicability of event-chain algorithms considerably.

Simultaneous learning of instantaneous and time-delayed genetic interactions using novel information theoretic scoring technique

PubMed Central

2012-01-01

Background Understanding gene interactions is a fundamental question in systems biology. Currently, modeling of gene regulations using the Bayesian Network (BN) formalism assumes that genes interact either instantaneously or with a certain amount of time delay. However in reality, biological regulations, both instantaneous and time-delayed, occur simultaneously. A framework that can detect and model both these two types of interactions simultaneously would represent gene regulatory networks more accurately. Results In this paper, we introduce a framework based on the Bayesian Network (BN) formalism that can represent both instantaneous and time-delayed interactions between genes simultaneously. A novel scoring metric having firm mathematical underpinnings is also proposed that, unlike other recent methods, can score both interactions concurrently and takes into account the reality that multiple regulators can regulate a gene jointly, rather than in an isolated pair-wise manner. Further, a gene regulatory network (GRN) inference method employing an evolutionary search that makes use of the framework and the scoring metric is also presented. Conclusion By taking into consideration the biological fact that both instantaneous and time-delayed regulations can occur among genes, our approach models gene interactions with greater accuracy. The proposed framework is efficient and can be used to infer gene networks having multiple orders of instantaneous and time-delayed regulations simultaneously. Experiments are carried out using three different synthetic networks (with three different mechanisms for generating synthetic data) as well as real life networks of Saccharomyces cerevisiae, E. coli and cyanobacteria gene expression data. The results show the effectiveness of our approach. PMID:22691450

Assessing the performance of dispersionless and dispersion-accounting methods: helium interaction with cluster models of the TiO2(110) surface.

PubMed

de Lara-Castells, María Pilar; Stoll, Hermann; Mitrushchenkov, Alexander O

2014-08-21

As a prototypical dispersion-dominated physisorption problem, we analyze here the performance of dispersionless and dispersion-accounting methodologies on the helium interaction with cluster models of the TiO2(110) surface. A special focus has been given to the dispersionless density functional dlDF and the dlDF+Das construction for the total interaction energy (K. Pernal, R. Podeswa, K. Patkowski, and K. Szalewicz, Phys. Rev. Lett. 2009, 109, 263201), where Das is an effective interatomic pairwise functional form for the dispersion. Likewise, the performance of symmetry-adapted perturbation theory (SAPT) method is evaluated, where the interacting monomers are described by density functional theory (DFT) with the dlDF, PBE, and PBE0 functionals. Our benchmarks include CCSD(T)-F12b calculations and comparative analysis on the nuclear bound states supported by the He-cluster potentials. Moreover, intra- and intermonomer correlation contributions to the physisorption interaction are analyzed through the method of increments (H. Stoll, J. Chem. Phys. 1992, 97, 8449) at the CCSD(T) level of theory. This method is further applied in conjunction with a partitioning of the Hartree-Fock interaction energy to estimate individual interaction energy components, comparing them with those obtained using the different SAPT(DFT) approaches. The cluster size evolution of dispersionless and dispersion-accounting energy components is then discussed, revealing the reduced role of the dispersionless interaction and intramonomer correlation when the extended nature of the surface is better accounted for. On the contrary, both post-Hartree-Fock and SAPT(DFT) results clearly demonstrate the high-transferability character of the effective pairwise dispersion interaction whatever the cluster model is. Our contribution also illustrates how the method of increments can be used as a valuable tool not only to achieve the accuracy of CCSD(T) calculations using large cluster models but also to evaluate the performance of SAPT(DFT) methods for the physically well-defined contributions to the total interaction energy. Overall, our work indicates the excellent performance of a dlDF+Das approach in which the parameters are optimized using the smallest cluster model of the target surface to treat van der Waals adsorbate-surface interactions.

KECSA-Movable Type Implicit Solvation Model (KMTISM)

PubMed Central

2015-01-01

Computation of the solvation free energy for chemical and biological processes has long been of significant interest. The key challenges to effective solvation modeling center on the choice of potential function and configurational sampling. Herein, an energy sampling approach termed the “Movable Type” (MT) method, and a statistical energy function for solvation modeling, “Knowledge-based and Empirical Combined Scoring Algorithm” (KECSA) are developed and utilized to create an implicit solvation model: KECSA-Movable Type Implicit Solvation Model (KMTISM) suitable for the study of chemical and biological systems. KMTISM is an implicit solvation model, but the MT method performs energy sampling at the atom pairwise level. For a specific molecular system, the MT method collects energies from prebuilt databases for the requisite atom pairs at all relevant distance ranges, which by its very construction encodes all possible molecular configurations simultaneously. Unlike traditional statistical energy functions, KECSA converts structural statistical information into categorized atom pairwise interaction energies as a function of the radial distance instead of a mean force energy function. Within the implicit solvent model approximation, aqueous solvation free energies are then obtained from the NVT ensemble partition function generated by the MT method. Validation is performed against several subsets selected from the Minnesota Solvation Database v2012. Results are compared with several solvation free energy calculation methods, including a one-to-one comparison against two commonly used classical implicit solvation models: MM-GBSA and MM-PBSA. Comparison against a quantum mechanics based polarizable continuum model is also discussed (Cramer and Truhlar’s Solvation Model 12). PMID:25691832

Colony-level variation in pollen collection and foraging preferences among wild-caught bumble bees (Hymenoptera: Apidae).

PubMed

Saifuddin, Mustafa; Jha, Shalene

2014-04-01

Given that many pollinators have exhibited dramatic declines related to habitat destruction, an improved understanding of pollinator resource collection across human-altered landscapes is essential to conservation efforts. Despite the importance of bumble bees (Bombus spp.) as global pollinators, little is known regarding how pollen collection patterns vary between individuals, colonies, and landscapes. In this study, Vosnesensky bumble bees (Bombus vosnesenskii Radoszkowski) were collected from a range of human-altered and natural landscapes in northern California. Extensive vegetation surveys and Geographic Information System (GIS)-based habitat classifications were conducted at each site, bees were genotyped to identify colony mates, and pollen loads were examined to identify visited plants. In contrast to predictions based on strong competitive interactions, pollen load composition was significantly more similar for bees captured in a shared study region compared with bees throughout the research area but was not significantly more similar for colony mates. Preference analyses revealed that pollen loads were not composed of the most abundant plant species per study region. The majority of ranked pollen preference lists were significantly correlated for pairwise comparisons of colony mates and individuals within a study region, whereas the majority of pairwise comparisons of ranked pollen preference lists between individuals located at separate study regions were uncorrelated. Results suggest that pollen load composition and foraging preferences are similar for bees throughout a shared landscape regardless of colony membership. The importance of native plant species in pollen collection is illustrated through preference analyses, and we suggest prioritization of specific rare native plant species for enhanced bumble bee pollen collection.

Structure based alignment and clustering of proteins (STRALCP)

DOEpatents

Zemla, Adam T.; Zhou, Carol E.; Smith, Jason R.; Lam, Marisa W.

2013-06-18

Disclosed are computational methods of clustering a set of protein structures based on local and pair-wise global similarity values. Pair-wise local and global similarity values are generated based on pair-wise structural alignments for each protein in the set of protein structures. Initially, the protein structures are clustered based on pair-wise local similarity values. The protein structures are then clustered based on pair-wise global similarity values. For each given cluster both a representative structure and spans of conserved residues are identified. The representative protein structure is used to assign newly-solved protein structures to a group. The spans are used to characterize conservation and assign a "structural footprint" to the cluster.

A water market simulator considering pair-wise trades between agents

NASA Astrophysics Data System (ADS)

Huskova, I.; Erfani, T.; Harou, J. J.

2012-04-01

In many basins in England no further water abstraction licences are available. Trading water between water rights holders has been recognized as a potentially effective and economically efficient strategy to mitigate increasing scarcity. A screening tool that could assess the potential for trade through realistic simulation of individual water rights holders would help assess the solution's potential contribution to local water management. We propose an optimisation-driven water market simulator that predicts pair-wise trade in a catchment and represents its interaction with natural hydrology and engineered infrastructure. A model is used to emulate licence-holders' willingness to engage in short-term trade transactions. In their simplest form agents are represented using an economic benefit function. The working hypothesis is that trading behaviour can be partially predicted based on differences in marginal values of water over space and time and estimates of transaction costs on pair-wise trades. We discuss the further possibility of embedding rules, norms and preferences of the different water user sectors to more realistically represent the behaviours, motives and constraints of individual licence holders. The potential benefits and limitations of such a social simulation (agent-based) approach is contrasted with our simulator where agents are driven by economic optimization. A case study based on the Dove River Basin (UK) demonstrates model inputs and outputs. The ability of the model to suggest impacts of water rights policy reforms on trading is discussed.

Strategies for informed sample size reduction in adaptive controlled clinical trials

NASA Astrophysics Data System (ADS)

Arandjelović, Ognjen

2017-12-01

Clinical trial adaptation refers to any adjustment of the trial protocol after the onset of the trial. The main goal is to make the process of introducing new medical interventions to patients more efficient. The principal challenge, which is an outstanding research problem, is to be found in the question of how adaptation should be performed so as to minimize the chance of distorting the outcome of the trial. In this paper, we propose a novel method for achieving this. Unlike most of the previously published work, our approach focuses on trial adaptation by sample size adjustment, i.e. by reducing the number of trial participants in a statistically informed manner. Our key idea is to select the sample subset for removal in a manner which minimizes the associated loss of information. We formalize this notion and describe three algorithms which approach the problem in different ways, respectively, using (i) repeated random draws, (ii) a genetic algorithm, and (iii) what we term pair-wise sample compatibilities. Experiments on simulated data demonstrate the effectiveness of all three approaches, with a consistently superior performance exhibited by the pair-wise sample compatibilities-based method.

Automatic Plagiarism Detection with PAIRwise 2.0

ERIC Educational Resources Information Center

Knight, Allan; Almeroth, Kevin

2011-01-01

As part of the research carried out at the University of California, Santa Barbara's Center for Information Technology and Society (CITS), the Paper Authentication and Integrity Research (PAIR) project was launched. We began by investigating how one recent technology affected student learning outcomes. One aspect of this research was to study the…

Crystal structure of 2-oxopyrrolidin-3-yl 4-(2-phenyl-diazen-1-yl)benzoate.

PubMed

Elkin, Igor; Maris, Thierry; Melkoumov, Alexandre; Hildgen, Patrice; Banquy, Xavier; Leclair, Grégoire; Barrett, Christopher

2018-04-01

In the title compound, C 17 H 15 N 3 O 3 , the plane of the pyrrolidone ring is inclined at an angle of 59.791 (2)° to that of the azo-benzene segment, which adopts a configuration close to planar. In the crystal, mol-ecules are oriented pairwise by (2-oxopyrrolidin-3-yl)-oxy moieties at an angle of 76.257 (3)°, linked by hydrogen bonds and π-stacking inter-actions, forming zigzag supra-molecular chains parallel to [010] further linked via additional C-H⋯π inter-actions.

CREDO: a structural interactomics database for drug discovery

PubMed Central

Schreyer, Adrian M.; Blundell, Tom L.

2013-01-01

CREDO is a unique relational database storing all pairwise atomic interactions of inter- as well as intra-molecular contacts between small molecules and macromolecules found in experimentally determined structures from the Protein Data Bank. These interactions are integrated with further chemical and biological data. The database implements useful data structures and algorithms such as cheminformatics routines to create a comprehensive analysis platform for drug discovery. The database can be accessed through a web-based interface, downloads of data sets and web services at http://www-cryst.bioc.cam.ac.uk/credo. Database URL: http://www-cryst.bioc.cam.ac.uk/credo PMID:23868908

Inverse Ising Inference Using All the Data

NASA Astrophysics Data System (ADS)

Aurell, Erik; Ekeberg, Magnus

2012-03-01

We show that a method based on logistic regression, using all the data, solves the inverse Ising problem far better than mean-field calculations relying only on sample pairwise correlation functions, while still computationally feasible for hundreds of nodes. The largest improvement in reconstruction occurs for strong interactions. Using two examples, a diluted Sherrington-Kirkpatrick model and a two-dimensional lattice, we also show that interaction topologies can be recovered from few samples with good accuracy and that the use of l1 regularization is beneficial in this process, pushing inference abilities further into low-temperature regimes.

Cell and Particle Interactions and Aggregation During Electrophoretic Motion

NASA Technical Reports Server (NTRS)

Wang, Hua; Zeng, Shulin; Loewenberg, Michael; Todd, Paul; Davis, Robert H.

1996-01-01

The stability and pairwise aggregation rates of small spherical particles under the collective effects of buoyancy-driven motion and electrophoretic migration are analyzed. The particles are assumed to be non-Brownian, with thin double-layers and different zeta potentials. The particle aggregation rates may be enhanced or reduced, respectively, by parallel and antiparallel alignments of the buoyancy-driven and electrophoretic velocities. For antiparallel alignments, with the buoyancy-driven relative velocity exceeding the electrophoretic relative velocity between two widely-separated particles, there is a 'collision-forbidden region' in parameter space due to hydrodynamic interactions; thus, the suspension becomes stable against aggregation.

Allowing for uncertainty due to missing continuous outcome data in pairwise and network meta-analysis.

PubMed

Mavridis, Dimitris; White, Ian R; Higgins, Julian P T; Cipriani, Andrea; Salanti, Georgia

2015-02-28

Missing outcome data are commonly encountered in randomized controlled trials and hence may need to be addressed in a meta-analysis of multiple trials. A common and simple approach to deal with missing data is to restrict analysis to individuals for whom the outcome was obtained (complete case analysis). However, estimated treatment effects from complete case analyses are potentially biased if informative missing data are ignored. We develop methods for estimating meta-analytic summary treatment effects for continuous outcomes in the presence of missing data for some of the individuals within the trials. We build on a method previously developed for binary outcomes, which quantifies the degree of departure from a missing at random assumption via the informative missingness odds ratio. Our new model quantifies the degree of departure from missing at random using either an informative missingness difference of means or an informative missingness ratio of means, both of which relate the mean value of the missing outcome data to that of the observed data. We propose estimating the treatment effects, adjusted for informative missingness, and their standard errors by a Taylor series approximation and by a Monte Carlo method. We apply the methodology to examples of both pairwise and network meta-analysis with multi-arm trials. © 2014 The Authors. Statistics in Medicine Published by John Wiley & Sons Ltd.

A discrete model of Ostwald ripening based on multiple pairwise interactions

NASA Astrophysics Data System (ADS)

Di Nunzio, Paolo Emilio

2018-06-01

A discrete multi-particle model of Ostwald ripening based on direct pairwise interactions is developed for particles with incoherent interfaces as an alternative to the classical LSW mean field theory. The rate of matter exchange depends on the average surface-to-surface interparticle distance, a characteristic feature of the system which naturally incorporates the effect of volume fraction of second phase. The multi-particle diffusion is described through the definition of an interaction volume containing all the particles involved in the exchange of solute. At small volume fractions this is proportional to the size of the central particle, at higher volume fractions it gradually reduces as a consequence of diffusion screening described on a geometrical basis. The topological noise present in real systems is also included. For volume fractions below about 0.1 the model predicts broad and right-skewed stationary size distributions resembling a lognormal function. Above this value, a transition to sharper, more symmetrical but still right-skewed shapes occurs. An excellent agreement with experiments is obtained for 3D particle size distributions of solid-solid and solid-liquid systems with volume fraction 0.07, 0.30, 0.52 and 0.74. The kinetic constant of the model depends on the cube root of volume fraction up to about 0.1, then increases rapidly with an upward concavity. It is in good agreement with the available literature data on solid-liquid mixtures in the volume fraction range from 0.20 to about 0.75.

Improving pairwise comparison of protein sequences with domain co-occurrence

PubMed Central

Gascuel, Olivier

2018-01-01

Comparing and aligning protein sequences is an essential task in bioinformatics. More specifically, local alignment tools like BLAST are widely used for identifying conserved protein sub-sequences, which likely correspond to protein domains or functional motifs. However, to limit the number of false positives, these tools are used with stringent sequence-similarity thresholds and hence can miss several hits, especially for species that are phylogenetically distant from reference organisms. A solution to this problem is then to integrate additional contextual information to the procedure. Here, we propose to use domain co-occurrence to increase the sensitivity of pairwise sequence comparisons. Domain co-occurrence is a strong feature of proteins, since most protein domains tend to appear with a limited number of other domains on the same protein. We propose a method to take this information into account in a typical BLAST analysis and to construct new domain families on the basis of these results. We used Plasmodium falciparum as a case study to evaluate our method. The experimental findings showed an increase of 14% of the number of significant BLAST hits and an increase of 25% of the proteome area that can be covered with a domain. Our method identified 2240 new domains for which, in most cases, no model of the Pfam database could be linked. Moreover, our study of the quality of the new domains in terms of alignment and physicochemical properties show that they are close to that of standard Pfam domains. Source code of the proposed approach and supplementary data are available at: https://gite.lirmm.fr/menichelli/pairwise-comparison-with-cooccurrence PMID:29293498

Delineating slowly and rapidly evolving fractions of the Drosophila genome.

PubMed

Keith, Jonathan M; Adams, Peter; Stephen, Stuart; Mattick, John S

2008-05-01

Evolutionary conservation is an important indicator of function and a major component of bioinformatic methods to identify non-protein-coding genes. We present a new Bayesian method for segmenting pairwise alignments of eukaryotic genomes while simultaneously classifying segments into slowly and rapidly evolving fractions. We also describe an information criterion similar to the Akaike Information Criterion (AIC) for determining the number of classes. Working with pairwise alignments enables detection of differences in conservation patterns among closely related species. We analyzed three whole-genome and three partial-genome pairwise alignments among eight Drosophila species. Three distinct classes of conservation level were detected. Sequences comprising the most slowly evolving component were consistent across a range of species pairs, and constituted approximately 62-66% of the D. melanogaster genome. Almost all (>90%) of the aligned protein-coding sequence is in this fraction, suggesting much of it (comprising the majority of the Drosophila genome, including approximately 56% of non-protein-coding sequences) is functional. The size and content of the most rapidly evolving component was species dependent, and varied from 1.6% to 4.8%. This fraction is also enriched for protein-coding sequence (while containing significant amounts of non-protein-coding sequence), suggesting it is under positive selection. We also classified segments according to conservation and GC content simultaneously. This analysis identified numerous sub-classes of those identified on the basis of conservation alone, but was nevertheless consistent with that classification. Software, data, and results available at www.maths.qut.edu.au/-keithj/. Genomic segments comprising the conservation classes available in BED format.

Emory University: High-Throughput Protein-Protein Interaction Dataset for Lung Cancer-Associated Genes | Office of Cancer Genomics

Cancer.gov

To discover novel PPI signaling hubs for lung cancer, CTD2 Center at Emory utilized large-scale genomics datasets and literature to compile a set of lung cancer-associated genes. A library of expression vectors were generated for these genes and utilized for detecting pairwise PPIs with cell lysate-based TR-FRET assays in high-throughput screening format. Read the abstract.

A Hybrid Physics-Based Data-Driven Approach for Point-Particle Force Modeling

NASA Astrophysics Data System (ADS)

Moore, Chandler; Akiki, Georges; Balachandar, S.

2017-11-01

This study improves upon the physics-based pairwise interaction extended point-particle (PIEP) model. The PIEP model leverages a physical framework to predict fluid mediated interactions between solid particles. While the PIEP model is a powerful tool, its pairwise assumption leads to increased error in flows with high particle volume fractions. To reduce this error, a regression algorithm is used to model the differences between the current PIEP model's predictions and the results of direct numerical simulations (DNS) for an array of monodisperse solid particles subjected to various flow conditions. The resulting statistical model and the physical PIEP model are superimposed to construct a hybrid, physics-based data-driven PIEP model. It must be noted that the performance of a pure data-driven approach without the model-form provided by the physical PIEP model is substantially inferior. The hybrid model's predictive capabilities are analyzed using more DNS. In every case tested, the hybrid PIEP model's prediction are more accurate than those of physical PIEP model. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship Program under Grant No. DGE-1315138 and the U.S. DOE, NNSA, ASC Program, as a Cooperative Agreement under Contract No. DE-NA0002378.

Multidrug resistant pathogens respond differently to the presence of co-pathogen, commensal, probiotic and host cells.

PubMed

Chan, Agnes P; Choi, Yongwook; Brinkac, Lauren M; Krishnakumar, Radha; DePew, Jessica; Kim, Maria; Hinkle, Mary K; Lesho, Emil P; Fouts, Derrick E

2018-06-05

In light of the ongoing antimicrobial resistance crisis, there is a need to understand the role of co-pathogens, commensals, and the local microbiome in modulating virulence and antibiotic resistance. To identify possible interactions that influence the expression of virulence or survival mechanisms in both the multidrug-resistant organisms (MDROs) and human host cells, unique cohorts of clinical isolates were selected for whole genome sequencing with enhanced assembly and full annotation, pairwise co-culturing, and transcriptome profiling. The MDROs were co-cultured in pairwise combinations either with: (1) another MDRO, (2) skin commensals (Staphylococcus epidermidis and Corynebacterium jeikeium), (3) the common probiotic Lactobacillus reuteri, and (4) human fibroblasts. RNA-Seq analysis showed distinct regulation of virulence and antimicrobial resistance gene responses across different combinations of MDROs, commensals, and human cells. Co-culture assays demonstrated that microbial interactions can modulate gene responses of both the target and pathogen/commensal species, and that the responses are specific to the identity of the pathogen/commensal species. In summary, bacteria have mechanisms to distinguish between friends, foe and host cells. These results provide foundational data and insight into the possibility of manipulating the local microbiome when treating complicated polymicrobial wound, intra-abdominal, or respiratory infections.

Enabling task-based information prioritization via semantic web encodings

NASA Astrophysics Data System (ADS)

Michaelis, James R.

2016-05-01

Modern Soldiers rely upon accurate and actionable information technology to achieve mission objectives. While increasingly rich sensor networks for Areas of Operation (AO) can offer many directions for aiding Soldiers, limitations are imposed by current tactical edge systems on the rate that content can be transmitted. Furthermore, mission tasks will often require very specific sets of information which may easily be drowned out by other content sources. Prior research on Quality and Value of Information (QoI/VoI) has aimed to define ways to prioritize information objects based on their intrinsic attributes (QoI) and perceived value to a consumer (VoI). As part of this effort, established ranking approaches for obtaining Subject Matter Expert (SME) recommendations, such as the Analytic Hierarchy Process (AHP) have been considered. However, limited work has been done to tie Soldier context - such as descriptions of their mission and tasks - back to intrinsic attributes of information objects. As a first step toward addressing the above challenges, this work introduces an ontology-backed approach - rooted in Semantic Web publication practices - for expressing both AHP decision hierarchies and corresponding SME feedback. Following a short discussion on related QoI/VoI research, an ontology-based data structure is introduced for supporting evaluation of Information Objects, using AHP rankings designed to facilitate information object prioritization. Consistent with alternate AHP approaches, prioritization in this approach is based on pairwise comparisons between Information Objects with respect to established criteria, as well as on pairwise comparison of the criteria to assess their relative importance. The paper concludes with a discussion of both ongoing and future work.

Pairwise-interaction extended point-particle model for particle-laden flows

NASA Astrophysics Data System (ADS)

Akiki, G.; Moore, W. C.; Balachandar, S.

2017-12-01

In this work we consider the pairwise interaction extended point-particle (PIEP) model for Euler-Lagrange simulations of particle-laden flows. By accounting for the precise location of neighbors the PIEP model goes beyond local particle volume fraction, and distinguishes the influence of upstream, downstream and laterally located neighbors. The two main ingredients of the PIEP model are (i) the undisturbed flow at any particle is evaluated as a superposition of the macroscale flow and a microscale flow that is approximated as a pairwise superposition of perturbation fields induced by each of the neighboring particles, and (ii) the forces and torque on the particle are then calculated from the undisturbed flow using the Faxén form of the force relation. The computational efficiency of the standard Euler-Lagrange approach is retained, since the microscale perturbation fields induced by a neighbor are pre-computed and stored as PIEP maps. Here we extend the PIEP force model of Akiki et al. [3] with a corresponding torque model to systematically include the effect of perturbation fields induced by the neighbors in evaluating the net torque. Also, we use DNS results from a uniform flow over two stationary spheres to further improve the PIEP force and torque models. We then test the PIEP model in three different sedimentation problems and compare the results against corresponding DNS to assess the accuracy of the PIEP model and improvement over the standard point-particle approach. In the case of two sedimenting spheres in a quiescent ambient the PIEP model is shown to capture the drafting-kissing-tumbling process. In cases of 5 and 80 sedimenting spheres a good agreement is obtained between the PIEP simulation and the DNS. For all three simulations, the DEM-PIEP was able to recreate, to a good extent, the results from the DNS, while requiring only a negligible fraction of the numerical resources required by the fully-resolved DNS.

Detection of protein-protein interactions by ribosome display and protein in situ immobilisation.

PubMed

He, Mingyue; Liu, Hong; Turner, Martin; Taussig, Michael J

2009-12-31

We describe a method for identification of protein-protein interactions by combining two cell-free protein technologies, namely ribosome display and protein in situ immobilisation. The method requires only PCR fragments as the starting material, the target proteins being made through cell-free protein synthesis, either associated with their encoding mRNA as ribosome complexes or immobilised on a solid surface. The use of ribosome complexes allows identification of interacting protein partners from their attached coding mRNA. To demonstrate the procedures, we have employed the lymphocyte signalling proteins Vav1 and Grb2 and confirmed the interaction between Grb2 and the N-terminal SH3 domain of Vav1. The method has promise for library screening of pairwise protein interactions, down to the analytical level of individual domain or motif mapping.

Using the principle of entropy maximization to infer genetic interaction networks from gene expression patterns

PubMed Central

Lezon, Timothy R.; Banavar, Jayanth R.; Cieplak, Marek; Maritan, Amos; Fedoroff, Nina V.

2006-01-01

We describe a method based on the principle of entropy maximization to identify the gene interaction network with the highest probability of giving rise to experimentally observed transcript profiles. In its simplest form, the method yields the pairwise gene interaction network, but it can also be extended to deduce higher-order interactions. Analysis of microarray data from genes in Saccharomyces cerevisiae chemostat cultures exhibiting energy metabolic oscillations identifies a gene interaction network that reflects the intracellular communication pathways that adjust cellular metabolic activity and cell division to the limiting nutrient conditions that trigger metabolic oscillations. The success of the present approach in extracting meaningful genetic connections suggests that the maximum entropy principle is a useful concept for understanding living systems, as it is for other complex, nonequilibrium systems. PMID:17138668

Structure of the N-terminal domain of human thioredoxin-interacting protein.

PubMed

Polekhina, Galina; Ascher, David Benjamin; Kok, Shie Foong; Beckham, Simone; Wilce, Matthew; Waltham, Mark

2013-03-01

Thioredoxin-interacting protein (TXNIP) is one of the six known α-arrestins and has recently received considerable attention owing to its involvement in redox signalling and metabolism. Various stress stimuli such as high glucose, heat shock, UV, H2O2 and mechanical stress among others robustly induce the expression of TXNIP, resulting in the sequestration and inactivation of thioredoxin, which in turn leads to cellular oxidative stress. While TXNIP is the only α-arrestin known to bind thioredoxin, TXNIP and two other α-arrestins, Arrdc4 and Arrdc3, have been implicated in metabolism. Furthermore, owing to its roles in the pathologies of diabetes and cardiovascular disease, TXNIP is considered to be a promising drug target. Based on their amino-acid sequences, TXNIP and the other α-arrestins are remotely related to β-arrestins. Here, the crystal structure of the N-terminal domain of TXNIP is reported. It provides the first structural information on any of the α-arrestins and reveals that although TXNIP adopts a β-arrestin fold as predicted, it is structurally more similar to Vps26 proteins than to β-arrestins, while sharing below 15% pairwise sequence identity with either.

Acetaldehyde dissociates the PTP1B–E-cadherin–β-catenin complex in Caco-2 cell monolayers by a phosphorylation-dependent mechanism

PubMed Central

Sheth, Parimal; Seth, Ankur; Atkinson, Katherine J.; Gheyi, Tarun; Kale, Gautam; Giorgianni, Francesco; Desiderio, Dominic M.; Li, Chunying; Naren, Anjaparavanda; Rao, Radhakrishna

2006-01-01

Interactions between E-cadherin, β-catenin and PTP1B (protein tyrosine phosphatase 1B) are crucial for the organization of AJs (adherens junctions) and epithelial cell–cell adhesion. In the present study, the effect of acetaldehyde on the AJs and on the interactions between E-cadherin, β-catenin and PTP1B was determined in Caco-2 cell monolayers. Treatment of cell monolayers with acetaldehyde induced redistribution of E-cadherin and β-catenin from the intercellular junctions by a tyrosine phosphorylation-dependent mechanism. The PTPase activity associated with E-cadherin and β-catenin was significantly reduced and the interaction of PTP1B with E-cadherin and β-catenin was attenuated by acetaldehyde. Acetaldehyde treatment resulted in phosphorylation of β-catenin on tyrosine residues, and abolished the interaction of β-catenin with E-cadherin by a tyrosine kinase-dependent mechanism. Protein binding studies showed that the treatment of cells with acetaldehyde reduced the binding of β-catenin to the C-terminal region of E-cadherin. Pairwise binding studies using purified proteins indicated that the direct interaction between E-cadherin and β-catenin was reduced by tyrosine phosphorylation of β-catenin, but was unaffected by tyrosine phosphorylation of E-cadherin-C. Treatment of cells with acetaldehyde also reduced the binding of E-cadherin to GST (glutathione S-transferase)–PTP1B. The pairwise binding study showed that GST–E-cadherin-C binds to recombinant PTP1B, but this binding was significantly reduced by tyrosine phosphorylation of E-cadherin. Acetaldehyde increased the phosphorylation of β-catenin on Tyr-331, Tyr-333, Tyr-654 and Tyr-670. These results show that acetaldehyde induces disruption of interactions between E-cadherin, β-catenin and PTP1B by a phosphorylation-dependent mechanism. PMID:17087658

A Study of the Use of Pairwise Comparison in the Context of Social Online Moderation

ERIC Educational Resources Information Center

Tarricone, Pina; Newhouse, C. Paul

2016-01-01

Traditional moderation of student assessments is often carried out with groups of teachers working face-to-face in a specified location making judgements concerning the quality of representations of achievement. This traditional model has relied little on modern information communications technologies and has been logistically challenging. We…

Non-random Patterns in the Distribution of NOR-bearing Chromosome Territories in Human Fibroblasts: A Network Model of Interactions

PubMed Central

Pliss, Artem; Fritz, Andrew J.; Stojkovic, Branislav; Ding, Hu; Mukherjee, Lopamudra; Bhattacharya, Sambit; Xu, Jinhui; Berezney, Ronald

2017-01-01

We present a 3-D mapping in WI38 human diploid fibroblast cells of chromosome territories (CT) 13,14,15,21, and 22, which contain the nucleolar organizing regions (NOR) and participate in the formation of nucleoli. The nuclear radial positioning of NOR-CT correlated with the size of chromosomes with smaller CT more interior. A high frequency of pairwise associations between NOR-CT ranging from 52% (CT13-21) to 82% (CT15-21) was detected as well as a triplet arrangement of CT15-21-22 (72%). The associations of homologous CT were significantly lower (24–36%). The arrangements of each pairwise CT varied from CT13-14 and CT13-22, which had a majority of cells with single associations, to CT13-15 and CT13-21 where a majority of cells had multiple interactions. In cells with multiple nucleoli, one of the nucleoli (termed “dominant”) always associated with a higher number of CT. Moreover, certain CT pairs more frequently contributed to the same nucleolus than to others. This nonrandom pattern suggests that a large number of the NOR-chromsomes are poised in close proximity during the postmitotic nucleolar recovery and through their NORs may contribute to the formation of the same nucleolus. A global data mining program termed the chromatic median determined the most probable interchromosomal arrangement of the entire NOR-CT population. This interactive network model was significantly above randomized simulation and was composed of 13 connections among the NOR-CT. We conclude that the NOR-CT form a global interactive network in the cell nucleus that may be a fundamental feature for the regulation of nucleolar and other genomic functions. PMID:25077974

Communicating with sentences: A multi-word naming game model

NASA Astrophysics Data System (ADS)

Lou, Yang; Chen, Guanrong; Hu, Jianwei

2018-01-01

Naming game simulates the process of naming an object by a single word, in which a population of communicating agents can reach global consensus asymptotically through iteratively pair-wise conversations. We propose an extension of the single-word model to a multi-word naming game (MWNG), simulating the case of describing a complex object by a sentence (multiple words). Words are defined in categories, and then organized as sentences by combining them from different categories. We refer to a formatted combination of several words as a pattern. In such an MWNG, through a pair-wise conversation, it requires the hearer to achieve consensus with the speaker with respect to both every single word in the sentence as well as the sentence pattern, so as to guarantee the correct meaning of the saying; otherwise, they fail reaching consensus in the interaction. We validate the model in three typical topologies as the underlying communication network, and employ both conventional and man-designed patterns in performing the MWNG.

Pairwise frictional profile between particles determines discontinuous shear thickening transition in non-colloidal suspensions

PubMed Central

Comtet, Jean; Chatté, Guillaume; Niguès, Antoine; Bocquet, Lydéric; Siria, Alessandro; Colin, Annie

2017-01-01

The process by which sheared suspensions go through a dramatic change in viscosity is known as discontinuous shear thickening. Although well-characterized on the macroscale, the microscopic mechanisms at play in this transition are still poorly understood. Here, by developing new experimental procedures based on quartz-tuning fork atomic force microscopy, we measure the pairwise frictional profile between approaching pairs of polyvinyl chloride and cornstarch particles in solvent. We report a clear transition from a low-friction regime, where pairs of particles support a finite normal load, while interacting purely hydrodynamically, to a high-friction regime characterized by hard repulsive contact between the particles and sliding friction. Critically, we show that the normal stress needed to enter the frictional regime at nanoscale matches the critical stress at which shear thickening occurs for macroscopic suspensions. Our experiments bridge nano and macroscales and provide long needed demonstration of the role of frictional forces in discontinuous shear thickening. PMID:28561032

Pairwise frictional profile between particles determines discontinuous shear thickening transition in non-colloidal suspensions.

PubMed

Comtet, Jean; Chatté, Guillaume; Niguès, Antoine; Bocquet, Lydéric; Siria, Alessandro; Colin, Annie

2017-05-31

The process by which sheared suspensions go through a dramatic change in viscosity is known as discontinuous shear thickening. Although well-characterized on the macroscale, the microscopic mechanisms at play in this transition are still poorly understood. Here, by developing new experimental procedures based on quartz-tuning fork atomic force microscopy, we measure the pairwise frictional profile between approaching pairs of polyvinyl chloride and cornstarch particles in solvent. We report a clear transition from a low-friction regime, where pairs of particles support a finite normal load, while interacting purely hydrodynamically, to a high-friction regime characterized by hard repulsive contact between the particles and sliding friction. Critically, we show that the normal stress needed to enter the frictional regime at nanoscale matches the critical stress at which shear thickening occurs for macroscopic suspensions. Our experiments bridge nano and macroscales and provide long needed demonstration of the role of frictional forces in discontinuous shear thickening.

InterEvDock2: an expanded server for protein docking using evolutionary and biological information from homology models and multimeric inputs.

PubMed

Quignot, Chloé; Rey, Julien; Yu, Jinchao; Tufféry, Pierre; Guerois, Raphaël; Andreani, Jessica

2018-05-08

Computational protein docking is a powerful strategy to predict structures of protein-protein interactions and provides crucial insights for the functional characterization of macromolecular cross-talks. We previously developed InterEvDock, a server for ab initio protein docking based on rigid-body sampling followed by consensus scoring using physics-based and statistical potentials, including the InterEvScore function specifically developed to incorporate co-evolutionary information in docking. InterEvDock2 is a major evolution of InterEvDock which allows users to submit input sequences - not only structures - and multimeric inputs and to specify constraints for the pairwise docking process based on previous knowledge about the interaction. For this purpose, we added modules in InterEvDock2 for automatic template search and comparative modeling of the input proteins. The InterEvDock2 pipeline was benchmarked on 812 complexes for which unbound homology models of the two partners and co-evolutionary information are available in the PPI4DOCK database. InterEvDock2 identified a correct model among the top 10 consensus in 29% of these cases (compared to 15-24% for individual scoring functions) and at least one correct interface residue among 10 predicted in 91% of these cases. InterEvDock2 is thus a unique protein docking server, designed to be useful for the experimental biology community. The InterEvDock2 web interface is available at http://bioserv.rpbs.univ-paris-diderot.fr/services/InterEvDock2/.

Discovery and identification of a series of alkyl decalin isomers in petroleum geological samples.

PubMed

Wang, Huitong; Zhang, Shuichang; Weng, Na; Zhang, Bin; Zhu, Guangyou; Liu, Lingyan

2015-07-07

The comprehensive two-dimensional gas chromatography/time-of-flight mass spectrometry (GC × GC/TOFMS) has been used to characterize a crude oil and a source rock extract sample. During the process, a series of pairwise components between monocyclic alkanes and mono-aromatics have been discovered. After tentative assignments of decahydronaphthalene isomers, a series of alkyl decalin isomers have been synthesized and used for identification and validation of these petroleum compounds. From both the MS and chromatography information, these pairwise compounds were identified as 2-alkyl-decahydronaphthalenes and 1-alkyl-decahydronaphthalenes. The polarity of 1-alkyl-decahydronaphthalenes was stronger. Their long chain alkyl substituent groups may be due to bacterial transformation or different oil cracking events. This systematic profiling of alkyl-decahydronaphthalene isomers provides further understanding and recognition of these potential petroleum biomarkers.

Preliminary Classification of Novel Hemorrhagic Fever-Causing Viruses Using Sequence-Based PAirwise Sequence Comparison (PASC) Analysis.

PubMed

Bào, Yīmíng; Kuhn, Jens H

2018-01-01

During the last decade, genome sequence-based classification of viruses has become increasingly prominent. Viruses can be even classified based on coding-complete genome sequence data alone. Nevertheless, classification remains arduous as experts are required to establish phylogenetic trees to depict the evolutionary relationships of such sequences for preliminary taxonomic placement. Pairwise sequence comparison (PASC) of genomes is one of several novel methods for establishing relationships among viruses. This method, provided by the US National Center for Biotechnology Information as an open-access tool, circumvents phylogenetics, and yet PASC results are often in agreement with those of phylogenetic analyses. Computationally inexpensive, PASC can be easily performed by non-taxonomists. Here we describe how to use the PASC tool for the preliminary classification of novel viral hemorrhagic fever-causing viruses.

Non-rigid multi-frame registration of cell nuclei in live cell fluorescence microscopy image data.

PubMed

Tektonidis, Marco; Kim, Il-Han; Chen, Yi-Chun M; Eils, Roland; Spector, David L; Rohr, Karl

2015-01-01

The analysis of the motion of subcellular particles in live cell microscopy images is essential for understanding biological processes within cells. For accurate quantification of the particle motion, compensation of the motion and deformation of the cell nucleus is required. We introduce a non-rigid multi-frame registration approach for live cell fluorescence microscopy image data. Compared to existing approaches using pairwise registration, our approach exploits information from multiple consecutive images simultaneously to improve the registration accuracy. We present three intensity-based variants of the multi-frame registration approach and we investigate two different temporal weighting schemes. The approach has been successfully applied to synthetic and live cell microscopy image sequences, and an experimental comparison with non-rigid pairwise registration has been carried out. Copyright © 2014 Elsevier B.V. All rights reserved.

Including non-additive genetic effects in Bayesian methods for the prediction of genetic values based on genome-wide markers

PubMed Central

2011-01-01

Background Molecular marker information is a common source to draw inferences about the relationship between genetic and phenotypic variation. Genetic effects are often modelled as additively acting marker allele effects. The true mode of biological action can, of course, be different from this plain assumption. One possibility to better understand the genetic architecture of complex traits is to include intra-locus (dominance) and inter-locus (epistasis) interaction of alleles as well as the additive genetic effects when fitting a model to a trait. Several Bayesian MCMC approaches exist for the genome-wide estimation of genetic effects with high accuracy of genetic value prediction. Including pairwise interaction for thousands of loci would probably go beyond the scope of such a sampling algorithm because then millions of effects are to be estimated simultaneously leading to months of computation time. Alternative solving strategies are required when epistasis is studied. Methods We extended a fast Bayesian method (fBayesB), which was previously proposed for a purely additive model, to include non-additive effects. The fBayesB approach was used to estimate genetic effects on the basis of simulated datasets. Different scenarios were simulated to study the loss of accuracy of prediction, if epistatic effects were not simulated but modelled and vice versa. Results If 23 QTL were simulated to cause additive and dominance effects, both fBayesB and a conventional MCMC sampler BayesB yielded similar results in terms of accuracy of genetic value prediction and bias of variance component estimation based on a model including additive and dominance effects. Applying fBayesB to data with epistasis, accuracy could be improved by 5% when all pairwise interactions were modelled as well. The accuracy decreased more than 20% if genetic variation was spread over 230 QTL. In this scenario, accuracy based on modelling only additive and dominance effects was generally superior to that of the complex model including epistatic effects. Conclusions This simulation study showed that the fBayesB approach is convenient for genetic value prediction. Jointly estimating additive and non-additive effects (especially dominance) has reasonable impact on the accuracy of prediction and the proportion of genetic variation assigned to the additive genetic source. PMID:21867519

Will kinematic Sunyaev-Zel'dovich measurements enhance the science return from galaxy redshift surveys?

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sugiyama, Naonori S.; Okumura, Teppei; Spergel, David N., E-mail: nao.s.sugiyama@gmail.com, E-mail: tokumura@asiaa.sinica.edu.tw, E-mail: dns@astro.princeton.edu

2017-01-01

Yes. Future CMB experiments such as Advanced ACTPol and CMB-S4 should achieve measurements with S/N of > 0.1 for the typical host halo of galaxies in redshift surveys. These measurements will provide complementary measurements of the growth rate of large scale structure f and the expansion rate of the Universe H to galaxy clustering measurements. This paper emphasizes that there is significant information in the anisotropy of the relative pairwise kSZ measurements. We expand the relative pairwise kSZ power spectrum in Legendre polynomials and consider up to its octopole. Assuming that the noise in the filtered maps is uncorrelated betweenmore » the positions of galaxies in the survey, we derive a simple analytic form for the power spectrum covariance of the relative pairwise kSZ temperature in redshift space. While many previous studies have assumed optimistically that the optical depth of the galaxies τ{sub T} in the survey is known, we marginalize over τ{sub T}, to compute constraints on the growth rate f and the expansion rate H . For realistic survey parameters, we find that combining kSZ and galaxy redshift survey data reduces the marginalized 1-σ errors on H and f to ∼50-70% compared to the galaxy-only analysis.« less

Will kinematic Sunyaev-Zel'dovich measurements enhance the science return from galaxy redshift surveys?

NASA Astrophysics Data System (ADS)

Sugiyama, Naonori S.; Okumura, Teppei; Spergel, David N.

2017-01-01

Yes. Future CMB experiments such as Advanced ACTPol and CMB-S4 should achieve measurements with S/N of > 0.1 for the typical host halo of galaxies in redshift surveys. These measurements will provide complementary measurements of the growth rate of large scale structure f and the expansion rate of the Universe H to galaxy clustering measurements. This paper emphasizes that there is significant information in the anisotropy of the relative pairwise kSZ measurements. We expand the relative pairwise kSZ power spectrum in Legendre polynomials and consider up to its octopole. Assuming that the noise in the filtered maps is uncorrelated between the positions of galaxies in the survey, we derive a simple analytic form for the power spectrum covariance of the relative pairwise kSZ temperature in redshift space. While many previous studies have assumed optimistically that the optical depth of the galaxies τT in the survey is known, we marginalize over τT, to compute constraints on the growth rate f and the expansion rate H. For realistic survey parameters, we find that combining kSZ and galaxy redshift survey data reduces the marginalized 1-σ errors on H and f to ~50-70% compared to the galaxy-only analysis.

Pairwise-Comparison Software

NASA Technical Reports Server (NTRS)

Ricks, Wendell R.

1995-01-01

Pairwise comparison (PWC) is computer program that collects data for psychometric scaling techniques now used in cognitive research. It applies technique of pairwise comparisons, which is one of many techniques commonly used to acquire the data necessary for analyses. PWC administers task, collects data from test subject, and formats data for analysis. Written in Turbo Pascal v6.0.

A systems biology model of the regulatory network in Populus leaves reveals interacting regulators and conserved regulation

PubMed Central

2011-01-01

Background Green plant leaves have always fascinated biologists as hosts for photosynthesis and providers of basic energy to many food webs. Today, comprehensive databases of gene expression data enable us to apply increasingly more advanced computational methods for reverse-engineering the regulatory network of leaves, and to begin to understand the gene interactions underlying complex emergent properties related to stress-response and development. These new systems biology methods are now also being applied to organisms such as Populus, a woody perennial tree, in order to understand the specific characteristics of these species. Results We present a systems biology model of the regulatory network of Populus leaves. The network is reverse-engineered from promoter information and expression profiles of leaf-specific genes measured over a large set of conditions related to stress and developmental. The network model incorporates interactions between regulators, such as synergistic and competitive relationships, by evaluating increasingly more complex regulatory mechanisms, and is therefore able to identify new regulators of leaf development not found by traditional genomics methods based on pair-wise expression similarity. The approach is shown to explain available gene function information and to provide robust prediction of expression levels in new data. We also use the predictive capability of the model to identify condition-specific regulation as well as conserved regulation between Populus and Arabidopsis. Conclusions We outline a computationally inferred model of the regulatory network of Populus leaves, and show how treating genes as interacting, rather than individual, entities identifies new regulators compared to traditional genomics analysis. Although systems biology models should be used with care considering the complexity of regulatory programs and the limitations of current genomics data, methods describing interactions can provide hypotheses about the underlying cause of emergent properties and are needed if we are to identify target genes other than those constituting the "low hanging fruit" of genomic analysis. PMID:21232107

Synergetic and Redundant Information Flow Detected by Unnormalized Granger Causality: Application to Resting State fMRI.

PubMed

Stramaglia, Sebastiano; Angelini, Leonardo; Wu, Guorong; Cortes, Jesus M; Faes, Luca; Marinazzo, Daniele

2016-12-01

We develop a framework for the analysis of synergy and redundancy in the pattern of information flow between subsystems of a complex network. The presence of redundancy and/or synergy in multivariate time series data renders difficulty to estimate the neat flow of information from each driver variable to a given target. We show that adopting an unnormalized definition of Granger causality, one may put in evidence redundant multiplets of variables influencing the target by maximizing the total Granger causality to a given target, over all the possible partitions of the set of driving variables. Consequently, we introduce a pairwise index of synergy which is zero when two independent sources additively influence the future state of the system, differently from previous definitions of synergy. We report the application of the proposed approach to resting state functional magnetic resonance imaging data from the Human Connectome Project showing that redundant pairs of regions arise mainly due to space contiguity and interhemispheric symmetry, while synergy occurs mainly between nonhomologous pairs of regions in opposite hemispheres. Redundancy and synergy, in healthy resting brains, display characteristic patterns, revealed by the proposed approach. The pairwise synergy index, here introduced, maps the informational character of the system at hand into a weighted complex network: the same approach can be applied to other complex systems whose normal state corresponds to a balance between redundant and synergetic circuits.

Computer-Based Readability Testing of Information Booklets for German Cancer Patients.

PubMed

Keinki, Christian; Zowalla, Richard; Pobiruchin, Monika; Huebner, Jutta; Wiesner, Martin

2018-04-12

Understandable health information is essential for treatment adherence and improved health outcomes. For readability testing, several instruments analyze the complexity of sentence structures, e.g., Flesch-Reading Ease (FRE) or Vienna-Formula (WSTF). Moreover, the vocabulary is of high relevance for readers. The aim of this study is to investigate the agreement of sentence structure and vocabulary-based (SVM) instruments. A total of 52 freely available German patient information booklets on cancer were collected from the Internet. The mean understandability level L was computed for 51 booklets. The resulting values of FRE, WSTF, and SVM were assessed pairwise for agreement with Bland-Altman plots and two-sided, paired t tests. For the pairwise comparison, the mean L values are L FRE  = 6.81, L WSTF  = 7.39, L SVM  = 5.09. The sentence structure-based metrics gave significantly different scores (P < 0.001) for all assessed booklets, confirmed by the Bland-Altman analysis. The study findings suggest that vocabulary-based instruments cannot be interchanged with FRE/WSTF. However, both analytical aspects should be considered and checked by authors to linguistically refine texts with respect to the individual target group. Authors of health information can be supported by automated readability analysis. Health professionals can benefit by direct booklet comparisons allowing for time-effective selection of suitable booklets for patients.

Evolution of biological sequences implies an extreme value distribution of type I for both global and local pairwise alignment scores.

PubMed

Bastien, Olivier; Maréchal, Eric

2008-08-07

Confidence in pairwise alignments of biological sequences, obtained by various methods such as Blast or Smith-Waterman, is critical for automatic analyses of genomic data. Two statistical models have been proposed. In the asymptotic limit of long sequences, the Karlin-Altschul model is based on the computation of a P-value, assuming that the number of high scoring matching regions above a threshold is Poisson distributed. Alternatively, the Lipman-Pearson model is based on the computation of a Z-value from a random score distribution obtained by a Monte-Carlo simulation. Z-values allow the deduction of an upper bound of the P-value (1/Z-value2) following the TULIP theorem. Simulations of Z-value distribution is known to fit with a Gumbel law. This remarkable property was not demonstrated and had no obvious biological support. We built a model of evolution of sequences based on aging, as meant in Reliability Theory, using the fact that the amount of information shared between an initial sequence and the sequences in its lineage (i.e., mutual information in Information Theory) is a decreasing function of time. This quantity is simply measured by a sequence alignment score. In systems aging, the failure rate is related to the systems longevity. The system can be a machine with structured components, or a living entity or population. "Reliability" refers to the ability to operate properly according to a standard. Here, the "reliability" of a sequence refers to the ability to conserve a sufficient functional level at the folded and maturated protein level (positive selection pressure). Homologous sequences were considered as systems 1) having a high redundancy of information reflected by the magnitude of their alignment scores, 2) which components are the amino acids that can independently be damaged by random DNA mutations. From these assumptions, we deduced that information shared at each amino acid position evolved with a constant rate, corresponding to the information hazard rate, and that pairwise sequence alignment scores should follow a Gumbel distribution, which parameters could find some theoretical rationale. In particular, one parameter corresponds to the information hazard rate. Extreme value distribution of alignment scores, assessed from high scoring segments pairs following the Karlin-Altschul model, can also be deduced from the Reliability Theory applied to molecular sequences. It reflects the redundancy of information between homologous sequences, under functional conservative pressure. This model also provides a link between concepts of biological sequence analysis and of systems biology.

Spatial interactions of yarded White-tailed Deer, Odocoileus virginianus

USGS Publications Warehouse

Nelson, M.E.; Sargeant, G.A.

2008-01-01

We examined the spatial interactions of nine female White-tailed Deer (Odocoileus virginianus) in two deeryards (winter aggregations) in northeastern Minnesota during February-April 1999. Global positioning system (GPS) collars yielded seven pair-wise comparisons of deer that were located at the same time (???1 minute apart) and mat used overlapping areas. Deer traveled separately and did not associate with one another. Within overlapping areas, comparisons of distances between deer and distances between random locations indicated deer moved without regard to each other. Similarly, comparisons of observed and expected probabilities of deer using areas overlapping those of other deer also evinced that deer moved independently.

Towards photonic quantum simulation of ground states of frustrated Heisenberg spin systems

PubMed Central

Ma, Xiao-song; Dakić, Borivoje; Kropatschek, Sebastian; Naylor, William; Chan, Yang-hao; Gong, Zhe-xuan; Duan, Lu-ming; Zeilinger, Anton; Walther, Philip

2014-01-01

Photonic quantum simulators are promising candidates for providing insight into other small- to medium-sized quantum systems. Recent experiments have shown that photonic quantum systems have the advantage to exploit quantum interference for the quantum simulation of the ground state of Heisenberg spin systems. Here we experimentally characterize this quantum interference at a tuneable beam splitter and further investigate the measurement-induced interactions of a simulated four-spin system by comparing the entanglement dynamics using pairwise concurrence. We also study theoretically a four-site square lattice with next-nearest neighbor interactions and a six-site checkerboard lattice, which might be in reach of current technology. PMID:24394808

Generalized thermodynamics of phase equilibria in scalar active matter

NASA Astrophysics Data System (ADS)

Solon, Alexandre P.; Stenhammar, Joakim; Cates, Michael E.; Kafri, Yariv; Tailleur, Julien

2018-02-01

Motility-induced phase separation (MIPS) arises generically in fluids of self-propelled particles when interactions lead to a kinetic slowdown at high densities. Starting from a continuum description of scalar active matter akin to a generalized Cahn-Hilliard equation, we give a general prescription for the mean densities of coexisting phases in flux-free steady states that amounts, at a hydrodynamics scale, to extremizing an effective free energy. We illustrate our approach on two well-known models: self-propelled particles interacting either through a density-dependent propulsion speed or via direct pairwise forces. Our theory accounts quantitatively for their phase diagrams, providing a unified description of MIPS.

Probabilistic visual and electromagnetic data fusion for robust drift-free sequential mosaicking: application to fetoscopy

PubMed Central

Tella-Amo, Marcel; Peter, Loic; Shakir, Dzhoshkun I.; Deprest, Jan; Iglesias, Juan Eugenio; Ourselin, Sebastien

2018-01-01

Abstract. The most effective treatment for twin-to-twin transfusion syndrome is laser photocoagulation of the shared vascular anastomoses in the placenta. Vascular connections are extremely challenging to locate due to their caliber and the reduced field-of-view of the fetoscope. Therefore, mosaicking techniques are beneficial to expand the scene, facilitate navigation, and allow vessel photocoagulation decision-making. Local vision-based mosaicking algorithms inherently drift over time due to the use of pairwise transformations. We propose the use of an electromagnetic tracker (EMT) sensor mounted at the tip of the fetoscope to obtain camera pose measurements, which we incorporate into a probabilistic framework with frame-to-frame visual information to achieve globally consistent sequential mosaics. We parametrize the problem in terms of plane and camera poses constrained by EMT measurements to enforce global consistency while leveraging pairwise image relationships in a sequential fashion through the use of local bundle adjustment. We show that our approach is drift-free and performs similarly to state-of-the-art global alignment techniques like bundle adjustment albeit with much less computational burden. Additionally, we propose a version of bundle adjustment that uses EMT information. We demonstrate the robustness to EMT noise and loss of visual information and evaluate mosaics for synthetic, phantom-based and ex vivo datasets. PMID:29487889

SANSparallel: interactive homology search against Uniprot.

PubMed

Somervuo, Panu; Holm, Liisa

2015-07-01

Proteins evolve by mutations and natural selection. The network of sequence similarities is a rich source for mining homologous relationships that inform on protein structure and function. There are many servers available to browse the network of homology relationships but one has to wait up to a minute for results. The SANSparallel webserver provides protein sequence database searches with immediate response and professional alignment visualization by third-party software. The output is a list, pairwise alignment or stacked alignment of sequence-similar proteins from Uniprot, UniRef90/50, Swissprot or Protein Data Bank. The stacked alignments are viewed in Jalview or as sequence logos. The database search uses the suffix array neighborhood search (SANS) method, which has been re-implemented as a client-server, improved and parallelized. The method is extremely fast and as sensitive as BLAST above 50% sequence identity. Benchmarks show that the method is highly competitive compared to previously published fast database search programs: UBLAST, DIAMOND, LAST, LAMBDA, RAPSEARCH2 and BLAT. The web server can be accessed interactively or programmatically at http://ekhidna2.biocenter.helsinki.fi/cgi-bin/sans/sans.cgi. It can be used to make protein functional annotation pipelines more efficient, and it is useful in interactive exploration of the detailed evidence supporting the annotation of particular proteins of interest. © The Author(s) 2015. Published by Oxford University Press on behalf of Nucleic Acids Research.

Trapped atoms along nanophotonic resonators

NASA Astrophysics Data System (ADS)

Fields, Brian; Kim, May; Chang, Tzu-Han; Hung, Chen-Lung

2017-04-01

Many-body systems subject to long-range interactions have remained a very challenging topic experimentally. Ultracold atoms trapped in extreme proximity to the surface of nanophotonic structures provides a dynamic system combining the strong atom-atom interactions mediated by guided mode photons with the exquisite control implemented with trapped atom systems. The hybrid system promises pair-wise tunability of long-range interactions between atomic pseudo spins, allowing studies of quantum magnetism extending far beyond nearest neighbor interactions. In this talk, we will discuss our current status developing high quality nanophotonic ring resonators, engineered on CMOS compatible optical chips with integrated nanostructures that, in combination with a side illuminating beam, can realize stable atom traps approximately 100nm above the surface. We will report on our progress towards loading arrays of cold atoms near the surface of these structures and studying atom-atom interaction mediated by photons with high cooperativity.

Scaling laws for van der Waals interactions in nanostructured materials.

PubMed

Gobre, Vivekanand V; Tkatchenko, Alexandre

2013-01-01

Van der Waals interactions have a fundamental role in biology, physics and chemistry, in particular in the self-assembly and the ensuing function of nanostructured materials. Here we utilize an efficient microscopic method to demonstrate that van der Waals interactions in nanomaterials act at distances greater than typically assumed, and can be characterized by different scaling laws depending on the dimensionality and size of the system. Specifically, we study the behaviour of van der Waals interactions in single-layer and multilayer graphene, fullerenes of varying size, single-wall carbon nanotubes and graphene nanoribbons. As a function of nanostructure size, the van der Waals coefficients follow unusual trends for all of the considered systems, and deviate significantly from the conventionally employed pairwise-additive picture. We propose that the peculiar van der Waals interactions in nanostructured materials could be exploited to control their self-assembly.

On stable Pareto laws in a hierarchical model of economy

NASA Astrophysics Data System (ADS)

Chebotarev, A. M.

2007-01-01

This study considers a model of the income distribution of agents whose pairwise interaction is asymmetric and price-invariant. Asymmetric transactions are typical for chain-trading groups who arrange their business such that commodities move from senior to junior partners and money moves in the opposite direction. The price-invariance of transactions means that the probability of a pairwise interaction is a function of the ratio of incomes, which is independent of the price scale or absolute income level. These two features characterize the hierarchical model. The income distribution in this class of models is a well-defined double-Pareto function, which possesses Pareto tails for the upper and lower incomes. For gross and net upper incomes, the model predicts definite values of the Pareto exponents, agross and anet, which are stable with respect to quantitative variation of the pair-interaction. The Pareto exponents are also stable with respect to the choice of a demand function within two classes of status-dependent behavior of agents: linear demand ( agross=1, anet=2) and unlimited slowly varying demand ( agross=anet=1). For the sigmoidal demand that describes limited returns, agross=anet=1+α, with some α>0 satisfying a transcendental equation. The low-income distribution may be singular or vanishing in the neighborhood of the minimal income; in any case, it is L1-integrable and its Pareto exponent is given explicitly. The theory used in the present study is based on a simple balance equation and new results from multiplicative Markov chains and exponential moments of random geometric progressions.

The influence of arene-ring size on stacking interaction with canonical base pairs

NASA Astrophysics Data System (ADS)

Formánek, Martin; Burda, Jaroslav V.

2014-04-01

Stacking interactions between aromatic molecules (benzene, p-cymene, biphenyl, and di- and tetra-hydrogen anthracene) and G.C and A.T canonical Watson-Crick (WC) base pairs are explored. Two functionals with dispersion corrections: ω-B97XD and B3LYP-D3 are used. For a comparison also the MP2 and B3LYP-D3/PCM methods were used for the most stable p-cymene…WC geometries. It was found that the stacking interaction increases with the size of π-conjugation system. Its extent is in agreement with experimental finding on anticancer activity of Ru(II) piano-stool complexes where intercalation of these aromatic molecules should play an important role. The explored structures are considered as ternary system so that decomposition of the interaction energy to pairwise and non-additivity contributions is also examined.

Energetics of the molecular interactions of L-alanine and L-serine with xylitol, D-sorbitol, and D-mannitol in aqueous solutions at 298.15 K

NASA Astrophysics Data System (ADS)

Mezhevoi, I. N.; Badelin, V. G.

2013-04-01

Integral enthalpies of dissolution Δsol H m of L-alanine and L-serine are measured via the calorimetry of dissolution in aqueous solutions of xylitol, D-sorbitol, and D-mannitol. Standard enthalpies of dissolution (Δsol H ○) and the transfer (Δtr H ○) of amino acids from water to binary solvent are calculated from the experimental data. Using the McMillan-Mayer theory, enthalpy coefficients of pairwise interactions h xy of amino acids with molecules of polyols are calculated that are negative. The obtained results are discussed within the theory of the prevalence of different types of interactions in mixed solutions and the effect of the structural features of interacting biomolecules on the thermochemical parameters of dissolution of amino acids.

How direct competition shapes coexistence and vaccine effects in multi-strain pathogen systems.

PubMed

Gjini, Erida; Valente, Carina; Sá-Leão, Raquel; Gomes, M Gabriela M

2016-01-07

We describe an integrated modeling framework for understanding strain coexistence in polymorphic pathogen systems. Previous studies have debated the utility of neutral formulations and focused on cross-immunity between strains as a major stabilizing mechanism. Here we convey that direct competition for colonization mediates stable coexistence only when competitive abilities amongst pathogen clones satisfy certain pairwise asymmetries. We illustrate our ideas with nested SIS models of single and dual colonization, applied to polymorphic pneumococcal bacteria. By fitting the models to cross-sectional prevalence data from Portugal (before and after the introduction of a seven-valent pneumococcal conjugate vaccine), we are able to not only statistically compare neutral and non-neutral epidemiological formulations, but also estimate vaccine efficacy, transmission and competition parameters simultaneously. Our study highlights that the response of polymorphic pathogen populations to interventions holds crucial information about strain interactions, which can be extracted by suitable nested modeling. Copyright © 2015 Elsevier Ltd. All rights reserved.

Exploring blocking assays using Octet, ProteOn, and Biacore biosensors.

PubMed

Abdiche, Yasmina N; Malashock, Dan S; Pinkerton, Alanna; Pons, Jaume

2009-03-15

We demonstrate the use of label-free real-time optical biosensors in competitive binding assays by epitope binning a panel of antibodies. We describe three assay orientations that we term in tandem, premix, and classical sandwich blocking, and we perform each of them on three platforms: ForteBio's Octet QK, Bio-Rad's ProteOn XPR36, and GE Healthcare's Biacore 3000. By testing whether antibodies block one another's binding to their antigen in a pairwise fashion, we establish a blocking profile for each antibody relative to the others in the panel. The blocking information is then used to create "bins" of antibodies with similar epitopes. The advantages and disadvantages of each biosensor, factors to consider when deciding on the most appropriate blocking assay orientation for a particular interaction system, and tips for dealing with ambiguous data are discussed. The data from our different assay orientations and biosensors agree very well, establishing these machines as valuable tools for characterizing antibody epitopes and multiprotein complexes of biological significance.

A Case Study in Integrating Multiple E-commerce Standards via Semantic Web Technology

NASA Astrophysics Data System (ADS)

Yu, Yang; Hillman, Donald; Setio, Basuki; Heflin, Jeff

Internet business-to-business transactions present great challenges in merging information from different sources. In this paper we describe a project to integrate four representative commercial classification systems with the Federal Cataloging System (FCS). The FCS is used by the US Defense Logistics Agency to name, describe and classify all items under inventory control by the DoD. Our approach uses the ECCMA Open Technical Dictionary (eOTD) as a common vocabulary to accommodate all different classifications. We create a semantic bridging ontology between each classification and the eOTD to describe their logical relationships in OWL DL. The essential idea is that since each classification has formal definitions in a common vocabulary, we can use subsumption to automatically integrate them, thus mitigating the need for pairwise mappings. Furthermore our system provides an interactive interface to let users choose and browse the results and more importantly it can translate catalogs that commit to these classifications using compiled mapping results.

Development of a methodology to compute solvation free energies on the basis of the theory of energy representation for solutions represented with a polarizable force field.

PubMed

Suzuoka, Daiki; Takahashi, Hideaki; Ishiyama, Tatsuya; Morita, Akihiro

2012-12-07

We have developed a method of molecular simulations utilizing a polarizable force field in combination with the theory of energy representation (ER) for the purpose of establishing an efficient and accurate methodology to compute solvation free energies. The standard version of the ER method is, however, based on the assumption that the solute-solvent interaction is pairwise additive for its construction. A crucial step in the present method is to introduce an intermediate state in the solvation process to treat separately the many-body interaction associated with the polarizable model. The intermediate state is chosen so that the solute-solvent interaction can be formally written in the pairwise form, though the solvent molecules are interacting with each other with polarizable charges dependent on the solvent configuration. It is, then, possible to extract the free energy contribution δμ due to the many-body interaction between solute and solvent from the total solvation free energy Δμ. It is shown that the free energy δμ can be computed by an extension of the recent development implemented in quantum mechanical∕molecular mechanical simulations. To assess the numerical robustness of the approach, we computed the solvation free energies of a water and a methanol molecule in water solvent, where two paths for the solvation processes were examined by introducing different intermediate states. The solvation free energies of a water molecule associated with the two paths were obtained as -5.3 and -5.8 kcal∕mol. Those of a methanol molecule were determined as -3.5 and -3.7 kcal∕mol. These results of the ER simulations were also compared with those computed by a numerically exact approach. It was demonstrated that the present approach produces the solvation free energies in comparable accuracies to simulations of thermodynamic integration (TI) method within a tenth of computational time used for the TI simulations.

Epistatic SNP interaction of ERCC6 with ERCC8 and their joint protein expression contribute to gastric cancer/atrophic gastritis risk.

PubMed

Jing, Jing-Jing; Lu, You-Zhu; Sun, Li-Ping; Liu, Jing-Wei; Gong, Yue-Hua; Xu, Qian; Dong, Nan-Nan; Yuan, Yuan

2017-06-27

Excision repair cross-complementing group 6 and 8 (ERCC6 and ERCC8) are two indispensable genes for the initiation of transcription-coupled nucleotide excision repair pathway. This study aimed to evaluate the interactions between single nucleotide polymorphisms of ERCC6 (rs1917799) and ERCC8 (rs158572 and rs158916) in gastric cancer and its precancerous diseases. Besides, protein level analysis were performed to compare ERCC6 and ERCC8 expression in different stages of gastric diseases, and to correlate SNPs jointly with gene expression. Sequenom MassARRAY platform method was used to detect polymorphisms of ERCC6 and ERCC8 in 1916 subjects. In situ ERCC6 and ERCC8 protein expression were detected by immunohistochemistry in 109 chronic superficial gastritis, 109 chronic atrophic gastritis and 109 gastric cancer cases. Our results demonstrated pairwise epistatic interactions between ERCC6 and ERCC8 SNPs that ERCC6 rs1917799-ERCC8 rs158572 combination was associated with decreased risk of chronic atrophic gastritis and increased risk of gastric cancer. ERCC6 rs1917799 also showed a significant interaction with ERCC8 rs158916 to reduce gastric cancer risk. The expressions of ERCC6, ERCC8 and ERCC6-ERCC8 combination have similarities that higher positivity was observed in chronic superficial gastritis compared with chronic atrophic gastritis and gastric cancer. As for the effects of ERCC6 and ERCC8 SNPs on the protein expression, single SNP had no correlation with corresponding gene expression, whereas the ERCC6 rs1917799-ERCC8 rs158572 pair had significant influence on ERCC6 and ERCC6-ERCC8 expression. In conclusion, ERCC6 rs1917799, ERCC8 rs158572 and rs158916 demonstrated pairwise epistatic interactions to associate with chronic atrophic gastritis and gastric cancer risk. The ERCC6 rs1917799-ERCC8 rs158572 pair significantly influence ERCC6 and ERCC6-ERCC8 expression.

Polar Desolvation and Position 226 of Pancreatic and Neutrophil Elastases Are Crucial to their Affinity for the Kunitz-Type Inhibitors ShPI-1 and ShPI-1/K13L.

PubMed

Hernández González, Jorge Enrique; García-Fernández, Rossana; Valiente, Pedro Alberto

2015-01-01

The Kunitz-type protease inhibitor ShPI-1 inhibits human neutrophil elastase (HNE, Ki = 2.35·10-8 M) but does not interact with the porcine pancreatic elastase (PPE); whereas its P1 site variant, ShPI-1/K13L, inhibits both HNE and PPE (Ki = 1.3·10-9 M, and Ki = 1.2·10-8 M, respectively). By employing a combination of molecular modeling tools, e.g., structural alignment, molecular dynamics simulations and Molecular Mechanics Generalized-Born/Poisson-Boltzmann Surface Area free energy calculations, we showed that D226 of HNE plays a critical role in the interaction of this enzyme with ShPI-1 through the formation of a strong salt bridge and hydrogen bonds with K13 at the inhibitor's P1 site, which compensate the unfavorable polar-desolvation penalty of the latter residue. Conversely, T226 of PPE is unable to establish strong interactions with K13, thereby precluding the insertion of K13 side-chain into the S1 subsite of this enzyme. An alternative conformation of K13 site-chain placed at the entrance of the S1 subsite of PPE, similar to that observed in the crystal structure of ShPI-1 in complex with chymotrypsin (PDB: 3T62), is also unfavorable due to the lack of stabilizing pair-wise interactions. In addition, our results suggest that the higher affinity of ShPI-1/K13L for both elastases mainly arises from the lower polar-desolvation penalty of L13 compared to that of K13, and not from stronger pair-wise interactions of the former residue with those of each enzyme. These results provide insights into the PPE and HNE inhibition and may contribute to the design of more potent and/or specific inhibitors toward one of these proteases.

Multiple Pairwise Analysis of Non-homologous Centromere Coupling Reveals Preferential Chromosome Size-Dependent Interactions and a Role for Bouquet Formation in Establishing the Interaction Pattern

PubMed Central

Lefrançois, Philippe; Rockmill, Beth; Xie, Pingxing; Roeder, G. Shirleen; Snyder, Michael

2016-01-01

During meiosis, chromosomes undergo a homology search in order to locate their homolog to form stable pairs and exchange genetic material. Early in prophase, chromosomes associate in mostly non-homologous pairs, tethered only at their centromeres. This phenomenon, conserved through higher eukaryotes, is termed centromere coupling in budding yeast. Both initiation of recombination and the presence of homologs are dispensable for centromere coupling (occurring in spo11 mutants and haploids induced to undergo meiosis) but the presence of the synaptonemal complex (SC) protein Zip1 is required. The nature and mechanism of coupling have yet to be elucidated. Here we present the first pairwise analysis of centromere coupling in an effort to uncover underlying rules that may exist within these non-homologous interactions. We designed a novel chromosome conformation capture (3C)-based assay to detect all possible interactions between non-homologous yeast centromeres during early meiosis. Using this variant of 3C-qPCR, we found a size-dependent interaction pattern, in which chromosomes assort preferentially with chromosomes of similar sizes, in haploid and diploid spo11 cells, but not in a coupling-defective mutant (spo11 zip1 haploid and diploid yeast). This pattern is also observed in wild-type diploids early in meiosis but disappears as meiosis progresses and homologous chromosomes pair. We found no evidence to support the notion that ancestral centromere homology plays a role in pattern establishment in S. cerevisiae post-genome duplication. Moreover, we found a role for the meiotic bouquet in establishing the size dependence of centromere coupling, as abolishing bouquet (using the bouquet-defective spo11 ndj1 mutant) reduces it. Coupling in spo11 ndj1 rather follows telomere clustering preferences. We propose that a chromosome size preference for centromere coupling helps establish efficient homolog recognition. PMID:27768699

Single-cell Hi-C for genome-wide detection of chromatin interactions that occur simultaneously in a single cell.

PubMed

Nagano, Takashi; Lubling, Yaniv; Yaffe, Eitan; Wingett, Steven W; Dean, Wendy; Tanay, Amos; Fraser, Peter

2015-12-01

Hi-C is a powerful method that provides pairwise information on genomic regions in spatial proximity in the nucleus. Hi-C requires millions of cells as input and, as genome organization varies from cell to cell, a limitation of Hi-C is that it only provides a population average of genome conformations. We developed single-cell Hi-C to create snapshots of thousands of chromatin interactions that occur simultaneously in a single cell. To adapt Hi-C to single-cell analysis, we modified the protocol to include in-nucleus ligation. This enables the isolation of single nuclei carrying Hi-C-ligated DNA into separate tubes, followed by reversal of cross-links, capture of biotinylated ligation junctions on streptavidin-coated magnetic beads and PCR amplification of single-cell Hi-C libraries. The entire laboratory protocol can be carried out in 1 week, and although we have demonstrated its use in mouse T helper (TH1) cells, it should be applicable to any cell type or species for which standard Hi-C has been successful. We also developed an analysis pipeline to filter noise and assess the quality of data sets in a few hours. Although the interactome maps produced by single-cell Hi-C are sparse, the data provide useful information to understand cellular variability in nuclear genome organization and chromosome structure. Standard wet and dry laboratory skills in molecular biology and computational analysis are required.

Differential Item Functioning Detection across Two Methods of Defining Group Comparisons: Pairwise and Composite Group Comparisons

ERIC Educational Resources Information Center

Sari, Halil Ibrahim; Huggins, Anne Corinne

2015-01-01

This study compares two methods of defining groups for the detection of differential item functioning (DIF): (a) pairwise comparisons and (b) composite group comparisons. We aim to emphasize and empirically support the notion that the choice of pairwise versus composite group definitions in DIF is a reflection of how one defines fairness in DIF…

Flocculation of copper(II) and tetracycline from water using a novel pH- and temperature-responsive flocculants.

PubMed

Yang, Zhen; Jia, Shuying; Zhuo, Ning; Yang, Weiben; Wang, Yuping

2015-12-01

Insufficient research is available on flocculation of combined pollutants of heavy metals and antibiotics, which widely exist in livestock wastewaters. Aiming at solving difficulties in flocculation of this sort of combined pollution, a novel pH- and temperature-responsive biomass-based flocculant, carboxymethyl chitosan-graft-poly(N-isoproyl acrylamide-co-diallyl dimethyl ammonium chloride) (denoted as CND) with two responsive switches [lower critical solution temperature (LCST) and isoelectric point (IEP)], was designed and synthesized. Its flocculation performance at different temperatures and pHs was evaluated using copper(II) and tetracycline (TC) as model contaminants. CND exhibited high efficiency for coremoval of both contaminants, whereas two commercial flocculants (polyaluminum chloride and polyacrylamide) did not. Especially, flocculation performance of the dual-responsive flocculant under conditions of temperature>LCST and IEP(contaminants)

Female elk contacts are neither frequency nor density dependent

USGS Publications Warehouse

Cross, Paul C.; Creech, Tyler G.; Ebinger, Michael R.; Manlove, Kezia R.; Irvine, Kathryn M.; Henningsen, John C.; Rogerson, Jared D.; Scurlock, Brandon M.; Creely, Scott

2013-01-01

Identifying drivers of contact rates among individuals is critical to understanding disease dynamics and implementing targeted control measures. We studied the interaction patterns of 149 female elk (Cervus canadensis) distributed across five different regions of western Wyoming over three years, defining a contact as an approach within one body length (∼2 m). Using hierarchical models that account for correlations within individuals, pairs, and groups, we found that pairwise contact rates within a group declined by a factor of three as group sizes increased 33-fold. Per capita contact rates, however, increased with group size according to a power function, such that female elk contact rates fell in between the predictions of density- or frequency-dependent disease models. We found similar patterns for the duration of contacts. Our results suggest that larger elk groups are likely to play a disproportionate role in the disease dynamics of directly transmitted infections in elk. Supplemental feeding of elk had a limited impact on pairwise interaction rates and durations, but per capita rates were more than two times higher on feeding grounds. Our statistical approach decomposes the variation in contact rate into individual, dyadic, and environmental effects, and provides insight into factors that may be targeted by disease control programs. In particular, female elk contact patterns were driven more by environmental factors such as group size than by either individual or dyad effects.

Net2Align: An Algorithm For Pairwise Global Alignment of Biological Networks

PubMed Central

Wadhwab, Gulshan; Upadhyayaa, K. C.

2016-01-01

The amount of data on molecular interactions is growing at an enormous pace, whereas the progress of methods for analysing this data is still lacking behind. Particularly, in the area of comparative analysis of biological networks, where one wishes to explore the similarity between two biological networks, this holds a potential problem. In consideration that the functionality primarily runs at the network level, it advocates the need for robust comparison methods. In this paper, we describe Net2Align, an algorithm for pairwise global alignment that can perform node-to-node correspondences as well as edge-to-edge correspondences into consideration. The uniqueness of our algorithm is in the fact that it is also able to detect the type of interaction, which is essential in case of directed graphs. The existing algorithm is only able to identify the common nodes but not the common edges. Another striking feature of the algorithm is that it is able to remove duplicate entries in case of variable datasets being aligned. This is achieved through creation of a local database which helps exclude duplicate links. In a pervasive computational study on gene regulatory network, we establish that our algorithm surpasses its counterparts in its results. Net2Align has been implemented in Java 7 and the source code is available as supplementary files. PMID:28356678

Probing long-range interactions by extracting free energies from genome-wide chromosome conformation capture data.

PubMed

Saberi, Saeed; Farré, Pau; Cuvier, Olivier; Emberly, Eldon

2015-05-23

A variety of DNA binding proteins are involved in regulating and shaping the packing of chromatin. They aid the formation of loops in the DNA that function to isolate different structural domains. A recent experimental technique, Hi-C, provides a method for determining the frequency of such looping between all distant parts of the genome. Given that the binding locations of many chromatin associated proteins have also been measured, it has been possible to make estimates for their influence on the long-range interactions as measured by Hi-C. However, a challenge in this analysis is the predominance of non-specific contacts that mask out the specific interactions of interest. We show that transforming the Hi-C contact frequencies into free energies gives a natural method for separating out the distance dependent non-specific interactions. In particular we apply Principal Component Analysis (PCA) to the transformed free energy matrix to identify the dominant modes of interaction. PCA identifies systematic effects as well as high frequency spatial noise in the Hi-C data which can be filtered out. Thus it can be used as a data driven approach for normalizing Hi-C data. We assess this PCA based normalization approach, along with several other normalization schemes, by fitting the transformed Hi-C data using a pairwise interaction model that takes as input the known locations of bound chromatin factors. The result of fitting is a set of predictions for the coupling energies between the various chromatin factors and their effect on the energetics of looping. We show that the quality of the fit can be used as a means to determine how much PCA filtering should be applied to the Hi-C data. We find that the different normalizations of the Hi-C data vary in the quality of fit to the pairwise interaction model. PCA filtering can improve the fit, and the predicted coupling energies lead to biologically meaningful insights for how various chromatin bound factors influence the stability of DNA loops in chromatin.

Calibration of Smartphone-Based Weather Measurements Using Pairwise Gossip.

PubMed

Zamora, Jane Louie Fresco; Kashihara, Shigeru; Yamaguchi, Suguru

2015-01-01

Accurate and reliable daily global weather reports are necessary for weather forecasting and climate analysis. However, the availability of these reports continues to decline due to the lack of economic support and policies in maintaining ground weather measurement systems from where these reports are obtained. Thus, to mitigate data scarcity, it is required to utilize weather information from existing sensors and built-in smartphone sensors. However, as smartphone usage often varies according to human activity, it is difficult to obtain accurate measurement data. In this paper, we present a heuristic-based pairwise gossip algorithm that will calibrate smartphone-based pressure sensors with respect to fixed weather stations as our referential ground truth. Based on actual measurements, we have verified that smartphone-based readings are unstable when observed during movement. Using our calibration algorithm on actual smartphone-based pressure readings, the updated values were significantly closer to the ground truth values.

Intrasubject multimodal groupwise registration with the conditional template entropy.

PubMed

Polfliet, Mathias; Klein, Stefan; Huizinga, Wyke; Paulides, Margarethus M; Niessen, Wiro J; Vandemeulebroucke, Jef

2018-05-01

Image registration is an important task in medical image analysis. Whereas most methods are designed for the registration of two images (pairwise registration), there is an increasing interest in simultaneously aligning more than two images using groupwise registration. Multimodal registration in a groupwise setting remains difficult, due to the lack of generally applicable similarity metrics. In this work, a novel similarity metric for such groupwise registration problems is proposed. The metric calculates the sum of the conditional entropy between each image in the group and a representative template image constructed iteratively using principal component analysis. The proposed metric is validated in extensive experiments on synthetic and intrasubject clinical image data. These experiments showed equivalent or improved registration accuracy compared to other state-of-the-art (dis)similarity metrics and improved transformation consistency compared to pairwise mutual information. Copyright © 2018 The Authors. Published by Elsevier B.V. All rights reserved.

Calibration of Smartphone-Based Weather Measurements Using Pairwise Gossip

PubMed Central

Yamaguchi, Suguru

2015-01-01

Accurate and reliable daily global weather reports are necessary for weather forecasting and climate analysis. However, the availability of these reports continues to decline due to the lack of economic support and policies in maintaining ground weather measurement systems from where these reports are obtained. Thus, to mitigate data scarcity, it is required to utilize weather information from existing sensors and built-in smartphone sensors. However, as smartphone usage often varies according to human activity, it is difficult to obtain accurate measurement data. In this paper, we present a heuristic-based pairwise gossip algorithm that will calibrate smartphone-based pressure sensors with respect to fixed weather stations as our referential ground truth. Based on actual measurements, we have verified that smartphone-based readings are unstable when observed during movement. Using our calibration algorithm on actual smartphone-based pressure readings, the updated values were significantly closer to the ground truth values. PMID:26421312

Analysis of laser printer and photocopier toners by spectral properties and chemometrics

NASA Astrophysics Data System (ADS)

Verma, Neha; Kumar, Raj; Sharma, Vishal

2018-05-01

The use of printers to generate falsified documents has become a common practice in today's world. The examination and identification of the printed matter in the suspected documents (civil or criminal cases) may provide important information about the authenticity of the document. In the present study, a total number of 100 black toner samples both from laser printers and photocopiers were examined using diffuse reflectance UV-Vis Spectroscopy. The present research is divided into two parts; visual discrimination and discrimination by using multivariate analysis. A comparison between qualitative and quantitative analysis showed that multivariate analysis (Principal component analysis) provides 99.59%pair-wise discriminating power for laser printer toners while 99.84% pair-wise discriminating power for photocopier toners. The overall results obtained confirm the applicability of UV-Vis spectroscopy and chemometrics, in the nondestructive analysis of toner printed documents while enhancing their evidential value for forensic applications.

Independence and totalness of subspaces in phase space methods

NASA Astrophysics Data System (ADS)

Vourdas, A.

2018-04-01

The concepts of independence and totalness of subspaces are introduced in the context of quasi-probability distributions in phase space, for quantum systems with finite-dimensional Hilbert space. It is shown that due to the non-distributivity of the lattice of subspaces, there are various levels of independence, from pairwise independence up to (full) independence. Pairwise totalness, totalness and other intermediate concepts are also introduced, which roughly express that the subspaces overlap strongly among themselves, and they cover the full Hilbert space. A duality between independence and totalness, that involves orthocomplementation (logical NOT operation), is discussed. Another approach to independence is also studied, using Rota's formalism on independent partitions of the Hilbert space. This is used to define informational independence, which is proved to be equivalent to independence. As an application, the pentagram (used in discussions on contextuality) is analysed using these concepts.

Short-term dynamics of causal information transfer in thalamocortical networks during natural inputs and microstimulation for somatosensory neuroprosthesis

PubMed Central

Semework, Mulugeta; DiStasio, Marcello

2014-01-01

Recording the activity of large populations of neurons requires new methods to analyze and use the large volumes of time series data thus created. Fast and clear methods for finding functional connectivity are an important step toward the goal of understanding neural processing. This problem presents itself readily in somatosensory neuroprosthesis (SSNP) research, which uses microstimulation (MiSt) to activate neural tissue to mimic natural stimuli, and has the capacity to potentiate, depotentiate, or even destroy functional connections. As the aim of SSNP engineering is artificially creating neural responses that resemble those observed during natural inputs, a central goal is describing the influence of MiSt on activity structure among groups of neurons, and how this structure may be altered to affect perception or behavior. In this paper, we demonstrate the concept of Granger causality, combined with maximum likelihood methods, applied to neural signals recorded before, during, and after natural and electrical stimulation. We show how these analyses can be used to evaluate the changing interactions in the thalamocortical somatosensory system in response to repeated perturbation. Using LFPs recorded from the ventral posterolateral thalamus (VPL) and somatosensory cortex (S1) in anesthetized rats, we estimated pair-wise functional interactions between functional microdomains. The preliminary results demonstrate input-dependent modulations in the direction and strength of information flow during and after application of MiSt. Cortico-cortical interactions during cortical MiSt and baseline conditions showed the largest causal influence differences, while there was no statistically significant difference between pre- and post-stimulation baseline causal activities. These functional connectivity changes agree with physiologically accepted communication patterns through the network, and their particular parameters have implications for both rehabilitation and brain—machine interface SSNP applications. PMID:25249973

fRMSDPred: Predicting Local RMSD Between Structural Fragments Using Sequence Information

DTIC Science & Technology

2007-04-04

machine learning approaches for estimating the RMSD value of a pair of protein fragments. These estimated fragment-level RMSD values can be used to construct the alignment, assess the quality of an alignment, and identify high-quality alignment segments. We present algorithms to solve this fragment-level RMSD prediction problem using a supervised learning framework based on support vector regression and classification that incorporates protein profiles, predicted secondary structure, effective information encoding schemes, and novel second-order pairwise exponential kernel

Power and sample-size estimation for microbiome studies using pairwise distances and PERMANOVA.

PubMed

Kelly, Brendan J; Gross, Robert; Bittinger, Kyle; Sherrill-Mix, Scott; Lewis, James D; Collman, Ronald G; Bushman, Frederic D; Li, Hongzhe

2015-08-01

The variation in community composition between microbiome samples, termed beta diversity, can be measured by pairwise distance based on either presence-absence or quantitative species abundance data. PERMANOVA, a permutation-based extension of multivariate analysis of variance to a matrix of pairwise distances, partitions within-group and between-group distances to permit assessment of the effect of an exposure or intervention (grouping factor) upon the sampled microbiome. Within-group distance and exposure/intervention effect size must be accurately modeled to estimate statistical power for a microbiome study that will be analyzed with pairwise distances and PERMANOVA. We present a framework for PERMANOVA power estimation tailored to marker-gene microbiome studies that will be analyzed by pairwise distances, which includes: (i) a novel method for distance matrix simulation that permits modeling of within-group pairwise distances according to pre-specified population parameters; (ii) a method to incorporate effects of different sizes within the simulated distance matrix; (iii) a simulation-based method for estimating PERMANOVA power from simulated distance matrices; and (iv) an R statistical software package that implements the above. Matrices of pairwise distances can be efficiently simulated to satisfy the triangle inequality and incorporate group-level effects, which are quantified by the adjusted coefficient of determination, omega-squared (ω2). From simulated distance matrices, available PERMANOVA power or necessary sample size can be estimated for a planned microbiome study. © The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.

Parallel kinetic Monte Carlo simulation framework incorporating accurate models of adsorbate lateral interactions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Nielsen, Jens; D’Avezac, Mayeul; Hetherington, James

2013-12-14

Ab initio kinetic Monte Carlo (KMC) simulations have been successfully applied for over two decades to elucidate the underlying physico-chemical phenomena on the surfaces of heterogeneous catalysts. These simulations necessitate detailed knowledge of the kinetics of elementary reactions constituting the reaction mechanism, and the energetics of the species participating in the chemistry. The information about the energetics is encoded in the formation energies of gas and surface-bound species, and the lateral interactions between adsorbates on the catalytic surface, which can be modeled at different levels of detail. The majority of previous works accounted for only pairwise-additive first nearest-neighbor interactions. Moremore » recently, cluster-expansion Hamiltonians incorporating long-range interactions and many-body terms have been used for detailed estimations of catalytic rate [C. Wu, D. J. Schmidt, C. Wolverton, and W. F. Schneider, J. Catal. 286, 88 (2012)]. In view of the increasing interest in accurate predictions of catalytic performance, there is a need for general-purpose KMC approaches incorporating detailed cluster expansion models for the adlayer energetics. We have addressed this need by building on the previously introduced graph-theoretical KMC framework, and we have developed Zacros, a FORTRAN2003 KMC package for simulating catalytic chemistries. To tackle the high computational cost in the presence of long-range interactions we introduce parallelization with OpenMP. We further benchmark our framework by simulating a KMC analogue of the NO oxidation system established by Schneider and co-workers [J. Catal. 286, 88 (2012)]. We show that taking into account only first nearest-neighbor interactions may lead to large errors in the prediction of the catalytic rate, whereas for accurate estimates thereof, one needs to include long-range terms in the cluster expansion.« less

Determination of ensemble-average pairwise root mean-square deviation from experimental B-factors.

PubMed

Kuzmanic, Antonija; Zagrovic, Bojan

2010-03-03

Root mean-square deviation (RMSD) after roto-translational least-squares fitting is a measure of global structural similarity of macromolecules used commonly. On the other hand, experimental x-ray B-factors are used frequently to study local structural heterogeneity and dynamics in macromolecules by providing direct information about root mean-square fluctuations (RMSF) that can also be calculated from molecular dynamics simulations. We provide a mathematical derivation showing that, given a set of conservative assumptions, a root mean-square ensemble-average of an all-against-all distribution of pairwise RMSD for a single molecular species, (1/2), is directly related to average B-factors () and (1/2). We show this relationship and explore its limits of validity on a heterogeneous ensemble of structures taken from molecular dynamics simulations of villin headpiece generated using distributed-computing techniques and the Folding@Home cluster. Our results provide a basis for quantifying global structural diversity of macromolecules in crystals directly from x-ray experiments, and we show this on a large set of structures taken from the Protein Data Bank. In particular, we show that the ensemble-average pairwise backbone RMSD for a microscopic ensemble underlying a typical protein x-ray structure is approximately 1.1 A, under the assumption that the principal contribution to experimental B-factors is conformational variability. 2010 Biophysical Society. Published by Elsevier Inc. All rights reserved.