BioModels Database logo

BioModels Database

spacer

BIOMD0000000178 - Lebeda2008_BoTN_Paralysis_4stepModel

 

 |   |   |  Send feedback
Reference Publication
Publication ID: 18551355
Lebeda FJ, Adler M, Erickson K, Chushak Y.
Onset dynamics of type A botulinum neurotoxin-induced paralysis.
J Pharmacokinet Pharmacodyn 2008 Jun; 35(3): 251-267
Integrated Toxicology Division, US Army Medical Research Institute of Infectious Diseases, Fort Detrick, MD 21702-5011, USA. Frank.Lebeda@amedd.army.mil  [more]
Model
Original Model: BIOMD0000000178.xml.origin
Submitter: Lukas Endler
Submission ID: MODEL8472829123
Submission Date: 11 Aug 2008 16:44:46 UTC
Last Modification Date: 01 Apr 2014 21:21:29 UTC
Creation Date: 08 Aug 2008 15:57:53 UTC
Encoders:  Lukas Endler
   Frank Lebeda
set #1
bqmodel:isDerivedFrom BioModels Database Lebeda2008_BoNT_Paralysis_3stepModel
PubMed 6243359
set #2
bqbiol:occursIn Taxonomy Mus musculus
Taxonomy Rattus norvegicus
Taxonomy Homo sapiens
set #3
bqbiol:isPartOf ICD A05.1
bqbiol:isVersionOf Reactome REACT_11184.1
Gene Ontology bontoxilysin activity
Gene Ontology negative regulation of synaptic transmission, cholinergic
Gene Ontology receptor-mediated endocytosis
Gene Ontology pathogenesis
Notes

This model is the 4-step model from the article:
Onset dynamics of type A botulinum neurotoxin-induced paralysis.
Lebeda FJ, Adler M, Erickson K, Chushak Y J Pharmacokinet Pharmacodyn2008 Jun; 35(3): 251-67 18551355,
Abstract:
Experimental studies have demonstrated that botulinum neurotoxin serotype A (BoNT/A) causes flaccid paralysis by a multi-step mechanism. Following its binding to specific receptors at peripheral cholinergic nerve endings, BoNT/A is internalized by receptor-mediated endocytosis. Subsequently its zinc-dependent catalytic domain translocates into the neuroplasm where it cleaves a vesicle-docking protein, SNAP-25, to block neurally evoked cholinergic neurotransmission. We tested the hypothesis that mathematical models having a minimal number of reactions and reactants can simulate published data concerning the onset of paralysis of skeletal muscles induced by BoNT/A at the isolated rat neuromuscular junction (NMJ) and in other systems. Experimental data from several laboratories were simulated with two different models that were represented by sets of coupled, first-order differential equations. In this study, the 3-step sequential model developed by Simpson (J Pharmacol Exp Ther 212:16-21,1980) was used to estimate upper limits of the times during which anti-toxins and other impermeable inhibitors of BoNT/A can exert an effect. The experimentally determined binding reaction rate was verified to be consistent with published estimates for the rate constants for BoNT/A binding to and dissociating from its receptors. Because this 3-step model was not designed to reproduce temporal changes in paralysis with different toxin concentrations, a new BoNT/A species and rate (k(S)) were added at the beginning of the reaction sequence to create a 4-step scheme. This unbound initial species is transformed at a rate determined by k(S) to a free species that is capable of binding. By systematically adjusting the values of k(S), the 4-step model simulated the rapid decline in NMJ function (k(S) >or= 0.01), the less rapid onset of paralysis in mice following i.m. injections (k (S) = 0.001), and the slow onset of the therapeutic effects of BoNT/A (k(S) < 0.001) in man. This minimal modeling approach was not only verified by simulating experimental results, it helped to quantitatively define the time available for an inhibitor to have some effect (t(inhib)) and the relation between this time and the rate of paralysis onset. The 4-step model predicted that as the rate of paralysis becomes slower, the estimated upper limits of (t(inhib)) for impermeable inhibitors become longer. More generally, this modeling approach may be useful in studying the kinetics of other toxins or viruses that invade host cells by similar mechanisms, e.g., receptor-mediated endocytosis.

Model updated by Viji on 07/09/2010.

This model is the extended model of BIOMD0000000267, which itself is the reduced form of the model developed by Simpson 1980; PMID 6243359

.

This model originates from BioModels Database: A Database of Annotated Published Models. It is copyright (c) 2005-2010 The BioModels Team.
For more information see the terms of use.
To cite BioModels Database, please use Le Novère N., Bornstein B., Broicher A., Courtot M., Donizelli M., Dharuri H., Li L., Sauro H., Schilstra M., Shapiro B., Snoep J.L., Hucka M. (2006) BioModels Database: A Free, Centralized Database of Curated, Published, Quantitative Kinetic Models of Biochemical and Cellular Systems Nucleic Acids Res., 34: D689-D691.

Model
Publication ID: 18551355 Submission Date: 11 Aug 2008 16:44:46 UTC Last Modification Date: 01 Apr 2014 21:21:29 UTC Creation Date: 08 Aug 2008 15:57:53 UTC
Mathematical expressions
Reactions
endocytosis translocation binding bulk_movement
Rules
Assignment Rule (variable: tension) Assignment Rule (variable: tot_free_BoNT/A)    
Physical entities
Compartments Species
extracellular tot_free_BoNT/A bulk_BoNT/A free_BoNT/A
bound_BoNT/A    
endosome transloc_BoNT/A    
neuroplasm lytic_BoNT/A    
Global parameters
tension      
Reactions (4)
 
 endocytosis [bound_BoNT/A] → [transloc_BoNT/A];  
 
 translocation [transloc_BoNT/A] → [lytic_BoNT/A];  
 
 binding [free_BoNT/A] → [bound_BoNT/A];  
 
 bulk_movement [bulk_BoNT/A] → [free_BoNT/A];  
 
Rules (2)
 
 Assignment Rule (name: tension) tension = 1-lytic
 
 Assignment Rule (name: BoNT) tot_free_BoNT/A = bulk+free
 
 extracellular Spatial dimensions: 3.0  Compartment size: 1.0
 
  tot_free_BoNT/A
Compartment: extracellular
 
 bulk_BoNT/A
Compartment: extracellular
Initial concentration: 1.0
 
 free_BoNT/A
Compartment: extracellular
Initial concentration: 0.0
 
 bound_BoNT/A
Compartment: extracellular
Initial concentration: 0.0
 
 endosome Spatial dimensions: 3.0  Compartment size: 1.0
 
 transloc_BoNT/A
Compartment: endosome
Initial concentration: 0.0
 
 neuroplasm Spatial dimensions: 3.0  Compartment size: 1.0
 
 lytic_BoNT/A
Compartment: neuroplasm
Initial concentration: 0.0
 
Global Parameters (1)
 
   tension
Value: NaN
 
endocytosis (1)
 
   kT
Value: 0.141   (Units: perminute)
Constant
 
translocation (1)
 
   kL
Value: 0.013   (Units: perminute)
Constant
 
binding (1)
 
   kB
Value: 0.058   (Units: perminute)
Constant
 
bulk_movement (1)
 
   kS
Value: 1.5E-4   (Units: perminute)
Constant
 
Representative curation result(s)
Representative curation result(s) of BIOMD0000000178

Curator's comment: (updated: 07 Sep 2010 15:24:42 BST)

Figure 3 and Figure 4 of the reference publication has been reproduced. The model was integrated and simulated using Copasi v4.5.31.

spacer
spacer