Effect of Binding Pose and Modeled Structures on SVMGen and GlideScore Enrichment of Chemical Libraries

Virtual screening consists of docking libraries of small molecules to a target protein followed by rank-ordering of the resulting structures using scoring functions. The ability of scoring methods to distinguish between actives and inactives depends on several factors that include the accuracy of the binding pose during the docking step and the quality of the three-dimensional structure of the target. Here, we build on our previous work to introduce a new scoring approach (SVMGen) that uses machine learning trained with features from statistical pair potentials obtained from three-dimensional crystal structures. We use SVMGen and GlideScore to explore how enrichment or rank-ordering is affected by binding pose accuracy. To that end, we create a validation set that consists strictly of proteins whose crystal structure was solved in complex with their inhibitors. For the rank-ordering studies, we use crystal structures from PDBbind along with corresponding binding affinity data provided in the database. In addition to binding pose, we investigate the effect of using modeled structures for the target on the enrichment performance of SVMGen and GlideScore. To accomplish this, we generated homology models for protein kinases in DUD-E for which crystal structures are available to enable comparison of enrichment between modeled and crystal structure. We also generate homology models for kinases in SARfari for which there are many known small-molecule inhibitors but no known crystal structure. These models are used to assess the ability of SVMGen and GlideScore to distinguish between actives and decoys. We focus our work on protein kinases considering the wealth of structural and binding affinity data that exists for this family of proteins

BIOCHEMICAL CHARACTERIZATION OF SMALL MOLECULES TARGETING RAL GTPASE

poster abstractThe Ral subfamily of GTPases consists of highly similar RalA and RalB isoforms that participate in diverse cellular functions including endocytosis, exocytosis, actin cytoskeletal dynamics, and transcription. A large body of evidence has implicated Ral GTPases with tumor cell growth, migration, and angiogenesis in bladder, prostate, lung, and pancreatic cancer. The purpose of this project was to target the activity of Ral GTPases and their association with effector proteins through the identification of small molecule inhibitors that block this interaction. In order to accomplish this, both direct binding to RalB as well as disruption of protein-protein interaction were investigated. The top 200 compounds from a larger computational library of 500,000 compounds targeting the RalBP1 binding site on RalB were tested. Differen-tial scanning fluorimetry (DSF) was used to measure the degree of direct binding between compound and protein through thermal melting shift. To measure disruption between RalB and RalBP1 by small molecules, a novel enzyme-linked immunosorbent assay (ELISA) was developed. Identification of a few key compounds binding to RalB as well as optimization of an ELISA assay for RalBP1 was accomplished. Further direction of this project would be to utilize the ELISA assay to test inhibition of the protein-protein interac-tion between RalB and RalBP1 using the top compounds from the DSF trials

Molecular Recognition in a Diverse Set of Protein-Ligand Interactions Studied with Molecular Dynamics Simulations and End-Point Free Energy Calculations

End-point free energy calculations using MM-GBSA and MM-PBSA provide a detailed understanding of molecular recognition in protein-ligand interactions. The binding free energy can be used to rank-order protein-ligand structures in virtual screening for compound or target identification. Here, we carry out free energy calculations for a diverse set of 11 proteins bound to 14 small molecules using extensive explicit-solvent MD simulations. The structure of these complexes was previously solved by crystallography and their binding studied with isothermal titration calorimetry (ITC) data enabling direct comparison to the MM-GBSA and MM-PBSA calculations. Four MM-GBSA and three MM-PBSA calculations reproduced the ITC free energy within 1 kcal•mol−1 highlighting the challenges in reproducing the absolute free energy from end-point free energy calculations. MM-GBSA exhibited better rank-ordering with a Spearman ρ of 0.68 compared to 0.40 for MM-PBSA with dielectric constant (ε = 1). An increase in ε resulted in significantly better rank-ordering for MM-PBSA (ρ = 0.91 for ε = 10). But larger ε significantly reduced the contributions of electrostatics, suggesting that the improvement is due to the non-polar and entropy components, rather than a better representation of the electrostatics. SVRKB scoring function applied to MD snapshots resulted in excellent rank-ordering (ρ = 0.81). Calculations of the configurational entropy using normal mode analysis led to free energies that correlated significantly better to the ITC free energy than the MD-based quasi-harmonic approach, but the computed entropies showed no correlation with the ITC entropy. When the adaptation energy is taken into consideration by running separate simulations for complex, apo and ligand (MM-PBSAADAPT), there is less agreement with the ITC data for the individual free energies, but remarkably good rank-ordering is observed (ρ = 0.89). Interestingly, filtering MD snapshots by pre-scoring protein-ligand complexes with a machine learning-based approach (SVMSP) resulted in a significant improvement in the MM-PBSA results (ε = 1) from ρ = 0.40 to ρ = 0.81. Finally, the non-polar components of MM-GBSA and MM-PBSA, but not the electrostatic components, showed strong correlation to the ITC free energy; the computed entropies did not correlate with the ITC entropy

