The JWS online simulation database

Summary: JWS Online is a web-based platform for construction, simulation and exchange of models in standard formats. We have extended the platform with a database for curated simulation experiments that can be accessed directly via a URL, allowing one-click reproduction of published results. Users can modify the simulation experiments and export them in standard formats. The Simulation database thus lowers the bar on exploring computational models, helps users create valid simulation descriptions and improves the reproducibility of published simulation experiments. Availability and Implementation: The Simulation Database is available on line at https://jjj.bio.vu. nl/models/experiments/

HMG-CoA reductase inhibitors, other lipid-lowering medication, antiplatelet therapy, and the risk of venous thrombosis

Background: Statins [3-hydroxymethyl-3-methylglutaryl coenzyme A (HMG-CoA) reductase inhibitors] and antiplatelet therapy reduce the risk of atherosclerotic disease. Besides a reduction of lipid levels, statins might also have antithrombotic and anti-inflammatory properties, and anti-platelet therapy reduces clot formation. We have studied the risk of venous thrombosis with use of statins, other lipid-lowering medication, and antiplatelet therapy. Materials and methods: Patients with a first episode of deep vein thrombosis in the leg or pulmonary embolism between March 1999 and September 2004 were included in a large population-based case–control study (MEGA study). Control subjects were partners of patients (53%) or recruited via a random-digit-dialing method (47%). Participants reported different all-medication use in a questionnaire. Results: Of 4538 patients, 154 used statins (3.3%), as did 354 of 5914 control subjects (5.7%). The use of statins [odds ratio (OR) 0.45; 95% confidence interval (CI) 0.36–0.56] but not other lipid-lowering medications (OR 1.22; 95% CI 0.62–2.43), was associated with a reduced venous thrombosis risk as compared with individuals who did not use any lipid-lowering medication, after adjustment for age, sex, body mass index, atherosclerotic disease, antiplatelet therapy and use of vitamin K antagonists. Different types and various durations of statin therapy were all associated with a decreased venous thrombosis risk. Antiplatelet therapy also reduced venous thrombosis risk (OR 0.56; 95% CI 0.42–0.74). However, sensitivity analyses suggested that this effect is most likely explained by a so-called ‘healthy user effect’. Simultaneous use of medication most strongly reduced venous thrombosis risk. Conclusion: These results suggest that the use of various types of statins is associated with a reduced risk of venous thrombosis, whereas antiplatelet therapy and other lipid-lowering medications are not.\u

Bringing metabolic networks to life: convenience rate law and thermodynamic constraints

BACKGROUND: Translating a known metabolic network into a dynamic model requires rate laws for all chemical reactions. The mathematical expressions depend on the underlying enzymatic mechanism; they can become quite involved and may contain a large number of parameters. Rate laws and enzyme parameters are still unknown for most enzymes. RESULTS: We introduce a simple and general rate law called "convenience kinetics". It can be derived from a simple random-order enzyme mechanism. Thermodynamic laws can impose dependencies on the kinetic parameters. Hence, to facilitate model fitting and parameter optimisation for large networks, we introduce thermodynamically independent system parameters: their values can be varied independently, without violating thermodynamical constraints. We achieve this by expressing the equilibrium constants either by Gibbs free energies of formation or by a set of independent equilibrium constants. The remaining system parameters are mean turnover rates, generalised Michaelis-Menten constants, and constants for inhibition and activation. All parameters correspond to molecular energies, for instance, binding energies between reactants and enzyme. CONCLUSION: Convenience kinetics can be used to translate a biochemical network – manually or automatically - into a dynamical model with plausible biological properties. It implements enzyme saturation and regulation by activators and inhibitors, covers all possible reaction stoichiometries, and can be specified by a small number of parameters. Its mathematical form makes it especially suitable for parameter estimation and optimisation. Parameter estimates can be easily computed from a least-squares fit to Michaelis-Menten values, turnover rates, equilibrium constants, and other quantities that are routinely measured in enzyme assays and stored in kinetic databases

Microbial catabolic activities are naturally selected by metabolic energy harvest rate

