LIBSA--a method for the determination of ligand-binding preference to allosteric sites on receptor ensembles.

Incorporation of receptor flexibility into computational drug discovery through the relaxed complex scheme is well suited for screening against a single binding site. In the absence of a known pocket or if there are multiple potential binding sites, it may be necessary to do docking against the entire surface of the target (global docking). However no suitable and easy-to-use tool is currently available to rank global docking results based on the preference of a ligand for a given binding site. We have developed a protocol, termed LIBSA for LIgand Binding Specificity Analysis, that analyzes multiple docked poses against a single or ensemble of receptor conformations and returns a metric for the relative binding to a specific region of interest. By using novel filtering algorithms and the signal-to-noise ratio (SNR), the relative ligand-binding frequency at different pockets can be calculated and compared quantitatively. Ligands can then be triaged by their tendency to bind to a site instead of ranking by affinity alone. The method thus facilitates screening libraries of ligand cores against a large library of receptor conformations without prior knowledge of specific pockets, which is especially useful to search for hits that selectively target a particular site. We demonstrate the utility of LIBSA by showing that it correctly identifies known ligand binding sites and predicts the relative preference of a set of related ligands for different pockets on the same receptor.

Isoform selectivity of adenylyl cyclase inhibitors: characterization of known and novel compounds.

Nine membrane-bound adenylyl cyclase (AC) isoforms catalyze the production of the second messenger cyclic AMP (cAMP) in response to various stimuli. Reduction of AC activity has well documented benefits, including benefits for heart disease and pain. These roles have inspired development of isoform-selective AC inhibitors, a lack of which currently limits exploration of functions and/or treatment of dysfunctions involving AC/cAMP signaling. However, inhibitors described as AC5- or AC1-selective have not been screened against the full panel of AC isoforms. We have measured pharmacological inhibitor profiles for all transmembrane AC isoforms. We found that 9-(tetrahydro-2-furanyl)-9H-purin-6-amine (SQ22,536), 2-amino-7-(furanyl)-7,8-dihydro-5(6H)-quinazolinone (NKY80), and adenine 9-ß-d-arabinofuranoside (Ara-A), described as supposedly AC5-selective, do not discriminate between AC5 and AC6, whereas the putative AC1-selective inhibitor 5-[[2-(6-amino-9H-purin-9-yl)ethyl]amino]-1-pentanol (NB001) does not directly target AC1 to reduce cAMP levels. A structure-based virtual screen targeting the ATP binding site of AC was used to identify novel chemical structures that show some preference for AC1 or AC2. Mutation of the AC2 forskolin binding pocket does not interfere with inhibition by SQ22,536 or the novel AC2 inhibitor, suggesting binding to the catalytic site. Thus, we show that compounds lacking the adenine chemical signature and targeting the ATP binding site can potentially be used to develop AC isoform-specific inhibitors, and discuss the need to reinterpret literature using AC5/6-selective molecules SQ22,536, NKY80, and Ara-A.

Andrographolide derivatives inhibit guanine nucleotide exchange and abrogate oncogenic Ras function.

Aberrant signaling by oncogenic mutant rat sarcoma (Ras) proteins occurs in ∼15% of all human tumors, yet direct inhibition of Ras by small molecules has remained elusive. Recently, several small-molecule ligands have been discovered that directly bind Ras and inhibit its function by interfering with exchange factor binding. However, it is unclear whether, or how, these ligands could lead to drugs that act against constitutively active oncogenic mutant Ras. Using a dynamics-based pocket identification scheme, ensemble docking, and innovative cell-based assays, here we show that andrographolide (AGP)--a bicyclic diterpenoid lactone isolated from Andrographis paniculata--and its benzylidene derivatives bind to transient pockets on Kirsten-Ras (K-Ras) and inhibit GDP-GTP exchange. As expected for inhibitors of exchange factor binding, AGP derivatives reduced GTP loading of wild-type K-Ras in response to acute EGF stimulation with a concomitant reduction in MAPK activation. Remarkably, however, prolonged treatment with AGP derivatives also reduced GTP loading of, and signal transmission by, oncogenic mutant K-RasG12V. In sum, the combined analysis of our computational and cell biology results show that AGP derivatives directly bind Ras, block GDP-GTP exchange, and inhibit both wild-type and oncogenic K-Ras signaling. Importantly, our findings not only show that nucleotide exchange factors are required for oncogenic Ras signaling but also demonstrate that inhibiting nucleotide exchange is a valid approach to abrogating the function of oncogenic mutant Ras.

Novel allosteric sites on Ras for lead generation.

Aberrant Ras activity is a hallmark of diverse cancers and developmental diseases. Unfortunately, conventional efforts to develop effective small molecule Ras inhibitors have met with limited success. We have developed a novel multi-level computational approach to discover potential inhibitors of previously uncharacterized allosteric sites. Our approach couples bioinformatics analysis, advanced molecular simulations, ensemble docking and initial experimental testing of potential inhibitors. Molecular dynamics simulation highlighted conserved allosteric coupling of the nucleotide-binding switch region with distal regions, including loop 7 and helix 5. Bioinformatics methods identified novel transient small molecule binding pockets close to these regions and in the vicinity of the conformationally responsive switch region. Candidate binders for these pockets were selected through ensemble docking of ZINC and NCI compound libraries. Finally, cell-based assays confirmed our hypothesis that the chosen binders can inhibit the downstream signaling activity of Ras. We thus propose that the predicted allosteric sites are viable targets for the development and optimization of new drugs.

Effect of normalization methods on biomarker identification utilizing mass spectrometry.

In biomarker identification using mass spectrometry, normalization makes it possible to compare mass spectra obtained from different samples. However, the relative influence of different normalization methods is an unexplored topic. In this study, we compared the most widely used normalization methods in a systemic manner to investigate impact of normalization. According to our experimental results, different normalization methods result in the selection of different features of different discriminatory power.