Onset dynamics of type A botulinum neurotoxin-induced paralysis

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 (kS) 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 kS to a free species that is capable of binding. By systematically adjusting the values of kS, the 4-step model simulated the rapid decline in NMJ function (kS ≥0.01), the less rapid onset of paralysis in mice following i.m. injections (kS = 0.001), and the slow onset of the therapeutic effects of BoNT/A (kS < 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 (tinhib) 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 (tinhib) 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

Light Chain Separated from the Rest of the Type A Botulinum Neurotoxin Molecule Is the Most Catalytically Active Form

Botulinum neurotoxins (BoNT) are the most potent of all toxins. The 50 kDa N-terminal endopeptidase catalytic light chain (LC) of BoNT is located next to its central, putative translocation domain. After binding to the peripheral neurons, the central domain of BoNT helps the LC translocate into cytosol where its proteolytic action on SNARE (soluble NSF attachment protein receptor) proteins blocks exocytosis of acetyl choline leading to muscle paralysis and eventual death. The translocation domain also contains 105 Å -long stretch of ∼100 residues, known as “belt,” that crosses over and wraps around the LC to shield the active site from solvent. It is not known if the LC gets dissociated from the rest of the molecule in the cytosol before catalysis. To investigate the structural identity of the protease, we prepared four variants of type A BoNT (BoNT/A) LC, and compared their catalytic parameters with those of BoNT/A whole toxin. The four variants were LC + translocation domain, a trypsin-nicked LC + translocation domain, LC + belt, and a free LC. Our results showed that Km for a 17-residue SNAP-25 (synaptosomal associated protein of 25 kDa) peptide for these constructs was not very different, but the turnover number (kcat) for the free LC was 6-100-fold higher than those of its four variants. Moreover, none of the four variants of the LC was prone to autocatalysis. Our results clearly demonstrated that in vitro, the LC minus the rest of the molecule is the most catalytically active form. The results may have implication as to the identity of the active, toxic moiety of BoNT/A in vivo