Stochastic simulation GUI for biochemical networks

MOTIVATION: This article describes the development of a useful graphical user interface for stochastic simulation of biochemical networks which allows model builders to run stochastic simulations of their models and perform statistical analysis on the results. These include the construction of correlations, power-spectral densities and transfer functions between selected inputs and outputs. AVAILABILITY: The software is licensed under the BSD open source license and is available at http://sourceforge.net/projects/jdesigner. In addition, a more detailed account of the algorithms employed in the tool can be found at the Wiki at http://www.sys-bio.org/sbwWiki. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Condor-COPASI: high-throughput computing for biochemical networks

ABSTRACT: BACKGROUND: Mathematical modelling has become a standard technique to improve our understanding of complex biological systems. As models become larger and more complex, simulations and analyses require increasing amounts of computational power. Clusters of computers in a high-throughput computing environment can help to provide the resources required for computationally expensive model analysis. However, exploiting such a system can be difficult for users without the necessary expertise. RESULTS: We present Condor-COPASI, a server-based software tool that integrates COPASI, a biological pathway simulation tool, with Condor, a high-throughput computing environment. Condor-COPASI provides a web-based interface, which makes it extremely easy for a user to run a number of model simulation and analysis tasks in parallel. Tasks are transparently split into smaller parts, and and submitted for execution on a Condor pool. Result output is presented to the user in a number of formats, including tables and interactive graphical displays. CONCLUSIONS: Condor-COPASI can effectively use a Condor high-throughput computing environment to provide significant gains in performance for a number of model simulation and analysis tasks. Condor-COPASI is free, open source software, released under the Artistic License 2.0, and is suitable for use by any institution with access to a Condor pool. Source code is freely available for download at http://code.google.com/p/condor-copasi/, along with full instructions on deployment and usage.     

Native structure-based modeling and simulation of biomolecular systems per mouse click

Background

Molecular dynamics (MD) simulations provide valuable insight into biomolecular systems at the atomic level. Notwithstanding the ever-increasing power of high performance computers current MD simulations face several challenges: the fastest atomic movements require time steps of a few femtoseconds which are small compared to biomolecular relevant timescales of milliseconds or even seconds for large conformational motions. At the same time, scalability to a large number of cores is limited mostly due to long-range interactions. An appealing alternative to atomic-level simulations is coarse-graining the resolution of the system or reducing the complexity of the Hamiltonian to improve sampling while decreasing computational costs. Native structure-based models, also called Gō-type models, are based on energy landscape theory and the principle of minimal frustration. They have been tremendously successful in explaining fundamental questions of, e.g., protein folding, RNA folding or protein function. At the same time, they are computationally sufficiently inexpensive to run complex simulations on smaller computing systems or even commodity hardware. Still, their setup and evaluation is quite complex even though sophisticated software packages support their realization.

Results

Here, we establish an efficient infrastructure for native structure-based models to support the community and enable high-throughput simulations on remote computing resources via GridBeans and UNICORE middleware. This infrastructure organizes the setup of such simulations resulting in increased comparability of simulation results. At the same time, complete workflows for advanced simulation protocols can be established and managed on remote resources by a graphical interface which increases reusability of protocols and additionally lowers the entry barrier into such simulations for, e.g., experimental scientists who want to compare their results against simulations. We demonstrate the power of this approach by illustrating it for protein folding simulations for a range of proteins.

Conclusions

We present software enhancing the entire workflow for native structure-based simulations including exception-handling and evaluations. Extending the capability and improving the accessibility of existing simulation packages the software goes beyond the state of the art in the domain of biomolecular simulations. Thus we expect that it will stimulate more individuals from the community to employ more confidently modeling in their research.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-292) contains supplementary material, which is available to authorized users.     

Improvements in GROMACS plugin for PyMOL including implicit solvent simulations and displaying results of PCA analysis

In order to get the dynamic molecule model from the static one, the molecular dynamics (MD) simulation needs to be performed. Some software sets such as GROMACS are used for that purpose. Unfortunately they lack GUI. The Dynamics PyMOL plugin allows researcher to perform MD simulations directly from the PyMOL software by GUI-based interface of GROMACS tools. This paper describes many improvements introduced into the Dynamics PyMOL plugin 2.0 including: an integration with ProDy library, possibility to use the implicit solvents, an ability to interpret the MD simulations, and implementation of some more GROMACS functionality.     

Java Web Simulation (JWS); a web based database of kinetic models

