Electrostatic potentials and electrostatic interaction energies of rat cytochrome b5 and a simulated anion-exchange adsorbent surface.

Electrostatic potentials were determined for the soluble tryptic core of rat cytochrome b5 (using a structure derived from homology modeling) and a simulated anion-exchange surface through application of the linearized finite-difference Poisson-Boltzmann equation with the simulation code UHBD. Objectives of this work included determination of the contributions of the various charged groups on the protein surface to electrostatic interactions with a simulated anion-exchange surface as a function of orientation, separation distance, and ionic strength, as well as examining the potential existence of a preferred contact orientation. Electrostatic interaction free energies for the complex of the model protein and the simulated surface were computed using the electrostatics section of UHBD employing a 110(3) grid. An initial coarse grid spacing of 2.0 A was required to obtain correct boundary conditions. The boundary conditions of the coarse grid were used in subsequent focusing steps until the electrostatic interaction free energies were relatively independent of grid spacing (at approximately 0.5 A). Explicit error analyses were performed to determine the effects of grid spacing and other model assumptions on the electrostatic interaction free energies. The computational results reveal the presence of a preferred interaction orientation; the interaction energy between these two entities, of opposite net charge, is repulsive over a range of orientations. The electrostatic interaction free energies appear to be the summation of multiple fractional interactions between the protein and the anion-exchange surface. The simulation results are compared with those of ion-exchange adsorption experiments with site-directed mutants of the recombinant protein. Comparisons of the results from the computational and experimental studies should lead to a better understanding of electrostatic interactions of proteins and charged surfaces.     

Identifying miRNA-disease association based on integrating miRNA topological similarity and functional similarity

Background: MicroRNAs (miRNAs) are a significant type of non-coding RNAs, which usually were encoded by endogenous genes with about ~22 nt nucleotides. Accumulating biological experiments have shown that miRNAs have close associations with various human diseases. Although traditional experimental methods achieve great successes in miRNA-disease interaction identification, these methods also have some limitations. Therefore, it is necessary to develop computational method to predict miRNA-disease interactions. Methods: Here, we propose a computational framework (MDVSI) to predict interactions between miRNAs and diseases by integrating miRNA topological similarity and functional similarity. Firstly, the CosRA index is utilized to measure miRNA similarity based on network topological feature. Then, in order to enhance the reliability of miRNA similarity, the functional similarity and CosRA similarity are integrated based on linear weight method. Further, the potential miRNA-disease associations are predicted by using recommendation method. In addition, in order to overcome limitation of recommendation method, for new disease, a new strategy is proposed to predict potential interactions between miRNAs and new disease based on disease functional similarity. Results: To evaluate the performance of different methods, we conduct ten-fold cross validation and de novo test in experiment and compare MDVSI with two the-state-of-art methods. The experimental result shows that MDVSI achieves an AUC of 0.91, which is at least 0.012 higher than other compared methods. Conclusions: In summary, we propose a computational framework (MDSVI) for miRNA-disease interaction prediction. The experiment results demonstrate that it outperforms other the-state-of-the-art methods. Case study shows that it can effectively identify potential miRNA-disease interactions.     

Using binding profiles to predict binding sites of target RNAs

Prediction of RNA-RNA interaction is a key to elucidating possible functions of small non-coding RNAs, and a number of computational methods have been proposed to analyze interacting RNA secondary structures. In this article, we focus on predicting binding sites of target RNAs that are expected to interact with regulatory antisense RNAs in a general form of interaction. For this purpose, we propose bistaRNA, a novel method for predicting multiple binding sites of target RNAs. bistaRNA employs binding profiles that represent scores for hybridized structures, leading to reducing the computational cost for interaction prediction. bistaRNA considers an ensemble of equilibrium interacting structures and seeks to maximize expected accuracy using dynamic programming. Experimental results on real interaction data validate good accuracy and fast computation time of bistaRNA as compared with several competitive methods. Moreover, we aim to find new targets given specific antisense RNAs, which provides interesting insights into antisense RNA regulation. bistaRNA is implemented in C++. The program and Supplementary Material are available at http://rna.naist.jp/program/bistarna/.     

Constraining ribosomal RNA conformational space

