Hydrogen Dynamics in Superprotonic CsHSO4

We present a detailed study of proton dynamics in the hydrogen-bonded superprotonic conductor CsHSO4 from first-principles molecular dynamics simulations, isolating the subtle interplay between the dynamics of the O--H chemical bonds, the O...H hydrogen bonds, and the SO4 tetrahedra in promoting proton diffusion. We find that the Grotthus mechanism of proton transport is primarily responsible for the dynamics of the chemical bonds, whereas the reorganization of the hydrogen-bond network is dominated by rapid angular hops in concert with small reorientations of the SO4 tetrahedra. Frequent proton jumping across the O--H...O complex is countered by a high rate of jump reversal, which we show is connected to the dynamics of the SO4 tetrahedra, resulting in a diminished CsHSO4/CsDSO4 isotope effect. We also find evidence of multiple timescales for SO4 reorientation events, leading to distinct diffusion mechanisms along the different crystal lattice directions. Finally, we employ graph-theoretic techniques to characterize the topology of the hydrogen-bond network and demonstrate a clear relationship between certain connectivity configurations and the likelihood for diffusive jump events.Comment: 12 pages, 10 figure

An interacting spin flip model for one-dimensional proton conduction

A discrete asymmetric exclusion process (ASEP) is developed to model proton conduction along one-dimensional water wires. Each lattice site represents a water molecule that can be in only one of three states; protonated, left-pointing, and right-pointing. Only a right(left)-pointing water can accept a proton from its left(right). Results of asymptotic mean field analysis and Monte-Carlo simulations for the three-species, open boundary exclusion model are presented and compared. The mean field results for the steady-state proton current suggest a number of regimes analogous to the low and maximal current phases found in the single species ASEP [B. Derrida, Physics Reports, {\bf 301}, 65-83, (1998)]. We find that the mean field results are accurate (compared with lattice Monte-Carlo simulations) only in the certain regimes. Refinements and extensions including more elaborate forces and pore defects are also discussed.Comment: 13pp, 6 fig

Arm-specific dynamics of chromosome evolution in malaria mosquitoes

<p>Abstract</p> <p>Background</p> <p>The malaria mosquito species of subgenus <it>Cellia </it>have rich inversion polymorphisms that correlate with environmental variables. Polymorphic inversions tend to cluster on the chromosomal arms 2R and 2L but not on X, 3R and 3L in <it>Anopheles gambiae </it>and homologous arms in other species. However, it is unknown whether polymorphic inversions on homologous chromosomal arms of distantly related species from subgenus <it>Cellia </it>nonrandomly share similar sets of genes. It is also unclear if the evolutionary breakage of inversion-poor chromosomal arms is under constraints.</p> <p>Results</p> <p>To gain a better understanding of the arm-specific differences in the rates of genome rearrangements, we compared gene orders and established syntenic relationships among <it>Anopheles gambiae, Anopheles funestus</it>, and <it>Anopheles stephensi</it>. We provided evidence that polymorphic inversions on the 2R arms in these three species nonrandomly captured similar sets of genes. This nonrandom distribution of genes was not only a result of preservation of ancestral gene order but also an outcome of extensive reshuffling of gene orders that created new combinations of homologous genes within independently originated polymorphic inversions. The statistical analysis of distribution of conserved gene orders demonstrated that the autosomal arms differ in their tolerance to generating evolutionary breakpoints. The fastest evolving 2R autosomal arm was enriched with gene blocks conserved between only a pair of species. In contrast, all identified syntenic blocks were preserved on the slowly evolving 3R arm of <it>An. gambiae </it>and on the homologous arms of <it>An. funestus </it>and <it>An. stephensi</it>.</p> <p>Conclusions</p> <p>Our results suggest that natural selection favors specific gene combinations within polymorphic inversions when distant species are exposed to similar environmental pressures. This knowledge could be useful for the discovery of genes responsible for an association of inversion polymorphisms with phenotypic variations in multiple species. Our data support the chromosomal arm specificity in rates of gene order disruption during mosquito evolution. We conclude that the distribution of breakpoint regions is evolutionary conserved on slowly evolving arms and tends to be lineage-specific on rapidly evolving arms.</p

CHAMPION: Chalmers Hierarchical Atomic, Molecular, Polymeric & Ionic Analysis Toolkit

We present CHAMPION: a software developed to automatically detect time-dependent bonds between atoms based on their dynamics, classify the local graph topology around them, and analyze the physicochemical properties of these topologies by statistical physics. In stark contrast to methodologies where bonds are detected based on static conditions such as cut-off distances, CHAMPION considers pairs of atoms to be bound only if they move together and act as a bound pair over time. Furthermore, the time-dependent global bond graph is possible to split into dynamically shifting connected components or subgraphs around a certain chemical motif and thereby allow the physicochemical properties of each such topology to be analyzed by statistical physics. Applicable to condensed matter and liquids in general, and electrolytes in particular, this allows both quantitative and qualitative descriptions of local structure, as well as dynamical processes such as speciation and diffusion. We present here a detailed overview of CHAMPION, including its underlying methodology, implementation and capabilities.Comment: 11 pages, 8 figure

Tyrosine Kinase Syk Non-Enzymatic Inhibitors and Potential Anti-Allergic Drug-Like Compounds Discovered by Virtual and In Vitro Screening