Software to make a database of kinetic models accessible via the internet has been developed and a core database has been set up at http://jjj.biochem.sun.ac.za/. This repository of models, available to everyone with internet access, opens a whole new way in which we can make our models public. Via the database, a user can change enzyme parameters and run time simulations or steady state analyses. The interface is user friendly and no additional software is necessary. The database currently contains 10 models, but since the generation of the program code to include new models has largely been automated the addition of new models is straightforward and people are invited to submit their models to be included in the database.     

EpiFire: An open source C++ library and application for contact network epidemiology

ABSTRACT: BACKGROUND: Contact network models have become increasingly common in epidemiology, but we lack a flexible programming framework for the generation and analysis of epidemiological contact networks and for the simulation of disease transmission through such networks. RESULTS: Here we present EpiFire, an applications programming interface and graphical user interface implemented in C++, which includes a fast and efficient library for generating, analyzing and manipulating networks. Network-based percolation and chain-binomial simulations of susceptible-infected-recovered disease transmission, as well as traditional non-network mass-action simulations, can be performed using EpiFire. CONCLUSIONS: EpiFire provides an open-source programming interface for the rapid development of network models with a focus in contact network epidemiology. EpiFire also provides a point-and-click interface for generating networks, conducting epidemic simulations, and creating figures. This interface is particularly useful as a pedagogical tool.     

AWTY (are we there yet?): a system for graphical exploration of MCMC convergence in Bayesian phylogenetics

A key element to a successful Markov chain Monte Carlo (MCMC) inference is the programming and run performance of the Markov chain. However, the explicit use of quality assessments of the MCMC simulations-convergence diagnostics-in phylogenetics is still uncommon. Here, we present a simple tool that uses the output from MCMC simulations and visualizes a number of properties of primary interest in a Bayesian phylogenetic analysis, such as convergence rates of posterior split probabilities and branch lengths. Graphical exploration of the output from phylogenetic MCMC simulations gives intuitive and often crucial information on the success and reliability of the analysis. The tool presented here complements convergence diagnostics already available in other software packages primarily designed for other applications of MCMC. Importantly, the common practice of using trace-plots of a single parameter or summary statistic, such as the likelihood score of sampled trees, can be misleading for assessing the success of a phylogenetic MCMC simulation. AVAILABILITY: The program is available as source under the GNU General Public License and as a web application at http://ceb.scs.fsu.edu/awty.     

IQuant: An automated pipeline for quantitative proteomics based upon isobaric tags

Quantitative proteomics technology based on isobaric tags is playing an important role in proteomic investigation. In this paper, we present an automated software, named IQuant, which integrates a postprocessing tool of protein identification and advanced statistical algorithms to process the MS/MS signals generated from the peptides labeled by isobaric tags and aims at proteomics quantification. The software of IQuant, which is freely downloaded at http://sourceforge.net/projects/iquant/ , can run from a graphical user interface and a command‐line interface, and can work on both Windows and Linux systems.     

Viola: A New Visual Programming Language Designed for the Rapid Development of Interacting Agent Systems

The construction of complex simulation models and the application of new computer hardware to ecological problems has resulted in the need for many ecologists to rely on computer programmers to develop their modelling software. However, this can lead to a lack of flexibility and understanding in model implementation and in resource problems for researchers. This paper presents a new programming language, Viola, based on a simple organisational concept which can be used by most researchers to develop complex simulations much more easily than could be achieved with standard programming languages such as C++. The language is object oriented and implemented through a visual interface. It is specifically designed to cope with complicated individual based behavioural simulations and comes with embedded concurrency handling abilities.     

QSim, a program for NMR simulations

We present QSim, a program for simulation of NMR experiments. Pulse sequences are implemented and analyzed in QSim using a mouse driven interface. QSim can handle almost any modern NMR experiment, using multiple channels, shaped pulses, mixing, decoupling, phase-cycling and pulsed field gradients. Any number of spins with any spin quantum number can, in theory, be used in simulations. Relaxation is accounted for during all steps of pulse sequences and relaxation interference effects are supported. Chemical kinetics between any numbers of states can be simulated. Both classical and quantum mechanical calculations can be performed. The result of a simulation can be presented either as magnetization as a function of time or as a processed spectrum.     

Predictive modelling of cervical disc implant wear

