Molecular dynamics simulations of complexes between wild-type and mutant anthrax protective antigen variants and a model anthrax toxin receptor

Bacillus anthracis, a spore-forming infectious bacterium, produces a toxin consisting of three proteins: lethal factor (LF), edema factor (EF), and protective antigen (PA). LF and EF possess intracellular enzymatic functions, the net effect of which is to severely compromise host innate immunity. During an anthrax infection PA plays the critical role of facilitating entry of both EF and LF toxins into host cell cytoplasm. Crystal structures of all three of the anthrax toxins have been determined, as well as the crystal structure of the (human) von Willebrand factor A (integrin VWA/I domain) - an anthrax toxin receptor. A theoretical structure of the complex between VWA/I and PA has also been reported. Here we report on the results of 1,000 psec molecular dynamics (MD) simulations carried out on complexes between the Anthrax Protective Antigen Domain 4 (PA-D4) and the von Willebrand Factor A (VWA/I). MD simulations (using Insight II software) were carried out for complexes containing wildtype (WT) PA-D4, as well as for complexes containing three different mutants of PA-D4, one containing three substitutions in the PA-D4 "small loop" (residues 679-693) (D683A/L685E/Y688C), one containing a single substitution at a key site at the PA-D4 - receptor interface (K679A) and another containing a deletion of eleven residues at the C-terminus of PA (triangle724-735). All three sets of PA mutations have been shown experimentally to result in serious deficiencies in PA function. Our MD results are consistent with these findings. Major disruptions in interactions were observed between the mutant PA-D4 domains and the anthrax receptor during the MD simulations. Many secondary structural features in PA-D4 are also severely compromised when VWA complexes with mutant variants of PA-D4 are subjected to MD simulations. These MD simulation results clearly indicate the importance of the mutated PA-D4 residues in both the "small loop" and at the carboxyl terminus in maintaining a PA conformation that is capable of effective interaction with the anthrax toxin receptor.