The fundamental trade-off between yield and rate of energy harvest per unit of substrate has been largely discussed as a main characteristic for microbial established cooperation or competition. In this study, this point is addressed by developing a generalized model that simulates competition between existing and not experimentally reported microbial catabolic activities defined only based on well-known biochemical pathways. No specific microbial physiological adaptations are considered, growth yield is calculated coupled to catabolism energetics and a common maximum biomass-specific catabolism rate (expressed as electron transfer rate) is assumed for all microbial groups. Under this approach, successful microbial metabolisms are predicted in line with experimental observations under the hypothesis of maximum energy harvest rate. Two microbial ecosystems, typically found in wastewater treatment plants, are simulated, namely: (i) the anaerobic fermentation of glucose and (ii) the oxidation and reduction of nitrogen under aerobic autotrophic (nitrification) and anoxic heterotrophic and autotrophic (denitrification) conditions. The experimentally observed cross feeding in glucose fermentation, through multiple intermediate fermentation pathways, towards ultimately methane and carbon dioxide is predicted. Analogously, two-stage nitrification (by ammonium and nitrite oxidizers) is predicted as prevailing over nitrification in one stage. Conversely, denitrification is predicted in one stage (by denitrifiers) as well as anammox (anaerobic ammonium oxidation). The model results suggest that these observations are a direct consequence of the different energy yields per electron transferred at the different steps of the pathways. Overall, our results theoretically support the hypothesis that successful microbial catabolic activities are selected by an overall maximum energy harvest rate

Integrating Quantitative Knowledge into a Qualitative Gene Regulatory Network

Despite recent improvements in molecular techniques, biological knowledge remains incomplete. Any theorizing about living systems is therefore necessarily based on the use of heterogeneous and partial information. Much current research has focused successfully on the qualitative behaviors of macromolecular networks. Nonetheless, it is not capable of taking into account available quantitative information such as time-series protein concentration variations. The present work proposes a probabilistic modeling framework that integrates both kinds of information. Average case analysis methods are used in combination with Markov chains to link qualitative information about transcriptional regulations to quantitative information about protein concentrations. The approach is illustrated by modeling the carbon starvation response in Escherichia coli. It accurately predicts the quantitative time-series evolution of several protein concentrations using only knowledge of discrete gene interactions and a small number of quantitative observations on a single protein concentration. From this, the modeling technique also derives a ranking of interactions with respect to their importance during the experiment considered. Such a classification is confirmed by the literature. Therefore, our method is principally novel in that it allows (i) a hybrid model that integrates both qualitative discrete model and quantities to be built, even using a small amount of quantitative information, (ii) new quantitative predictions to be derived, (iii) the robustness and relevance of interactions with respect to phenotypic criteria to be precisely quantified, and (iv) the key features of the model to be extracted that can be used as a guidance to design future experiments

The development of computational biology in South Africa: successes achieved and lessons learnt

Bioinformatics is now a critical skill in many research and commercial environments as biological data are increasing in both size and complexity. South African researchers recognized this need in the mid-1990s and responded by working with the government as well as international bodies to develop initiatives to build bioinformatics capacity in the country. Significant injections of support from these bodies provided a springboard for the establishment of computational biology units at multiple universities throughout the country, which took on teaching, basic research and support roles. Several challenges were encountered, for example with unreliability of funding, lack of skills, and lack of infrastructure. However, the bioinformatics community worked together to overcome these, and South Africa is now arguably the leading country in bioinformatics on the African continent. Here we discuss how the discipline developed in the country, highlighting the challenges, successes, and lessons learnt

Minimum Information About a Simulation Experiment (MIASE)

The original publication is available at www.ploscompbiol.orgReproducibility of experiments is a basic requirement for science. Minimum Information (MI) guidelines have proved a helpful means of enabling reuse of existing work in modern biology. The Minimum Information Required in the Annotation of Models (MIRIAM) guidelines promote the exchange and reuse of biochemical computational models. However, information about a model alone is not sufficient to enable its efficient reuse in a computational setting. Advanced numerical algorithms and complex modeling workflows used in modern computational biology make reproduction of simulations difficult. It is therefore essential to define the core information necessary to perform simulations of those models. The Minimum Information About a Simulation Experiment describes the minimal set of information that must be provided to make the description of a simulation experiment available to others. It includes the list of models to use and their modifications, all the simulation procedures to apply and in which order, the processing of the raw numerical results, and the description of the final output. MIASE allows for the reproduction of any simulation experiment. The provision of this information, along with a set of required models, guarantees that the simulation experiment represents the intention of the original authors. Following MIASE guidelines will thus improve the quality of scientific reporting, and will also allow collaborative, more distributed efforts in computational modeling and simulation of biological processes.The discussions that led to the definition of MIASE benefited from the support of a Japan Partnering Award by the UK Biotechnology and Biological Sciences Research Council. DW was supported by the Marie Curie program and by the German Research Association (DFG Research Training School ‘‘dIEM oSiRiS’’ 1387/1). This publication is based on work (EJC) supported in part by Award No KUK-C1-013-04, made by King Abdullah University of Science and Technology (KAUST). FTB acknowledges support by the NIH (grant 1R01GM081070- 01). JC is supported by the European Commission, DG Information Society, through the Seventh Framework Programme of Information and Communication Technologies, under the VPH NoE project (grant number 223920). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.Publishers versio