Small-Molecule Binding Sites to Explore New Targets in the Cancer Proteome

The Cancer Genome Atlas (TCGA) offers an unprecedented opportunity to identify small-molecule binding sites on proteins with overexpressed mRNA levels that correlate with poor survival. Here, we analyze RNA-seq and clinical data for 10 tumor types to identify genes that are both overexpressed and correlate with patient survival. Protein products of these genes were scanned for binding sites that possess shape and physicochemical properties that can accommodate small-molecule probes or therapeutic agents (druggable). These binding sites were classified as enzyme active sites (ENZ), protein-protein interaction sites (PPI), or other sites whose function is unknown (OTH). Interestingly, the overwhelming majority of binding sites were classified as OTH. We find that ENZ, PPI, and OTH binding sites often occurred on the same structure suggesting that many of these OTH cavities can be used for allosteric modulation of enzyme activity or protein-protein interactions with small molecules. We discovered several ENZ (PYCR1, QPRT, and HSPA6) and PPI (CASC5, ZBTB32, and CSAD) binding sites on proteins that have been seldom explored in cancer. We also found proteins that have been extensively studied in cancer that have not been previously explored with small molecules that harbor ENZ (PKMYT1, STEAP3, and NNMT) and PPI (HNF4A, MEF2B, and CBX2) binding sites. All binding sites were classified by the signaling pathways to which the protein that harbors them belongs using KEGG. In addition, binding sites were mapped onto structural protein-protein interaction networks to identify promising sites for drug discovery. Finally, we identify pockets that harbor missense mutations previously identified from analysis of the TCGA data. The occurrence of mutations in these binding sites provides new opportunities to develop small-molecule probes to explore their function in cancer.

A new class of orthosteric uPAR·uPA small-molecule antagonists are allosteric inhibitors of the uPAR·vitronectin interaction

The urokinase receptor (uPAR) is a GPI-anchored cell surface receptor that is at the center of an intricate network of protein-protein interactions. Its immediate binding partners are the serine proteinase urokinase (uPA), and vitronectin (VTN), a component of the extracellular matrix. uPA and VTN bind at distinct sites on uPAR to promote extracellular matrix degradation and integrin signaling, respectively. Here, we report the discovery of a new class of pyrrolone small-molecule inhibitors of the tight ∼1 nM uPAR·uPA protein-protein interaction. These compounds were designed to bind to the uPA pocket on uPAR. The highest affinity compound, namely 7, displaced a fluorescently labeled α-helical peptide (AE147-FAM) with an inhibition constant Ki of 0.7 μM and inhibited the tight uPAR·uPAATF interaction with an IC50 of 18 μM. Biophysical studies with surface plasmon resonance showed that VTN binding is highly dependent on uPA. This cooperative binding was confirmed as 7, which binds at the uPAR·uPA interface, also inhibited the distal VTN·uPAR interaction. In cell culture, 7 blocked the uPAR·uPA interaction in uPAR-expressing human embryonic kidney (HEK-293) cells and impaired cell adhesion to VTN, a process that is mediated by integrins. As a result, 7 inhibited integrin signaling in MDA-MB-231 cancer cells as evidenced by a decrease in focal adhesion kinase (FAK) phosphorylation and Rac1 GTPase activation. Consistent with these results, 7 blocked breast MDA-MB-231 cancer cell invasion with IC50 values similar to those observed in ELISA and surface plasmon resonance competition studies. Explicit-solvent molecular dynamics simulations show that the cooperativity between uPA and VTN is attributed to stabilization of uPAR motion by uPA. In addition, free energy calculations revealed that uPA stabilizes the VTNSMB·uPAR interaction through more favorable electrostatics and entropy. Disruption of the uPAR·VTNSMB interaction by 7 is consistent with the cooperative binding to uPAR by uPA and VTN. Interestingly, the VTNSMB·uPAR interaction was less favorable in the VTNSMB·uPAR·7 complex suggesting potential cooperativity between 7 and VTN. Compound 7 provides an excellent starting point for the development of more potent derivatives to explore uPAR biology