This study presents a chain of simulations aimed at estimating the wear in a cervical disc implant and providing insight into the in vivo biomechanical performance of the implant. The simulation chain can start with determining a representative maximum range of motion (ROM) of a person's head. The ROM is used as motion input to a kinematic simulation of the cervical spine containing a disc implant. The cervical spine geometry is obtained from computed tomography (CT) scans and converted to STL format using reverse engineering software. The time histories of the loads imposed by the adjacent vertebrae on the implant, as well as the vertebral relative rotations can be extracted from the kinematic simulation. Alternatively, force and motion profiles prescribed by wear test protocols (e.g. ISO 18192-1 and ASTM F2423-05) can be used. The force and motion profiles are applied as boundary conditions to a non-linear finite element model (FEM) of the implant to determine the time-varying contact stress and slip velocity distributions at the interface between the two halves of the implant. The stresses and slip velocities are used in a linear wear model to estimate the wear rate distribution at the FEM's nodal points where contact occurs. Reverse engineering software is used to triangulate the contact surface so that the total wear volume can be calculated. The simulation chain's predicted wear rate shows good agreement with in vitro results in the literature. The simulation chain is thereby demonstrated to be suitable for comparative pre-experimental studies of spinal implant designs.     

A new JAVA interface implementation of THESIAS: testing haplotype effects in association studies

THESIAS (Testing Haplotype EffectS In Association Studies) is a popular software for carrying haplotype association analysis in unrelated individuals. In addition to the command line interface, a graphical JAVA interface is now proposed allowing one to run THESIAS in a user-friendly manner. Besides, new functionalities have been added to THESIAS including the possibility to analyze polychotomous phenotype and X-linked polymorphisms. AVAILABILITY: The software package including documentation and example data files is freely available at http://genecanvas.ecgene.net. The source codes are also available upon request.     

QTLNetwork: mapping and visualizing genetic architecture of complex traits in experimental populations

SUMMARY: QTLNetwork is a software package for mapping and visualizing the genetic architecture underlying complex traits for experimental populations derived from a cross between two inbred lines. It can simultaneously map quantitative trait loci (QTL) with individual effects, epistasis and QTL-environment interaction. Currently, it is able to handle data from F(2), backcross, recombinant inbred lines and double-haploid populations, as well as populations from specific mating designs (immortalized F(2) and BC(n)F(n) populations). The Windows version of QTLNetwork was developed with a graphical user interface. Alternatively, the command-line versions have the facility to be run in other prevalent operating systems, such as Linux, Unix and MacOS. AVAILABILITY: http://ibi.zju.edu.cn/software/qtlnetwork.     

Stochastic simulation of biological reactions, and its applications for studying actin polymerization

Molecular events in biological cells occur in local subregions, where the molecules tend to be small in number. The cytoskeleton, which is important for both the structural changes of cells and their functions, is also a countable entity because of its long fibrous shape. To simulate the local environment using a computer, stochastic simulations should be run. We herein report a new method of stochastic simulation based on random walk and reaction by the collision of all molecules. The microscopic reaction rate P(r) is calculated from the macroscopic rate constant k. The formula involves only local parameters embedded for each molecule. The results of the stochastic simulations of simple second-order, polymerization, Michaelis-Menten-type and other reactions agreed quite well with those of deterministic simulations when the number of molecules was sufficiently large. An analysis of the theory indicated a relationship between variance and the number of molecules in the system, and results of multiple stochastic simulation runs confirmed this relationship. We simulated Ca2(+) dynamics in a cell by inward flow from a point on the cell surface and the polymerization of G-actin forming F-actin. Our results showed that this theory and method can be used to simulate spatially inhomogeneous events.     

Experimental verification of lipid bilayer structure through multi-scale modeling

Great progress has been made in applying coarse-grain molecular dynamics (CGMD) simulations to the investigation of membrane biophysics. In order to validate the accuracy of CGMD simulations of membranes, atomistic scale detail is necessary for direct comparison to structural experiments. Here, we present our strategy for verifying CGMD lipid bilayer simulations. Through reverse coarse graining and subsequent calculation of the bilayer electron density profile, we are able to compare the simulations to our experimental low angle X-ray scattering (LAXS) data. In order to determine the best match to the experimental data, atomistic simulations are run at a range of areas (in the NP_NAT ensemble), starting from distinct configurations extracted from the CGMD simulation (run in the NPT ensemble). We demonstrate the effectiveness of this procedure with two small, single-component bilayers, and suggest that the greater utility of our algorithm will be for CGMD simulations of more complex structures.