Simulating botulinum neurotoxin with constant pH molecular dynamics in Generalized Born implicit solvent
Introduction
Botulinum neurotoxins (BoNTs) are among the most potent toxins to human beings. There are seven different serotypes (A-G) of BoNT [1]. They were believed causative agents of botulism, a potentially fatal condition of neuromuscular paralysis and researched by many groups [2], [3], [4], [5], [6]. Currently no effective chemical antidote is available against botulism principally because of the lack of knowledge of the molecular structures and the mechanism of toxicity [1].
In the whole three-step toxic action, translocation domain across membranes attracts great attention since there are interesting phenomena in channel forming and enormous size transport. An acidic environment is believed to induce conformational changes in the translocation domain. Unfortunately, there was no experimental approach available for the range pH 4.0 to 5.0 [3]. For this reason, the property of BoNT/A in this range is of greater interests and an efficient constant pH MD modeling technique is highly desirable.
Section snippets
Constant pH MD
It has been long and well known that protein structures are strongly dependent on solvent pH [7], [8]. On the other hand, it is known that a structural change also affects the values of titratable residues. Therefore, there is a tight coupling between protein conformation and protonation state [8]. Modeling such correlation requires a reliable constant pH molecular dynamics method and such a method was indeed proposed and implemented in AMBER 8 package by Mongan et al. [9]. This method
Numerical experiment
Our MD experiments were performed by using AMBER 8. The ff99 force field was employed. A refined GB model () was used for solvation. Salt concentration was set at 0.2 M. The cutoff for nonbonded interactions and computation of effective Born radii was 30 Å. The time step was 2 fs. Sander [9] was used as MD engine and ptraj [9] was used to analyze the trajectories and calculate the RMSD values in each case.
Results
A valid constant pH MD simulation method should yield predictions consistent with experimental values. It can rapidly converge to those predictions and maintain a stable trajectory. Furthermore, the method should be computationally efficient [8].
Analysis and conclusions
To explain these results, the following reasons are considered.
(1) Sufficiently conformational sampling is infeasible for the system of this size.
The instantaneous is strongly dependent on conformation, so it is reasonable that they would have different protonation state population if two simulations sample conformation space differently. The random error due to incomplete conformational sampling may produce larger effects than those caused by a small change in pH [8]. Of all 1277 residues
Acknowledgement
We thank Brookhaven National Laboratory (BNL) for their support and assistance.
References (11)
- et al.
Recombinant expression and purification of the botulinum neurotoxin type A translocation domain
Protein Expression and Purification
(1997) - et al.
A common subunit structure in Clostridium botulinum type A, B and E toxins
Biochem. Biophys. Res. Commun.
(1972) Intimate details of the most poisonous poison
Nat. Struct. Biol.
(2000)Clostridium botulinum neurotoxin
Micro-Biol. Rev.
(1980)- et al.
Spectroscopic analysis of pH-induced changes in the molecular features of type a botulinum neurotoxin light chain
Biochemistry
(2000)