Despite the potential for many possible secondary-structure conformations, the native sequence of ribosomal RNA (rRNA) is able to find the correct and universally conserved core fold. This study reports a computational analysis investigating two mechanisms that appear to constrain rRNA secondary-structure conformational space: ribosomal proteins and rRNA sequence composition. The analysis was carried out by using rRNA–ribosomal protein interaction data for the Escherichia coli 16S rRNA and free energy minimization software for secondary-structure prediction. The results indicate that selection pressures on rRNA sequence composition and ribosomal protein–rRNA interaction play a key role in constraining the rRNA secondary structure to a single stable form.     

An Evaluation of Alternative Designs for a Grid Information Service

Computational grids consisting of large and diverse sets of distributed resources have recently been adopted by organizations such as NASA and the NSF. One key component of a computational grid is an information services that provides information about resources, services, and applications to users and their tools. This information is required to use a computational grid and therefore should be available in a timely and reliable manner. In this work, we describe the Globus information service, describe how this service is used, analyze its current performance, and perform trace-driven simulations to evaluate alternative implementations of this grid information service. We find that the majority of the transactions with the information service are changes to the data maintained by the service. We also find that of the three servers we evaluate, one of the commercial products provides the best performance for our workload and that the response time of the information service was not improved during the single experiment we performed with data distributed across two servers.     

Unstructured hexahedral mesh generation of complex vascular trees using a multi-block grid-based approach

The trend towards realistic numerical models of (pathologic) patient-specific vascular structures brings along larger computational domains and more complex geometries, increasing both the computation time and the operator time. Hexahedral grids effectively lower the computational run time and the required computational infrastructure, but at high cost in terms of operator time and minimal cell quality, especially when the computational analyses are targeting complex geometries such as aneurysm necks, severe stenoses and bifurcations. Moreover, such grids generally do not allow local refinements. As an attempt to overcome these limitations, a novel approach to hexahedral meshing is proposed in this paper, which combines the automated generation of multi-block structures with a grid-based method. The robustness of the novel approach is tested on common complex geometries, such as tree-like structures (including trifurcations), stenoses, and aneurysms. Additionally, the performance of the generated grid is assessed using two numerical examples. In the first example, a grid sensitivity analysis is performed for blood flow simulated in an abdominal mouse aorta and compared to tetrahedral grids with a prismatic boundary layer. In the second example, the fluid–structure interaction in a model of an aorta with aortic coarctation is simulated and the effect of local grid refinement is analyzed.     

A fault detection service for wide area distributed computations

The potential for faults in distributed computing systems is a significant complicating factor for application developers. While a variety of techniques exist for detecting and correcting faults, the implementation of these techniques in a particular context can be difficult. Hence, we propose a fault detection service designed to be incorporated, in a modular fashion, into distributed computing systems, tools, or applications. This service uses well-known techniques based on unreliable fault detectors to detect and report component failure, while allowing the user to trade off timeliness of reporting against false positive rates. We describe the architecture of this service, report on experimental results that quantify its cost and accuracy, and describe its use in two applications, monitoring the status of system components of the GUSTO computational grid testbed and as part of the NetSolve network-enabled numerical solver. This revised version was published online in July 2006 with corrections to the Cover Date.     

FastContact: rapid estimate of contact and binding free energies

Interaction free energies are crucial for analyzing binding propensities in proteins. Although the problem of computing binding free energies remains open, approximate estimates have become very useful for filtering potential binding complexes. We report on the implementation of a fast computational estimate of the binding free energy based on a statistically determined desolvation contact potential and Coulomb electrostatics with a distance-dependent dielectric constant, and validated in the Critical Assessment of PRotein Interactions experiment. The application also reports residue contact free energies that rapidly highlight the hotspots of the interaction. AVAILABILITY: The program was written in Fortran. The executable and full documentation is freely available at http://structure.pitt.edu/software/FastContact     

AGScan: a pluggable microarray image quantification software based on the ImageJ library

