Empirical valence bond simulations of the chemical mechanism of ATP to cAMP conversion by anthrax edema factor.

The two-metal catalysis by the adenylyl cyclase domain of the anthrax edema factor toxin was simulated using the empirical valence bond (EVB) quantum mechanical/molecular mechanical approach. These calculations considered the energetics of the nucleophile deprotonation and the formation of a new P-O bond in aqueous solution and in the enzyme-substrate complex present in the crystal structure models of the reactant and product states of the reaction. Our calculations support a reaction pathway that involves metal-assisted transfer of a proton from the nucleophile to the bulk aqueous solution followed by subsequent formation of an unstable pentavalent intermediate that decomposes into cAMP and pyrophosphate (PPi). This pathway involves ligand exchange in the first solvation sphere of the catalytic metal. At 12.9 kcal/mol, the barrier for the last step of the reaction, the cleavage of the P-O bond to PPi, corresponds to the highest point on the free energy profile for this reaction pathway. However, this energy is too close to the value of 11.4 kcal/mol calculated for the barrier of the nucleophilic attack step to reach a definitive conclusion about the rate-limiting step. The calculated reaction mechanism is supported by reasonable agreement between the experimental and calculated catalytic rate constant decrease caused by the mutation of the active site lysine 346 to arginine.