Tackling Exascale Software Challenges in Molecular Dynamics Simulations with GROMACS

GROMACS is a widely used package for biomolecular simulation, and over the last two decades it has evolved from small-scale efficiency to advanced heterogeneous acceleration and multi-level parallelism targeting some of the largest supercomputers in the world. Here, we describe some of the ways we have been able to realize this through the use of parallelization on all levels, combined with a constant focus on absolute performance. Release 4.6 of GROMACS uses SIMD acceleration on a wide range of architectures, GPU offloading acceleration, and both OpenMP and MPI parallelism within and between nodes, respectively. The recent work on acceleration made it necessary to revisit the fundamental algorithms of molecular simulation, including the concept of neighborsearching, and we discuss the present and future challenges we see for exascale simulation - in particular a very fine-grained task parallelism. We also discuss the software management, code peer review and continuous integration testing required for a project of this complexity.Comment: EASC 2014 conference proceedin

The genome and transcriptome of Trichormus sp NMC-1: insights into adaptation to extreme environments on the Qinghai-Tibet Plateau

The Qinghai-Tibet Plateau (QTP) has the highest biodiversity for an extreme environment worldwide, and provides an ideal natural laboratory to study adaptive evolution. In this study, we generated a draft genome sequence of cyanobacteria Trichormus sp. NMC-1 in the QTP and performed whole transcriptome sequencing under low temperature to investigate the genetic mechanism by which T. sp. NMC-1 adapted to the specific environment. Its genome sequence was 5.9 Mb with a G+C content of 39.2% and encompassed a total of 5362 CDS. A phylogenomic tree indicated that this strain belongs to the Trichormus and Anabaena cluster. Genome comparison between T. sp. NMC-1 and six relatives showed that functionally unknown genes occupied a much higher proportion (28.12%) of the T. sp. NMC-1 genome. In addition, functions of specific, significant positively selected, expanded orthogroups, and differentially expressed genes involved in signal transduction, cell wall/membrane biogenesis, secondary metabolite biosynthesis, and energy production and conversion were analyzed to elucidate specific adaptation traits. Further analyses showed that the CheY-like genes, extracellular polysaccharide and mycosporine-like amino acids might play major roles in adaptation to harsh environments. Our findings indicate that sophisticated genetic mechanisms are involved in cyanobacterial adaptation to the extreme environment of the QTP

Nucleotide Binding Switches the Information Flow in Ras GTPases

The Ras superfamily comprises many guanine nucleotide-binding proteins (G proteins) that are essential to intracellular signal transduction. The guanine nucleotide-dependent intrinsic flexibility patterns of five G proteins were investigated in atomic detail through Molecular Dynamics simulations of the GDP- and GTP-bound states (SGDP and SGTP, respectively). For all the considered systems, the intrinsic flexibility of SGDP was higher than that of SGTP, suggesting that Guanine Exchange Factor (GEF) recognition and nucleotide switch require higher amplitude motions than effector recognition or GTP hydrolysis. Functional mode, dynamic domain, and interaction energy correlation analyses highlighted significant differences in the dynamics of small G proteins and Gα proteins, especially in the inactive state. Indeed, SGDP of Gαt, is characterized by a more extensive energy coupling between nucleotide binding site and distal regions involved in GEF recognition compared to small G proteins, which attenuates in the active state. Moreover, mechanically distinct domains implicated in nucleotide switch could be detected in the presence of GDP but not in the presence of GTP. Finally, in small G proteins, functional modes are more detectable in the inactive state than in the active one and involve changes in solvent exposure of two highly conserved amino acids in switches I and II involved in GEF recognition. The average solvent exposure of these amino acids correlates in turn with the rate of GDP release, suggesting for them either direct or indirect roles in the process of nucleotide switch. Collectively, nucleotide binding changes the information flow through the conserved Ras-like domain, where GDP enhances the flexibility of mechanically distinct portions involved in nucleotide switch, and favors long distance allosteric communication (in Gα proteins), compared to GTP

International Veterinary Epilepsy Task Force consensus proposal: Medical treatment of canine epilepsy in Europe

In Europe, the number of antiepileptic drugs (AEDs) licensed for dogs has grown considerably over the last years. Nevertheless, the same questions remain, which include, 1) when to start treatment, 2) which drug is best used initially, 3) which adjunctive AED can be advised if treatment with the initial drug is unsatisfactory, and 4) when treatment changes should be considered. In this consensus proposal, an overview is given on the aim of AED treatment, when to start long-term treatment in canine epilepsy and which veterinary AEDs are currently in use for dogs. The consensus proposal for drug treatment protocols, 1) is based on current published evidence-based literature, 2) considers the current legal framework of the cascade regulation for the prescription of veterinary drugs in Europe, and 3) reflects the authors’ experience. With this paper it is aimed to provide a consensus for the management of canine idiopathic epilepsy. Furthermore, for the management of structural epilepsy AEDs are inevitable in addition to treating the underlying cause, if possible

The Spread of Fecally Transmitted Parasites in Socially-Structured Populations