Many different programs are available to analyze microarray images. Most programs are commercial packages, some are free. In the latter group only few propose automatic grid alignment and batch mode. More often than not a program implements only one quantification algorithm. AGScan is an open source program that works on all major platforms. It is based on the ImageJ library [Rasband (1997-2006)] and offers a plug-in extension system to add new functions to manipulate images, align grid and quantify spots. It is appropriate for daily laboratory use and also as a framework for new algorithms. AVAILABILITY: The program is freely distributed under X11 Licence. The install instructions can be found in the user manual. The software can be downloaded from http://mulcyber.toulouse.inra.fr/projects/agscan/. The questions and plug-ins can be sent to the contact listed below.     

The molecular basis of IGF-II/IGF2R recognition: a combined molecular dynamics simulation,free-energy calculation and computational alanine scanning study

Insulin-like growth factor-II (IGF-II) is a key regulator of cell growth, survival, migration and differentiation, and is thus pivotal in many cancers. An individual with a high IGF-II level is at high risk of developing cancer, whereas IGF2R is implicated as being important in tumor suppression. Thus, uncovering the essence of the IGF-II/IGF2R interaction is very important to understanding the origin of the tumor-suppressing effect of IGF2R. In this study, in order to investigate the interaction of the IGF-II/IGF2R complex and to characterize the binding hot spots of this interaction, a 10 ns molecular dynamics simulation combined with MM-PBSA/MM-GBSA computations and computational alanine scanning was performed on the IGF-II/IGF2R complex. From the results of the free-energy decomposition and the computational alanine scanning calculation, we identified the key residues in the IGF-II/IGF-2R interaction. The results from the calculation were consistent with reported experimental mutagenesis studies. The information on the interaction of IGF-II and IGF2R obtained is vital for understanding how the structure of IGF2R influences the function of IGF-II in growth and development. This study will also lead to new opportunities to develop molecular probes that can assist in diagnostic screening, and even novel approaches to controlling tumor development.     

外来物种入侵大数据采集方法的建立与应用

为准确、快速地获取入侵生物野外调查数据, 我们基于全球卫星导航系统、地理信息系统、移动互联网等现代信息技术提出了外来物种入侵大数据采集方法, 设计并研发了数据表单可自定义的野外调查工具软件——云采集。该系统以Android手机为数据采集终端, 采用C#和Java语言设计开发, 运用卫星导航定位技术实现野外调查发生位置的快速采集, 通过定义9种调查指标的数据类型及指标(列值)默认值、图像拍摄、语音录入、排序等4个辅助属性, 建立调查指标与手机客户端数据录入界面的关联, 实现用户界面可定制的数据录入模式。该系统在国家重点研发项目、福建省科技重大专项及福建省红火蚁(Solenopsis invicta)疫情普查等项目的调查任务中予以应用。实践检验表明: 该系统实现了野外调查数据的离线采集、数据同步、数据查询与输出管理, 将移动智能终端采集取代传统的纸笔记录, 简化了野外调查的流程, 提高了入侵生物野外调查的数据质量, 为外来生物入侵野外调查大数据采集提供了信息化支持。     

Pressure-induced metallization of condensed phase β-HMX under shock loadings via molecular dynamics simulations in conjunction with multi-scale shock technique

Thanks to the advances in grid technologies, we are able to propose here an evolution of our molecular simulator that, when moving to larger systems, instead of reducing the granularity of the dynamical treatment (as is often done in molecular dynamics studies of such systems) exploits the extra power of the grid approach to the end of preserving the detailed nature of theatomistic formulation of the interaction. Key steps of such evolution are: (1) the assemblage of the interaction based on a composition of the ab initio intramolecular data and a portable parameterization of the intermolecular potential linking ab initio evaluation of intramolecular potentials and the partitioning of molecular polarizability; (2) the exploitation of an efficient coordinated porting and running of molecular dynamics codes on the European grid distributed computing infrastructure. As a prototype case study, the N-methylacetamide dimer in vacuo has been considered and the formation of possible conformers is analyzed.     

Empirical Study of an Adaptive Multiscale Model for Simulating Cardiac Conduction

