Botulinum neurotoxins have remarkable persistence (∼weeks to months in cells), outlasting the small-molecule inhibitors designed to target them. To address this disconnect, inhibitors bearing two pharmacophores-a zinc binding group and a Cys-reactive warhead-were designed to leverage both affinity and reactivity. A series of first-generation bifunctional inhibitors was achieved through structure-based inhibitor design. Through X-ray crystallography, engagement of both the catalytic Zn2+ and Cys165 was confirmed. A second-generation series improved on affinity by incorporating known reversible inhibitor pharmacophores; the mechanism was confirmed by exhaustive dialysis, mass spectrometry, and in vitro evaluation against the C165S mutant. Finally, a third-generation inhibitor was shown to have good cellular activity and low toxicity. In addition to our findings, an alternative method of modeling time-dependent inhibition that simplifies assay setup and allows comparison of inhibition models is discussed.

The lethal factor (LF) enzyme secreted by Bacillus anthracis is a zinc hydrolase that is chiefly responsible for anthrax-related cell death. Although many studies of the design of small molecule LF inhibitors have been conducted, no LF inhibitor is yet available as a therapeutic agent. Inhibitors with considerable chemical diversity have been developed and investigated; however, the LF S2' subsite has not yet been systematically explored as a potential target for lead optimization. Here we present synthesis, experimental evaluation, modeling, and structural biology for a novel series of sulfonamide hydroxamate LF inhibitor analogues specifically designed to extend into, and probe chemical preferences of, this S2' subsite. We discovered that this region accommodates a wide variety of chemical functionalities and that a broad selection of ligand structural modifications directed to this area can be incorporated without significant deleterious alterations in biological activity. We also identified key residues in this subsite that can potentially be targeted to improve inhibitor binding.

The Bacillus anthracis lethal factor (LF) is one component of a tripartite exotoxin partly responsible for persistent anthrax cytotoxicity after initial bacterial infection. Inhibitors of the zinc metalloproteinase have been investigated as potential therapeutic agents, but LF is a challenging target because inhibitors lack sufficient selectivity or possess poor pharmaceutical properties. These structural studies reveal an alternate conformation of the enzyme, induced upon binding of specific inhibitors, that opens a previously unobserved deep pocket termed S1'(∗) which might afford new opportunities to design selective inhibitors that target this subsite.

The secreted anthrax toxin consists of three components: the protective antigen (PA), edema factor (EF) and lethal factor (LF). LF, a zinc metalloproteinase, compromises the host immune system primarily by targeting mitogen-activated protein kinase kinases in macrophages. Peptide substrates and small-molecule inhibitors bind LF in the space between domains 3 and 4 of the hydrolase. Domain 3 is attached on a hinge to domain 2 via residues Ile300 and Pro385, and can move through an angular arc of greater than 35° in response to the binding of different ligands. Here, multiple LF structures including five new complexes with co-crystallized inhibitors are compared and three frequently populated LF conformational states termed `bioactive', `open' and `tight' are identified. The bioactive position is observed with large substrate peptides and leaves all peptide-recognition subsites open and accessible. The tight state is seen in unliganded and small-molecule complex structures. In this state, domain 3 is clamped over certain substrate subsites, blocking access. The open position appears to be an intermediate state between these extremes and is observed owing to steric constraints imposed by specific bound ligands. The tight conformation may be the lowest-energy conformation among the reported structures, as it is the position observed with no bound ligand, while the open and bioactive conformations are likely to be ligand-induced.

In the present paper, we apply the electrostatically embedded many-body expansion of the correlation energy (EE-MB-CE) to the calculation of zinc-ligand and cadmium-ligand bond dissociation energies, and we analyze the errors due to various fragmentation schemes in a variety of neutral, positively charged, and negatively charged Zn2+ and Cd2+ coordination complexes. As a result of the analysis, we are able to present a new, simple, and unambiguous fragmentation strategy. Following this strategy, we show that both methods perform well for zinc-ligand and cadmium-ligand bond dissociation energies for all systems studied in the paper, including a model of the catalytic site of the zinc-bearing anthrax toxin lethal factor (LF), which has garnered substantial attention as a target for drug development. To draw general conclusions we consider ten pentacoordinate and hexacoordinate zinc and cadmium containing coordination complexes, each with 10 or 15 different fragmentation schemes. By analyzing errors, we developed a prescription for the optimal fragmentation strategy. With this scheme, and using MP2 correlation energies as a test, we find that the electrostatically embedded three-body expansion of the correlation energy (EE-3B-CE) method is able to reproduce all 53 conventionally calculated bond energies with an average absolute error of only 0.59 kcal/mol. The paper also presents EE-MB-CE calculations using the CCSD(T) level of theory on an LF model system. With CCSD(T), EE-3B-CE has an average error of 0.30 kcal/mol.

The electrostatically embedded many-body (EE-MB) method has proven accurate for calculating cohesive and conformational energies in clusters, and it has recently been extended to obtain bond dissociation energies for metal-ligand bonds in positively charged inorganic coordination complexes. In the present paper, we present four key guidelines that maximize the accuracy and efficiency of EE-MB calculations for metal centers. Then, following these guidelines, we show that the EE-MB method can also perform well for bond dissociation energies in a variety of neutral and negatively charged inorganic coordination systems representing metalloenzyme active sites, including a model of the catalytic site of the zinc-bearing anthrax toxin lethal factor, a popular target for drug development. In particular, we find that the electrostatically embedded three-body (EE-3B) method is able to reproduce conventionally calculated bond-breaking energies in a series of pentacoordinate and hexacoordinate zinc-containing systems with an average absolute error (averaged over 25 cases) of only 0.98 kcal/mol.