Mammals are infected by a wide array of gastrointestinal parasites, including parasites that also infect humans and domesticated animals. Many of these parasites are acquired through contact with infectious stages present in soil, feces or vegetation, suggesting that ranging behavior will have a major impact on their spread. We developed an individual-based spatial simulation model to investigate how range use intensity, home range overlap, and defecation rate impact the spread of fecally transmitted parasites in a population composed of social groups (i.e., a socially structured population). We also investigated the effects of epidemiological parameters involving host and parasite mortality rates, transmissibility, disease–related mortality, and group size. The model was spatially explicit and involved the spillover of a gastrointestinal parasite from a reservoir population along the edge of a simulated reserve, which was designed to mimic the introduction pathogens into protected areas. Animals ranged randomly within a “core” area, with biased movement toward the range center when outside the core. We systematically varied model parameters using a Latin hypercube sampling design. Analyses of simulation output revealed a strong positive association between range use intensity and the prevalence of infection. Moreover, the effects of range use intensity were similar in magnitude to effects of group size, mortality rates, and the per-contact probability of transmission. Defecation rate covaried positively with gastrointestinal parasite prevalence. Greater home range overlap had no positive effects on prevalence, with a smaller core resulting in less range overlap yet more intensive use of the home range and higher prevalence. Collectively, our results reveal that parasites with fecal-oral transmission spread effectively in socially structured populations. Future application should focus on parameterizing the model with empirically derived ranging behavior for different species or populations and data on transmission characteristics of different infectious organisms

Econometric Modelling: Basics

AbstractThis chapter addresses basic topics related to choice data analysis. It starts by describing the coding of attribute levels and choosing the functional form of the attributes in the utility function. Next, it focuses on econometric models with special attention devoted to the random parameter mixed logit model. In this context, the chapter compares different coefficient distributions to be used, addresses specifics of the cost attribute coefficient and it pays attention to potential correlations between random coefficients. Finally, topics related to the estimation procedure such as assuring its convergence or random draws are discussed

Role of structural dynamics at the receptor G protein interface for signal transduction

GPCRs catalyze GDP/GTP exchange in the α-subunit of heterotrimeric G proteins (Gαßγ) through displacement of the Gα C-terminal α5 helix, which directly connects the interface of the active receptor (R*) to the nucleotide binding pocket of G. Hydrogen-deuterium exchange mass spectrometry and kinetic analysis of R* catalysed G protein activation have suggested that displacement of α5 starts from an intermediate GDP bound complex (R*•GGDP). To elucidate the structural basis of receptor-catalysed displacement of α5, we modelled the structure of R*•GGDP. A flexible docking protocol yielded an intermediate R*•GGDP complex, with a similar overall arrangement as in the X-ray structure of the nucleotide free complex (R*•Gempty), however with the α5 C-terminus (GαCT) forming different polar contacts with R*. Starting molecular dynamics simulations of GαCT bound to R* in the intermediate position, we observe a screw-like motion, which restores the specific interactions of α5 with R* in R*•Gempty. The observed rotation of α5 by 60° is in line with experimental data. Reformation of hydrogen bonds, water expulsion and formation of hydrophobic interactions are driving forces of the α5 displacement. We conclude that the identified interactions between R* and G protein define a structural framework in which the α5 displacement promotes direct transmission of the signal from R* to the GDP binding pocket

MSOAR 2.0: Incorporating tandem duplications into ortholog assignment based on genome rearrangement

<p>Abstract</p> <p>Background</p> <p>Ortholog assignment is a critical and fundamental problem in comparative genomics, since orthologs are considered to be functional counterparts in different species and can be used to infer molecular functions of one species from those of other species. MSOAR is a recently developed high-throughput system for assigning one-to-one orthologs between closely related species on a genome scale. It attempts to reconstruct the evolutionary history of input genomes in terms of genome rearrangement and gene duplication events. It assumes that a gene duplication event inserts a duplicated gene into the genome of interest at a random location (<it>i.e.</it>, the random duplication model). However, in practice, biologists believe that genes are often duplicated by tandem duplications, where a duplicated gene is located next to the original copy (<it>i.e.</it>, the tandem duplication model).</p> <p>Results</p> <p>In this paper, we develop MSOAR 2.0, an improved system for one-to-one ortholog assignment. For a pair of input genomes, the system first focuses on the tandemly duplicated genes of each genome and tries to identify among them those that were duplicated after the speciation (<it>i.e.</it>, the so-called inparalogs), using a simple phylogenetic tree reconciliation method. For each such set of tandemly duplicated inparalogs, all but one gene will be deleted from the concerned genome (because they cannot possibly appear in any one-to-one ortholog pairs), and MSOAR is invoked. Using both simulated and real data experiments, we show that MSOAR 2.0 is able to achieve a better sensitivity and specificity than MSOAR. In comparison with the well-known genome-scale ortholog assignment tool InParanoid, Ensembl ortholog database, and the orthology information extracted from the well-known whole-genome multiple alignment program MultiZ, MSOAR 2.0 shows the highest sensitivity. Although the specificity of MSOAR 2.0 is slightly worse than that of InParanoid in the real data experiments, it is actually better than that of InParanoid in the simulation tests.</p> <p>Conclusions</p> <p>Our preliminary experimental results demonstrate that MSOAR 2.0 is a highly accurate tool for one-to-one ortholog assignment between closely related genomes. The software is available to the public for free and included as online supplementary material.</p