We modify and empirically study an adaptive multiscale model for simulating cardiac action potential propagation along a strand of cardiomyocytes. The model involves microscale partial differential equations posed over cells near the action potential upstroke and macroscale partial differential equations posed over the remainder of the tissue. An important advantage of the modified model of this paper is that, unlike our original model, it does not require perfect alignment between myocytes and the macroscale computational grid. We study the effects of gap-junctional coupling, ephaptic coupling, and macroscale grid spacing on the accuracy of the multiscale model. Our simulations reveal that the multiscale method accurately reproduces both the wavespeed and the waveform, including both upstroke and recovery, of fully microscale models. They also reveal that perfect alignment between myocytes and the macroscale grid is not necessary to reproduce the dynamics of a traveling action potential. Further, our simulations suggest that the macroscale grid spacing used in an adaptive multiscale model need not be much finer than the spatial width of an action potential. These results are demonstrated to hold under high, low, and zero gap-junctional coupling regimes.     

Integrating ARC grid middleware with Taverna workflows

Summary: This work presents two independent approaches for aseamless integration of computational grids with the bioinformaticsworkflow suite Taverna. These are supported by a unique relationaldatabase to link applications with grid resources and presentsthose as workflow elements. A web portal facilitates its collaborativemaintenance. The first approach implements a gateway serviceto handle authentication certificates and all communicationwith the grid. It reads the database to spawn web services forworkflow elements which are in turn used by Taverna. The secondapproach lets Taverna communicate with the grid on its own,by means of a newly developed plug-in. It reads the databaseand executes the needed tasks directly on the grid. While thegateway service is non-intrusive, the plug-in has technicaladvantages, e.g. by allowing data to remain on the grid whilebeing passed between workflow elements. Availability: http://grid.inb.uni-luebeck.de/ Contact: bayer{at}inb.uni-luebeck.de Associate Editor: Alfonso Valencia     

Detecting Protein Complexes in Protein Interaction Networks Modeled as Gene Expression Biclusters

Developing suitable methods for the detection of protein complexes in protein interaction networks continues to be an intriguing area of research. The importance of this objective originates from the fact that protein complexes are key players in most cellular processes. The more complexes we identify, the better we can understand normal as well as abnormal molecular events. Up till now, various computational methods were designed for this purpose. However, despite their notable performance, questions arise regarding potential ways to improve them, in addition to ameliorative guidelines to introduce novel approaches. A close interpretation leads to the assent that the way in which protein interaction networks are initially viewed should be adjusted. These networks are dynamic in reality and it is necessary to consider this fact to enhance the detection of protein complexes. In this paper, we present “DyCluster”, a framework to model the dynamic aspect of protein interaction networks by incorporating gene expression data, through biclustering techniques, prior to applying complex-detection algorithms. The experimental results show that DyCluster leads to higher numbers of correctly-detected complexes with better evaluation scores. The high accuracy achieved by DyCluster in detecting protein complexes is a valid argument in favor of the proposed method. DyCluster is also able to detect biologically meaningful protein groups. The code and datasets used in the study are downloadable from https://github.com/emhanna/DyCluster.     

Optimized detection of transcription factor-binding sites in ChIP-seq experiments

We developed a computational procedure for optimizing the binding site detections in a given ChIP-seq experiment by maximizing their reproducibility under bootstrap sampling. We demonstrate how the procedure can improve the detection accuracies beyond those obtained with the default settings of popular peak calling software, or inform the user whether the peak detection results are compromised, circumventing the need for arbitrary re-iterative peak calling under varying parameter settings. The generic, open-source implementation is easily extendable to accommodate additional features and to promote its widespread application in future ChIP-seq studies. The peakROTS R-package and user guide are freely available at http://www.nic.funet.fi/pub/sci/molbio/peakROTS.     

OBIYagns: a grid-based biochemical simulator with a parameter estimator

SUMMARY: OBIYagns (yet another gene network simulator) is a biochemical system simulator that comprises a multiple-user Web-based graphical interface, an ordinary differential equation solver and a parameter estimators distributed over an open bioinformatics grid (OBIGrid). This grid-based biochemical simulation system can achieve high performance and provide a secure simulation environment for estimating kinetic parameters in an acceptable time period. OBIYagns can be applied to larger system biology-oriented simulation projects. AVAILABILITY: OBIYagns example models, methods and user guide are available at https://access.obigrid.org/yagns/ SUPPLEMENTARY INFORMATION: Please refer to Bioinformatics online.