Purpose
The Virtual Cell, designed to be accessible to the experimental biologist, is a fully modular computational framework that provides a general approach to modeling the spatially organized and interdependent chemical events that underlie dynamic cellular processes. A Java biological interface allows experimentalists to input specific geometry, in the form of experimental images, specify compartment topology, define and assign species, chemical reactions and transport kinetics, and computational mesh.
The user interface for the Virtual Cell allows full interaction from the WWW. This interface facilitates the creation and analysis of abstract models, generation of specific simulations, control and monitoring of simulations, and data analysis. This interface aims to present the physiology of chemical and structural phenomena and mechanisms in a natural and consistent manner.
| Current Implementation The current Virtual Cell application utilizes a newly designed system-level architecture based on unified modeling language (UML). The system is decomposed into an application framework and system services. The application framework is the Modeling Framework. The system services are the Database Service, the Simulation Control Service (which encapsulates the Simulation Library), and the Simulation Data Service. The architecture is designed such that the location of the user interface and the corresponding back end services are transparent to the majority of the application. The typical configuration is a Java applet running in a WWW browser, with the Database, Simulation Control, and Simulation Data services executing on a remote machine (WWW server). Alternatively, the software may be executed as a standalone application on a local machine with the requirement that the Java Runtime Environment and a C++ compiler are installed. |
(Click on image for larger view)
![]() |
| Physiological Models The physiology section allows for the specification of cellular structures (topology) and molecular species (with default initial concentrations and diffusion rates) within the cellular model under study. The model is further defined by a collection of biochemical reactions acting on molecular species within the cellular structures. Biochemical reactions are objects that represent complete descriptions of reaction stoichiometry and kinetics. Reactions can be either a collection of related reaction steps occurring at or near a single structure, or trans-membrane fluxes. Membranes are automatically generated when a feature is enclosed within another feature. |
![]() |
| Geometric Model geometric model represents the morphometry of a particular cell or a portion of a cell. Such a model is required to represent the cellular system. The geometry may be captured by a variety of imaging modalities, or it may be analytically defined for very regular structures or symmetric cells. Internal structures that are included in the physiological model, but are not spatially resolved, may be mapped continuously as a volume fraction of the enclosing compartment. The geometry may be specified as compartmental, 1-D, 2-D or 3-D segmented images. Segmentation is necessary to define discrete regions that represent the anatomical features of interest. Membranes are implicitly represented as the boundary between dissimilar compartments. A compartmental model represents a single point simulation based on the defined physiological model, the surface to volume ratio and volume fractions. Compartmental models are useful in calculating quick approximate solutions. |
![]() |
| Applications The regions defined in the geometry model are mapped to structures defined in the physiology model. Compartmental models represent a single point simulation and multiple regions within such a model are still mapped to a single compartment. The application allows for the specification of spatially accurate initial conditions and boundary conditions for each species of each structure. It is also a means of isolating simulation specific assumptions from the rest of the model. Electrical mapping is also defined within the application component of the software. |
![]() |
| Simulations The Equation Viewer displays the equations generated as a result of mapping the physiological model to either a cellular geometry model (spatial simulation) or a single point approximation (compartmental model). The parameter values may be substituted (and the expression simplified), or left in their symbolic representation. |
![]() |
| Compartmental Deterministic Simulation The Compartmental Deterministic Simulation component executes a compartmental (single point) simulation based on the defined physiological model and the geometric assumptions entered in the Structure Mapping (surface to volume ratios and volume fractions). The simulation is based on a set of nonlinear ordinary differential equations (ODE) that typically are solved in seconds. This allows an interactive, though manual, modification of parameters and a quick determination of the effect over time. Once the simulation is complete, the results, either for a single species or for multiple species, can be viewed easily. The simulation tool includes a local sensitivity analysis option. The analysis is performed by direct solution for variable sensitivities to a small variation of a single simulation parameter. Equations for the sensitivities are generated by symbolic differentiation. Normalized (log) sensitivities are displayed as functions of time. For a quick preview, the perturbed state trajectories are obtained by extrapolation. |
![]() |
| Spatial Simulation Partial differential equations that correspond to diffusive species and ordinary differential equations for non-diffusive species are generated for the solution of a complete spatial simulation. The equations are sent to a remote server where an executable is generated, run and the results are collected and stored. The simulation data server coordinates the client access to the server-sided simulation for display and analysis. The analysis capability available includes graphing the spatial distribution of a species as a line scan and graphing a time series at a single point. |
(Click on Image for Larger View)![]() |
| Storage The current implementation allows whole physiological models, geometric models, and applications to be stored and retrieved. The application is stored such that it encapsulates all of the information to reproduce and describe a particular spatial simulation. There is, however, no ability to querying the stored models for specific attributes. The models are currently stored intact using Java's Object Serialization capability. Experimental segmented images are uploaded to the image database. Within the database the user can assign spatial scales and identify regions based on anatomical features. |
Virtual Cell Hardware
|
| Data Export, Analysis, and Interoperability The simulation results are in the form of multivariate, time series data that can have up to 3 spatial dimensions. To allow flexibility in analysis, exporting facilities have been introduced. Data may be exported as spreadsheet data (comma separated value .csv), image files (GIF, NRRD), or movies (Quicktime or Animated GIFs). The capability for exporting and importing models between the Virtual Cell format and SMBL Levels 1 and 2 and limited capability to import from CellML to a math model is also currently available. Equations from the Virtual Cell can be exported into MATLAB (version 5 or 6) and reports can be generated via .pdf, .rtf or .html. |
![]() |

| Technology | Modeling Process | ||
| Software Architecture | Modeling Framework | ||
| Math Framework | User Interface | ||
| Testing Framework | VCML Specification |