A threonine turnstile defines a dynamic amphiphilic binding motif in the AAA ATPase p97 allosteric binding site.

The turnstile motion of two neighboring threonines sets up a dynamic side chain interplay that can accommodate both polar and apolar ligands in a small molecule allosteric protein binding site. A computational model based on SAR data and both X-ray and cryo-EM structures of the AAA ATPase p97 was used to analyze the effects of paired threonines at the inhibitor site. Specifically, the Thr side chain hydroxyl groups form a hydrogen bonding network that readily accommodates small, highly polar ligand substituents. Conversely, diametric rotation of the χ1 torsion by 150-180° orients the side chain ß-methyl groups into the binding cleft, creating a hydrophobic pocket that can accommodate small, apolar substituents. This motif was found to be critical for rationalizing the affinities of a structurally focused set of inhibitors of p97 covering a > 2000-fold variation in potencies, with a preference for either small-highly polar or small-apolar groups. The threonine turnstile motif was further validated by a PDB search that identified analogous binding modes in ligand interactions in PKB, as well as by an analysis of NMR structures demonstrating additional gear-like interactions between adjacent Thr pairs. Combined, these data suggest that the threonine turnstile motif may be a general feature of interest in protein binding pockets.

Calculation methods for the enhancement of pharmaceutical properties in small molecules: estimating the cationic pKa.

In this review, a summary of methodologies is covered to enable medicinal chemists to access an overview of pK(a) estimation devices. In order to stave overutilization of costly synthetic resources, the chemist requires an accurate and computationally tractable solution for estimating a pK(a) of a candidate molecule. We focus on the cationic moieties, since they are so fundamentally important in the chemistry of drugs, and possess unique requirements to obtain a reasonably reliable pK(a) estimation.

Three-dimensional database mining identifies a unique chemotype that unites structurally diverse botulinum neurotoxin serotype A inhibitors in a three-zone pharmacophore.

A search query consisting of two aromatic centers and two cationic centers was defined based on previously identified small molecule inhibitors of the botulinum neurotoxin serotype A light chain (BoNT/A LC) and used to mine the National Cancer Institute Open Repository. Ten small molecule hits were identified, and upon testing, three demonstrated inhibitory activity. Of these, one was structurally unique, possessing a rigid diazachrysene scaffold. The steric limitations of the diazachrysene imposed a separation between the overlaps of previously identified inhibitors, revealing an extended binding mode. As a result, the pharmacophore for BoNT/A LC inhibition has been modified to encompass three zones. To demonstrate the utility of this model, a novel three-zone inhibitor was mined and its activity was confirmed.

A refined pharmacophore identifies potent 4-amino-7-chloroquinoline-based inhibitors of the botulinum neurotoxin serotype A metalloprotease.

