Edema Factor (EF) is one of three major toxins of anthrax. EF is an adenylyl cyclase that disrupts cell signaling by accelerating the conversion of ATP into cyclic-AMP. EF has a much higher catalytic rate than that of mammalian adenylyl cyclases (mACs). Crystal structures were obtained for mACs and EF, but the molecular basis for different catalytic activities remained poorly understood. In particular, the arrangement of the active site in EF is unclear in what concerns the number of ions present and the conformation of the substrate. Here, we use quantum mechanics-molecular mechanics simulations to estimate the free-energy profiles for the reaction catalyzed by EF and a mAC. We found that EF catalysis is possible, and faster than that of mACs, in both one and two Mg(2+)-ion-binding modes, providing adaptive plasticity to host-cell environments. In both enzymes, the reaction mechanisms are highly associative. However, mechanistic differences exist. In the mAC, the nucleophile oxygen (ATP-O3') is consistently coordinated to one of the Mg(2+) ions, increasing its acidity. In EF, on the other hand, this coordination is eventual and not essential for the reaction to proceed. The persistent coordination of O3' to the ion is favored in mACs by a greater ion partial charge. In EF, the reduced acidity of the O3' oxygen is compensated by the presence of the His351 residue for proton abstraction. As proton transfer in EF does not require persistent attachment of the substrate to an ion, the substrate (ATP) and transition state display greater conformational flexibilities. These greater flexibilities allow the sampling of lower-energy conformations and might represent an entropic advantage for catalytic efficiency.

No abstract

Human Mn-containing superoxide dismutase (hMnSOD) is a mitochondrial enzyme that metabolizes superoxide radical (O2(•-)). O2(•-) reacts at diffusional rates with nitric oxide to yield a potent nitrating species, peroxynitrite anion (ONOO(-)). MnSOD is nitrated and inactivated in vivo, with active site Tyr34 as the key oxidatively modified residue. We previously reported a k of ∼1.0 × 10(5) M(-1) s(-1) for the reaction of hMnSOD with ONOO(-) by direct stopped-flow spectroscopy and the critical role of Mn in the nitration process. In this study, we further established the mechanism of the reaction of hMnSOD with ONOO(-), including the necessary re-examination of the second-order rate constant by an independent method and the delineation of the microscopic steps that lead to the regio-specific nitration of Tyr34. The redetermination of k was performed by competition kinetics utilizing coumarin boronic acid, which reacts with ONOO(-) at a rate of ∼1 × 10(6) M(-1) s(-1) to yield the fluorescence product, 7-hydroxycoumarin. Time-resolved fluorescence studies in the presence of increasing concentrations of hMnSOD provided a k of ∼1.0 × 10(5) M(-1) s(-1), fully consistent with the direct method. Proteomic analysis indicated that ONOO(-), but not other nitrating agents, mediates the selective modification of active site Tyr34. Hybrid quantum-classical (quantum mechanics/molecular mechanics) simulations supported a series of steps that involve the initial reaction of ONOO(-) with Mn(III) to yield Mn(IV) and intermediates that ultimately culminate in 3-nitroTyr34. The data reported herein provide a kinetic and mechanistic basis for rationalizing how MnSOD constitutes an intramitochondrial target for ONOO(-) and the microscopic events, with atomic level resolution, that lead to selective and efficient nitration of critical Tyr34.

The fundamental object for studying a (bio)chemical reaction obtained from simulations is the free energy profile, which can be directly related to experimentally determined properties. Although quite accurate hybrid quantum (DFT based)-classical methods are available, achieving statistically accurate and well converged results at a moderate computational cost is still an open challenge. Here, we present and thoroughly test a hybrid differential relaxation algorithm (HyDRA), which allows faster equilibration of the classical environment during the nonequilibrium steering of a (bio)chemical reaction. We show and discuss why (in the context of Jarzynski's Relationship) this method allows obtaining accurate free energy profiles with smaller number of independent trajectories and/or faster pulling speeds, thus reducing the overall computational cost. Moreover, due to the availability and straightforward implementation of the method, we expect that it will foster theoretical studies of key enzymatic processes.