Linked Lists and the Method of Lights in Molecular Dynamics Simulation - Search for the Best Method of Forces Evaluation in Sequential MD Codes

Abstract

In this paper two methods of forces evaluation used in the MD codes are presented and compared against the classical linked lists algorithm [3,4] and its modified version [5]. The first algorithm, so called the method of lights is a sequential version of the CYBER 205 vector oriented code [6]. A new algorithm of forces evaluation is also proposed, which incorporates advantages of the method of lights and the linked lists technique.     

Integration of multiple types of genetic markers for neuroblastoma may contribute to improved prediction of the overall survival

Background

Modern experimental techniques deliver data sets containing profiles of tens of thousands of potential molecular and genetic markers that can be used to improve medical diagnostics. Previous studies performed with three different experimental methods for the same set of neuroblastoma patients create opportunity to examine whether augmenting gene expression profiles with information on copy number variation can lead to improved predictions of patients survival. We propose methodology based on comprehensive cross-validation protocol, that includes feature selection within cross-validation loop and classification using machine learning. We also test dependence of results on the feature selection process using four different feature selection methods.

Results

The models utilising features selected based on information entropy are slightly, but significantly, better than those using features obtained with t-test. The synergy between data on genetic variation and gene expression is possible, but not confirmed. A slight, but statistically significant, increase of the predictive power of machine learning models has been observed for models built on combined data sets. It was found while using both out of bag estimate and in cross-validation performed on a single set of variables. However, the improvement was smaller and non-significant when models were built within full cross-validation procedure that included feature selection within cross-validation loop. Good correlation between performance of the models in the internal and external cross-validation was observed, confirming the robustness of the proposed protocol and results.

Conclusions

We have developed a protocol for building predictive machine learning models. The protocol can provide robust estimates of the model performance on unseen data. It is particularly well-suited for small data sets. We have applied this protocol to develop prognostic models for neuroblastoma, using data on copy number variation and gene expression. We have shown that combining these two sources of information may increase the quality of the models. Nevertheless, the increase is small and larger samples are required to reduce noise and bias arising due to overfitting.

Reviewers

This article was reviewed by Lan Hu, Tim Beissbarth and Dimitar Vassilev.

     

Integer Interparticle Distances in Molecular Dynamics Simulation

The paper describes a sequential MD algorithm in which distances between particles are evaluated using fixed point arithmetics. Errors introduced by the method are estimated. Some simulation timings as well as fluctuations of the total energy are compared to the results obtained using floating point arithmetics.     

Parallel splitting solvers for the isogeometric analysis of the Cahn-Hilliard equation

Abstract

Modeling tumor growth in biological systems is a challenging problem with important consequences for diagnosis and treatment of various forms of cancer. This growth process requires large simulation complexity due to evolving biological and chemical processes in living tissue and interactions of cellular and vascular constituents in living organisms. Herein, we describe with a phase-field model, namely the Cahn-Hilliard equation the intricate interactions between the tumors and their host tissue. The spatial discretization uses highly-continuous isogeometric elements. For fast simulation of the time-dependent Cahn-Hilliard equation, we employ an alternating directions implicit methodology. Thus, we reduce the original problems to Kronecker products of 1?D matrices that can be factorized in a linear computational cost. The implementation enables parallel multi-core simulations and shows good scalability on shared-memory multi-core machines. Combined with the high accuracy of isogeometric elements, our method shows high efficiency in solving the Cahn-Hilliard equation on tensor-product meshes.     

Structural insight into an ankyrin-sensitive lipid-binding site of erythroid beta-spectrin

It was recently shown that the region within beta-spectrin responsible for interactions with ankyrin includes a lipid-binding site which displayed sensitivity to inhibition by ankyrin. We studied its structure by constructing a series of single and double spin-labeled beta-spectrin-derived peptides and analyzing their spin-spin distances via electron paramagnetic resonance spectroscopy and the Fourier deconvolution method. The results indicate that the whole ankyrin-sensitive lipid-binding site of beta-spectrin exhibits a helical conformation revealing a distinct 3(10)-helix contribution at its N-terminus. The start of the helix was located five residues upstream along the sequence compared to the theoretical predictions. A model based on the obtained data provides direct evidence that the examined lipid-binding site is a highly amphipathic helix, which is correlated with the specific conformation of its N-terminal fragment.     

A Metastable Contact and Structural Disorder in the Estrogen Receptor Transactivation Domain

1. Download high-res image (178KB)
2. Download full-size image

     

Fecal filtrate transplantation protects against necrotizing enterocolitis

Necrotizing enterocolitis (NEC) is a life-threatening gastrointestinal disorder afflicting preterm infants, which is currently unpreventable. Fecal microbiota transplantation (FMT) is a promising preventive therapy, but the transfer of pathogenic microbes or toxic compounds raise concern. Removal of bacteria from donor feces by micropore filtering may reduce this risk of bacterial infection, while residual bacteriophages could maintain the NEC-preventive effects. We aimed to assess preclinical efficacy and safety of fecal filtrate transplantation (FFT). Using fecal material from healthy suckling piglets, we compared rectal FMT administration (FMT, n = 16) with cognate FFT by either rectal (FFTr, n = 14) or oro-gastric administration (FFTo, n = 13) and saline (CON, n = 16) in preterm, cesarean-delivered piglets as models for preterm infants. We assessed gut pathology and analyzed mucosal and luminal bacterial and viral composition using 16S rRNA gene amplicon and meta-virome sequencing. Finally, we used isolated ileal mucosa, coupled with RNA-Seq, to gauge the host response to the different treatments. Oro-gastric FFT completely prevented NEC, which was confirmed by microscopy, whereas FMT did not perform better than control. Oro-gastric FFT increased viral diversity and reduced Proteobacteria relative abundance in the ileal mucosa relative to control. An induction of mucosal immunity was observed in response to FMT but not FFT. As preterm infants are extremely vulnerable to infections, rational NEC-preventive strategies need incontestable safety profiles. We show in a clinically relevant animal model that FFT, as opposed to FMT, efficiently prevents NEC without any recognizable side effects.Subject terms: Bacteriophages, Microbial ecology, Inflammatory bowel disease     

Mapping of branchpoint nucleotides in mutant pre-mRNAs expressed in plant cells

Pre-mRNA encoding rubisco activase in the Arabidopsis thaliana mutant rca contains a GU to AU change at the 5' splice site of intron 3 and this mutation results in accumulation of splicing intermediates bearing an incompletely processed intron. It has been demonstrated that one of the intermediates contains intron 3 in the form of a lariat and the branchpoint nucleotide has been mapped to the A residue at position −32 forming part of the sequence UUG A U. Analysis of a similar GU to AU 5' splice site mutation, present in a synthetic pre-mRNA context expressed in transfected protoplasts of Nicotiana plumbaginifolia , also suggests formation of lariats with branching occurring at A₋₃₁. A small fraction (approximately 10%) of this mutant pre-mRNA also underwent the second step of splicing. In addition to the consensus AG, an AU dinucleotide was used as splicing acceptor.