We previously identified structurally diverse small molecule (non-peptidic) inhibitors (SMNPIs) of the botulinum neurotoxin serotype A (BoNT/A) light chain (LC). Of these, several (including antimalarial drugs) contained a 4-amino-7-chloroquinoline (ACQ) substructure and a separate positive ionizable amine component. The same antimalarials have also been found to interfere with BoNT/A translocation into neurons, via pH elevation of the toxin-mediated endosome. Thus, this structural class of small molecules may serve as dual-function BoNT/A inhibitors. In this study, we used a refined pharmacophore for BoNT/A LC inhibition to identify four new, potent inhibitors of this structural class (IC50's ranged from 3.2 to 17 muM). Molecular docking indicated that the binding modes for the new SMNPIs are consistent with those of other inhibitors that we have identified, further supporting our structure-based pharmacophore. Finally, structural motifs of the new SMNPIs, as well as two structure-based derivatives, were examined for activity, providing valuable information about pharmacophore component contributions to inhibition.

Inhibition of metalloprotease botulinum serotype A from a pseudo-peptide binding mode to a small molecule that is active in primary neurons.

An efficient research strategy integrating empirically guided, structure-based modeling and chemoinformatics was used to discover potent small molecule inhibitors of the botulinum neurotoxin serotype A light chain. First, a modeled binding mode for inhibitor 2-mercapto-3-phenylpropionyl-RATKML (K(i) = 330 nM) was generated, and required the use of a molecular dynamic conformer of the enzyme displaying the reorientation of surface loops bordering the substrate binding cleft. These flexible loops are conformationally variable in x-ray crystal structures, and the model predicted that they were pivotal for providing complementary binding surfaces and solvent shielding for the pseudo-peptide. The docked conformation of 2-mercapto-3-phenylpropionyl-RATKML was then used to refine our pharmacophore for botulinum serotype A light chain inhibition. Data base search queries derived from the pharmacophore were employed to mine small molecule (non-peptidic) inhibitors from the National Cancer Institute's Open Repository. Four of the inhibitors possess K(i) values ranging from 3.0 to 10.0 microM. Of these, NSC 240898 is a promising lead for therapeutic development, as it readily enters neurons, exhibits no neuronal toxicity, and elicits dose-dependent protection of synaptosomal-associated protein (of 25 kDa) in a primary culture of embryonic chicken neurons. Isothermal titration calorimetry showed that the interaction between NSC 240898 and the botulinum A light chain is largely entropy-driven, and occurs with a 1:1 stoichiometry and a dissociation constant of 4.6 microM.

A common pharmacophore for a diverse set of colchicine site inhibitors using a structure-based approach.

Modulating the structure and function of tubulin and microtubules is an important route to anticancer therapeutics, and therefore, small molecules that bind to tubulin and cause mitotic arrest are of immense interest. A large number of synthetic and natural compounds with diverse structures have been shown to bind at the colchicine site, one of the major binding sites on tubulin, and inhibit tubulin assembly. Using the recently determined X-ray structure of the tubulin:colchicinoid complex as the template, we employed docking studies to determine the binding modes of a set of structurally diverse colchicine site inhibitors. These binding models were subsequently used to construct a comprehensive, structure-based pharmacophore that in combination with molecular dynamics simulations confirms and extends our understanding of binding interactions at the colchicine site.

An all-atom model of the pore-like structure of hexameric VP40 from Ebola: structural insights into the monomer-hexamer transition.

The matrix protein VP40 is an indispensable component of viral assembly and budding by the Ebola virus. VP40 is a monomer in solution, but can fold into hexameric and octameric states, two oligomeric conformations that play central roles in the Ebola viral life cycle. While the X-ray structures of monomeric and octameric VP40 have been determined, the structure of hexameric VP40 has only been solved by three-dimensional electron microscopy (EM) to a resolution of approximately 30 A. In this paper, we present the refinement of the EM reconstruction of truncated hexameric VP40 to approximately 20 A and the construction of an all-atom model (residues 44-212) using the EM model at approximately 20 A and the X-ray structure of monomeric VP40 as templates. The hexamer model suggests that the monomer-hexamer transition involves a conformational change in the N-terminal domain that is not evident during octamerization and therefore, may provide the basis for elucidating the biological function of VP40.

Conformational sampling of the botulinum neurotoxin serotype A light chain: implications for inhibitor binding.

Botulinum neurotoxins (BoNTs) are the most potent of the known biological toxins, and consequently are listed as category A biowarfare agents. Currently, the only treatments against BoNTs include preventative antitoxins and long-term supportive care. Consequently, there is an urgent need for therapeutics to counter these enzymes--post exposure. In a previous study, we identified a number of small, nonpeptidic lead inhibitors of BoNT serotype A light chain (BoNT/A LC) metalloprotease activity, and we identified a common pharmacophore for these molecules. In this study, we have focused on how the dynamic movement of amino acid residues in and surrounding the substrate binding cleft of the BoNT/A LC might affect inhibitor binding modes. The X-ray crystal structures of two BoNT/A LCs (PDB refcodes=3BTA and 1E1H) were examined. Results from these analyses indicate that the core structural features of the examined BoNT/A LCs, including alpha-helices and beta-sheets, remained relatively unchanged during 1 ns dynamics trajectories. However, conformational flexibility was observed in surface loops bordering the substrate binding clefts in both examined structures. Our analyses indicate that these loops may possess the ability to decrease the solvent accessibility of the substrate binding cleft, while at the same time creating new residue contacts for the inhibitors. Loop movements and conformational/positional analyses of residues within the substrate binding cleft are discussed with respect to BoNT/A LC inhibitor binding and our common pharmacophore for inhibition. The results from these studies may aid in the future identification/development of more potent small molecule inhibitors that take advantage of new binding contacts in the BoNT/A LC.

Identification of small molecule inhibitors of anthrax lethal factor.

The virulent spore-forming bacterium Bacillus anthracis secretes anthrax toxin composed of protective antigen (PA), lethal factor (LF) and edema factor (EF). LF is a Zn-dependent metalloprotease that inactivates key signaling molecules, such as mitogen-activated protein kinase kinases (MAPKK), to ultimately cause cell death. We report here the identification of small molecule (nonpeptidic) inhibitors of LF. Using a two-stage screening assay, we determined the LF inhibitory properties of 19 compounds. Here, we describe six inhibitors on the basis of a pharmacophoric relationship determined using X-ray crystallographic data, molecular docking studies and three-dimensional (3D) database mining from the US National Cancer Institute (NCI) chemical repository. Three of these compounds have K(i) values in the 0.5-5 microM range and show competitive inhibition. These molecular scaffolds may be used to develop therapeutically viable inhibitors of LF.

Novel small molecule inhibitors of botulinum neurotoxin A metalloprotease activity.

Botulinum neurotoxins (BoNTs) are among the most lethal biological substances to have been weaponized and are listed as biodefense category A agents. Currently, no small molecule (non-peptidic) therapeutics exist to counter this threat; hence, identifying and developing compounds that inhibit BoNTs is a high priority. In the present study, a high-throughput assay was used to identify small molecules that inhibit the metalloprotease activity of BoNT serotype A light chain (BoNT/A LC). All inhibitors were further verified using a HPLC-based assay. Conformational analyses of these compounds, in conjunction with molecular docking studies, were used to predict structural features that contribute to inhibitor binding and potency. Based on these results, a common pharmacophore for BoNT/A LC inhibitors is proposed. This is the first study to report small molecules (non-peptidics) that inhibit BoNT/A LC metalloprotease activity in the low microM range.