Mimicking Intermolecular Interactions of Tight Protein–Protein Complexes for Small-Molecule Antagonists

Tight protein–protein interactions (Kd1000 Å2) are highly challenging to disrupt with small molecules. Historically, the design of small molecules to inhibit protein–protein interactions has focused on mimicking the position of interface protein ligand side chains. Here, we explore mimicry of the pairwise intermolecular interactions of the native protein ligand with residues of the protein receptor to enrich commercial libraries for small-molecule inhibitors of tight protein–protein interactions. We use the high-affinity interaction (Kd=1 nm) between the urokinase receptor (uPAR) and its ligand urokinase (uPA) to test our methods. We introduce three methods for rank-ordering small molecules docked to uPAR: 1) a new fingerprint approach that represents uPA′s pairwise interaction energies with uPAR residues; 2) a pharmacophore approach to identify small molecules that mimic the position of uPA interface residues; and 3) a combined fingerprint and pharmacophore approach. Our work led to small molecules with novel chemotypes that inhibited a tight uPAR⋅uPA protein–protein interaction with single-digit micromolar IC50 values. We also report the extensive work that identified several of the hits as either lacking stability, thiol reactive, or redox active. This work suggests that mimicking the binding profile of the native ligand and the position of interface residues can be an effective strategy to enrich commercial libraries for small-molecule inhibitors of tight protein–protein interactions

Enrichment of Chemical Libraries Docked to Protein Conformational Ensembles and Application to Aldehyde Dehydrogenase 2

Molecular recognition is a complex process that involves a large ensemble of structures of the receptor and ligand. Yet, most structure-based virtual screening is carried out on a single structure typically from X-ray crystallography. Explicit-solvent molecular dynamics (MD) simulations offer an opportunity to sample multiple conformational states of a protein. Here we evaluate our recently developed scoring method SVMSP in its ability to enrich chemical libraries docked to MD structures of seven proteins from the Directory of Useful Decoys (DUD). SVMSP is a target-specific rescoring method that combines machine learning with statistical potentials. We find that enrichment power as measured by the area under the ROC curve (ROC-AUC) is not affected by increasing the number of MD structures. Among individual MD snapshots, many exhibited enrichment that was significantly better than the crystal structure, but no correlation between enrichment and structural deviation from crystal structure was found. We followed an innovative approach by training SVMSP scoring models using MD structures (SVMSPMD). The resulting models were applied to two difficult cases (p38 and CDK2) for which enrichment was not better than random. We found remarkable increase in enrichment power, particularly for p38, where the ROC-AUC increased by 0.30 to 0.85. Finally, we explored approaches for a priori identification of MD snapshots with high enrichment power from an MD simulation in the absence of active compounds. We found that the use of randomly selected compounds docked to the target of interest using SVMSP led to notable enrichment for EGFR and Src MD snapshots. SVMSP rescoring of protein–compound MD structures was applied for the search of small-molecule inhibitors of the mitochondrial enzyme aldehyde dehydrogenase 2 (ALDH2). Rank-ordering of a commercial library of 50 000 compounds docked to MD structures of ALDH2 led to five small-molecule inhibitors. Four compounds had IC50s below 5 μM. These compounds serve as leads for the design and synthesis of more potent and selective ALDH2 inhibitors

Structure-Based Target-Specific Screening Leads to Small-Molecule CaMKII Inhibitors

