Motivation: Simulation and modeling is becoming a standard approach to understand complex biochemical processes. Therefore, there is a big need for software tools that allow access to diverse simulation and modeling methods as well as support for the usage of these methods. Results: Here, we present COPASI, a platform-independent and userfriendly biochemical simulator that offers several unique features. We discuss numerical issues with these features; in particular, the criteria to switch between stochastic and deterministic simulation methods, hybrid deterministic-stochastic methods, and the importance of random number generator numerical resolution in stochastic simulation. Availability: The complete software is available in binary (executable) for MS Windows, OS X, Linux (Intel) and Sun Solaris (SPARC), as well as the full source code under an open source license from http://www.

735hardly imagine today's electronics industry, with its powerful, visually oriented design and automation tools, without having first established standard notations for circuit diagrams. Such was not the case in biology 2 . Despite the visual nature of much of the information exchange, the field was permeated with ad hoc graphical notations having little in common between different researchers, publications, textbooks and software tools. No standard visual language existed for describing biochemical interaction networks, inter-and intracellular signaling gene regulation-concepts at the core of much of today's research in molecular, systems and synthetic biology. The closest to a standard is the notation long used in many metabolic and signaling pathway maps, but in reality, even that lacks uniformity between sources and suffers from undesirable ambiguities (Fig. 1). Moreover, the existing tentative representations, however well crafted, were ambiguous, and only suitable for specific needs, such as representing metabolic networks or signaling pathways or gene regulation.The molecular biology era, and more recently the rise of genomics and other high-throughput technologies, have brought a staggering increase in data to be interpreted. It also favored the routine use of software to help formulate hypotheses, design experiments and interpret results. As a group of biochemists, modelers and computer scientists working in systems biology, we believe establishing standard graphical notations is an important step toward more efficient and accurate transmission of biological knowledge among our different communities. Toward this goal, we initiated the SBGN project in 2005, with the aim of developing and standardizing a systematic and unambiguous graphical notation for applications in molecular and systems biology. Historical antecedentsGraphical representation of biochemical and cellular processes has been used in biochemical textbooks as far back as sixty years ago 3 , reaching an apex in the wall charts hand drawn by Nicholson 4 and Michal 5 . Those graphs describe the processes that transform a set of inputs into a set of outputs, in effect being process, or state transition, diagrams. This style was emulated in the first database systems that depicted metabolic networks, including EMP 6 , EcoCyc 7 and KEGG 8 . More notations have been 'defined' by virtue of their implementation in specialized software tools such as pathway and network designers (e.g., NetBuilder 9 , Patika 10 , JDesigner 11 , CellDesigner 12 ). Those "Un bon croquis vaut mieux qu'un long discours" ("A good sketch is better than a long speech"), said Napoleon Bonaparte. This claim is nowhere as true as for technical illustrations. Diagrams naturally engage innate cognitive faculties 1 that humans have possessed since before the time of our cave-drawing ancestors. Little wonder that we find ourselves turning to them in every field of endeavor. Just as with written human languages, communication involving diagrams requires that authors and readers agr...

Understanding the mechanism and specificity of substrate binding in the cytochrome P450 (P450) superfamily is an important step toward explaining its key role in drug metabolism, toxicity, xenobiotic degradation, and several biosynthetic pathways. Here we investigate the ligand exit pathways and mechanisms of P450cam (CYP101), P450BM-3 (CYP102), and P450eryF (CYP107A1) by using random expulsion molecular dynamics and classical molecular dynamics simulations. Although several different pathways are found for each protein, one pathway is common to all three. The mechanism of ligand exit along this pathway is, however, quite different in the three different proteins. For P450cam, small backbone conformational changes, in combination with aromatic side chain rotation, allow for the passage of the rather rigid, compact, and hydrophobic substrate, camphor. In P450BM-3, larger transient backbone changes are observed on ligand exit. R47, situated at the entrance to the channel, appears important in guiding negatively charged fatty acid substrates in and out of the active site. In P450eryF, an isolated buried arginine, R185, stabilized by four hydrogen bonds to backbone carbonyl oxygen atoms, is located in the exit channel and is identified as having a particularly unusual functionality, dynamically gating channel opening. The results for these three P450s suggest that the channel opening mechanisms are adjusted to the physico-chemical properties of the substrate and can kinetically modulate protein-substrate specificity.

Computational modeling and simulation of biochemical networks is at the core of systems biology and this includes many types of analyses that can aid understanding of how these systems work. COPASI is a generic software package for modeling and simulation of biochemical networks which provides many of these analyses in convenient ways that do not require the user to program or to have deep knowledge of the numerical algorithms. Here we provide a description of how these modeling techniques can be applied to biochemical models using COPASI. The focus is both on practical aspects of software usage as well as on the utility of these analyses in aiding biological understanding. Practical examples are described for steady-state and time-course simulations, stoichiometric analyses, parameter scanning, sensitivity analysis (including metabolic control analysis), global optimization, parameter estimation, and stochastic simulation. The examples used are all published models that are available in the BioModels database in SBML format.