In the past decade, the spleen tyrosine kinase (Syk) has shown a high potential for the discovery of new treatments for inflammatory and autoimmune disorders. Pharmacological inhibitors of Syk catalytic site bearing therapeutic potential have been developed, with however limited specificity towards Syk. To address this topic, we opted for the design of drug-like compounds that could impede the interaction of Syk with its cellular partners while maintaining an active kinase protein. To achieve this challenging task, we used the powerful potential of intracellular antibodies for the modulation of cellular functions in vivo, combined to structure-based in silico screening. In our previous studies, we reported the anti-allergic properties of the intracellular antibody G4G11. With the aim of finding functional mimics of G4G11, we developed an Antibody Displacement Assay and we isolated the drug-like compound C-13, with promising in vivo anti-allergic activity. The likely binding cavity of this compound is located at the close vicinity of G4G11 epitope, far away from the catalytic site of Syk. Here we report the virtual screen of a collection of 500,000 molecules against this new cavity, which led to the isolation of 1000 compounds subsequently evaluated for their in vitro inhibitory effects using the Antibody Displacement Assay. Eighty five compounds were selected and evaluated for their ability to inhibit the liberation of allergic mediators from mast cells. Among them, 10 compounds inhibited degranulation with IC50 values ≤10 µM. The most bioactive compounds combine biological activity, significant inhibition of antibody binding and strong affinity for Syk. Moreover, these molecules show a good potential for oral bioavailability and are not kinase catalytic site inhibitors. These bioactive compounds could be used as starting points for the development of new classes of non-enzymatic inhibitors of Syk and for drug discovery endeavour in the field of inflammation related disorders

The Type 2 Diabetes Knowledge Portal: an open access genetic resource dedicated to type 2 diabetes and related traits

Associations between human genetic variation and clinical phenotypes have become a foundation of biomedical research. Most repositories of these data seek to be disease-agnostic and therefore lack disease-focused views. The Type 2 Diabetes Knowledge Portal (T2DKP) is a public resource of genetic datasets and genomic annotations dedicated to type 2 diabetes (T2D) and related traits. Here, we seek to make the T2DKP more accessible to prospective users and more useful to existing users. First, we evaluate the T2DKP's comprehensiveness by comparing its datasets with those of other repositories. Second, we describe how researchers unfamiliar with human genetic data can begin using and correctly interpreting them via the T2DKP. Third, we describe how existing users can extend their current workflows to use the full suite of tools offered by the T2DKP. We finally discuss the lessons offered by the T2DKP toward the goal of democratizing access to complex disease genetic results

Non-equilibrium statistical mechanics: From a paradigmatic model to biological transport

Unlike equilibrium statistical mechanics, with its well-established foundations, a similar widely-accepted framework for non-equilibrium statistical mechanics (NESM) remains elusive. Here, we review some of the many recent activities on NESM, focusing on some of the fundamental issues and general aspects. Using the language of stochastic Markov processes, we emphasize general properties of the evolution of configurational probabilities, as described by master equations. Of particular interest are systems in which the dynamics violate detailed balance, since such systems serve to model a wide variety of phenomena in nature. We next review two distinct approaches for investigating such problems. One approach focuses on models sufficiently simple to allow us to find exact, analytic, non-trivial results. We provide detailed mathematical analyses of a one-dimensional continuous-time lattice gas, the totally asymmetric exclusion process (TASEP). It is regarded as a paradigmatic model for NESM, much like the role the Ising model played for equilibrium statistical mechanics. It is also the starting point for the second approach, which attempts to include more realistic ingredients in order to be more applicable to systems in nature. Restricting ourselves to the area of biophysics and cellular biology, we review a number of models that are relevant for transport phenomena. Successes and limitations of these simple models are also highlighted.Comment: 72 pages, 18 figures, Accepted to: Reports on Progress in Physic

Integrating sequence and array data to create an improved 1000 Genomes Project haplotype reference panel

A major use of the 1000 Genomes Project (1000GP) data is genotype imputation in genome-wide association studies (GWAS). Here we develop a method to estimate haplotypes from low-coverage sequencing data that can take advantage of single-nucleotide polymorphism (SNP) microarray genotypes on the same samples. First the SNP array data are phased to build a backbone (or 'scaffold') of haplotypes across each chromosome. We then phase the sequence data 'onto' this haplotype scaffold. This approach can take advantage of relatedness between sequenced and non-sequenced samples to improve accuracy. We use this method to create a new 1000GP haplotype reference set for use by the human genetic community. Using a set of validation genotypes at SNP and bi-allelic indels we show that these haplotypes have lower genotype discordance and improved imputation performance into downstream GWAS samples, especially at low-frequency variants. © 2014 Macmillan Publishers Limited. All rights reserved

Application of 3D Zernike descriptors to shape-based ligand similarity searching

Background: The identification of promising drug leads from a large database of compounds is an important step in the preliminary stages of drug design. Although shape is known to play a key role in the molecular recognition process, its application to virtual screening poses significant hurdles both in terms of the encoding scheme and speed. Results: In this study, we have examined the efficacy of the alignment independent three-dimensional Zernike descriptor (3DZD) for fast shape based similarity searching. Performance of this approach was compared with several other methods including the statistical moments based ultrafast shape recognition scheme (USR) and SIMCOMP, a graph matching algorithm that compares atom environments. Three benchmark datasets are used to thoroughly test the methods in terms of their ability for molecular classification, retrieval rate, and performance under the situation that simulates actual virtual screening tasks over a large pharmaceutical database. The 3DZD performed better than or comparable to the other methods examined, depending on the datasets and evaluation metrics used. Reasons for the success and the failure of the shape based methods for specific cases are investigated. Based on the results for the three datasets, general conclusions are drawn with regard to their efficiency and applicability