Target-specific scoring methods are more commonly used to identify small-molecule inhibitors among compounds docked to a target of interest. Top candidates that emerge from these methods have rarely been tested for activity and specificity across a family of proteins. In this study we docked a chemical library into CaMKIIδ, a member of the Ca2+ /calmodulin (CaM)-dependent protein kinase (CaMK) family, and re-scored the resulting protein-compound structures using Support Vector Machine SPecific (SVMSP), a target-specific method that we developed previously. Among the 35 selected candidates, three hits were identified, such as quinazoline compound 1 (KIN-1; N4-[7-chloro-2-[(E)-styryl]quinazolin-4-yl]-N1,N1-diethylpentane-1,4-diamine), which was found to inhibit CaMKIIδ kinase activity at single-digit micromolar IC50 . Activity across the kinome was assessed by profiling analogues of 1, namely 6 (KIN-236; N4-[7-chloro-2-[(E)-2-(2-chloro-4,5-dimethoxyphenyl)vinyl]quinazolin-4-yl]-N1,N1-diethylpentane-1,4-diamine), and an analogue of hit compound 2 (KIN-15; 2-[4-[(E)-[(5-bromobenzofuran-2-carbonyl)hydrazono]methyl]-2-chloro-6-methoxyphenoxy]acetic acid), namely 14 (KIN-332; N-[(E)-[4-(2-anilino-2-oxoethoxy)-3-chlorophenyl]methyleneamino]benzofuran-2-carboxamide), against 337 kinases. Interestingly, for compound 6, CaMKIIδ and homologue CaMKIIγ were among the top ten targets. Among the top 25 targets of 6, IC50 values ranged from 5 to 22 μm. Compound 14 was found to be not specific toward CaMKII kinases, but it does inhibit two kinases with sub-micromolar IC50 values among the top 25. Derivatives of 1 were tested against several kinases including several members of the CaMK family. These data afforded a limited structure-activity relationship study. Molecular dynamics simulations with explicit solvent followed by end-point MM-GBSA free-energy calculations revealed strong engagement of specific residues within the ATP binding pocket, and also changes in the dynamics as a result of binding. This work suggests that target-specific scoring approaches such as SVMSP may hold promise for the identification of small-molecule kinase inhibitors that exhibit some level of specificity toward the target of interest across a large number of proteins

Discovery of Novel Regulators of Aldehyde Dehydrogenase Isoenzymes

Over the past three years we have been involved in high-throughput screening in an effort to discover novel small molecular modulators of aldehyde dehydrogenase (ALDH) activity. In particular, we have been interested in both the activation and inhibition of the three commonly studied isoenzymes, ALDH1A1, ALDH2 and ALDH3A1, as their distinct, yet overlapping substrate specificities, present a particularly difficult challenge for inhibitor discovery and design. Activation of ALDH2 has been shown to benefit cardiovascular outcome following periods of ischemia and renewed interest in specific inhibition of ALDH2 has application for alcohol aversion therapy, and more recently, in cocaine addiction. In contrast, inhibition of either ALDH1A1 or ALDH3A1 has application in cancer treatments where the isoenzymes are commonly over-expressed and serve as markers for cancer stem cells. We are taking two distinct approaches for these screens: in vitro enzyme activity screens using chemical libraries and virtual computational screens using the structures of the target enzymes as filters for identifying potential inhibitors, followed by in vitro testing of their ability to inhibit their intended targets. We have identified selective inhibitors of each of these three isoenzymes with inhibition constants in the high nanomolar to low micromolar range from these screening procedures. Together, these inhibitors provide proof for concept that selective inhibition of these broad specificity general detoxication enzymes through small molecule discovery and design is possible

A Computational Investigation of Small-Molecule Engagement of Hot Spots at Protein–Protein Interaction Interfaces

The binding affinity of a protein–protein interaction is concentrated at amino acids known as hot spots. It has been suggested that small molecules disrupt protein–protein interactions by either (i) engaging receptor protein hot spots or (ii) mimicking hot spots of the protein ligand. Yet, no systematic studies have been done to explore how effectively existing small-molecule protein–protein interaction inhibitors mimic or engage hot spots at protein interfaces. Here, we employ explicit-solvent molecular dynamics simulations and end-point MM-GBSA free energy calculations to explore this question. We select 36 compounds for which high-quality binding affinity and cocrystal structures are available. Five complexes that belong to three classes of protein–protein interactions (primary, secondary, and tertiary) were considered, namely, BRD4•H4, XIAP•Smac, MDM2•p53, Bcl-xL•Bak, and IL-2•IL-2Rα. Computational alanine scanning using MM-GBSA identified hot-spot residues at the interface of these protein interactions. Decomposition energies compared the interaction of small molecules with individual receptor hot spots to those of the native protein ligand. Pharmacophore analysis was used to investigate how effectively small molecules mimic the position of hot spots of the protein ligand. Finally, we study whether small molecules mimic the effects of the native protein ligand on the receptor dynamics. Our results show that, in general, existing small-molecule inhibitors of protein–protein interactions do not optimally mimic protein–ligand hot spots, nor do they effectively engage protein receptor hot spots. The more effective use of hot spots in future drug design efforts may result in smaller compounds with higher ligand efficiencies that may lead to greater success in clinical trials