Previous topic

Biochemical Reaction Network Modeling

Next topic

Simulation Results

Modeling Framework StructureΒΆ

The figure below shows the conceptual flow of information through BioTRaNS. A Chemical Mixture (concentration of a single or multiple chemicals) and Reaction Environment (types and amounts of enzymes and other agents) are specified by the user. The Chemical Mixture is converted to a set of Virtual Molecules (vMols), each with its own specifications (i.e., a canonical SMILES and list of computed geometric, energetic, and physicochemical properties). SMILES (Simplified Molecular Input Line Entry Specification) is a general-purpose chemical nomenclature and data-exchange format used in computational chemistry. In general, properties are retrieved from associated databases when available, or are computed as needed.

The Reaction Environment is translated into a set of Virtual Enzymes (vEnzs), representing the actual enzymes and Virtual Agents (vAgnts), representing nonenzymatic reactions. Each vEnz is endowed with an appropriate binding calculator and transforms governing the biotransformations that the vEnz can mediate. The binding calculator computes the feasibility of a particular vMol binding to a vEnz, based on quantitative structure-activity relationships (or decision tree) derived from properties of known substrates for each enzyme. Transformations are stored as a list of SMIRKS-based representations. SMIRKS (Smiles Reac[k]tion Specification) is a superset of the SMILES language that describes chemical transformations. Our methodology is unique in implementing transforms that describe steps of a reaction mechanism, including those for enzyme-mediated reactions. The use of mechanistic steps is key to automatically and accurately predicting metabolites (vide infra for an example). After the vMol is created, the binding calculator for each vEnz assesses binding feasibility; if feasible, the vMol becomes an eligible reactant. All eligible reactants then undergo appropriate virtual biotransformations, creating specific chemical reactions and converting the vMol into one or more metabolites, which also are checked for reaction feasibility. The information on the metabolites and the associated transform/agent interconnecting each substrate-product pair is converted to a Virtual Reaction (vRxn). Thus, each vRxn comprises the vMols and vEnzs from which it was derived, as well as appropriate kinetic properties.

Based on the aggregate information in the vRxns, a reaction network is generated in the form of a mathematical graph, wherein the metabolites constitute the vertices or nodes, and the reactions form the edges. Well-described in discrete mathematics, graphs can be efficiently analyzed using algorithms of computer science; graph theory allows efficient determination of characteristics such as common reactants between various pathways, the shortest path to a given metabolite, and common (or all) routes to given metabolites. These network graphs can be depicted at different levels of detail, from highly detailed, showing steps of a reaction mechanism, to overview, showing only key metabolites. Thus, highly complex pathways can be easily simplified. This ability to delineate the metabolites and visualize where pathways intersect is a key tool in gaining new insight into interrelationships among chemicals and their metabolic pathways.

Simultaneous to the graph creation, kinetic properties in each vRxn are used to create the appropriate reaction rate equations (ordinary differential equations, ODE).

These properties include

  1. rate constants (e.g., Michaelis constant, Km, and maximum velocity, Vmax, for enzyme-catalyzed reactions and k for nonenzymatic reactions),
  2. inhibitor constants, Ki, and
  3. modes of inhibition or allosterism.

The total set of rate equations and specified initial conditions forms an initial value problem that is solved by a stiff ODE equation solver for the concentrations of all species as a function of time. Virtual Enzymes and Other Virtual Agents. The constituent transforms for the each virtual enzyme are compiled by carefully culling the literature for data on enzymes known to act on the chemicals and chemical metabolites of interest.

_images/brnm_flow.png